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(221.102
+ ,220.892
+ ,225.537
+ ,219.334
+ ,216.126
+ ,224.573
+ ,213.991
+ ,221.865
+ ,215.382
+ ,213.962
+ ,217.009
+ ,318.110
+ ,219.662
+ ,212.650
+ ,209.307
+ ,210.541
+ ,210.609
+ ,208.910
+ ,207.541
+ ,207.699
+ ,205.005
+ ,205.747
+ ,203.644
+ ,229.937
+ ,214.446
+ ,210.194
+ ,206.535
+ ,216.524
+ ,198.243
+ ,208.274
+ ,207.493
+ ,215.525
+ ,207.562
+ ,213.355
+ ,209.048
+ ,220.497
+ ,214.563
+ ,211.571
+ ,216.385
+ ,211.496
+ ,209.683
+ ,206.304
+ ,213.925
+ ,204.829
+ ,205.729
+ ,200.296
+ ,207.960
+ ,207.729
+ ,208.327
+ ,207.794
+ ,213.700
+ ,213.136
+ ,210.354
+ ,202.205
+ ,220.425
+ ,210.038
+ ,205.201
+ ,197.441
+ ,200.021
+ ,203.165
+ ,201.699
+ ,201.003
+ ,232.308
+ ,211.412
+ ,209.661
+ ,213.587
+ ,193.611
+ ,192.543
+ ,196.040
+ ,191.407
+ ,189.726
+ ,191.272
+ ,187.924
+ ,198.213
+ ,199.352
+ ,199.289
+ ,195.475
+ ,198.045
+ ,197.615
+ ,189.015
+ ,189.668
+ ,189.120
+ ,194.168
+ ,192.304
+ ,185.913
+ ,197.599
+ ,186.085
+ ,190.566
+ ,187.054
+ ,193.222
+ ,189.856
+ ,190.608
+ ,190.588
+ ,186.773
+ ,183.510
+ ,180.106
+ ,180.150
+ ,178.412
+ ,182.353
+ ,203.805
+ ,186.054
+ ,184.290
+ ,187.483
+ ,187.111
+ ,189.561
+ ,184.439
+ ,182.985
+ ,183.828
+ ,184.036
+ ,183.214
+ ,183.464
+ ,173.718
+ ,180.210
+ ,171.252
+ ,172.705
+ ,174.006
+ ,172.043
+ ,169.445
+ ,169.449
+ ,177.073
+ ,170.799
+ ,171.648
+ ,172.220
+ ,165.795
+ ,167.466
+ ,165.528
+ ,162.851
+ ,165.864
+ ,162.094
+ ,162.385
+ ,164.293
+ ,165.983
+ ,159.680
+ ,161.739
+ ,159.302
+ ,167.795
+ ,164.242
+ ,159.743
+ ,160.887
+ ,163.844
+ ,161.172
+ ,159.330
+ ,155.570
+ ,156.749
+ ,155.012
+ ,163.419
+ ,153.630
+ ,154.535
+ ,151.543
+ ,152.955
+ ,150.166
+ ,151.416
+ ,150.332
+ ,152.196
+ ,153.422
+ ,147.435)
> par2 = '12'
> par1 = 'additive'
> par2 <- '12'
> par1 <- 'additive'
> #'GNU S' R Code compiled by R2WASP v. 1.2.327 ()
> #Author: root
> #To cite this work: Wessa P., (2013), Classical Decomposition (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_decompose.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> par2 <- as.numeric(par2)
> x <- ts(x,freq=par2)
> m <- decompose(x,type=par1)
> m$figure
 [1] -1.80767998 -0.07103414  1.83153183  2.60278183 -0.96346470 -0.67330498
 [7] -0.70541262 -0.72937442 -2.18863831 -3.90938137 -3.15322164  9.76719850
> postscript(file="/var/wessaorg/rcomp/tmp/1qivf1386347165.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(m)
> dev.off()
null device 
          1 
> mylagmax <- length(x)/2
> postscript(file="/var/wessaorg/rcomp/tmp/209ob1386347165.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow = c(2,2))
> acf(as.numeric(x),lag.max = mylagmax,main='Observed')
> acf(as.numeric(m$trend),na.action=na.pass,lag.max = mylagmax,main='Trend')
> acf(as.numeric(m$seasonal),na.action=na.pass,lag.max = mylagmax,main='Seasonal')
> acf(as.numeric(m$random),na.action=na.pass,lag.max = mylagmax,main='Random')
> par(op)
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/3apv81386347165.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow = c(2,2))
> spectrum(as.numeric(x),main='Observed')
> spectrum(as.numeric(m$trend[!is.na(m$trend)]),main='Trend')
> spectrum(as.numeric(m$seasonal[!is.na(m$seasonal)]),main='Seasonal')
> spectrum(as.numeric(m$random[!is.na(m$random)]),main='Random')
> par(op)
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/47pmd1386347165.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow = c(2,2))
> cpgram(as.numeric(x),main='Observed')
> cpgram(as.numeric(m$trend[!is.na(m$trend)]),main='Trend')
> cpgram(as.numeric(m$seasonal[!is.na(m$seasonal)]),main='Seasonal')
> cpgram(as.numeric(m$random[!is.na(m$random)]),main='Random')
> par(op)
> 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,'Classical Decomposition by Moving Averages',6,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'t',header=TRUE)
> a<-table.element(a,'Observations',header=TRUE)
> a<-table.element(a,'Fit',header=TRUE)
> a<-table.element(a,'Trend',header=TRUE)
> a<-table.element(a,'Seasonal',header=TRUE)
> a<-table.element(a,'Random',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:length(m$trend)) {
+ a<-table.row.start(a)
+ a<-table.element(a,i,header=TRUE)
+ a<-table.element(a,x[i])
+ if (par1 == 'additive') a<-table.element(a,signif(m$trend[i]+m$seasonal[i],6)) else a<-table.element(a,signif(m$trend[i]*m$seasonal[i],6))
+ a<-table.element(a,signif(m$trend[i],6))
+ a<-table.element(a,signif(m$seasonal[i],6))
+ a<-table.element(a,signif(m$random[i],6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/5ozgf1386347165.tab") 
> 
> try(system("convert tmp/1qivf1386347165.ps tmp/1qivf1386347165.png",intern=TRUE))
character(0)
> try(system("convert tmp/209ob1386347165.ps tmp/209ob1386347165.png",intern=TRUE))
character(0)
> try(system("convert tmp/3apv81386347165.ps tmp/3apv81386347165.png",intern=TRUE))
character(0)
> try(system("convert tmp/47pmd1386347165.ps tmp/47pmd1386347165.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  5.117   0.810   5.896