source('/home/pw/wessanet/cretab')



myrfcuid = 'r0678754'

x <- c(58.5,59.8,64.6,62.2,68,64.3,58.9,64.8,67.5,76.2,73.7,70.4,67.7,63.7,72.4,66,70.1,70.4,66.6,72.6,74,79,76.1,72.3,71.6,67.2,73.8,70.8,71.4,70.4,70.7,70.6,75.5,82.1,74.3,76.3,74.5,71.1,73.3,73.8,69,71.1,71.9,69,77.3,82.8,74,77.6,72.3,70.7,81,76.4,72.3,79.5,73.3,74.5,82.7,83.8,81.6,85.5,76.7,71.8,80.2,76.8,76.1,80.7,71.3,80.9,85,84.5,87.7,87.7,80.2,74.4,85.8,77,84.5,83.6,77.7,85.7,87.9,93.7,92.3,87,89.1,81.3,92.7,83.9,87.3,89.1,86.9,91.7,93,105.3,101.6,94.2,100.5,95.8,95.8,102.1,96,96.8,98.9,93.4,105.5,110.9,98.6,102.6,93.5,90.8,99.7,97.8,91.1,98.1,96,93.5,101.2,105.2,98.9,101.3,92.1,90.6,105.4,98.4,92.7,101.2,93.4,98.3,104.3,107,107.7,108.9,99.6,96.1,109,99.5,104.6,99.9,94.1,105.3,110.4,110.5,110,108.5,104.3,101.2,109.2,99.6,105.6,106.2,102.2,107.5,105.8,120.5,113.2,104.3,107.7,99.2,105.1,104.3,106.1,100.8,106.7,101.6,104.4,114.8,105.4,104,102,96.5,102.3,105.3,101.9,102.2,102.8,100.4,110.7,116.4,106,109.2,103,99.8,109.8,107.3,101.2,111.8,106.9,103.5,113.1,119.4,113.3,115,104.7,107.2,116.6,111.3,111.4,115,102.4,111.4,113.2,112.9,114.2,115.6,107.1,102.3,117.9,105.8,114.3,113.1,102.9,112.2)
par4 = '12'
par3 = '0'
par2 = '0'
par1 = '0.5'
par4 <- '12'
par3 <- '0'
par2 <- '0'
par1 <- '0.5'
#'GNU S' R Code compiled by R2WASP v. 1.2.327 (Fri, 21 Jul 2017 20:25:18 +0200)
#Author: root
#To cite this work: Wessa P., (2017), Spectral Analysis (v1.0.9) in Free Statistics Software (v$_version), Office for Research Development and Education, URL https://www.wessa.net/rwasp_spectrum.wasp/
#Source of accompanying publication: Office for Research, Development, and Education
#
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.numeric(par4)
if (par1 == 0) {
x <- log(x)
} else {
x <- (x ^ par1 - 1) / par1
}
if (par2 > 0) x <- diff(x,lag=1,difference=par2)
if (par3 > 0) x <- diff(x,lag=par4,difference=par3)
postscript(file="/home/pw/wessanet/rcomp/tmp/1p0lc1513855596.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
r <- spectrum(x,main='Raw Periodogram')
dev.off()
postscript(file="/home/pw/wessanet/rcomp/tmp/21qj01513855596.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
cpgram(x,main='Cumulative Periodogram')
dev.off()

a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Raw Periodogram',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Parameter',header=TRUE)
a<-table.element(a,'Value',header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Box-Cox transformation parameter (lambda)',header=TRUE)
a<-table.element(a,par1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of non-seasonal differencing (d)',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of seasonal differencing (D)',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Seasonal Period (s)',header=TRUE)
a<-table.element(a,par4)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Frequency (Period)',header=TRUE)
a<-table.element(a,'Spectrum',header=TRUE)
a<-table.row.end(a)
for (i in 1:length(r$freq)) {
a<-table.row.start(a)
mylab <- round(r$freq[i],4)
mylab <- paste(mylab,' (',sep='')
mylab <- paste(mylab,round(1/r$freq[i],4),sep='')
mylab <- paste(mylab,')',sep='')
a<-table.element(a,mylab,header=TRUE)
a<-table.element(a,round(r$spec[i],6))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file="/home/pw/wessanet/rcomp/tmp/3y9cl1513855596.tab") 

try(system("convert /home/pw/wessanet/rcomp/tmp/1p0lc1513855596.ps /home/pw/wessanet/rcomp/tmp/1p0lc1513855596.png",intern=TRUE))
try(system("convert /home/pw/wessanet/rcomp/tmp/21qj01513855596.ps /home/pw/wessanet/rcomp/tmp/21qj01513855596.png",intern=TRUE))