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(1045620,1050293,1030185,1015746,999514,990204,962206,947314,958768)
> par1 = '12'
> par1 <- '12'
> #'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)
> nx <- length(x)
> x <- ts(x,frequency=par1)
> m <- StructTS(x,type='BSM')
> m$coef
   level    slope     seas  epsilon 
       0 50611106        0 53630569 
> m$fitted
          level      slope       sea
Jan 1 1045620.0      0.000   0.00000
Feb 1 1048002.8   2167.292 107.77238
Mar 1 1035294.1  -7637.468  43.34047
Apr 1 1018528.2 -13200.151  81.87061
May 1 1000821.7 -15960.564 101.06127
Jun 1  988869.2 -13489.170  92.01850
Jul 1  965241.1 -19751.767  99.27754
Aug 1  946807.3 -18937.301  99.25933
Sep 1  951392.6  -4399.088 102.99029
> m$resid
         Jan        Feb        Mar        Apr        May        Jun        Jul
1  0.0000000  0.4360717 -1.3893756 -0.7989860 -0.3935342  0.3485216 -0.8805258
         Aug        Sep
1  0.1144863  2.0436141
> mylevel <- as.numeric(m$fitted[,'level'])
> myslope <- as.numeric(m$fitted[,'slope'])
> myseas <- as.numeric(m$fitted[,'sea'])
> myresid <- as.numeric(m$resid)
> myfit <- mylevel+myseas
> mylagmax <- nx/2
> postscript(file="/var/fisher/rcomp/tmp/1ct1y1384431570.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(mylevel,na.action=na.pass,lag.max = mylagmax,main='Level')
> acf(myseas,na.action=na.pass,lag.max = mylagmax,main='Seasonal')
> acf(myresid,na.action=na.pass,lag.max = mylagmax,main='Standardized Residals')
> par(op)
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/2u21s1384431570.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(mylevel,main='Level')
> spectrum(myseas,main='Seasonal')
> spectrum(myresid,main='Standardized Residals')
> par(op)
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/3yecm1384431570.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(mylevel,main='Level')
> cpgram(myseas,main='Seasonal')
> cpgram(myresid,main='Standardized Residals')
> par(op)
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/40b5i1384431570.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(as.numeric(m$resid),main='Standardized Residuals',ylab='Residuals',xlab='time',type='b')
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/5qkpv1384431570.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow = c(2,2))
> hist(m$resid,main='Residual Histogram')
> plot(density(m$resid),main='Residual Kernel Density')
> qqnorm(m$resid,main='Residual Normal QQ Plot')
> qqline(m$resid)
> plot(m$resid^2, myfit^2,main='Sq.Resid vs. Sq.Fit',xlab='Squared residuals',ylab='Squared Fit')
> par(op)
> 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,'Structural Time Series Model',6,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'t',header=TRUE)
> a<-table.element(a,'Observed',header=TRUE)
> a<-table.element(a,'Level',header=TRUE)
> a<-table.element(a,'Slope',header=TRUE)
> a<-table.element(a,'Seasonal',header=TRUE)
> a<-table.element(a,'Stand. Residuals',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:nx) {
+ a<-table.row.start(a)
+ a<-table.element(a,i,header=TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,mylevel[i])
+ a<-table.element(a,myslope[i])
+ a<-table.element(a,myseas[i])
+ a<-table.element(a,myresid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/fisher/rcomp/tmp/6lemx1384431570.tab") 
> 
> try(system("convert tmp/1ct1y1384431570.ps tmp/1ct1y1384431570.png",intern=TRUE))
character(0)
> try(system("convert tmp/2u21s1384431570.ps tmp/2u21s1384431570.png",intern=TRUE))
character(0)
> try(system("convert tmp/3yecm1384431570.ps tmp/3yecm1384431570.png",intern=TRUE))
character(0)
> try(system("convert tmp/40b5i1384431570.ps tmp/40b5i1384431570.png",intern=TRUE))
character(0)
> try(system("convert tmp/5qkpv1384431570.ps tmp/5qkpv1384431570.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.552   0.565   3.099