R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. Natural language support but running in an English locale 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(6 + ,73 + ,62 + ,66 + ,4 + ,58 + ,54 + ,54 + ,5 + ,68 + ,41 + ,82 + ,4 + ,62 + ,49 + ,61 + ,4 + ,65 + ,49 + ,65 + ,6 + ,81 + ,72 + ,77 + ,6 + ,73 + ,78 + ,66 + ,4 + ,64 + ,58 + ,66 + ,4 + ,68 + ,58 + ,66 + ,6 + ,51 + ,23 + ,48 + ,4 + ,68 + ,39 + ,57 + ,6 + ,61 + ,63 + ,80 + ,5 + ,69 + ,46 + ,60 + ,4 + ,73 + ,58 + ,70 + ,6 + ,61 + ,39 + ,85 + ,3 + ,62 + ,44 + ,59 + ,5 + ,63 + ,49 + ,72 + ,6 + ,69 + ,57 + ,70 + ,4 + ,47 + ,76 + ,74 + ,6 + ,66 + ,63 + ,70 + ,2 + ,58 + ,18 + ,51 + ,7 + ,63 + ,40 + ,70 + ,5 + ,69 + ,59 + ,71 + ,2 + ,59 + ,62 + ,72 + ,4 + ,59 + ,70 + ,50 + ,4 + ,63 + ,65 + ,69 + ,6 + ,65 + ,56 + ,73 + ,6 + ,65 + ,45 + ,66 + ,5 + ,71 + ,57 + ,73 + ,6 + ,60 + ,50 + ,58 + ,6 + ,81 + ,40 + ,78 + ,4 + ,67 + ,58 + ,83 + ,6 + ,66 + ,49 + ,76 + ,6 + ,62 + ,49 + ,77 + ,6 + ,63 + ,27 + ,79 + ,2 + ,73 + ,51 + ,71 + ,4 + ,55 + ,75 + ,79 + ,5 + ,59 + ,65 + ,60 + ,3 + ,64 + ,47 + ,73 + ,7 + ,63 + ,49 + ,70 + ,5 + ,64 + ,65 + ,42 + ,3 + ,73 + ,61 + ,74 + ,8 + ,54 + ,46 + ,68 + ,8 + ,76 + ,69 + ,83 + ,5 + ,74 + ,55 + ,62 + ,6 + ,63 + ,78 + ,79 + ,3 + ,73 + ,58 + ,61 + ,5 + ,67 + ,34 + ,86 + ,4 + ,68 + ,67 + ,64 + ,5 + ,66 + ,45 + ,75 + ,5 + ,62 + ,68 + ,59 + ,6 + ,71 + ,49 + ,82 + ,5 + ,63 + ,19 + ,61 + ,6 + ,75 + ,72 + ,69 + ,6 + ,77 + ,59 + ,60 + ,4 + ,62 + ,46 + ,59 + ,8 + ,74 + ,56 + ,81 + ,6 + ,67 + ,45 + ,65 + ,4 + ,56 + ,53 + ,60 + ,6 + ,60 + ,67 + ,60 + ,5 + ,58 + ,73 + ,45 + ,5 + ,65 + ,46 + ,75 + ,6 + ,49 + ,70 + ,84 + ,6 + ,61 + ,38 + ,77 + ,6 + ,66 + ,54 + ,64 + ,6 + ,64 + ,46 + ,54 + ,6 + ,65 + ,46 + ,72 + ,6 + ,46 + ,45 + ,56 + ,7 + ,65 + ,47 + ,67 + ,4 + ,81 + ,25 + ,81 + ,4 + ,72 + ,63 + ,73 + ,3 + ,65 + ,46 + ,67 + ,6 + ,74 + ,69 + ,72 + ,5 + ,59 + ,43 + ,69 + ,5 + ,69 + ,49 + ,71 + ,3 + ,58 + ,39 + ,77 + ,5 + ,71 + ,65 + ,63 + ,4 + ,79 + ,54 + ,49 + ,3 + ,68 + ,50 + ,74 + ,7 + ,66 + ,42 + ,76 + ,4 + ,62 + ,45 + ,65 + ,4 + ,69 + ,50 + ,65 + ,5 + ,63 + ,55 + ,69 + ,6 + ,62 + ,38 + ,71 + ,2 + ,61 + ,40 + ,68 + ,2 + ,65 + ,51 + ,49 + ,6 + ,64 + ,49 + ,86 + ,4 + ,56 + ,39 + ,63 + ,5 + ,56 + ,57 + ,77 + ,6 + ,48 + ,30 + ,52 + ,7 + ,74 + ,51 + ,73 + ,8 + ,69 + ,48 + ,63 + ,6 + ,62 + ,56 + ,54 + ,6 + ,73 + ,66 + ,56 + ,3 + ,64 + ,72 + ,54 + ,7 + ,57 + ,28 + ,61 + ,3 + ,57 + ,52 + ,70 + ,6 + ,60 + ,53 + ,68 + ,4 + ,61 + ,70 + ,63 + ,4 + ,72 + ,63 + ,76 + ,6 + ,57 + ,46 + ,69 + ,6 + ,51 + ,45 + ,71 + ,6 + ,63 + ,68 + ,39 + ,4 + ,54 + ,54 + ,54 + ,7 + ,72 + ,60 + ,64 + ,5 + ,62 + ,50 + ,70 + ,7 + ,68 + ,66 + ,76 + ,4 + ,62 + ,56 + ,71 + ,6 + ,63 + ,54 + ,73 + ,6 + ,77 + ,72 + ,81 + ,6 + ,57 + ,34 + ,50 + ,5 + ,57 + ,39 + ,42 + ,5 + ,61 + ,66 + ,66 + ,6 + ,65 + ,27 + ,77 + ,7 + ,63 + ,63 + ,62 + ,4 + ,66 + ,65 + ,66 + ,4 + ,68 + ,63 + ,69 + ,8 + ,72 + ,49 + ,72 + ,6 + ,68 + ,42 + ,67 + ,3 + ,59 + ,51 + ,59 + ,4 + ,56 + ,50 + ,66 + ,5 + ,62 + ,64 + ,68 + ,5 + ,72 + ,68 + ,72 + ,6 + ,68 + ,66 + ,73 + ,8 + ,67 + ,59 + ,69 + ,2 + ,54 + ,32 + ,57 + ,4 + ,69 + ,62 + ,55 + ,7 + ,61 + ,52 + ,72 + ,5 + ,55 + ,34 + ,68 + ,6 + ,75 + ,63 + ,83 + ,6 + ,55 + ,48 + ,74 + ,4 + ,49 + ,53 + ,72 + ,5 + ,54 + ,39 + ,66 + ,6 + ,66 + ,51 + ,61 + ,6 + ,73 + ,60 + ,86 + ,6 + ,63 + ,70 + ,81 + ,6 + ,61 + ,40 + ,79 + ,5 + ,74 + ,61 + ,73 + ,5 + ,81 + ,35 + ,59 + ,6 + ,62 + ,39 + ,64 + ,4 + ,64 + ,31 + ,75 + ,6 + ,62 + ,36 + ,68 + ,3 + ,85 + ,51 + ,84 + ,6 + ,74 + ,55 + ,68 + ,8 + ,51 + ,67 + ,68 + ,4 + ,66 + ,40 + ,69) + ,dim=c(4 + ,146) + ,dimnames=list(c('Celebrity' + ,'TotNV' + ,'TotAnx' + ,'TotGR') + ,1:146)) > y <- array(NA,dim=c(4,146),dimnames=list(c('Celebrity','TotNV','TotAnx','TotGR'),1:146)) > 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 = '3' > par2 = 'none' > 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(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] "Celebrity" > x[,par1] [1] 6 4 5 4 4 6 6 4 4 6 4 6 5 4 6 3 5 6 4 6 2 7 5 2 4 4 6 6 5 6 6 4 6 6 6 2 4 [38] 5 3 7 5 3 8 8 5 6 3 5 4 5 5 6 5 6 6 4 8 6 4 6 5 5 6 6 6 6 6 6 7 4 4 3 6 5 [75] 5 3 5 4 3 7 4 4 5 6 2 2 6 4 5 6 7 8 6 6 3 7 3 6 4 4 6 6 6 4 7 5 7 4 6 6 6 [112] 5 5 6 7 4 4 8 6 3 4 5 5 6 8 2 4 7 5 6 6 4 5 6 6 6 6 5 5 6 4 6 3 6 8 4 > 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]) 2 3 4 5 6 7 8 6 11 32 28 52 10 7 > colnames(x) [1] "Celebrity" "TotNV" "TotAnx" "TotGR" > colnames(x)[par1] [1] "Celebrity" > x[,par1] [1] 6 4 5 4 4 6 6 4 4 6 4 6 5 4 6 3 5 6 4 6 2 7 5 2 4 4 6 6 5 6 6 4 6 6 6 2 4 [38] 5 3 7 5 3 8 8 5 6 3 5 4 5 5 6 5 6 6 4 8 6 4 6 5 5 6 6 6 6 6 6 7 4 4 3 6 5 [75] 5 3 5 4 3 7 4 4 5 6 2 2 6 4 5 6 7 8 6 6 3 7 3 6 4 4 6 6 6 4 7 5 7 4 6 6 6 [112] 5 5 6 7 4 4 8 6 3 4 5 5 6 8 2 4 7 5 6 6 4 5 6 6 6 6 5 5 6 4 6 3 6 8 4 > if (par2 == 'none') { + m <- ctree(as.formula(paste(colnames(x)[par1],' ~ .',sep='')),data = x) + } > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/10h3s1292184433.tab") + } + } > m Conditional inference tree with 1 terminal nodes Response: Celebrity Inputs: TotNV, TotAnx, TotGR Number of observations: 146 1)* weights = 146 > postscript(file="/var/www/html/freestat/rcomp/tmp/2aqkd1292184433.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/freestat/rcomp/tmp/3aqkd1292184433.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 6 5.143836 0.8561644 2 4 5.143836 -1.1438356 3 5 5.143836 -0.1438356 4 4 5.143836 -1.1438356 5 4 5.143836 -1.1438356 6 6 5.143836 0.8561644 7 6 5.143836 0.8561644 8 4 5.143836 -1.1438356 9 4 5.143836 -1.1438356 10 6 5.143836 0.8561644 11 4 5.143836 -1.1438356 12 6 5.143836 0.8561644 13 5 5.143836 -0.1438356 14 4 5.143836 -1.1438356 15 6 5.143836 0.8561644 16 3 5.143836 -2.1438356 17 5 5.143836 -0.1438356 18 6 5.143836 0.8561644 19 4 5.143836 -1.1438356 20 6 5.143836 0.8561644 21 2 5.143836 -3.1438356 22 7 5.143836 1.8561644 23 5 5.143836 -0.1438356 24 2 5.143836 -3.1438356 25 4 5.143836 -1.1438356 26 4 5.143836 -1.1438356 27 6 5.143836 0.8561644 28 6 5.143836 0.8561644 29 5 5.143836 -0.1438356 30 6 5.143836 0.8561644 31 6 5.143836 0.8561644 32 4 5.143836 -1.1438356 33 6 5.143836 0.8561644 34 6 5.143836 0.8561644 35 6 5.143836 0.8561644 36 2 5.143836 -3.1438356 37 4 5.143836 -1.1438356 38 5 5.143836 -0.1438356 39 3 5.143836 -2.1438356 40 7 5.143836 1.8561644 41 5 5.143836 -0.1438356 42 3 5.143836 -2.1438356 43 8 5.143836 2.8561644 44 8 5.143836 2.8561644 45 5 5.143836 -0.1438356 46 6 5.143836 0.8561644 47 3 5.143836 -2.1438356 48 5 5.143836 -0.1438356 49 4 5.143836 -1.1438356 50 5 5.143836 -0.1438356 51 5 5.143836 -0.1438356 52 6 5.143836 0.8561644 53 5 5.143836 -0.1438356 54 6 5.143836 0.8561644 55 6 5.143836 0.8561644 56 4 5.143836 -1.1438356 57 8 5.143836 2.8561644 58 6 5.143836 0.8561644 59 4 5.143836 -1.1438356 60 6 5.143836 0.8561644 61 5 5.143836 -0.1438356 62 5 5.143836 -0.1438356 63 6 5.143836 0.8561644 64 6 5.143836 0.8561644 65 6 5.143836 0.8561644 66 6 5.143836 0.8561644 67 6 5.143836 0.8561644 68 6 5.143836 0.8561644 69 7 5.143836 1.8561644 70 4 5.143836 -1.1438356 71 4 5.143836 -1.1438356 72 3 5.143836 -2.1438356 73 6 5.143836 0.8561644 74 5 5.143836 -0.1438356 75 5 5.143836 -0.1438356 76 3 5.143836 -2.1438356 77 5 5.143836 -0.1438356 78 4 5.143836 -1.1438356 79 3 5.143836 -2.1438356 80 7 5.143836 1.8561644 81 4 5.143836 -1.1438356 82 4 5.143836 -1.1438356 83 5 5.143836 -0.1438356 84 6 5.143836 0.8561644 85 2 5.143836 -3.1438356 86 2 5.143836 -3.1438356 87 6 5.143836 0.8561644 88 4 5.143836 -1.1438356 89 5 5.143836 -0.1438356 90 6 5.143836 0.8561644 91 7 5.143836 1.8561644 92 8 5.143836 2.8561644 93 6 5.143836 0.8561644 94 6 5.143836 0.8561644 95 3 5.143836 -2.1438356 96 7 5.143836 1.8561644 97 3 5.143836 -2.1438356 98 6 5.143836 0.8561644 99 4 5.143836 -1.1438356 100 4 5.143836 -1.1438356 101 6 5.143836 0.8561644 102 6 5.143836 0.8561644 103 6 5.143836 0.8561644 104 4 5.143836 -1.1438356 105 7 5.143836 1.8561644 106 5 5.143836 -0.1438356 107 7 5.143836 1.8561644 108 4 5.143836 -1.1438356 109 6 5.143836 0.8561644 110 6 5.143836 0.8561644 111 6 5.143836 0.8561644 112 5 5.143836 -0.1438356 113 5 5.143836 -0.1438356 114 6 5.143836 0.8561644 115 7 5.143836 1.8561644 116 4 5.143836 -1.1438356 117 4 5.143836 -1.1438356 118 8 5.143836 2.8561644 119 6 5.143836 0.8561644 120 3 5.143836 -2.1438356 121 4 5.143836 -1.1438356 122 5 5.143836 -0.1438356 123 5 5.143836 -0.1438356 124 6 5.143836 0.8561644 125 8 5.143836 2.8561644 126 2 5.143836 -3.1438356 127 4 5.143836 -1.1438356 128 7 5.143836 1.8561644 129 5 5.143836 -0.1438356 130 6 5.143836 0.8561644 131 6 5.143836 0.8561644 132 4 5.143836 -1.1438356 133 5 5.143836 -0.1438356 134 6 5.143836 0.8561644 135 6 5.143836 0.8561644 136 6 5.143836 0.8561644 137 6 5.143836 0.8561644 138 5 5.143836 -0.1438356 139 5 5.143836 -0.1438356 140 6 5.143836 0.8561644 141 4 5.143836 -1.1438356 142 6 5.143836 0.8561644 143 3 5.143836 -2.1438356 144 6 5.143836 0.8561644 145 8 5.143836 2.8561644 146 4 5.143836 -1.1438356 > 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/freestat/rcomp/tmp/4lzjg1292184433.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/freestat/rcomp/tmp/5p0i41292184433.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/freestat/rcomp/tmp/6a0yr1292184433.tab") + } Warning message: In cor(result$Forecasts, result$Actuals) : the standard deviation is zero > 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/freestat/rcomp/tmp/7djff1292184433.tab") + } > > try(system("convert tmp/2aqkd1292184433.ps tmp/2aqkd1292184433.png",intern=TRUE)) character(0) > try(system("convert tmp/3aqkd1292184433.ps tmp/3aqkd1292184433.png",intern=TRUE)) character(0) > try(system("convert tmp/4lzjg1292184433.ps tmp/4lzjg1292184433.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 3.915 0.722 4.089