R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(168.67
+ ,164.83
+ ,184.38
+ ,180.81
+ ,190.54
+ ,181.41
+ ,155.67
+ ,135.99
+ ,125.88
+ ,126.09
+ ,114.86
+ ,127.98
+ ,127.98
+ ,125.11
+ ,125.93
+ ,128.2
+ ,125.93
+ ,111.94
+ ,120.01
+ ,124.09
+ ,126.02
+ ,136.41
+ ,143.79
+ ,141.67
+ ,143.9
+ ,155
+ ,144.83
+ ,141.4
+ ,137
+ ,141.02
+ ,131.11
+ ,132.83
+ ,136.73
+ ,141.18
+ ,137.86
+ ,133.79
+ ,128.53
+ ,125.87
+ ,124.27
+ ,123.96
+ ,128.15
+ ,126.4
+ ,127.86
+ ,129.31
+ ,132.56
+ ,141.28
+ ,145.55
+ ,146.54
+ ,143.14
+ ,145.72
+ ,148.21
+ ,150.4
+ ,149.94
+ ,146.66
+ ,143.37
+ ,145.29
+ ,140.24
+ ,136.12
+ ,140.25
+ ,140.64
+ ,145.58
+ ,143.73
+ ,141.27
+ ,140.66
+ ,141.94
+ ,141.16
+ ,134.31
+ ,132.93
+ ,133.07
+ ,140.48
+ ,154.85
+ ,196.77
+ ,235.3
+ ,226.52
+ ,237.62
+ ,224.07
+ ,208.74
+ ,174.54
+ ,170.63
+ ,172.23
+ ,198.36
+ ,175.91
+ ,154.63
+ ,134.31
+ ,121.75
+ ,119.6
+ ,102.04
+ ,106.3
+ ,116.38
+ ,103.72
+ ,98.56
+ ,100.9
+ ,110
+ ,118.26
+ ,124.77
+ ,125.22
+ ,126.38
+ ,137.14
+ ,134.74
+ ,134.3
+ ,136.39
+ ,141.83
+ ,139.24
+ ,128.89
+ ,134.83
+ ,130.43
+ ,132.09
+ ,144.95
+ ,149.5
+ ,137.57
+ ,139.38
+ ,143.06
+ ,138.65
+ ,123.21
+ ,85.91
+ ,77.4
+ ,77.84
+ ,67.76
+ ,70.72
+ ,72.55
+ ,75.83
+ ,84.01
+ ,93.96
+ ,93.73
+ ,92.02
+ ,88.26
+ ,86.48
+ ,94.42
+ ,94.92
+ ,91.41
+ ,84.84
+ ,89.89
+ ,86.32
+ ,89.57
+ ,93.72
+ ,92.27
+ ,87.59
+ ,85.5
+ ,82.81
+ ,81.62
+ ,87.45
+ ,79.86
+ ,78.52
+ ,75.1
+ ,72.99
+ ,67.88
+ ,70.14
+ ,65.43
+ ,60.26
+ ,58.38
+ ,57.68
+ ,52.42
+ ,52.73
+ ,61.4
+ ,67.13
+ ,77.46
+ ,68.66
+ ,67.46
+ ,62.77
+ ,56.88
+ ,61.48
+ ,61.99
+ ,71.56
+ ,76.56
+ ,79.82
+ ,75.05
+ ,77.07
+ ,80
+ ,77.21
+ ,82.16
+ ,85.57
+ ,89.23
+ ,121.98
+ ,142.56
+ ,217.67
+ ,198.07
+ ,220.1
+ ,198.68
+ ,181.64
+ ,167.47
+ ,172.33
+ ,168.71
+ ,178.22
+ ,172.81
+ ,168.83
+ ,152.25
+ ,143.83
+ ,151.41
+ ,131.87
+ ,125.38
+ ,123.23
+ ,103.99
+ ,109.38
+ ,123.79
+ ,119.05
+ ,122.01
+ ,128.56
+ ,127.91
+ ,120.47
+ ,122.49
+ ,114.05
+ ,120.62
+ ,119.61
+ ,115.01
+ ,131.83
+ ,167.2
+ ,193.82
+ ,204.43
+ ,264.5
+ ,212.55
+ ,186.52
+ ,185.17
+ ,184.38
+ ,161.45
+ ,154.15
+ ,174.25
+ ,175.04
+ ,175.87
+ ,154.82
+ ,147.08
+ ,134.35
+ ,121.56
+ ,113.86
+ ,119.89
+ ,108.07
+ ,107.07
+ ,115.14
+ ,116.03
+ ,111.48
+ ,103.24
+ ,103.23
+ ,99.69
+ ,108.91
+ ,104.21
+ ,90.85
+ ,87.64
+ ,81.06
+ ,92.2
+ ,114.02
+ ,123.56
+ ,109.17
+ ,101.65
+ ,97.95
+ ,92.56
+ ,91.76
+ ,84.1
+ ,84.67
+ ,74.52
+ ,73.83
+ ,75.37
+ ,70.47
+ ,64.5
+ ,64.98
+ ,66.94
+ ,65.93
+ ,65.51
+ ,68.94
+ ,63.67
+ ,58.47
+ ,59.68
+ ,57.71
+ ,56.53
+ ,58.96
+ ,55.6
+ ,57.34
+ ,60.51
+ ,66.38
+ ,65.78
+ ,58.43
+ ,55.16
+ ,53.09
+ ,52.02
+ ,57.58
+ ,64.05
+ ,70.18
+ ,63.86
+ ,65.22
+ ,67.6
+ ,61.66
+ ,65.32
+ ,66.18
+ ,61.34
+ ,62.29
+ ,63.6
+ ,65.51
+ ,62.58
+ ,62.36
+ ,64.88
+ ,73.73
+ ,77.51
+ ,77.47
+ ,74.34
+ ,75.81
+ ,82.16
+ ,73.96
+ ,73.17
+ ,80.99
+ ,79.81
+ ,89.51
+ ,102.57
+ ,107.11
+ ,122.23
+ ,134.69
+ ,128.79
+ ,126.16
+ ,119.98
+ ,108.45
+ ,108.43
+ ,98.17
+ ,106.09
+ ,108.81
+ ,103.03
+ ,124.36
+ ,118.52
+ ,112.2
+ ,114.71
+ ,107.96
+ ,101.21
+ ,102.77
+ ,112.13
+ ,109.36
+ ,110.91
+ ,123.57
+ ,129.95
+ ,124.46
+ ,122.34
+ ,116.61
+ ,114.59
+ ,112.52
+ ,118.67
+ ,116.8
+ ,123.63
+ ,128.04
+ ,134.57
+ ,130.33
+ ,136.47
+ ,139.05
+ ,158.21
+ ,148.07
+ ,137.74
+ ,139.74
+ ,144.08
+ ,145.35
+ ,145.77
+ ,140.56
+ ,121.41
+ ,120.44
+ ,116.97
+ ,128.03
+ ,128.51
+ ,127.76
+ ,134.58
+ ,147.64
+ ,144.46
+ ,137.6
+ ,146.87
+ ,145.67
+ ,151.95
+ ,150.23
+ ,155.86
+ ,154.4
+ ,156.36
+ ,162.13
+ ,171.06
+ ,174.01
+ ,193.52
+ ,205.26
+ ,212.8
+ ,222.1)
> par4 = '12'
> par3 = '1'
> par2 = '1'
> 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/12tjp1291716631.ps",horizontal=F,onefile=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/2v2ja1291716631.ps",horizontal=F,onefile=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/3exnp1291716631.tab") 
> try(system("convert tmp/12tjp1291716631.ps tmp/12tjp1291716631.png",intern=TRUE))
character(0)
> try(system("convert tmp/2v2ja1291716631.ps tmp/2v2ja1291716631.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.724   0.378   1.776