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(13,96,60.1,281.25,548.59,470.325,456.25,456.25,483.5,498.25,457.5,490.7,303.075,510,417.7,631.4,595.625,531.15,581.725,581.725,544.75,582,622.225,638.666,469.9,532.78,538.5,884.65,559.08,1045.51,1319.34,0,1300,1322.6,962.57,1030.715,1038.31,967.87,1228.2,1300.285,1095.995,1181.955,1499.46,1523.38,1237.765,1119.545,1136.165,1081.63,1060.065,1216.67,1186.17,1217.475,1096.95,1685.6,1758.5,1786.6,2049.895,1845.895,2015.02,1609.63,918.725,1240.96,1671.785)
> par10 = '0.1'
> par9 = '1'
> par8 = 'dumresult'
> par7 = 'dum'
> par6 = '12'
> par5 = 'ZZZ'
> par4 = 'NA'
> par3 = 'NA'
> par2 = 'Croston'
> par1 = 'Input box'
> par10 <- '0.1'
> par9 <- '1'
> par8 <- 'dumresult'
> 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
63       1438.466 1137.521 1739.411 978.2107 1898.722
64       1438.466 1136.020 1740.912 975.9151 1901.018
65       1438.466 1134.527 1742.406 973.6309 1903.302
66       1438.466 1133.041 1743.892 971.3579 1905.575
67       1438.466 1131.561 1745.371 969.0958 1907.837
68       1438.466 1130.090 1746.843 966.8446 1910.088
69       1438.466 1128.625 1748.308 964.6042 1912.329
70       1438.466 1127.166 1749.766 962.3742 1914.559
71       1438.466 1125.715 1751.218 960.1546 1916.778
72       1438.466 1124.271 1752.662 957.9453 1918.987

$period
   Point Forecast     Lo 80    Hi 80     Lo 95    Hi 95
63       1.004239 0.8345035 1.173975 0.7446508 1.263827
64       1.004239 0.8336569 1.174821 0.7433561 1.265122
65       1.004239 0.8328145 1.175664 0.7420678 1.266410
66       1.004239 0.8319762 1.176502 0.7407858 1.267692
67       1.004239 0.8311420 1.177336 0.7395100 1.268968
68       1.004239 0.8303118 1.178166 0.7382403 1.270238
69       1.004239 0.8294856 1.178993 0.7369766 1.271502
70       1.004239 0.8286632 1.179815 0.7357189 1.272759
71       1.004239 0.8278447 1.180634 0.7344671 1.274011
72       1.004239 0.8270299 1.181448 0.7332210 1.275257


In-sample error measures:
        ME       RMSE        MAE        MPE       MAPE       MASE 
233.018350 349.012743 279.974981       -Inf        Inf   1.497651 

Forecasts:
   Point Forecast
64       1432.394
65       1432.394
66       1432.394
67       1432.394
68       1432.394
69       1432.394
70       1432.394
71       1432.394
72       1432.394
73       1432.394
> postscript(file="/var/www/html/rcomp/tmp/1cwjp1271584395.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/2cwjp1271584395.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/38ozg1271584395.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/41xyj1271584395.tab") 
> try(system("convert tmp/1cwjp1271584395.ps tmp/1cwjp1271584395.png",intern=TRUE))
character(0)
> try(system("convert tmp/2cwjp1271584395.ps tmp/2cwjp1271584395.png",intern=TRUE))
character(0)
> 
> #-SERVER-wessa.org
> 
> 
> 
> proc.time()
   user  system elapsed 
  3.035   0.344   4.162