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



myrfcuid = 'r0370864'

x <- c(1932.8,1861.4,2170.2,1999.6,2225.5,2195.7,2713.1,2412,2568.3,2623.7,3185.5,2722.6,3046.3,2854.2,3337.6,2920.3,3058.3,2933.7,3773.4,3193.5,3472.2,3345.5,4028.4,3463.1,3675.4,3500.8,4142.1,3598,3765.3,3557.7,4303.6,3620.1,3691.1,3678.1,4505.8,3695,3894.1,3718.9,4749.8,3855.9,4011.7,3907.6,4812.5,4071.3,4163.4,4077.6,5109.2,4207.6,4320.8,4396.9,5358.8)
par4 = '1'
par3 = '0'
par2 = '1'
par1 = '1'
par4 <- '1'
par3 <- '0'
par2 <- '1'
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/1fmax1482260656.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/2wczu1482260656.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/3gvv81482260656.tab") 

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