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(358.59
+ ,362.96
+ ,362.42
+ ,364.97
+ ,364.04
+ ,361.06
+ ,358.48
+ ,352.96
+ ,359.59
+ ,360.39
+ ,357.40
+ ,362.93
+ ,364.55
+ ,365.73
+ ,364.70
+ ,364.65
+ ,359.43
+ ,362.14
+ ,356.97
+ ,354.82
+ ,353.17
+ ,357.06
+ ,356.18
+ ,355.01
+ ,355.65
+ ,357.31
+ ,357.07
+ ,357.91
+ ,358.48
+ ,358.97
+ ,351.77
+ ,352.16
+ ,359.08
+ ,360.35
+ ,359.53
+ ,359.30
+ ,358.41
+ ,359.68
+ ,355.31
+ ,357.08
+ ,349.71
+ ,354.13
+ ,345.49
+ ,341.69
+ ,344.25
+ ,340.17
+ ,342.47
+ ,344.43
+ ,333.23
+ ,339.72
+ ,342.61
+ ,346.36
+ ,339.09
+ ,339.73
+ ,341.12
+ ,335.94
+ ,333.46
+ ,335.66
+ ,341.12
+ ,342.21
+ ,342.62
+ ,346.06
+ ,344.43
+ ,346.65
+ ,343.74
+ ,335.67
+ ,342.75
+ ,341.77
+ ,345.84
+ ,346.52
+ ,350.79
+ ,345.44
+ ,345.87
+ ,338.48
+ ,337.21
+ ,340.81
+ ,339.86
+ ,342.86
+ ,343.33
+ ,341.73
+ ,351.38
+ ,351.13
+ ,345.99
+ ,347.55
+ ,346.02
+ ,345.29
+ ,347.03
+ ,348.01
+ ,345.48
+ ,349.40
+ ,351.05
+ ,349.70
+ ,350.86
+ ,354.45
+ ,355.30
+ ,357.48
+ ,355.24
+ ,351.79
+ ,355.22
+ ,351.02
+ ,350.28
+ ,350.17
+ ,348.16
+ ,340.30
+ ,343.75
+ ,344.71
+ ,344.13
+ ,342.14
+ ,345.04
+ ,346.02
+ ,346.43
+ ,347.07
+ ,339.33
+ ,339.10
+ ,337.19
+ ,339.58
+ ,327.85
+ ,326.81
+ ,321.73
+ ,320.45
+ ,327.69
+ ,323.95
+ ,320.47
+ ,322.13
+ ,316.34
+ ,314.78
+ ,308.90
+ ,308.62
+ ,314.41
+ ,306.88
+ ,310.60
+ ,321.60
+ ,321.50
+ ,325.68
+ ,324.35
+ ,320.01
+ ,326.88
+ ,332.39
+ ,331.48
+ ,332.62
+ ,324.79
+ ,327.12
+ ,328.91
+ ,328.37
+ ,324.83
+ ,325.90
+ ,326.18
+ ,328.94
+ ,333.78
+ ,328.06
+ ,325.87
+ ,325.41
+ ,318.86
+ ,319.13
+ ,310.16
+ ,311.73
+ ,306.54
+ ,311.16
+ ,311.98
+ ,306.72
+ ,308.05
+ ,300.76
+ ,301.90
+ ,293.09
+ ,292.76
+ ,294.58
+ ,289.90
+ ,296.69
+ ,297.21
+ ,293.31
+ ,296.25
+ ,298.60
+ ,296.87
+ ,301.02
+ ,304.73
+ ,301.92
+ ,295.72
+ ,293.18
+ ,298.35
+ ,297.99
+ ,299.85
+ ,299.85
+ ,304.45
+ ,299.45
+ ,298.14
+ ,298.78
+ ,297.02
+ ,301.33
+ ,294.96
+ ,296.69
+ ,300.73
+ ,301.96
+ ,297.38
+ ,293.87
+ ,285.96
+ ,285.41
+ ,283.70
+ ,284.76
+ ,277.11
+ ,274.73
+ ,274.73
+ ,274.73
+ ,274.73
+ ,274.69
+ ,275.42
+ ,264.15
+ ,276.24
+ ,268.88
+ ,277.97
+ ,280.49
+ ,281.09
+ ,276.16
+ ,272.58
+ ,270.94
+ ,284.31
+ ,283.94
+ ,284.18
+ ,282.83
+ ,283.84
+ ,282.71
+ ,279.29
+ ,280.70
+ ,274.47
+ ,273.44
+ ,275.49
+ ,279.46
+ ,280.19
+ ,288.21
+ ,284.80
+ ,281.41
+ ,283.39
+ ,287.97
+ ,290.77
+ ,290.60
+ ,289.67
+ ,289.84
+ ,298.55
+ ,296.07
+ ,297.14
+ ,295.34
+ ,296.25
+ ,294.30
+ ,296.15
+ ,296.49
+ ,298.05
+ ,301.03
+ ,300.52
+ ,301.50
+ ,296.93
+ ,289.84
+ ,291.44
+ ,286.88
+ ,286.74
+ ,288.93
+ ,292.19
+ ,295.39
+ ,295.86
+ ,293.36
+ ,292.86
+ ,292.73
+ ,296.73
+ ,285.02
+ ,285.24
+ ,288.62
+ ,283.36
+ ,285.84
+ ,291.48
+ ,291.41
+ ,287.77
+ ,284.97
+ ,286.05
+ ,278.19
+ ,281.21
+ ,277.92
+ ,280.08
+ ,269.24
+ ,268.48
+ ,268.83
+ ,269.54
+ ,262.37
+ ,265.12
+ ,265.34
+ ,263.32
+ ,267.18
+ ,260.75
+ ,261.78
+ ,257.27
+ ,255.63
+ ,251.39
+ ,259.49
+ ,261.18
+ ,261.65
+ ,262.01
+ ,265.23
+ ,268.10
+ ,262.27
+ ,263.59
+ ,257.85
+ ,265.69
+ ,271.15
+ ,266.69
+ ,265.77
+ ,262.32
+ ,270.48
+ ,273.03
+ ,269.13
+ ,280.65
+ ,282.75
+ ,281.44
+ ,281.99
+ ,282.86
+ ,287.21
+ ,283.11
+ ,280.66
+ ,282.39
+ ,280.83
+ ,284.71
+ ,279.99
+ ,283.50
+ ,284.88
+ ,288.60
+ ,284.80
+ ,287.20
+ ,286.22
+ ,286.54
+ ,279.58
+ ,283.08
+ ,288.88
+ ,280.18
+ ,284.16
+ ,290.57
+ ,286.82
+ ,273.00
+ ,278.69
+ ,264.54
+ ,271.92
+ ,283.60
+ ,269.25
+ ,263.58
+ ,264.16
+ ,268.85
+ ,269.67
+ ,249.41
+ ,268.99
+ ,268.65
+ ,260.16
+ ,256.55
+ ,251.47
+ ,234.93
+ ,232.96
+ ,215.49
+ ,213.68
+ ,236.07
+ ,235.41
+ ,214.77
+ ,225.85
+ ,224.64
+ ,238.26
+ ,232.44
+ ,222.50
+ ,225.28
+ ,220.49
+ ,216.86
+ ,234.70
+ ,230.06
+ ,238.27
+ ,238.56
+ ,242.70
+ ,249.14
+ ,234.89
+ ,227.78
+ ,234.04
+ ,230.70
+ ,230.17
+ ,218.23
+ ,232.20
+ ,220.76
+ ,215.60
+ ,217.69
+ ,204.35
+ ,191.44
+ ,203.84
+ ,211.86
+ ,210.57
+ ,219.57
+ ,219.98
+ ,226.01
+ ,207.04
+ ,212.52
+ ,217.92
+ ,210.45
+ ,218.53
+ ,223.32
+ ,218.76
+ ,217.63)
> par4 = '1'
> par3 = '0'
> 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/1agnu1260968074.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/2i15i1260968074.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/3wtye1260968075.tab") 
> 
> try(system("convert tmp/1agnu1260968074.ps tmp/1agnu1260968074.png",intern=TRUE))
character(0)
> try(system("convert tmp/2i15i1260968074.ps tmp/2i15i1260968074.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.725   0.336   0.911