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



myrfcuid = 'r0708166'

x <- c(78,100.1,113.2,93.1,115.4,103.3,45.1,104.7,111.3,111.5,100.9,82.1,85.4,97.7,106.6,92.6,109.2,110,52.5,105.3,102.3,118.5,100,74.4,89.2,91.9,107,103.6,101.8,105.1,55.5,92.1,109.8,112.7,98.5,70.3,84.5,91.1,107.6,102.2,96,107.3,59.9,90.2,116.3,115.6,92,76.5,87.9,95.8,116.9,102.9,95.8,117.3,52.8,100.1,116.3,111.8,98.5,86.2,79.9,92.3,100.5,112.5,101.1,121.5,49.6,104.8,120.4,108.3,105.2,85.7,86.8,95.1,117,100.1,112.3,119.6,51.8,105.5,119.9,115.4,112.8,85.1,96.2,103.6,119.9,103.7,109,119.6,57,109.2,112.6,126,109.7,80.1,105.8,114.1,98.3,125.3,111.6,119.7,65,99,124.5,119,98.8,81.8,90.3,102,119.3,104.3,102.8,118.8,60.9,101,122.6,122.2,95,75.6,83.1,89.8,126.1,108.6,98.9,124.3,56.8,102.7,121.7,118.2,101,69,88.6,109.6,128.2,102,122.7,110.5,54,108.1,125,114.1,112.4,87.3,95.4,96.9,125.8,102,112.5,118.9,62.7,110,114.7,124.4,111.9,77,84.1,96.5,106.8,107.9,107.5,114.3,66.6,97.9,117.8,123.8,103.3,84.2,103.6,103.6,112.2,102.7,100.8,109.4,63.5,92.3,119.2,121.5,97.6,78.3,95.6,97.9,114.4,100.9,94.4,117.2,61,95.8,116.2,118.5,94.3,74.4,94.9,102,102.9,109.5,99.7,118.3,56.2,100.3,116.9,108.7,93.9,85.3,85.3,102.4,121.6,91.4,110.2,112.7,55.7,100.1)
par4 = '12'
par3 = '1'
par2 = '0'
par1 = '1'
par4 <- '12'
par3 <- '1'
par2 <- '0'
par1 <- '1'
#'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/1fm7c1513086527.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/2tf5l1513086527.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/3oiki1513086527.tab") 

try(system("convert /home/pw/wessanet/rcomp/tmp/1fm7c1513086527.ps /home/pw/wessanet/rcomp/tmp/1fm7c1513086527.png",intern=TRUE))
try(system("convert /home/pw/wessanet/rcomp/tmp/2tf5l1513086527.ps /home/pw/wessanet/rcomp/tmp/2tf5l1513086527.png",intern=TRUE))