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(353.4
+ ,329.08
+ ,331.89
+ ,339.94
+ ,330.8
+ ,361.26
+ ,358.02
+ ,356.15
+ ,322.56
+ ,306.1
+ ,303.99
+ ,322.23
+ ,330.2
+ ,343.91
+ ,367.07
+ ,375.22
+ ,375.35
+ ,389.81
+ ,371.18
+ ,387.18
+ ,395.43
+ ,387.86
+ ,392.46
+ ,375.11
+ ,417.03
+ ,408.79
+ ,412.68
+ ,403.67
+ ,414.95
+ ,415.35
+ ,408.2
+ ,424.19
+ ,414.03
+ ,417.8
+ ,418.66
+ ,431.35
+ ,435.7
+ ,438.78
+ ,443.38
+ ,451.67
+ ,440.19
+ ,450.23
+ ,450.54
+ ,448.13
+ ,463.55
+ ,458.93
+ ,467.83
+ ,461.93
+ ,466.51
+ ,481.6
+ ,467.19
+ ,445.66
+ ,450.91
+ ,456.5
+ ,444.27
+ ,458.28
+ ,475.49
+ ,462.69
+ ,472.26
+ ,453.55
+ ,459.21
+ ,470.42
+ ,487.39
+ ,500.7
+ ,514.76
+ ,533.4
+ ,544.75
+ ,562.06
+ ,561.88
+ ,584.41
+ ,581.5
+ ,605.37
+ ,615.93
+ ,636.02
+ ,640.43
+ ,645.5
+ ,654.17
+ ,669.12
+ ,670.63
+ ,639.95
+ ,651.99
+ ,687.31
+ ,705.27
+ ,757.02
+ ,740.74
+ ,786.16
+ ,790.82
+ ,757.12
+ ,801.34
+ ,848.28
+ ,885.14
+ ,954.29
+ ,899.47
+ ,947.28
+ ,914.62
+ ,955.4
+ ,970.43
+ ,980.28
+ ,1049.34
+ ,1101.75
+ ,1111.75
+ ,1090.82
+ ,1133.84
+ ,1120.67
+ ,957.28
+ ,1017.01
+ ,1098.67
+ ,1163.63
+ ,1129.23
+ ,1279.64
+ ,1238.33
+ ,1286.37
+ ,1335.18
+ ,1301.84
+ ,1372.71
+ ,1328.72
+ ,1320.41
+ ,1282.71
+ ,1362.93
+ ,1388.91
+ ,1469.25
+ ,1394.46
+ ,1366.42
+ ,1498.58
+ ,1452.43
+ ,1420.6
+ ,1454.6
+ ,1430.83
+ ,1517.68
+ ,1436.52
+ ,1429.4
+ ,1314.95
+ ,1320.28
+ ,1366.01
+ ,1239.94
+ ,1160.33
+ ,1249.46
+ ,1255.82
+ ,1224.42
+ ,1211.23
+ ,1133.58
+ ,1040.94
+ ,1059.78
+ ,1139.45
+ ,1148.08
+ ,1130.2
+ ,1106.73
+ ,1147.39
+ ,1076.92
+ ,1067.14
+ ,989.82
+ ,911.62
+ ,916.07
+ ,815.28
+ ,885.76
+ ,936.31
+ ,879.82
+ ,855.7
+ ,841.15
+ ,848.18
+ ,916.92
+ ,963.59
+ ,974.5
+ ,990.31
+ ,1008.01
+ ,995.97
+ ,1050.71
+ ,1058.2
+ ,1111.92
+ ,1131.13
+ ,1144.94
+ ,1113.89
+ ,1107.3
+ ,1120.68
+ ,1140.84
+ ,1101.72
+ ,1104.24
+ ,1114.58
+ ,1130.2
+ ,1173.78
+ ,1211.92
+ ,1181.27
+ ,1203.6
+ ,1180.59
+ ,1156.85
+ ,1191.5
+ ,1191.33
+ ,1234.18
+ ,1220.33
+ ,1228.81
+ ,1207.01
+ ,1249.48
+ ,1248.29
+ ,1280.08
+ ,1280.66
+ ,1302.88
+ ,1310.61
+ ,1270.05
+ ,1270.06
+ ,1278.53
+ ,1303.8
+ ,1335.83
+ ,1377.76
+ ,1400.63
+ ,1418.03
+ ,1437.9
+ ,1406.8
+ ,1420.83
+ ,1482.37
+ ,1530.63
+ ,1504.66
+ ,1455.18
+ ,1473.96
+ ,1527.29
+ ,1545.79
+ ,1479.63
+ ,1467.97
+ ,1378.6
+ ,1330.45
+ ,1326.41
+ ,1385.97
+ ,1399.62
+ ,1276.69
+ ,1269.42
+ ,1287.83
+ ,1164.17
+ ,968.67
+ ,888.61
+ ,902.99
+ ,823.09
+ ,729.57
+ ,793.59
+ ,872.74
+ ,923.26
+ ,920.82
+ ,990.22
+ ,1019.52
+ ,1054.91
+ ,1036.18
+ ,1098.89)
> par4 = '12'
> par3 = '0'
> par2 = '1'
> par1 = '0.0'
> #'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/18ofn1261230795.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/2p10v1261230795.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/3iek01261230795.tab") 
> 
> try(system("convert tmp/18ofn1261230795.ps tmp/18ofn1261230795.png",intern=TRUE))
character(0)
> try(system("convert tmp/2p10v1261230795.ps tmp/2p10v1261230795.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.616   0.326   0.749