R version 2.6.0 (2007-10-03)
Copyright (C) 2007 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> x <- c(122302.01
+ ,109264.65
+ ,103674.75
+ ,103890.3
+ ,75512.66
+ ,83121.3
+ ,125096.81
+ ,74206.73
+ ,88481.63
+ ,111598.17
+ ,146919.48
+ ,150790.85
+ ,113780.5
+ ,110870.76
+ ,118785.32
+ ,112820.5
+ ,102188.92
+ ,97092.73
+ ,114067.82
+ ,89690.15
+ ,89267.9
+ ,96198.64
+ ,129599.75
+ ,169424.7
+ ,152510.91
+ ,121850.2
+ ,144737.64
+ ,121381.88
+ ,106894.86
+ ,94305.06
+ ,116800.42
+ ,77584.28
+ ,100680.88
+ ,106634.05
+ ,168390.77
+ ,211971.89
+ ,136163.28
+ ,168950.25
+ ,89816.88
+ ,85406.93
+ ,66055.52
+ ,73311.68
+ ,85674.51
+ ,82822.59
+ ,94277.63
+ ,100991.65
+ ,149245.88
+ ,208517.17
+ ,40733.51
+ ,121352.23
+ ,104020.11
+ ,99566.82
+ ,101352.17
+ ,106628.41
+ ,109696.95
+ ,248696.37
+ ,105628.33
+ ,120449.17
+ ,136547.7
+ ,140896.42
+ ,131509.91
+ ,95450.31
+ ,133592.64
+ ,110332.9
+ ,88110.54
+ ,64931.25
+ ,98446.22
+ ,84212.38
+ ,77519.55
+ ,124806.02
+ ,102185.94
+ ,151348.79
+ ,124378.28
+ ,101433.13
+ ,126724.22
+ ,87461.88
+ ,95288.27
+ ,129055.33
+ ,107753.06
+ ,96364.03
+ ,71662.75
+ ,125666.24
+ ,456841.51
+ ,167642.32
+ ,167154.73
+ ,139685.18
+ ,119275.2
+ ,122746.05
+ ,107337.43
+ ,112584.89
+ ,133183.08
+ ,121152.57
+ ,119815.6
+ ,122858.44
+ ,152077.17
+ ,157221.96
+ ,140435.08
+ ,101455.09
+ ,104791.29
+ ,77226.59
+ ,84477.43
+ ,66227.74
+ ,89076.23
+ ,108924.43
+ ,83926.11
+ ,91764.8
+ ,120892.76
+ ,129952.42
+ ,135865.14
+ ,105512.77
+ ,96486.62
+ ,78064.88
+ ,92370.22
+ ,98454.46
+ ,96703.93
+ ,83170.95)
> par4 = '12'
> par3 = '1'
> 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/html/rcomp/tmp/12vxn1200483312.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> r <- spectrum(x,main='Raw Periodogram')
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/2zj7y1200483312.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> cpgram(x,main='Cumulative Periodogram')
> dev.off()
null device 
          1 
> load(file='/var/www/html/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/rcomp/tmp/3t3no1200483312.tab") 
> 
> system("convert tmp/12vxn1200483312.ps tmp/12vxn1200483312.png")
> system("convert tmp/2zj7y1200483312.ps tmp/2zj7y1200483312.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.828   0.345   1.019