| 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)
 bitmap(file='test1.png')
 r <- spectrum(x,main='Raw Periodogram')
 dev.off()
 bitmap(file='test2.png')
 cpgram(x,main='Cumulative Periodogram')
 dev.off()
 load(file='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='mytable.tab')
 
 |