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 <- array(list(24 + ,14 + ,11 + ,12 + ,24 + ,26 + ,25 + ,11 + ,7 + ,8 + ,25 + ,23 + ,17 + ,6 + ,17 + ,8 + ,30 + ,25 + ,18 + ,12 + ,10 + ,8 + ,19 + ,23 + ,18 + ,8 + ,12 + ,9 + ,22 + ,19 + ,16 + ,10 + ,12 + ,7 + ,22 + ,29 + ,20 + ,10 + ,11 + ,4 + ,25 + ,25 + ,16 + ,11 + ,11 + ,11 + ,23 + ,21 + ,18 + ,16 + ,12 + ,7 + ,17 + ,22 + ,17 + ,11 + ,13 + ,7 + ,21 + ,25 + ,23 + ,13 + ,14 + ,12 + ,19 + ,24 + ,30 + ,12 + ,16 + ,10 + ,19 + ,18 + ,23 + ,8 + ,11 + ,10 + ,15 + ,22 + ,18 + ,12 + ,10 + ,8 + ,16 + ,15 + ,15 + ,11 + ,11 + ,8 + ,23 + ,22 + ,12 + ,4 + ,15 + ,4 + ,27 + ,28 + ,21 + ,9 + ,9 + ,9 + ,22 + ,20 + ,15 + ,8 + ,11 + ,8 + ,14 + ,12 + ,20 + ,8 + ,17 + ,7 + ,22 + ,24 + ,31 + ,14 + ,17 + ,11 + ,23 + ,20 + ,27 + ,15 + ,11 + ,9 + ,23 + ,21 + ,34 + ,16 + ,18 + ,11 + ,21 + ,20 + ,21 + ,9 + ,14 + ,13 + ,19 + ,21 + ,31 + ,14 + ,10 + ,8 + ,18 + ,23 + ,19 + ,11 + ,11 + ,8 + ,20 + ,28 + ,16 + ,8 + ,15 + ,9 + ,23 + ,24 + ,20 + ,9 + ,15 + ,6 + ,25 + ,24 + ,21 + ,9 + ,13 + ,9 + ,19 + ,24 + ,22 + ,9 + ,16 + ,9 + ,24 + ,23 + ,17 + ,9 + ,13 + ,6 + ,22 + ,23 + ,24 + ,10 + ,9 + ,6 + ,25 + ,29 + ,25 + ,16 + ,18 + ,16 + ,26 + ,24 + ,26 + ,11 + ,18 + ,5 + ,29 + ,18 + ,25 + ,8 + ,12 + ,7 + ,32 + ,25 + ,17 + ,9 + ,17 + ,9 + ,25 + ,21 + ,32 + ,16 + ,9 + ,6 + ,29 + ,26 + ,33 + ,11 + ,9 + ,6 + ,28 + ,22 + ,13 + ,16 + ,12 + ,5 + ,17 + ,22 + ,32 + ,12 + ,18 + ,12 + ,28 + ,22 + ,25 + ,12 + ,12 + ,7 + ,29 + ,23 + ,29 + ,14 + ,18 + ,10 + ,26 + ,30 + ,22 + ,9 + ,14 + ,9 + ,25 + ,23 + ,18 + ,10 + ,15 + ,8 + ,14 + ,17 + ,17 + ,9 + ,16 + ,5 + ,25 + ,23 + ,20 + ,10 + ,10 + ,8 + ,26 + ,23 + ,15 + ,12 + ,11 + ,8 + ,20 + ,25 + ,20 + ,14 + ,14 + ,10 + ,18 + ,24 + ,33 + ,14 + ,9 + ,6 + ,32 + ,24 + ,29 + ,10 + ,12 + ,8 + ,25 + ,23 + ,23 + ,14 + ,17 + ,7 + ,25 + ,21 + ,26 + ,16 + ,5 + ,4 + ,23 + ,24 + ,18 + ,9 + ,12 + ,8 + ,21 + ,24 + ,20 + ,10 + ,12 + ,8 + ,20 + ,28 + ,11 + ,6 + ,6 + ,4 + ,15 + ,16 + ,28 + ,8 + ,24 + ,20 + ,30 + ,20 + ,26 + ,13 + ,12 + ,8 + ,24 + ,29 + ,22 + ,10 + ,12 + ,8 + ,26 + ,27 + ,17 + ,8 + ,14 + ,6 + ,24 + ,22 + ,12 + ,7 + ,7 + ,4 + ,22 + ,28 + ,14 + ,15 + ,13 + ,8 + ,14 + ,16 + ,17 + ,9 + ,12 + ,9 + ,24 + ,25 + ,21 + ,10 + ,13 + ,6 + ,24 + ,24 + ,19 + ,12 + ,14 + ,7 + ,24 + ,28 + ,18 + ,13 + ,8 + ,9 + ,24 + ,24 + ,10 + ,10 + ,11 + ,5 + ,19 + ,23 + ,29 + ,11 + ,9 + ,5 + ,31 + ,30 + ,31 + ,8 + ,11 + ,8 + ,22 + ,24 + ,19 + ,9 + ,13 + ,8 + ,27 + ,21 + ,9 + ,13 + ,10 + ,6 + ,19 + ,25 + ,20 + ,11 + ,11 + ,8 + ,25 + ,25 + ,28 + ,8 + ,12 + ,7 + ,20 + ,22 + ,19 + ,9 + ,9 + ,7 + ,21 + ,23 + ,30 + ,9 + ,15 + ,9 + ,27 + ,26 + ,29 + ,15 + ,18 + ,11 + ,23 + ,23 + ,26 + ,9 + ,15 + ,6 + ,25 + ,25 + ,23 + ,10 + ,12 + ,8 + ,20 + ,21 + ,13 + ,14 + ,13 + ,6 + ,21 + ,25 + ,21 + ,12 + ,14 + ,9 + ,22 + ,24 + ,19 + ,12 + ,10 + ,8 + ,23 + ,29 + ,28 + ,11 + ,13 + ,6 + ,25 + ,22 + ,23 + ,14 + ,13 + ,10 + ,25 + ,27 + ,18 + ,6 + ,11 + ,8 + ,17 + ,26 + ,21 + ,12 + ,13 + ,8 + ,19 + ,22 + ,20 + ,8 + ,16 + ,10 + ,25 + ,24 + ,23 + ,14 + ,8 + ,5 + ,19 + ,27 + ,21 + ,11 + ,16 + ,7 + ,20 + ,24 + ,21 + ,10 + ,11 + ,5 + ,26 + ,24 + ,15 + ,14 + ,9 + ,8 + ,23 + ,29 + ,28 + ,12 + ,16 + ,14 + ,27 + ,22 + ,19 + ,10 + ,12 + ,7 + ,17 + ,21 + ,26 + ,14 + ,14 + ,8 + ,17 + ,24 + ,10 + ,5 + ,8 + ,6 + ,19 + ,24 + ,16 + ,11 + ,9 + ,5 + ,17 + ,23 + ,22 + ,10 + ,15 + ,6 + ,22 + ,20 + ,19 + ,9 + ,11 + ,10 + ,21 + ,27 + ,31 + ,10 + ,21 + ,12 + ,32 + ,26 + ,31 + ,16 + ,14 + ,9 + ,21 + ,25 + ,29 + ,13 + ,18 + ,12 + ,21 + ,21 + ,19 + ,9 + ,12 + ,7 + ,18 + ,21 + ,22 + ,10 + ,13 + ,8 + ,18 + ,19 + ,23 + ,10 + ,15 + ,10 + ,23 + ,21 + ,15 + ,7 + ,12 + ,6 + ,19 + ,21 + ,20 + ,9 + ,19 + ,10 + ,20 + ,16 + ,18 + ,8 + ,15 + ,10 + ,21 + ,22 + ,23 + ,14 + ,11 + ,10 + ,20 + ,29 + ,25 + ,14 + ,11 + ,5 + ,17 + ,15 + ,21 + ,8 + ,10 + ,7 + ,18 + ,17 + ,24 + ,9 + ,13 + ,10 + ,19 + ,15 + ,25 + ,14 + ,15 + ,11 + ,22 + ,21 + ,17 + ,14 + ,12 + ,6 + ,15 + ,21 + ,13 + ,8 + ,12 + ,7 + ,14 + ,19 + ,28 + ,8 + ,16 + ,12 + ,18 + ,24 + ,21 + ,8 + ,9 + ,11 + ,24 + ,20 + ,25 + ,7 + ,18 + ,11 + ,35 + ,17 + ,9 + ,6 + ,8 + ,11 + ,29 + ,23 + ,16 + ,8 + ,13 + ,5 + ,21 + ,24 + ,19 + ,6 + ,17 + ,8 + ,25 + ,14 + ,17 + ,11 + ,9 + ,6 + ,20 + ,19 + ,25 + ,14 + ,15 + ,9 + ,22 + ,24 + ,20 + ,11 + ,8 + ,4 + ,13 + ,13 + ,29 + ,11 + ,7 + ,4 + ,26 + ,22 + ,14 + ,11 + ,12 + ,7 + ,17 + ,16 + ,22 + ,14 + ,14 + ,11 + ,25 + ,19 + ,15 + ,8 + ,6 + ,6 + ,20 + ,25 + ,19 + ,20 + ,8 + ,7 + ,19 + ,25 + ,20 + ,11 + ,17 + ,8 + ,21 + ,23 + ,15 + ,8 + ,10 + ,4 + ,22 + ,24 + ,20 + ,11 + ,11 + ,8 + ,24 + ,26 + ,18 + ,10 + ,14 + ,9 + ,21 + ,26 + ,33 + ,14 + ,11 + ,8 + ,26 + ,25 + ,22 + ,11 + ,13 + ,11 + ,24 + ,18 + ,16 + ,9 + ,12 + ,8 + ,16 + ,21 + ,17 + ,9 + ,11 + ,5 + ,23 + ,26 + ,16 + ,8 + ,9 + ,4 + ,18 + ,23 + ,21 + ,10 + ,12 + ,8 + ,16 + ,23 + ,26 + ,13 + ,20 + ,10 + ,26 + ,22 + ,18 + ,13 + ,12 + ,6 + ,19 + ,20 + ,18 + ,12 + ,13 + ,9 + ,21 + ,13 + ,17 + ,8 + ,12 + ,9 + ,21 + ,24 + ,22 + ,13 + ,12 + ,13 + ,22 + ,15 + ,30 + ,14 + ,9 + ,9 + ,23 + ,14 + ,30 + ,12 + ,15 + ,10 + ,29 + ,22 + ,24 + ,14 + ,24 + ,20 + ,21 + ,10 + ,21 + ,15 + ,7 + ,5 + ,21 + ,24 + ,21 + ,13 + ,17 + ,11 + ,23 + ,22 + ,29 + ,16 + ,11 + ,6 + ,27 + ,24 + ,31 + ,9 + ,17 + ,9 + ,25 + ,19 + ,20 + ,9 + ,11 + ,7 + ,21 + ,20 + ,16 + ,9 + ,12 + ,9 + ,10 + ,13 + ,22 + ,8 + ,14 + ,10 + ,20 + ,20 + ,20 + ,7 + ,11 + ,9 + ,26 + ,22 + ,28 + ,16 + ,16 + ,8 + ,24 + ,24 + ,38 + ,11 + ,21 + ,7 + ,29 + ,29 + ,22 + ,9 + ,14 + ,6 + ,19 + ,12 + ,20 + ,11 + ,20 + ,13 + ,24 + ,20 + ,17 + ,9 + ,13 + ,6 + ,19 + ,21 + ,28 + ,14 + ,11 + ,8 + ,24 + ,24 + ,22 + ,13 + ,15 + ,10 + ,22 + ,22 + ,31 + ,16 + ,19 + ,16 + ,17 + ,20) + ,dim=c(6 + ,159) + ,dimnames=list(c('CM' + ,'D' + ,'PE' + ,'PC' + ,'PS' + ,'O') + ,1:159)) > y <- array(NA,dim=c(6,159),dimnames=list(c('CM','D','PE','PC','PS','O'),1:159)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par4 = 'no' > par3 = '' > par2 = 'none' > par1 = '5' > #'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(party) Loading required package: survival Loading required package: splines Loading required package: grid Loading required package: modeltools Loading required package: stats4 Loading required package: coin Loading required package: mvtnorm Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric Loading required package: sandwich Loading required package: strucchange Loading required package: vcd Loading required package: MASS Loading required package: colorspace > library(Hmisc) Attaching package: 'Hmisc' The following object(s) are masked from package:survival : untangle.specials The following object(s) are masked from package:base : format.pval, round.POSIXt, trunc.POSIXt, units > par1 <- as.numeric(par1) > par3 <- as.numeric(par3) > x <- data.frame(t(y)) > is.data.frame(x) [1] TRUE > x <- x[!is.na(x[,par1]),] > k <- length(x[1,]) > n <- length(x[,1]) > colnames(x)[par1] [1] "PS" > x[,par1] [1] 24 25 30 19 22 22 25 23 17 21 19 19 15 16 23 27 22 14 22 23 23 21 19 18 20 [26] 23 25 19 24 22 25 26 29 32 25 29 28 17 28 29 26 25 14 25 26 20 18 32 25 25 [51] 23 21 20 15 30 24 26 24 22 14 24 24 24 24 19 31 22 27 19 25 20 21 27 23 25 [76] 20 21 22 23 25 25 17 19 25 19 20 26 23 27 17 17 19 17 22 21 32 21 21 18 18 [101] 23 19 20 21 20 17 18 19 22 15 14 18 24 35 29 21 25 20 22 13 26 17 25 20 19 [126] 21 22 24 21 26 24 16 23 18 16 26 19 21 21 22 23 29 21 21 23 27 25 21 10 20 [151] 26 24 29 19 24 19 24 22 17 > if (par2 == 'kmeans') { + cl <- kmeans(x[,par1], par3) + print(cl) + clm <- matrix(cbind(cl$centers,1:par3),ncol=2) + clm <- clm[sort.list(clm[,1]),] + for (i in 1:par3) { + cl$cluster[cl$cluster==clm[i,2]] <- paste('C',i,sep='') + } + cl$cluster <- as.factor(cl$cluster) + print(cl$cluster) + x[,par1] <- cl$cluster + } > if (par2 == 'quantiles') { + x[,par1] <- cut2(x[,par1],g=par3) + } > if (par2 == 'hclust') { + hc <- hclust(dist(x[,par1])^2, 'cen') + print(hc) + memb <- cutree(hc, k = par3) + dum <- c(mean(x[memb==1,par1])) + for (i in 2:par3) { + dum <- c(dum, mean(x[memb==i,par1])) + } + hcm <- matrix(cbind(dum,1:par3),ncol=2) + hcm <- hcm[sort.list(hcm[,1]),] + for (i in 1:par3) { + memb[memb==hcm[i,2]] <- paste('C',i,sep='') + } + memb <- as.factor(memb) + print(memb) + x[,par1] <- memb + } > if (par2=='equal') { + ed <- cut(as.numeric(x[,par1]),par3,labels=paste('C',1:par3,sep='')) + x[,par1] <- as.factor(ed) + } > table(x[,par1]) 10 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 35 1 1 4 3 3 9 7 16 11 17 14 13 14 17 9 5 2 6 2 1 3 1 > colnames(x) [1] "CM" "D" "PE" "PC" "PS" "O" > colnames(x)[par1] [1] "PS" > x[,par1] [1] 24 25 30 19 22 22 25 23 17 21 19 19 15 16 23 27 22 14 22 23 23 21 19 18 20 [26] 23 25 19 24 22 25 26 29 32 25 29 28 17 28 29 26 25 14 25 26 20 18 32 25 25 [51] 23 21 20 15 30 24 26 24 22 14 24 24 24 24 19 31 22 27 19 25 20 21 27 23 25 [76] 20 21 22 23 25 25 17 19 25 19 20 26 23 27 17 17 19 17 22 21 32 21 21 18 18 [101] 23 19 20 21 20 17 18 19 22 15 14 18 24 35 29 21 25 20 22 13 26 17 25 20 19 [126] 21 22 24 21 26 24 16 23 18 16 26 19 21 21 22 23 29 21 21 23 27 25 21 10 20 [151] 26 24 29 19 24 19 24 22 17 > if (par2 == 'none') { + m <- ctree(as.formula(paste(colnames(x)[par1],' ~ .',sep='')),data = x) + } > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > if (par2 != 'none') { + m <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data = x) + if (par4=='yes') { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'10-Fold Cross Validation',3+2*par3,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'',1,TRUE) + a<-table.element(a,'Prediction (training)',par3+1,TRUE) + a<-table.element(a,'Prediction (testing)',par3+1,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Actual',1,TRUE) + for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE) + a<-table.element(a,'CV',1,TRUE) + for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE) + a<-table.element(a,'CV',1,TRUE) + a<-table.row.end(a) + for (i in 1:10) { + ind <- sample(2, nrow(x), replace=T, prob=c(0.9,0.1)) + m.ct <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data =x[ind==1,]) + if (i==1) { + m.ct.i.pred <- predict(m.ct, newdata=x[ind==1,]) + m.ct.i.actu <- x[ind==1,par1] + m.ct.x.pred <- predict(m.ct, newdata=x[ind==2,]) + m.ct.x.actu <- x[ind==2,par1] + } else { + m.ct.i.pred <- c(m.ct.i.pred,predict(m.ct, newdata=x[ind==1,])) + m.ct.i.actu <- c(m.ct.i.actu,x[ind==1,par1]) + m.ct.x.pred <- c(m.ct.x.pred,predict(m.ct, newdata=x[ind==2,])) + m.ct.x.actu <- c(m.ct.x.actu,x[ind==2,par1]) + } + } + print(m.ct.i.tab <- table(m.ct.i.actu,m.ct.i.pred)) + numer <- 0 + for (i in 1:par3) { + print(m.ct.i.tab[i,i] / sum(m.ct.i.tab[i,])) + numer <- numer + m.ct.i.tab[i,i] + } + print(m.ct.i.cp <- numer / sum(m.ct.i.tab)) + print(m.ct.x.tab <- table(m.ct.x.actu,m.ct.x.pred)) + numer <- 0 + for (i in 1:par3) { + print(m.ct.x.tab[i,i] / sum(m.ct.x.tab[i,])) + numer <- numer + m.ct.x.tab[i,i] + } + print(m.ct.x.cp <- numer / sum(m.ct.x.tab)) + for (i in 1:par3) { + a<-table.row.start(a) + a<-table.element(a,paste('C',i,sep=''),1,TRUE) + for (jjj in 1:par3) a<-table.element(a,m.ct.i.tab[i,jjj]) + a<-table.element(a,round(m.ct.i.tab[i,i]/sum(m.ct.i.tab[i,]),4)) + for (jjj in 1:par3) a<-table.element(a,m.ct.x.tab[i,jjj]) + a<-table.element(a,round(m.ct.x.tab[i,i]/sum(m.ct.x.tab[i,]),4)) + a<-table.row.end(a) + } + a<-table.row.start(a) + a<-table.element(a,'Overall',1,TRUE) + for (jjj in 1:par3) a<-table.element(a,'-') + a<-table.element(a,round(m.ct.i.cp,4)) + for (jjj in 1:par3) a<-table.element(a,'-') + a<-table.element(a,round(m.ct.x.cp,4)) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/1es0t1292573946.tab") + } + } > m Conditional inference tree with 3 terminal nodes Response: PS Inputs: CM, D, PE, PC, O Number of observations: 159 1) CM <= 24; criterion = 1, statistic = 28.54 2) O <= 17; criterion = 1, statistic = 18.463 3)* weights = 16 2) O > 17 4)* weights = 96 1) CM > 24 5)* weights = 47 > postscript(file="/var/www/html/rcomp/tmp/271he1292573946.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(m) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/371he1292573946.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,par1] ~ as.factor(where(m)),main='Response by Terminal Node',xlab='Terminal Node',ylab='Response') > dev.off() null device 1 > if (par2 == 'none') { + forec <- predict(m) + result <- as.data.frame(cbind(x[,par1],forec,x[,par1]-forec)) + colnames(result) <- c('Actuals','Forecasts','Residuals') + print(result) + } Actuals Forecasts Residuals 1 24 21.58333 2.41666667 2 25 24.91489 0.08510638 3 30 21.58333 8.41666667 4 19 21.58333 -2.58333333 5 22 21.58333 0.41666667 6 22 21.58333 0.41666667 7 25 21.58333 3.41666667 8 23 21.58333 1.41666667 9 17 21.58333 -4.58333333 10 21 21.58333 -0.58333333 11 19 21.58333 -2.58333333 12 19 24.91489 -5.91489362 13 15 21.58333 -6.58333333 14 16 17.37500 -1.37500000 15 23 21.58333 1.41666667 16 27 21.58333 5.41666667 17 22 21.58333 0.41666667 18 14 17.37500 -3.37500000 19 22 21.58333 0.41666667 20 23 24.91489 -1.91489362 21 23 24.91489 -1.91489362 22 21 24.91489 -3.91489362 23 19 21.58333 -2.58333333 24 18 24.91489 -6.91489362 25 20 21.58333 -1.58333333 26 23 21.58333 1.41666667 27 25 21.58333 3.41666667 28 19 21.58333 -2.58333333 29 24 21.58333 2.41666667 30 22 21.58333 0.41666667 31 25 21.58333 3.41666667 32 26 24.91489 1.08510638 33 29 24.91489 4.08510638 34 32 24.91489 7.08510638 35 25 21.58333 3.41666667 36 29 24.91489 4.08510638 37 28 24.91489 3.08510638 38 17 21.58333 -4.58333333 39 28 24.91489 3.08510638 40 29 24.91489 4.08510638 41 26 24.91489 1.08510638 42 25 21.58333 3.41666667 43 14 17.37500 -3.37500000 44 25 21.58333 3.41666667 45 26 21.58333 4.41666667 46 20 21.58333 -1.58333333 47 18 21.58333 -3.58333333 48 32 24.91489 7.08510638 49 25 24.91489 0.08510638 50 25 21.58333 3.41666667 51 23 24.91489 -1.91489362 52 21 21.58333 -0.58333333 53 20 21.58333 -1.58333333 54 15 17.37500 -2.37500000 55 30 24.91489 5.08510638 56 24 24.91489 -0.91489362 57 26 21.58333 4.41666667 58 24 21.58333 2.41666667 59 22 21.58333 0.41666667 60 14 17.37500 -3.37500000 61 24 21.58333 2.41666667 62 24 21.58333 2.41666667 63 24 21.58333 2.41666667 64 24 21.58333 2.41666667 65 19 21.58333 -2.58333333 66 31 24.91489 6.08510638 67 22 24.91489 -2.91489362 68 27 21.58333 5.41666667 69 19 21.58333 -2.58333333 70 25 21.58333 3.41666667 71 20 24.91489 -4.91489362 72 21 21.58333 -0.58333333 73 27 24.91489 2.08510638 74 23 24.91489 -1.91489362 75 25 24.91489 0.08510638 76 20 21.58333 -1.58333333 77 21 21.58333 -0.58333333 78 22 21.58333 0.41666667 79 23 21.58333 1.41666667 80 25 24.91489 0.08510638 81 25 21.58333 3.41666667 82 17 21.58333 -4.58333333 83 19 21.58333 -2.58333333 84 25 21.58333 3.41666667 85 19 21.58333 -2.58333333 86 20 21.58333 -1.58333333 87 26 21.58333 4.41666667 88 23 21.58333 1.41666667 89 27 24.91489 2.08510638 90 17 21.58333 -4.58333333 91 17 24.91489 -7.91489362 92 19 21.58333 -2.58333333 93 17 21.58333 -4.58333333 94 22 21.58333 0.41666667 95 21 21.58333 -0.58333333 96 32 24.91489 7.08510638 97 21 24.91489 -3.91489362 98 21 24.91489 -3.91489362 99 18 21.58333 -3.58333333 100 18 21.58333 -3.58333333 101 23 21.58333 1.41666667 102 19 21.58333 -2.58333333 103 20 17.37500 2.62500000 104 21 21.58333 -0.58333333 105 20 21.58333 -1.58333333 106 17 24.91489 -7.91489362 107 18 17.37500 0.62500000 108 19 17.37500 1.62500000 109 22 24.91489 -2.91489362 110 15 21.58333 -6.58333333 111 14 21.58333 -7.58333333 112 18 24.91489 -6.91489362 113 24 21.58333 2.41666667 114 35 24.91489 10.08510638 115 29 21.58333 7.41666667 116 21 21.58333 -0.58333333 117 25 17.37500 7.62500000 118 20 21.58333 -1.58333333 119 22 24.91489 -2.91489362 120 13 17.37500 -4.37500000 121 26 24.91489 1.08510638 122 17 17.37500 -0.37500000 123 25 21.58333 3.41666667 124 20 21.58333 -1.58333333 125 19 21.58333 -2.58333333 126 21 21.58333 -0.58333333 127 22 21.58333 0.41666667 128 24 21.58333 2.41666667 129 21 21.58333 -0.58333333 130 26 24.91489 1.08510638 131 24 21.58333 2.41666667 132 16 21.58333 -5.58333333 133 23 21.58333 1.41666667 134 18 21.58333 -3.58333333 135 16 21.58333 -5.58333333 136 26 24.91489 1.08510638 137 19 21.58333 -2.58333333 138 21 17.37500 3.62500000 139 21 21.58333 -0.58333333 140 22 17.37500 4.62500000 141 23 24.91489 -1.91489362 142 29 24.91489 4.08510638 143 21 17.37500 3.62500000 144 21 21.58333 -0.58333333 145 23 21.58333 1.41666667 146 27 24.91489 2.08510638 147 25 24.91489 0.08510638 148 21 21.58333 -0.58333333 149 10 17.37500 -7.37500000 150 20 21.58333 -1.58333333 151 26 21.58333 4.41666667 152 24 24.91489 -0.91489362 153 29 24.91489 4.08510638 154 19 17.37500 1.62500000 155 24 21.58333 2.41666667 156 19 21.58333 -2.58333333 157 24 24.91489 -0.91489362 158 22 21.58333 0.41666667 159 17 24.91489 -7.91489362 > if (par2 != 'none') { + print(cbind(as.factor(x[,par1]),predict(m))) + myt <- table(as.factor(x[,par1]),predict(m)) + print(myt) + } > postscript(file="/var/www/html/rcomp/tmp/4zshz1292573946.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > if(par2=='none') { + op <- par(mfrow=c(2,2)) + plot(density(result$Actuals),main='Kernel Density Plot of Actuals') + plot(density(result$Residuals),main='Kernel Density Plot of Residuals') + plot(result$Forecasts,result$Actuals,main='Actuals versus Predictions',xlab='Predictions',ylab='Actuals') + plot(density(result$Forecasts),main='Kernel Density Plot of Predictions') + par(op) + } > if(par2!='none') { + plot(myt,main='Confusion Matrix',xlab='Actual',ylab='Predicted') + } > dev.off() null device 1 > if (par2 == 'none') { + detcoef <- cor(result$Forecasts,result$Actuals) + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goodness of Fit',2,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Correlation',1,TRUE) + a<-table.element(a,round(detcoef,4)) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'R-squared',1,TRUE) + a<-table.element(a,round(detcoef*detcoef,4)) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'RMSE',1,TRUE) + a<-table.element(a,round(sqrt(mean((result$Residuals)^2)),4)) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/51n4w1292573946.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Actuals, Predictions, and Residuals',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'#',header=TRUE) + a<-table.element(a,'Actuals',header=TRUE) + a<-table.element(a,'Forecasts',header=TRUE) + a<-table.element(a,'Residuals',header=TRUE) + a<-table.row.end(a) + for (i in 1:length(result$Actuals)) { + a<-table.row.start(a) + a<-table.element(a,i,header=TRUE) + a<-table.element(a,result$Actuals[i]) + a<-table.element(a,result$Forecasts[i]) + a<-table.element(a,result$Residuals[i]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/6txlz1292573946.tab") + } > if (par2 != 'none') { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Confusion Matrix (predicted in columns / actuals in rows)',par3+1,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'',1,TRUE) + for (i in 1:par3) { + a<-table.element(a,paste('C',i,sep=''),1,TRUE) + } + a<-table.row.end(a) + for (i in 1:par3) { + a<-table.row.start(a) + a<-table.element(a,paste('C',i,sep=''),1,TRUE) + for (j in 1:par3) { + a<-table.element(a,myt[i,j]) + } + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/7suug1292573946.tab") + } > try(system("convert tmp/271he1292573946.ps tmp/271he1292573946.png",intern=TRUE)) character(0) > try(system("convert tmp/371he1292573946.ps tmp/371he1292573946.png",intern=TRUE)) character(0) > try(system("convert tmp/4zshz1292573946.ps tmp/4zshz1292573946.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.917 0.591 6.092