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



myrfcuid = 'r0382442'

x <- c(13663,11635,9606,8784.5,9415.5,10418,11344.5,11271,11895,12152.5,12731,12951,10692,8563.5,6217,5562,6294.5,7422,9254.5,10607,11268,12041,12962.5,12200.5,10400.5,8765,7000,6677,7318,7999,8762,9696,10373,10682.5,10935.5,10815.5,8669,7079.5,5640,5238.5,5777.5,6479,7290,7343,7810.5,8171.5,8532,8719,7281.5,5923.5,4837,4675.5,4585.5,5083,5766,6201,6778,7393.5,7849.5,8282.5,7610,6192.5,4693.5,4869,5149,5648.5,6230.5,7032,7727,8087.5,8443,9002,7717.5,6374.5,4995.5,4655,5198,5501,6119.5,6922,7390,7466.5,7773,7865,6567,5132.5,3656.5,3623,4045.5,4617,5374,6022.5,6464.5,7058,7484.5,7955,6801,5499,4179.5,4305.5,3304,5773.5,6419.5,6938,7760,8224,8381,8667,7304.5,5565.5,4023,3932.5,4508.5,5491,6284)
par4 = '12'
par3 = '0'
par2 = '0'
par1 = '1'
par4 <- '12'
par3 <- '0'
par2 <- '0'
par1 <- '1'
#'GNU S' R Code compiled by R2WASP v. 1.2.327 (Tue, 16 Aug 2016 15:29:54 +0200)
#Author: root
#To cite this work: Wessa P., (2012), Spectral Analysis (v1.0.8) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://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/1rn601481306683.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/2eejf1481306683.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/3yv571481306683.tab") 

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