R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-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(1657
+ ,1418
+ ,1501
+ ,1315
+ ,1621
+ ,2308
+ ,3554
+ ,3318
+ ,3252
+ ,2921
+ ,2133
+ ,2040
+ ,1858
+ ,1833
+ ,2094
+ ,2173
+ ,2366
+ ,2074
+ ,2522
+ ,1822
+ ,1952
+ ,2232
+ ,1755
+ ,1791
+ ,2075
+ ,1850
+ ,2137
+ ,2467
+ ,2154
+ ,2289
+ ,2628
+ ,2074
+ ,2798
+ ,2194
+ ,2442
+ ,2565
+ ,2063
+ ,2070
+ ,2539
+ ,1898
+ ,2139
+ ,2408
+ ,2725
+ ,2201
+ ,2311
+ ,2548
+ ,2276
+ ,2351
+ ,2280
+ ,2057
+ ,2479
+ ,2379
+ ,2295
+ ,2456
+ ,2546
+ ,2844
+ ,2260
+ ,2981
+ ,2678
+ ,3440
+ ,2842
+ ,2450
+ ,2669
+ ,2570
+ ,2540
+ ,2318
+ ,2930
+ ,2947
+ ,2799
+ ,2695
+ ,2498
+ ,2260
+ ,2160
+ ,2058
+ ,2533
+ ,2150
+ ,2172
+ ,2155
+ ,3016
+ ,2333
+ ,2355
+ ,2825
+ ,2214
+ ,2360
+ ,2299
+ ,1746
+ ,2069
+ ,2267
+ ,1878
+ ,2266
+ ,2282
+ ,2085
+ ,2277
+ ,2251
+ ,1828
+ ,1954
+ ,1851
+ ,1570
+ ,1852
+ ,2187
+ ,1855
+ ,2218
+ ,2253
+ ,2028
+ ,2169
+ ,1997
+ ,2034
+ ,1791
+ ,1627
+ ,1631
+ ,2319
+ ,1707
+ ,1747
+ ,2397
+ ,2059
+ ,2251
+ ,2558
+ ,2406
+ ,2049
+ ,2074
+ ,1734
+ ,1983
+ ,2121
+ ,1905
+ ,2126
+ ,2363
+ ,2173
+ ,2710
+ ,2137
+ ,2742
+ ,2419
+ ,2194
+ ,2660
+ ,2189
+ ,2310
+ ,2349
+ ,2540
+ ,2434
+ ,2916
+ ,2446
+ ,2375
+ ,3032
+ ,2218
+ ,1920
+ ,2039
+ ,1889
+ ,2014
+ ,2105
+ ,2153
+ ,2309
+ ,2955
+ ,2225
+ ,2160
+ ,2386
+ ,1653
+ ,1099
+ ,5010
+ ,2672
+ ,2729
+ ,2955
+ ,2409
+ ,3086
+ ,3384
+ ,2458
+ ,2913
+ ,2448
+ ,2215
+ ,2179
+ ,2461
+ ,2098
+ ,2621
+ ,2703
+ ,2388
+ ,3880
+ ,3310
+ ,3093
+ ,3237
+ ,3002
+ ,2670
+ ,2311
+ ,2062
+ ,2059
+ ,2465
+ ,2213
+ ,2028
+ ,2322
+ ,2825
+ ,2687
+ ,2373
+ ,2889
+ ,2708
+ ,2542
+ ,2477
+ ,2419
+ ,2977
+ ,3001
+ ,3075
+ ,2870
+ ,3756
+ ,3443
+ ,2948
+ ,3560
+ ,3257
+ ,2600
+ ,2741
+ ,2349
+ ,2783
+ ,2845
+ ,2987
+ ,2696
+ ,3874
+ ,2912
+ ,2743
+ ,3857
+ ,2660
+ ,2226
+ ,2942
+ ,2420
+ ,2516
+ ,2421
+ ,2631
+ ,2887
+ ,3328
+ ,2587
+ ,2695
+ ,3669
+ ,2773
+ ,2527
+ ,2750
+ ,2014
+ ,2763
+ ,2726
+ ,1826
+ ,2713
+ ,3040
+ ,2405
+ ,2526
+ ,2526
+ ,2529
+ ,2474
+ ,2576
+ ,2219
+ ,2900
+ ,2274
+ ,2184
+ ,2629
+ ,2739
+ ,2933
+ ,3144
+ ,3354
+ ,3357
+ ,3329)
> 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/wessaorg/rcomp/tmp/1wtfu1323124149.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/wessaorg/rcomp/tmp/2dkvm1323124149.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/35lwd1323124149.tab") 
> 
> try(system("convert tmp/1wtfu1323124149.ps tmp/1wtfu1323124149.png",intern=TRUE))
character(0)
> try(system("convert tmp/2dkvm1323124149.ps tmp/2dkvm1323124149.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.805   0.123   0.935