x <- c(13.193,15.234,14.718,16.961,13.945,15.876,16.226,18.316,16.748,17.904,17.209,18.950,17.225,18.710,17.236,18.687,17.580,19.568,17.381,19.580,17.260,18.661,15.658,18.674,15.908,17.475,17.725,19.562,16.368,19.555,17.743,19.867,15.703,19.324,18.162,19.074,15.323,19.704,18.375,18.352,13.927,17.795,16.761,18.902,16.239,19.158,18.279,15.698,16.239,18.431,18.414,19.801,14.995,18.706,18.232,19.409,16.263,19.017,20.298,19.891,15.203,17.845,17.502,18.532,15.737,17.770,17.224,17.601,14.940,18.507,17.635,19.392,15.699,17.661,18.243,19.643,15.770,17.344,17.229,17.322,16.152,17.919,16.918,18.114,16.308,17.759,16.021,17.952,15.954,17.762,16.610,17.751,15.458,18.106,15.990,15.349,13.185,15.409,16.007,16.633,14.800,15.974,15.693)
par4 = '4'
par3 = '0'
par2 = '0'
par1 = '1'
#'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
#Author: Prof. Dr. P. Wessa
#To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
#Source of accompanying publication: Office for Research, Development, and Education
#Technical description: Write here your technical program description (don't use hard returns!)
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="/var/www/rcomp/tmp/1vqzm1292928706.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="/var/www/rcomp/tmp/2ozgp1292928706.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
cpgram(x,main='Cumulative Periodogram')
dev.off()

#Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/www/rcomp/createtable")

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="/var/www/rcomp/tmp/3kref1292928706.tab") 

try(system("convert tmp/1vqzm1292928706.ps tmp/1vqzm1292928706.png",intern=TRUE))
try(system("convert tmp/2ozgp1292928706.ps tmp/2ozgp1292928706.png",intern=TRUE))