R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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(897262
+ ,1133132
+ ,1384548
+ ,2324057
+ ,2502808
+ ,2516762
+ ,5579822
+ ,4945991
+ ,2019915
+ ,1830905
+ ,1251016
+ ,949902
+ ,923000
+ ,1215747
+ ,1479112
+ ,2371781
+ ,2521576
+ ,2350559
+ ,5673323
+ ,4414295
+ ,2016902
+ ,1958302
+ ,1284086
+ ,1186305
+ ,957833
+ ,1255719
+ ,1482709
+ ,2361136
+ ,2508100
+ ,2254488
+ ,5669953
+ ,4227480
+ ,2067790
+ ,1958419
+ ,1318158
+ ,1287921
+ ,1076982
+ ,1293669
+ ,1582053
+ ,2393005
+ ,2310531
+ ,2597899
+ ,5507587
+ ,4194133
+ ,2185092
+ ,2122018
+ ,1413348
+ ,1338342
+ ,1052655
+ ,1370046
+ ,1887027
+ ,2448017
+ ,2550796
+ ,2655837
+ ,5269499
+ ,4247405
+ ,2109722
+ ,2143145
+ ,1582013
+ ,1413221
+ ,1118520
+ ,1478655
+ ,2000108
+ ,2085234
+ ,2651805
+ ,2522176
+ ,5170142
+ ,4150129
+ ,2104254
+ ,2211398
+ ,1505900
+ ,1524305
+ ,1093144
+ ,1449647
+ ,1771197
+ ,2445932
+ ,2678945
+ ,2400737
+ ,4796880
+ ,4118001
+ ,2125714
+ ,2125515
+ ,1508760
+ ,1508765
+ ,1091075
+ ,1514814
+ ,1748997
+ ,2424406
+ ,2747942
+ ,2377332
+ ,5210706
+ ,3882821
+ ,2197469
+ ,2271155
+ ,1618917
+ ,1391579
+ ,1143249
+ ,1445785
+ ,1870242
+ ,2597788
+ ,2436231
+ ,2684184
+ ,4705109
+ ,4331347
+ ,2369192
+ ,2283947
+ ,1749607
+ ,1598601
+ ,1221234
+ ,1497778
+ ,1823567
+ ,2489908
+ ,2532837
+ ,2456065
+ ,4627018
+ ,4276894
+ ,2314950
+ ,2238987
+ ,1652753
+ ,1561968
+ ,1115878
+ ,1596714
+ ,1910242
+ ,2286450
+ ,2772441
+ ,2394538
+ ,4715128
+ ,4402420
+ ,2325392
+ ,2306683
+ ,1725282
+ ,1541370
+ ,1168142
+ ,1457835
+ ,1816380
+ ,2446552
+ ,2575774
+ ,2537852
+ ,4728097
+ ,4372685
+ ,2302672
+ ,2346402
+ ,1689915
+ ,1576183)
> par4 = '12'
> 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/html/rcomp/tmp/1lr011229623199.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/2ohx61229623199.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> cpgram(x,main='Cumulative Periodogram')
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> 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/3iszf1229623199.tab") 
> 
> system("convert tmp/1lr011229623199.ps tmp/1lr011229623199.png")
> system("convert tmp/2ohx61229623199.ps tmp/2ohx61229623199.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.581   0.377   1.346