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(24 + ,26 + ,38 + ,23 + ,10 + ,11 + ,25 + ,23 + ,36 + ,15 + ,10 + ,11 + ,30 + ,25 + ,23 + ,25 + ,10 + ,11 + ,19 + ,23 + ,30 + ,18 + ,10 + ,11 + ,22 + ,19 + ,26 + ,21 + ,10 + ,11 + ,22 + ,29 + ,26 + ,19 + ,10 + ,11 + ,25 + ,25 + ,30 + ,15 + ,13 + ,12 + ,23 + ,21 + ,27 + ,22 + ,10 + ,11 + ,17 + ,22 + ,34 + ,19 + ,10 + ,11 + ,21 + ,25 + ,28 + ,20 + ,13 + ,9 + ,19 + ,24 + ,36 + ,26 + ,10 + ,11 + ,19 + ,18 + ,42 + ,26 + ,10 + ,11 + ,15 + ,22 + ,31 + ,21 + ,10 + ,11 + ,23 + ,22 + ,26 + ,19 + ,10 + ,11 + ,27 + ,28 + ,16 + ,19 + ,13 + ,12 + ,14 + ,12 + ,23 + ,19 + ,10 + ,11 + ,23 + ,20 + ,45 + ,28 + ,10 + ,11 + ,19 + ,21 + ,30 + ,27 + ,10 + ,11 + ,18 + ,23 + ,45 + ,18 + ,10 + ,11 + ,20 + ,28 + ,30 + ,19 + ,10 + ,11 + ,23 + ,24 + ,24 + ,24 + ,10 + ,11 + ,25 + ,24 + ,29 + ,21 + ,13 + ,12 + ,19 + ,24 + ,30 + ,22 + ,13 + ,9 + ,24 + ,23 + ,31 + ,25 + ,10 + ,11 + ,25 + ,29 + ,34 + ,15 + ,10 + ,11 + ,26 + ,24 + ,41 + ,34 + ,10 + ,11 + ,29 + ,18 + ,37 + ,23 + ,10 + ,11 + ,32 + ,25 + ,33 + ,19 + ,10 + ,11 + ,29 + ,26 + ,48 + ,15 + ,10 + ,11 + ,28 + ,22 + ,44 + ,15 + ,10 + ,11 + ,17 + ,22 + ,29 + ,17 + ,10 + ,11 + ,28 + ,22 + ,44 + ,30 + ,13 + ,9 + ,26 + ,30 + ,43 + ,28 + ,10 + ,11 + ,25 + ,23 + ,31 + ,23 + ,10 + ,11 + ,14 + ,17 + ,28 + ,23 + ,10 + ,11 + ,25 + ,23 + ,26 + ,21 + ,10 + ,11 + ,26 + ,23 + ,30 + ,18 + ,10 + ,11 + ,20 + ,25 + ,27 + ,19 + ,15 + ,11 + ,18 + ,24 + ,34 + ,24 + ,10 + ,11 + ,32 + ,24 + ,47 + ,15 + ,10 + ,11 + ,25 + ,21 + ,37 + ,24 + ,13 + ,16 + ,21 + ,24 + ,27 + ,20 + ,10 + ,11 + ,20 + ,28 + ,30 + ,20 + ,10 + ,11 + ,30 + ,20 + ,36 + ,44 + ,10 + ,11 + ,24 + ,29 + ,39 + ,20 + ,10 + ,11 + ,26 + ,27 + ,32 + ,20 + ,10 + ,11 + ,24 + ,22 + ,25 + ,20 + ,10 + ,11 + ,22 + ,28 + ,19 + ,11 + ,10 + ,11 + ,14 + ,16 + ,29 + ,21 + ,10 + ,11 + ,24 + ,25 + ,26 + ,21 + ,13 + ,9 + ,24 + ,24 + ,31 + ,19 + ,13 + ,12 + ,24 + ,28 + ,31 + ,21 + ,10 + ,11 + ,24 + ,24 + ,31 + ,17 + ,10 + ,11 + ,22 + ,24 + ,39 + ,19 + ,10 + ,11 + ,27 + ,21 + ,28 + ,21 + ,10 + ,11 + ,19 + ,25 + ,22 + ,16 + ,10 + ,11 + ,25 + ,25 + ,31 + ,19 + ,10 + ,11 + ,20 + ,22 + ,36 + ,19 + ,10 + ,11 + ,21 + ,23 + ,28 + ,16 + ,10 + ,11 + ,27 + ,26 + ,39 + ,24 + ,10 + ,11 + ,25 + ,25 + ,35 + ,21 + ,10 + ,11 + ,20 + ,21 + ,33 + ,20 + ,10 + ,11 + ,21 + ,25 + ,27 + ,19 + ,10 + ,11 + ,22 + ,24 + ,33 + ,23 + ,10 + ,11 + ,23 + ,29 + ,31 + ,18 + ,10 + ,11 + ,25 + ,22 + ,39 + ,19 + ,10 + ,11 + ,25 + ,27 + ,37 + ,23 + ,10 + ,11 + ,17 + ,26 + ,24 + ,19 + ,10 + ,11 + ,25 + ,24 + ,28 + ,26 + ,13 + ,12 + ,19 + ,27 + ,37 + ,13 + ,13 + ,12 + ,20 + ,24 + ,32 + ,23 + ,10 + ,11 + ,26 + ,24 + ,31 + ,16 + ,13 + ,12 + ,23 + ,29 + ,29 + ,17 + ,13 + ,12 + ,27 + ,22 + ,40 + ,30 + ,10 + ,11 + ,17 + ,24 + ,40 + ,22 + ,10 + ,11 + ,19 + ,24 + ,15 + ,14 + ,10 + ,11 + ,17 + ,23 + ,27 + ,14 + ,13 + ,9 + ,22 + ,20 + ,32 + ,21 + ,13 + ,9 + ,21 + ,27 + ,28 + ,21 + ,10 + ,11 + ,32 + ,26 + ,41 + ,33 + ,10 + ,11 + ,21 + ,25 + ,47 + ,23 + ,10 + ,11 + ,21 + ,21 + ,42 + ,30 + ,10 + ,11 + ,18 + ,19 + ,32 + ,21 + ,11 + ,17 + ,23 + ,21 + ,33 + ,25 + ,10 + ,11 + ,20 + ,16 + ,29 + ,29 + ,10 + ,11 + ,20 + ,29 + ,37 + ,21 + ,10 + ,11 + ,17 + ,15 + ,39 + ,16 + ,10 + ,11 + ,18 + ,17 + ,29 + ,17 + ,10 + ,11 + ,19 + ,15 + ,33 + ,23 + ,10 + ,11 + ,15 + ,21 + ,31 + ,18 + ,13 + ,9 + ,14 + ,19 + ,21 + ,19 + ,10 + ,11 + ,18 + ,24 + ,36 + ,28 + ,10 + ,11 + ,35 + ,17 + ,32 + ,29 + ,10 + ,11 + ,29 + ,23 + ,15 + ,19 + ,10 + ,11 + ,25 + ,14 + ,25 + ,25 + ,13 + ,9 + ,20 + ,19 + ,28 + ,15 + ,10 + ,11 + ,22 + ,24 + ,39 + ,24 + ,10 + ,11 + ,13 + ,13 + ,31 + ,12 + ,13 + ,9 + ,26 + ,22 + ,40 + ,11 + ,10 + ,11 + ,17 + ,16 + ,25 + ,19 + ,10 + ,11 + ,25 + ,19 + ,36 + ,25 + ,10 + ,11 + ,20 + ,25 + ,23 + ,12 + ,10 + ,11 + ,19 + ,25 + ,39 + ,15 + ,10 + ,11 + ,21 + ,23 + ,31 + ,25 + ,10 + ,11 + ,22 + ,24 + ,23 + ,14 + ,10 + ,11 + ,24 + ,26 + ,31 + ,19 + ,10 + ,11 + ,21 + ,26 + ,28 + ,23 + ,13 + ,9 + ,26 + ,25 + ,47 + ,19 + ,13 + ,9 + ,16 + ,21 + ,25 + ,20 + ,10 + ,11 + ,23 + ,26 + ,26 + ,16 + ,13 + ,9 + ,18 + ,23 + ,24 + ,13 + ,12 + ,18 + ,21 + ,13 + ,30 + ,22 + ,10 + ,11 + ,21 + ,24 + ,25 + ,21 + ,13 + ,16 + ,23 + ,14 + ,44 + ,18 + ,15 + ,13 + ,21 + ,10 + ,38 + ,44 + ,10 + ,11 + ,21 + ,24 + ,36 + ,12 + ,10 + ,11 + ,23 + ,22 + ,34 + ,28 + ,13 + ,12 + ,27 + ,24 + ,45 + ,17 + ,13 + ,16 + ,21 + ,20 + ,29 + ,18 + ,10 + ,11 + ,10 + ,13 + ,25 + ,21 + ,10 + ,11 + ,20 + ,20 + ,30 + ,24 + ,10 + ,11 + ,26 + ,22 + ,27 + ,20 + ,10 + ,11 + ,24 + ,24 + ,44 + ,24 + ,10 + ,11 + ,24 + ,20 + ,31 + ,33 + ,10 + ,11 + ,22 + ,22 + ,35 + ,25 + ,10 + ,11 + ,17 + ,20 + ,47 + ,35 + ,10 + ,11) + ,dim=c(6 + ,126) + ,dimnames=list(c('PS' + ,'O' + ,'CMD' + ,'PEC' + ,'happiness' + ,'depression') + ,1:126)) > y <- array(NA,dim=c(6,126),dimnames=list(c('PS','O','CMD','PEC','happiness','depression'),1:126)) > 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 = '2' > #'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] "O" > x[,par1] [1] 26 23 25 23 19 29 25 21 22 25 24 18 22 22 28 12 20 21 23 28 24 24 24 23 29 [26] 24 18 25 26 22 22 22 30 23 17 23 23 25 24 24 21 24 28 20 29 27 22 28 16 25 [51] 24 28 24 24 21 25 25 22 23 26 25 21 25 24 29 22 27 26 24 27 24 24 29 22 24 [76] 24 23 20 27 26 25 21 19 21 16 29 15 17 15 21 19 24 17 23 14 19 24 13 22 16 [101] 19 25 25 23 24 26 26 25 21 26 23 13 24 14 10 24 22 24 20 13 20 22 24 20 22 [126] 20 > 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 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 1 3 2 2 3 3 2 5 7 9 14 12 24 14 8 4 5 6 1 > colnames(x) [1] "PS" "O" "CMD" "PEC" "happiness" [6] "depression" > colnames(x)[par1] [1] "O" > x[,par1] [1] 26 23 25 23 19 29 25 21 22 25 24 18 22 22 28 12 20 21 23 28 24 24 24 23 29 [26] 24 18 25 26 22 22 22 30 23 17 23 23 25 24 24 21 24 28 20 29 27 22 28 16 25 [51] 24 28 24 24 21 25 25 22 23 26 25 21 25 24 29 22 27 26 24 27 24 24 29 22 24 [76] 24 23 20 27 26 25 21 19 21 16 29 15 17 15 21 19 24 17 23 14 19 24 13 22 16 [101] 19 25 25 23 24 26 26 25 21 26 23 13 24 14 10 24 22 24 20 13 20 22 24 20 22 [126] 20 > 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/1wwt41291994884.tab") + } + } > m Conditional inference tree with 3 terminal nodes Response: O Inputs: PS, CMD, PEC, happiness, depression Number of observations: 126 1) PS <= 15; criterion = 0.998, statistic = 12.504 2)* weights = 8 1) PS > 15 3) PEC <= 24; criterion = 0.998, statistic = 12.306 4)* weights = 92 3) PEC > 24 5)* weights = 26 > postscript(file="/var/www/html/freestat/rcomp/tmp/2wwt41291994884.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/3wwt41291994884.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 26 23.55435 2.44565217 2 23 23.55435 -0.55434783 3 25 21.07692 3.92307692 4 23 23.55435 -0.55434783 5 19 23.55435 -4.55434783 6 29 23.55435 5.44565217 7 25 23.55435 1.44565217 8 21 23.55435 -2.55434783 9 22 23.55435 -1.55434783 10 25 23.55435 1.44565217 11 24 21.07692 2.92307692 12 18 21.07692 -3.07692308 13 22 16.62500 5.37500000 14 22 23.55435 -1.55434783 15 28 23.55435 4.44565217 16 12 16.62500 -4.62500000 17 20 21.07692 -1.07692308 18 21 21.07692 -0.07692308 19 23 23.55435 -0.55434783 20 28 23.55435 4.44565217 21 24 23.55435 0.44565217 22 24 23.55435 0.44565217 23 24 23.55435 0.44565217 24 23 21.07692 1.92307692 25 29 23.55435 5.44565217 26 24 21.07692 2.92307692 27 18 23.55435 -5.55434783 28 25 23.55435 1.44565217 29 26 23.55435 2.44565217 30 22 23.55435 -1.55434783 31 22 23.55435 -1.55434783 32 22 21.07692 0.92307692 33 30 21.07692 8.92307692 34 23 23.55435 -0.55434783 35 17 16.62500 0.37500000 36 23 23.55435 -0.55434783 37 23 23.55435 -0.55434783 38 25 23.55435 1.44565217 39 24 23.55435 0.44565217 40 24 23.55435 0.44565217 41 21 23.55435 -2.55434783 42 24 23.55435 0.44565217 43 28 23.55435 4.44565217 44 20 21.07692 -1.07692308 45 29 23.55435 5.44565217 46 27 23.55435 3.44565217 47 22 23.55435 -1.55434783 48 28 23.55435 4.44565217 49 16 16.62500 -0.62500000 50 25 23.55435 1.44565217 51 24 23.55435 0.44565217 52 28 23.55435 4.44565217 53 24 23.55435 0.44565217 54 24 23.55435 0.44565217 55 21 23.55435 -2.55434783 56 25 23.55435 1.44565217 57 25 23.55435 1.44565217 58 22 23.55435 -1.55434783 59 23 23.55435 -0.55434783 60 26 23.55435 2.44565217 61 25 23.55435 1.44565217 62 21 23.55435 -2.55434783 63 25 23.55435 1.44565217 64 24 23.55435 0.44565217 65 29 23.55435 5.44565217 66 22 23.55435 -1.55434783 67 27 23.55435 3.44565217 68 26 23.55435 2.44565217 69 24 21.07692 2.92307692 70 27 23.55435 3.44565217 71 24 23.55435 0.44565217 72 24 23.55435 0.44565217 73 29 23.55435 5.44565217 74 22 21.07692 0.92307692 75 24 23.55435 0.44565217 76 24 23.55435 0.44565217 77 23 23.55435 -0.55434783 78 20 23.55435 -3.55434783 79 27 23.55435 3.44565217 80 26 21.07692 4.92307692 81 25 23.55435 1.44565217 82 21 21.07692 -0.07692308 83 19 23.55435 -4.55434783 84 21 21.07692 -0.07692308 85 16 21.07692 -5.07692308 86 29 23.55435 5.44565217 87 15 23.55435 -8.55434783 88 17 23.55435 -6.55434783 89 15 23.55435 -8.55434783 90 21 16.62500 4.37500000 91 19 16.62500 2.37500000 92 24 21.07692 2.92307692 93 17 21.07692 -4.07692308 94 23 23.55435 -0.55434783 95 14 21.07692 -7.07692308 96 19 23.55435 -4.55434783 97 24 23.55435 0.44565217 98 13 16.62500 -3.62500000 99 22 23.55435 -1.55434783 100 16 23.55435 -7.55434783 101 19 21.07692 -2.07692308 102 25 23.55435 1.44565217 103 25 23.55435 1.44565217 104 23 21.07692 1.92307692 105 24 23.55435 0.44565217 106 26 23.55435 2.44565217 107 26 23.55435 2.44565217 108 25 23.55435 1.44565217 109 21 23.55435 -2.55434783 110 26 23.55435 2.44565217 111 23 23.55435 -0.55434783 112 13 23.55435 -10.55434783 113 24 23.55435 0.44565217 114 14 23.55435 -9.55434783 115 10 21.07692 -11.07692308 116 24 23.55435 0.44565217 117 22 21.07692 0.92307692 118 24 23.55435 0.44565217 119 20 23.55435 -3.55434783 120 13 16.62500 -3.62500000 121 20 23.55435 -3.55434783 122 22 23.55435 -1.55434783 123 24 23.55435 0.44565217 124 20 21.07692 -1.07692308 125 22 21.07692 0.92307692 126 20 21.07692 -1.07692308 > 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/40era1291994884.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/5lf8y1291994884.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/66xom1291994884.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/freestat/rcomp/tmp/7hpop1291994884.tab") + } > > try(system("convert tmp/2wwt41291994884.ps tmp/2wwt41291994884.png",intern=TRUE)) character(0) > try(system("convert tmp/3wwt41291994884.ps tmp/3wwt41291994884.png",intern=TRUE)) character(0) > try(system("convert tmp/40era1291994884.ps tmp/40era1291994884.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.238 0.756 5.252