R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 The R Foundation for Statistical Computing
Platform: i686-pc-linux-gnu (32-bit)

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(1356.876334
+ ,1267.205322
+ ,1296.765977
+ ,1088.768988
+ ,974.9813884
+ ,901.3596576
+ ,697.4384802
+ ,619.518781
+ ,446.3076104
+ ,279.7073872
+ ,388.5496968
+ ,290.3616669
+ ,278.7900756
+ ,433.3678289
+ ,436.029714
+ ,642.0796877
+ ,800.8666721
+ ,1005.060868
+ ,986.1196671
+ ,1320.620732
+ ,1293.652614
+ ,1315.564345
+ ,1438.509653
+ ,1311.68492
+ ,1281.851742
+ ,1229.184547
+ ,1157.483378
+ ,1049.930474
+ ,720.8534348
+ ,756.8141819
+ ,587.8054359
+ ,440.0632504
+ ,476.8424862
+ ,376.9041092
+ ,408.9909215
+ ,335.6074609
+ ,510.7385144
+ ,784.0262709
+ ,868.9048671
+ ,930.5891366
+ ,1232.387738
+ ,1433.973937
+ ,1441.341759
+ ,1560.894767
+ ,1570.121112
+ ,1539.736335
+ ,1509.516021
+ ,1264.261299
+ ,1280.1511
+ ,1104.198071
+ ,867.0732077
+ ,832.9038416
+ ,551.1656185
+ ,578.9382082
+ ,548.7690841
+ ,526.7505509
+ ,543.7872669
+ ,387.4285031
+ ,583.0257696
+ ,858.5103593
+ ,835.2939603
+ ,1234.576103
+ ,1307.198241
+ ,1550.368145
+ ,1513.93417
+ ,1607.433912
+ ,1768.513076
+ ,1778.655313
+ ,1668.394241
+ ,1383.577733
+ ,1271.618447
+ ,1140.579064
+ ,998.7028247
+ ,804.5968421
+ ,743.1134557
+ ,661.133487
+ ,635.0605076
+ ,334.5931783
+ ,502.1494165
+ ,450.8321675
+ ,700.3681247
+ ,799.3587889
+ ,1068.293463
+ ,1229.161035
+ ,1369.983001
+ ,1586.670588
+ ,1649.664042
+ ,1871.680788
+ ,1816.035687
+ ,1919.694672
+ ,1817.473725
+ ,1707.303794
+ ,1606.456459
+ ,1446.737592
+ ,1005.045182
+ ,913.67958
+ ,694.7422417
+ ,647.2193014
+ ,603.2233095
+ ,561.3616255
+ ,477.9286605
+ ,648.8516322
+ ,676.4550378
+ ,1016.256499
+ ,1139.927185
+ ,1269.321991
+ ,1546.025319
+ ,1753.304467
+ ,1958.22963
+ ,2094.065774
+ ,1977.440699
+ ,2027.739055
+ ,1842.633768
+ ,1852.639287
+ ,1706.985147
+ ,1471.71997
+ ,1159.935346
+ ,888.3265125
+ ,747.4889418
+ ,573.6363
+ ,696.6030203
+ ,736.0360216
+ ,625.9859074
+ ,668.0129992
+ ,878.1341576
+ ,925.7355157
+ ,1384.068603
+ ,1504.264376
+ ,1677.525248
+ ,1935.129054
+ ,1995.385689
+ ,2158.635275
+ ,2278.786793
+ ,2004.468115
+ ,2055.40499
+ ,1967.431468
+ ,1770.334016
+ ,1699.341376
+ ,1476.810831
+ ,970.7915825
+ ,997.474492
+ ,652.0843291
+ ,708.677725
+ ,656.5517255
+ ,770.8461221
+ ,901.3881145
+ ,832.0016083
+ ,1028.128263
+ ,1296.071512
+ ,1785.233557
+ ,1804.294853
+ ,2165.319782
+ ,2274.425758
+ ,2444.192672
+ ,2254.896853
+ ,2475.21364
+ ,2096.602738
+ ,2180.683121
+ ,2069.506046
+ ,1776.663609
+ ,1612.449679
+ ,1322.927331
+ ,914.8881288
+ ,749.3597793
+ ,985.6467572
+ ,880.7232311
+ ,819.4581121
+ ,760.5542696
+ ,985.0475722
+ ,1295.758114
+ ,1683.545809
+ ,1787.214405
+ ,1972.790239
+ ,2134.412385
+ ,2263.189757
+ ,2432.266646
+ ,2369.774565
+ ,2360.715385
+ ,2432.992026
+ ,2437.706015
+ ,1921.308588
+ ,1960.743986
+ ,1481.140265
+ ,1436.91711
+ ,1299.84893
+ ,985.736701
+ ,805.8439046
+ ,736.9963963
+ ,898.3310247
+ ,1065.843227
+ ,1113.559476
+ ,1211.918524
+ ,1462.155718
+ ,1801.38398
+ ,2045.731198
+ ,2405.387519
+ ,2376.281401
+ ,2542.194732
+ ,2535.172687
+ ,2518.76254
+ ,2467.712934
+ ,2483.547793
+ ,2111.5245
+ ,1987.469819
+ ,1900.284914
+ ,1545.174509
+ ,1241.063215
+ ,1244.17204
+ ,1001.510714
+ ,1004.787213
+ ,1064.602774
+ ,1250.752996
+ ,1244.505238
+ ,1565.513983
+ ,1904.45973
+ ,2200.393469
+ ,2489.302736
+ ,2382.022655
+ ,2740.656964
+ ,2692.201467
+ ,2707.985974
+ ,2793.734301
+ ,2909.368555
+ ,2611.616903
+ ,2367.927913
+ ,2036.392713
+ ,1794.090508
+ ,1658.074271
+ ,1450.206111
+ ,1159.117531
+ ,1144.297116
+ ,1188.107289
+ ,1082.864685
+ ,1152.30001
+ ,1370.470201
+ ,1560.414285
+ ,1878.82661
+ ,2213.171169
+ ,2306.794427
+ ,2878.274067
+ ,3043.563683
+ ,3163.511145
+ ,3218.690493
+ ,3177.601319
+ ,2941.518453
+ ,2922.823356
+ ,2325.88093
+ ,2414.009426
+ ,1874.014382
+ ,1866.949762
+ ,1330.453111
+ ,1413.735476
+ ,988.2221384
+ ,1046.533458
+ ,1398.652585
+ ,1444.958119
+ ,1637.387943
+ ,1926.364707
+ ,2088.421239
+ ,2250.892557
+ ,2803.393628
+ ,2920.006864
+ ,3028.392467
+ ,3284.085805
+ ,3116.539112
+ ,3120.244711
+ ,2946.135624
+ ,2740.67997
+ ,2810.140837
+ ,2233.823039
+ ,2114.011418
+ ,2076.542437
+ ,1571.481405
+ ,1567.958655
+ ,1217.444043
+ ,1136.643661
+ ,1481.287919
+ ,1528.477337
+ ,1797.720002
+ ,1819.469046
+ ,2395.358987
+ ,2707.314708
+ ,2921.060577
+ ,3246.433366
+ ,3051.477548
+ ,3444.481052
+ ,3432.580153
+ ,3421.04804
+ ,3381.830601
+ ,3169.679136
+ ,2872.11508
+ ,2601.76563
+ ,2143.011809
+ ,2159.598897
+ ,1645.12941
+ ,1622.54215
+ ,1350.064447
+ ,1372.486802
+ ,1475.078831
+ ,1356.08523
+ ,1974.230192
+ ,1835.493872
+ ,2495.225532
+ ,2813.660064
+ ,3199.803837
+ ,3128.428211
+ ,3614.588977
+ ,3741.816099
+ ,3711.62971
+ ,3501.116488
+ ,3620.888016
+ ,3348.935195)
> par4 = '12'
> par3 = '0'
> par2 = '0'
> par1 = '1'
> 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/fisher/rcomp/tmp/11hpu1386179792.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/fisher/rcomp/tmp/23cbh1386179792.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/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/fisher/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/fisher/rcomp/tmp/3175c1386179792.tab") 
> 
> try(system("convert tmp/11hpu1386179792.ps tmp/11hpu1386179792.png",intern=TRUE))
character(0)
> try(system("convert tmp/23cbh1386179792.ps tmp/23cbh1386179792.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.004   0.907   3.875