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(0,0.5,0,0.5,0,0,0,4.5,4.5,10,3.75,58.45,7.45,22.275,22.21,5.8,36.43,31.65,50,43,0,68.5,33.5,23,0.5,69.25,32,39.213,46.426,46.855,153.135,64,31,2.25,2.25,2.3,22.6,1.5,10.65,34,81.75,106.5,0.525,24.025,5.25,9,12.8,25.05,0.3,75.75,54.75,1.526,102,3.752,17.25,9.2,50.25,2.25,3.95,60,55.8,6.75,61.95,7.025,85.75,18.525,6,25.35,46.775,51.025,30,3,30,44,80.75,27.5,39.725,29.25,32.725,56.25,28.65,51.75,32.26,72,65.4,33.75,77.85,10.875)
> par10 = '0.1'
> par9 = '3'
> par8 = 'B611crostonm'
> par7 = 'dum'
> par6 = '12'
> par5 = 'ZZZ'
> par4 = 'NA'
> par3 = 'NA'
> par2 = 'Croston'
> par1 = 'Input box'
> par10 <- '0.1'
> par9 <- '3'
> par8 <- 'B611crostonm'
> par7 <- 'dum'
> par6 <- '12'
> par5 <- 'ZZZ'
> par4 <- 'NA'
> par3 <- 'NA'
> par2 <- 'Croston'
> par1 <- 'Input box'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Dr. Ian E. Holliday
> #To cite this work: Ian E. Holliday, 2009, 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: 
> #Technical description: 
> if(par3!='NA') par3 <- as.numeric(par3) else par3 <- NA
> if(par4!='NA') par4 <- as.numeric(par4) else par4 <- NA
> par6 <- as.numeric(par6) #Seasonal Period
> par9 <- as.numeric(par9) #Forecast Horizon
> par10 <- as.numeric(par10) #Alpha
> library(forecast)
Loading required package: tseries
Loading required package: quadprog
Loading required package: zoo

Attaching package: 'zoo'


	The following object(s) are masked from package:base :

	 as.Date.numeric 

This is forecast 2.03 
> if (par1 == 'CSV') {
+ xarr <- read.csv(file=paste('tmp/',par7,'.csv',sep=''),header=T)
+ numseries <- length(xarr[1,])-1
+ n <- length(xarr[,1])
+ nmh <- n - par9
+ nmhp1 <- nmh + 1
+ rarr <- array(NA,dim=c(n,numseries))
+ farr <- array(NA,dim=c(n,numseries))
+ parr <- array(NA,dim=c(numseries,8))
+ colnames(parr) = list('ME','RMSE','MAE','MPE','MAPE','MASE','ACF1','TheilU')
+ for(i in 1:numseries) {
+ sindex <- i+1
+ x <- xarr[,sindex]
+ if(par2=='Croston') {
+ if (i==1) m <- croston(x,alpha=par10)
+ if (i==1) mydemand <- m$model$demand[]
+ fit <- croston(x[1:nmh],h=par9,alpha=par10)
+ }
+ if(par2=='ARIMA') {
+ m <- auto.arima(ts(x,freq=par6),d=par3,D=par4)
+ mydemand <- forecast(m)
+ fit <- auto.arima(ts(x[1:nmh],freq=par6),d=par3,D=par4)
+ }
+ if(par2=='ETS') {
+ m <- ets(ts(x,freq=par6),model=par5)
+ mydemand <- forecast(m)
+ fit <- ets(ts(x[1:nmh],freq=par6),model=par5)
+ }
+ try(rarr[,i] <- mydemand$resid,silent=T)
+ try(farr[,i] <- mydemand$mean,silent=T)
+ if (par2!='Croston') parr[i,] <- accuracy(forecast(fit,par9),x[nmhp1:n])
+ if (par2=='Croston') parr[i,] <- accuracy(fit,x[nmhp1:n])
+ }
+ write.csv(farr,file=paste('tmp/',par8,'_f.csv',sep=''))
+ write.csv(rarr,file=paste('tmp/',par8,'_r.csv',sep=''))
+ write.csv(parr,file=paste('tmp/',par8,'_p.csv',sep=''))
+ }
> if (par1 == 'Input box') {
+ numseries <- 1
+ n <- length(x)
+ if(par2=='Croston') {
+ m <- croston(x)
+ mydemand <- m$model$demand[]
+ }
+ if(par2=='ARIMA') {
+ m <- auto.arima(ts(x,freq=par6),d=par3,D=par4)
+ mydemand <- forecast(m)
+ }
+ if(par2=='ETS') {
+ m <- ets(ts(x,freq=par6),model=par5)
+ mydemand <- forecast(m)
+ }
+ summary(m)
+ }

Forecast method: Croston's method

Model Information:
$demand
   Point Forecast    Lo 80    Hi 80     Lo 95    Hi 95
83       42.44498 3.258099 81.63186 -17.48618 102.3761
84       42.44498 3.062652 81.82731 -17.78509 102.6751
85       42.44498 2.868170 82.02179 -18.08253 102.9725
86       42.44498 2.674640 82.21532 -18.37851 103.2685
87       42.44498 2.482046 82.40792 -18.67305 103.5630
88       42.44498 2.290376 82.59959 -18.96619 103.8562
89       42.44498 2.099617 82.79035 -19.25793 104.1479
90       42.44498 1.909755 82.98021 -19.54830 104.4383
91       42.44498 1.720779 83.16918 -19.83731 104.7273
92       42.44498 1.532675 83.35729 -20.12499 105.0150

$period
   Point Forecast     Lo 80    Hi 80     Lo 95    Hi 95
83       1.000387 0.5242854 1.476488 0.2722525 1.728521
84       1.000387 0.5219108 1.478863 0.2686209 1.732153
85       1.000387 0.5195479 1.481226 0.2650072 1.735766
86       1.000387 0.5171966 1.483577 0.2614112 1.739362
87       1.000387 0.5148567 1.485917 0.2578326 1.742941
88       1.000387 0.5125280 1.488246 0.2542712 1.746502
89       1.000387 0.5102104 1.490563 0.2507267 1.750047
90       1.000387 0.5079037 1.492870 0.2471989 1.753575
91       1.000387 0.5056077 1.495166 0.2436875 1.757086
92       1.000387 0.5033223 1.497451 0.2401923 1.760581


In-sample error measures:
       ME      RMSE       MAE       MPE      MAPE      MASE 
 6.374529 30.555575 23.128768      -Inf       Inf  0.809942 

Forecasts:
   Point Forecast
89       42.42857
90       42.42857
91       42.42857
92       42.42857
93       42.42857
94       42.42857
95       42.42857
96       42.42857
97       42.42857
98       42.42857
> postscript(file="/var/www/html/rcomp/tmp/1u2ly1272290205.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow=c(2,1))
> if (par2=='Croston') plot(m)
> if ((par2=='ARIMA') | par2=='ETS') plot(forecast(m))
> plot(mydemand$resid,type='l',main='Residuals', ylab='residual value', xlab='time')
> par(op)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/24bkj1272290205.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow=c(2,2))
> acf(mydemand$resid, lag.max=n/3, main='Residual ACF', ylab='autocorrelation', xlab='time lag')
> pacf(mydemand$resid,lag.max=n/3, main='Residual PACF', ylab='partial autocorrelation', xlab='time lag')
> cpgram(mydemand$resid, main='Cumulative Periodogram of Residuals')
> qqnorm(mydemand$resid); qqline(mydemand$resid, col=2)
> par(op)
> 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,'Demand Forecast',6,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Point',header=TRUE)
> a<-table.element(a,'Forecast',header=TRUE)
> a<-table.element(a,'95% LB',header=TRUE)
> a<-table.element(a,'80% LB',header=TRUE)
> a<-table.element(a,'80% UB',header=TRUE)
> a<-table.element(a,'95% UB',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:length(mydemand$mean)) {
+ a<-table.row.start(a)
+ a<-table.element(a,i+n,header=TRUE)
+ a<-table.element(a,as.numeric(mydemand$mean[i]))
+ a<-table.element(a,as.numeric(mydemand$lower[i,2]))
+ a<-table.element(a,as.numeric(mydemand$lower[i,1]))
+ a<-table.element(a,as.numeric(mydemand$upper[i,1]))
+ a<-table.element(a,as.numeric(mydemand$upper[i,2]))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/3qc1p1272290205.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Actuals and Interpolation',3,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Time',header=TRUE)
> a<-table.element(a,'Actual',header=TRUE)
> a<-table.element(a,'Forecast',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i,header=TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i] - as.numeric(m$resid[i]))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/4tchv1272290205.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'What is next?',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/Patrick.Wessa/rwasp_demand_forecasting_simulate.wasp',sep=''),'Simulate Time Series','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/Patrick.Wessa/rwasp_demand_forecasting_croston.wasp',sep=''),'Generate Forecasts','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/Patrick.Wessa/rwasp_demand_forecasting_analysis.wasp',sep=''),'Forecast Analysis','',target=''))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/543gf1272290205.tab") 
> try(system("convert tmp/1u2ly1272290205.ps tmp/1u2ly1272290205.png",intern=TRUE))
character(0)
> try(system("convert tmp/24bkj1272290205.ps tmp/24bkj1272290205.png",intern=TRUE))
character(0)
> 
> #-SERVER-wessa.org
> 
> 
> 
> proc.time()
   user  system elapsed 
  3.940   0.390   4.159