x <- c(5,4,5,6,6,6,7,8,7,8,7,8,8,9,9,8,9,9,10,11,12,13,13,13,14,14,15,15,16,16,17,18,19,20,22,20,22,25,24,25,28,26,27,26,25,27,28,30,31,32,34,34,33,32,34,36,37,40,38,38,36,40,40,42,44,45,47,49,47,49,52,50,50,57,58,58,58,61,61,64,68,40,34,46,36,34,45,55,50,56,72,76,78,77,90,88,97,93,84,67,72,75,71,75,90,78,73,62,65,61,58,33,39,56,79,82,79,73,87,85,83,82,83,92,95,97,87,84,84,89,103,106,109,106,105,115,120,124,121,131,139,133,119,123,120,128,134,126,115,106,99,100,99,99,100,100,108,109,115,114,108,113,118,122,118,121,118,121,121,112,119,116,110,111,106,108)
par4 = '1'
par3 = '0'
par2 = '0'
par1 = '1'
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/html/freestat/rcomp/tmp/15wpu1292743473.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/html/freestat/rcomp/tmp/25wpu1292743473.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/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/www/html/freestat/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/html/freestat/rcomp/tmp/3cf461292743473.tab") 

try(system("convert tmp/15wpu1292743473.ps tmp/15wpu1292743473.png",intern=TRUE))
try(system("convert tmp/25wpu1292743473.ps tmp/25wpu1292743473.png",intern=TRUE))