R version 2.10.1 (2009-12-14)
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 <- array(list(275.0916667
+ ,NA
+ ,275.182487301474
+ ,274.816575230748
+ ,341.6167
+ ,401.3423333
+ ,275.0916667
+ ,400.972242874856
+ ,295.968930016752
+ ,379.0167
+ ,260.6333333
+ ,287.71673336
+ ,260.845832809383
+ ,338.328662343278
+ ,346.35
+ ,286.1
+ ,285.008393354
+ ,286.165919817652
+ ,308.265916154109
+ ,302.25
+ ,368.0833333
+ ,285.1175540186
+ ,367.853139451386
+ ,301.071991695963
+ ,239.7667
+ ,385.1916667
+ ,293.41413194674
+ ,385.000486465639
+ ,324.590521416548
+ ,386.9333
+ ,181.5333333
+ ,302.591885422066
+ ,182.007237238174
+ ,345.060642650557
+ ,229.25
+ ,145.1
+ ,290.486030209860
+ ,145.595640852172
+ ,289.120656303052
+ ,182.2667
+ ,203.1833333
+ ,275.947427188874
+ ,203.328434462353
+ ,240.050659948121
+ ,272.3333
+ ,227.5833333
+ ,268.671017799986
+ ,227.574414189821
+ ,227.504420390666
+ ,234.4
+ ,239.0833333
+ ,264.562249349988
+ ,239.026738533428
+ ,227.532064816018
+ ,315.8
+ ,109.175
+ ,262.014357744989
+ ,109.701909322141
+ ,231.457107359672
+ ,118.75
+ ,231.0833333
+ ,246.73042197049
+ ,230.854906936773
+ ,189.871114035390
+ ,200.6667
+ ,216.9166667
+ ,245.165713103441
+ ,216.849097935047
+ ,203.886458556338
+ ,147.3333
+ ,229.8583333
+ ,242.340808463097
+ ,229.753265490343
+ ,208.317704231102
+ ,304.6167
+ ,272.7916667
+ ,241.092560946787
+ ,272.526813215067
+ ,215.643057115841
+ ,187.5167
+ ,155.0833333
+ ,244.262471522108
+ ,155.425389569831
+ ,235.07741198107
+ ,114.5167)
+ ,dim=c(5
+ ,17)
+ ,dimnames=list(c('Actuals'
+ ,'Croston'
+ ,'ETS'
+ ,'Arima'
+ ,'Kaneka')
+ ,1:17))
>  y <- array(NA,dim=c(5,17),dimnames=list(c('Actuals','Croston','ETS','Arima','Kaneka'),1:17))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
> par2 = 'No'
> par1 = '1'
> #'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: 
> library(rpart)
Warning message:
package 'rpart' was built under R version 2.8.1 and help may not work correctly 
> library(partykit)
Loading required package: grid
Loading required package: mvtnorm
> par1 <- as.numeric(par1)
> autoprune <- function ( tree, method='Minimum CV'){
+ xerr <- tree$cptable[,'xerror']
+ cpmin.id <- which.min(xerr)
+ if (method == 'Minimum CV Error plus 1 SD'){
+ xstd <- tree$cptable[,'xstd']
+ errt <- xerr[cpmin.id] + xstd[cpmin.id]
+ cpSE1.min <- which.min( errt < xerr )
+ mycp <- (tree$cptable[,'CP'])[cpSE1.min]
+ }
+ if (method == 'Minimum CV') {
+ mycp <- (tree$cptable[,'CP'])[cpmin.id]
+ }
+ return (mycp)
+ }
> conf.multi.mat <- function(true, new)
+ {
+ if ( all( is.na(match( levels(true),levels(new) ) )) )
+ stop ( 'conflict of vector levels')
+ multi.t <- list()
+ for (mylev in levels(true) ) {
+ true.tmp <- true
+ new.tmp <- new
+ left.lev <- levels (true.tmp)[- match(mylev,levels(true) ) ]
+ levels(true.tmp) <- list ( mylev = mylev, all = left.lev )
+ levels(new.tmp)  <- list ( mylev = mylev, all = left.lev )
+ curr.t <- conf.mat ( true.tmp , new.tmp )
+ multi.t[[mylev]] <- curr.t
+ multi.t[[mylev]]$precision <-
+ round( curr.t$conf[1,1] / sum( curr.t$conf[1,] ), 2 )
+ }
+ return (multi.t)
+ }
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> m <- rpart(as.data.frame(x1))
> par2
[1] "No"
> if (par2 != 'No') {
+ mincp <- autoprune(m,method=par2)
+ print(mincp)
+ m <- prune(m,cp=mincp)
+ }
> m$cptable
    CP nsplit rel error
1 0.01      0         1
> postscript(file="/var/www/rcomp/tmp/170081274874923.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(as.party(m),tp_args=list(id=FALSE))
Error: length(i) == 1 & is.numeric(i) is not TRUE
Execution halted