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(4 + ,1 + ,27 + ,5 + ,26 + ,49 + ,35 + ,4 + ,1 + ,36 + ,4 + ,25 + ,45 + ,34 + ,5 + ,1 + ,25 + ,4 + ,17 + ,54 + ,13 + ,2 + ,1 + ,27 + ,3 + ,37 + ,36 + ,35 + ,3 + ,2 + ,25 + ,3 + ,35 + ,36 + ,28 + ,5 + ,2 + ,44 + ,3 + ,15 + ,53 + ,32 + ,4 + ,1 + ,50 + ,4 + ,27 + ,46 + ,35 + ,4 + ,1 + ,41 + ,4 + ,36 + ,42 + ,36 + ,4 + ,1 + ,48 + ,5 + ,25 + ,41 + ,27 + ,4 + ,2 + ,43 + ,4 + ,30 + ,45 + ,29 + ,5 + ,2 + ,47 + ,2 + ,27 + ,47 + ,27 + ,4 + ,2 + ,41 + ,3 + ,33 + ,42 + ,28 + ,3 + ,1 + ,44 + ,2 + ,29 + ,45 + ,29 + ,4 + ,2 + ,47 + ,5 + ,30 + ,40 + ,28 + ,3 + ,2 + ,40 + ,3 + ,25 + ,45 + ,30 + ,3 + ,2 + ,46 + ,3 + ,23 + ,40 + ,25 + ,4 + ,1 + ,28 + ,3 + ,26 + ,42 + ,15 + ,3 + ,1 + ,56 + ,3 + ,24 + ,45 + ,33 + ,4 + ,2 + ,49 + ,4 + ,35 + ,47 + ,31 + ,2 + ,2 + ,25 + ,4 + ,39 + ,31 + ,37 + ,4 + ,2 + ,41 + ,4 + ,23 + ,46 + ,37 + ,3 + ,2 + ,26 + ,3 + ,32 + ,34 + ,34 + ,4 + ,1 + ,50 + ,5 + ,29 + ,43 + ,32 + ,4 + ,1 + ,47 + ,4 + ,26 + ,45 + ,21 + ,3 + ,1 + ,52 + ,2 + ,21 + ,42 + ,25 + ,3 + ,2 + ,37 + ,5 + ,35 + ,51 + ,32 + ,2 + ,2 + ,41 + ,3 + ,23 + ,44 + ,28 + ,4 + ,1 + ,45 + ,4 + ,21 + ,47 + ,22 + ,5 + ,2 + ,26 + ,4 + ,28 + ,47 + ,25 + ,4 + ,1 + ,NA + ,3 + ,30 + ,41 + ,26 + ,2 + ,1 + ,52 + ,4 + ,21 + ,44 + ,34 + ,5 + ,1 + ,46 + ,2 + ,29 + ,51 + ,34 + ,4 + ,1 + ,58 + ,3 + ,28 + ,46 + ,36 + ,3 + ,1 + ,54 + ,5 + ,19 + ,47 + ,36 + ,4 + ,1 + ,29 + ,3 + ,26 + ,46 + ,26 + ,2 + ,2 + ,50 + ,3 + ,33 + ,38 + ,26 + ,3 + ,1 + ,43 + ,2 + ,34 + ,50 + ,34 + ,3 + ,2 + ,30 + ,3 + ,33 + ,48 + ,33 + ,3 + ,2 + ,47 + ,2 + ,40 + ,36 + ,31 + ,5 + ,1 + ,45 + ,3 + ,24 + ,51 + ,33 + ,NA + ,2 + ,48 + ,1 + ,35 + ,35 + ,22 + ,4 + ,2 + ,48 + ,3 + ,35 + ,49 + ,29 + ,4 + ,2 + ,26 + ,4 + ,32 + ,38 + ,24 + ,4 + ,1 + ,46 + ,5 + ,20 + ,47 + ,37 + ,2 + ,2 + ,NA + ,3 + ,35 + ,36 + ,32 + ,4 + ,2 + ,50 + ,3 + ,35 + ,47 + ,23 + ,3 + ,1 + ,25 + ,4 + ,21 + ,46 + ,29 + ,4 + ,1 + ,47 + ,2 + ,33 + ,43 + ,35 + ,1 + ,2 + ,47 + ,2 + ,40 + ,53 + ,20 + ,2 + ,1 + ,41 + ,3 + ,22 + ,55 + ,28 + ,2 + ,2 + ,45 + ,2 + ,35 + ,39 + ,26 + ,4 + ,2 + ,41 + ,4 + ,20 + ,55 + ,36 + ,3 + ,2 + ,45 + ,5 + ,28 + ,41 + ,26 + ,4 + ,2 + ,40 + ,3 + ,46 + ,33 + ,33 + ,3 + ,1 + ,29 + ,4 + ,18 + ,52 + ,25 + ,3 + ,2 + ,34 + ,5 + ,22 + ,42 + ,29 + ,5 + ,1 + ,45 + ,5 + ,20 + ,56 + ,32 + ,3 + ,2 + ,52 + ,3 + ,25 + ,46 + ,35 + ,2 + ,2 + ,41 + ,4 + ,31 + ,33 + ,24 + ,1 + ,2 + ,48 + ,3 + ,21 + ,51 + ,31 + ,2 + ,2 + ,45 + ,3 + ,23 + ,46 + ,29 + ,5 + ,1 + ,54 + ,2 + ,26 + ,46 + ,27 + ,4 + ,2 + ,25 + ,3 + ,34 + ,50 + ,29 + ,4 + ,2 + ,26 + ,4 + ,31 + ,46 + ,29 + ,3 + ,1 + ,28 + ,4 + ,23 + ,51 + ,27 + ,4 + ,2 + ,50 + ,4 + ,31 + ,48 + ,34 + ,4 + ,2 + ,48 + ,4 + ,26 + ,44 + ,32 + ,2 + ,2 + ,51 + ,3 + ,36 + ,38 + ,31 + ,3 + ,2 + ,53 + ,3 + ,28 + ,42 + ,31 + ,4 + ,1 + ,37 + ,3 + ,34 + ,39 + ,31 + ,3 + ,1 + ,56 + ,2 + ,25 + ,45 + ,16 + ,2 + ,1 + ,43 + ,3 + ,33 + ,31 + ,25 + ,4 + ,1 + ,34 + ,3 + ,46 + ,29 + ,27 + ,4 + ,1 + ,42 + ,3 + ,24 + ,48 + ,32 + ,3 + ,2 + ,32 + ,3 + ,32 + ,38 + ,28 + ,5 + ,2 + ,31 + ,5 + ,33 + ,55 + ,25 + ,1 + ,1 + ,46 + ,3 + ,42 + ,32 + ,25 + ,3 + ,2 + ,30 + ,5 + ,17 + ,51 + ,36 + ,3 + ,2 + ,47 + ,4 + ,36 + ,53 + ,36 + ,5 + ,2 + ,33 + ,4 + ,40 + ,47 + ,36 + ,2 + ,1 + ,25 + ,4 + ,30 + ,45 + ,27 + ,3 + ,1 + ,25 + ,5 + ,19 + ,33 + ,29 + ,3 + ,2 + ,21 + ,4 + ,33 + ,49 + ,32 + ,4 + ,2 + ,36 + ,5 + ,35 + ,46 + ,29 + ,2 + ,2 + ,50 + ,3 + ,23 + ,42 + ,31 + ,4 + ,2 + ,48 + ,3 + ,15 + ,56 + ,34 + ,3 + ,2 + ,48 + ,2 + ,38 + ,35 + ,27 + ,3 + ,1 + ,25 + ,3 + ,37 + ,40 + ,28 + ,3 + ,1 + ,48 + ,4 + ,23 + ,44 + ,32 + ,2 + ,2 + ,49 + ,5 + ,41 + ,46 + ,33 + ,3 + ,1 + ,27 + ,5 + ,34 + ,46 + ,29 + ,2 + ,1 + ,28 + ,3 + ,38 + ,39 + ,32 + ,4 + ,2 + ,43 + ,2 + ,45 + ,35 + ,35 + ,4 + ,2 + ,48 + ,3 + ,27 + ,48 + ,33 + ,2 + ,2 + ,48 + ,4 + ,46 + ,42 + ,27 + ,1 + ,1 + ,25 + ,1 + ,26 + ,39 + ,16 + ,5 + ,2 + ,49 + ,4 + ,44 + ,39 + ,32 + ,4 + ,1 + ,26 + ,3 + ,36 + ,41 + ,26 + ,4 + ,1 + ,51 + ,3 + ,20 + ,52 + ,32 + ,4 + ,2 + ,25 + ,4 + ,44 + ,45 + ,38 + ,3 + ,1 + ,29 + ,3 + ,27 + ,42 + ,24 + ,3 + ,1 + ,29 + ,4 + ,27 + ,44 + ,26 + ,1 + ,1 + ,43 + ,2 + ,41 + ,33 + ,19 + ,5 + ,2 + ,46 + ,3 + ,30 + ,42 + ,37 + ,3 + ,1 + ,44 + ,3 + ,33 + ,46 + ,25 + ,3 + ,1 + ,25 + ,3 + ,37 + ,45 + ,24 + ,2 + ,1 + ,51 + ,2 + ,30 + ,40 + ,23 + ,4 + ,1 + ,42 + ,5 + ,20 + ,48 + ,28 + ,4 + ,2 + ,53 + ,5 + ,44 + ,32 + ,38 + ,3 + ,1 + ,25 + ,4 + ,20 + ,53 + ,28 + ,4 + ,2 + ,49 + ,2 + ,33 + ,39 + ,28 + ,4 + ,1 + ,51 + ,3 + ,31 + ,45 + ,26 + ,2 + ,2 + ,20 + ,3 + ,23 + ,36 + ,21 + ,3 + ,2 + ,44 + ,3 + ,33 + ,38 + ,35 + ,3 + ,2 + ,38 + ,4 + ,33 + ,49 + ,31 + ,3 + ,1 + ,46 + ,5 + ,32 + ,46 + ,34 + ,4 + ,2 + ,42 + ,4 + ,25 + ,43 + ,30 + ,5 + ,1 + ,29 + ,NA + ,22 + ,37 + ,30 + ,3 + ,2 + ,46 + ,4 + ,16 + ,48 + ,24 + ,3 + ,2 + ,49 + ,2 + ,36 + ,45 + ,27 + ,2 + ,2 + ,51 + ,3 + ,35 + ,32 + ,26 + ,3 + ,1 + ,38 + ,3 + ,25 + ,46 + ,30 + ,1 + ,1 + ,41 + ,1 + ,27 + ,20 + ,15 + ,4 + ,2 + ,47 + ,3 + ,32 + ,42 + ,28 + ,4 + ,2 + ,44 + ,3 + ,36 + ,45 + ,34 + ,4 + ,2 + ,47 + ,3 + ,51 + ,29 + ,29 + ,3 + ,2 + ,46 + ,3 + ,30 + ,51 + ,26 + ,5 + ,1 + ,44 + ,4 + ,20 + ,55 + ,31 + ,2 + ,2 + ,28 + ,3 + ,29 + ,50 + ,28 + ,2 + ,2 + ,47 + ,4 + ,26 + ,44 + ,33 + ,3 + ,2 + ,28 + ,4 + ,20 + ,41 + ,32 + ,3 + ,1 + ,41 + ,5 + ,40 + ,40 + ,33 + ,2 + ,2 + ,45 + ,4 + ,29 + ,47 + ,31 + ,1 + ,2 + ,46 + ,4 + ,32 + ,42 + ,37 + ,3 + ,1 + ,46 + ,4 + ,33 + ,40 + ,27 + ,5 + ,2 + ,22 + ,3 + ,32 + ,51 + ,19 + ,4 + ,2 + ,33 + ,3 + ,34 + ,43 + ,27 + ,4 + ,1 + ,41 + ,4 + ,24 + ,45 + ,31 + ,4 + ,2 + ,47 + ,5 + ,25 + ,41 + ,38 + ,3 + ,1 + ,25 + ,3 + ,41 + ,41 + ,22 + ,5 + ,2 + ,42 + ,3 + ,39 + ,37 + ,35 + ,3 + ,2 + ,47 + ,3 + ,21 + ,46 + ,35 + ,3 + ,2 + ,50 + ,3 + ,38 + ,38 + ,30 + ,3 + ,1 + ,55 + ,5 + ,28 + ,39 + ,41 + ,3 + ,1 + ,21 + ,3 + ,37 + ,45 + ,25 + ,4 + ,1 + ,NA + ,3 + ,26 + ,46 + ,28 + ,2 + ,1 + ,52 + ,3 + ,30 + ,39 + ,45 + ,2 + ,2 + ,49 + ,4 + ,25 + ,21 + ,21 + ,4 + ,2 + ,46 + ,4 + ,38 + ,31 + ,33 + ,3 + ,1 + ,NA + ,4 + ,31 + ,35 + ,25 + ,3 + ,2 + ,45 + ,3 + ,31 + ,49 + ,29 + ,2 + ,2 + ,52 + ,3 + ,27 + ,40 + ,31 + ,3 + ,1 + ,NA + ,3 + ,21 + ,45 + ,29 + ,3 + ,2 + ,40 + ,4 + ,26 + ,46 + ,31 + ,4 + ,2 + ,49 + ,4 + ,37 + ,45 + ,31 + ,1 + ,1 + ,38 + ,5 + ,28 + ,34 + ,25 + ,1 + ,1 + ,32 + ,5 + ,29 + ,41 + ,27 + ,5 + ,2 + ,46 + ,4 + ,33 + ,43 + ,26 + ,4 + ,2 + ,32 + ,3 + ,41 + ,45 + ,26 + ,3 + ,2 + ,41 + ,3 + ,19 + ,48 + ,23 + ,3 + ,2 + ,43 + ,3 + ,37 + ,43 + ,27 + ,4 + ,1 + ,44 + ,4 + ,36 + ,45 + ,24 + ,3 + ,1 + ,47 + ,5 + ,27 + ,45 + ,35 + ,2 + ,2 + ,28 + ,3 + ,33 + ,34 + ,24 + ,1 + ,1 + ,52 + ,1 + ,29 + ,40 + ,32 + ,1 + ,1 + ,27 + ,2 + ,42 + ,40 + ,24 + ,5 + ,2 + ,45 + ,5 + ,27 + ,55 + ,24 + ,4 + ,1 + ,27 + ,4 + ,47 + ,44 + ,38 + ,3 + ,1 + ,25 + ,4 + ,17 + ,44 + ,36 + ,4 + ,1 + ,28 + ,4 + ,34 + ,48 + ,24 + ,5 + ,1 + ,25 + ,3 + ,32 + ,51 + ,18 + ,4 + ,1 + ,52 + ,4 + ,25 + ,49 + ,34 + ,4 + ,1 + ,44 + ,3 + ,27 + ,33 + ,23 + ,2 + ,2 + ,43 + ,3 + ,37 + ,43 + ,35 + ,3 + ,2 + ,47 + ,4 + ,34 + ,44 + ,22 + ,4 + ,2 + ,52 + ,4 + ,27 + ,44 + ,34 + ,3 + ,2 + ,40 + ,2 + ,37 + ,41 + ,28 + ,4 + ,1 + ,42 + ,3 + ,32 + ,45 + ,34 + ,3 + ,1 + ,45 + ,5 + ,26 + ,44 + ,32 + ,4 + ,1 + ,45 + ,2 + ,29 + ,44 + ,24 + ,1 + ,1 + ,50 + ,5 + ,28 + ,40 + ,34 + ,2 + ,1 + ,49 + ,3 + ,19 + ,48 + ,33 + ,3 + ,1 + ,52 + ,2 + ,46 + ,49 + ,33 + ,3 + ,2 + ,48 + ,3 + ,31 + ,46 + ,29 + ,5 + ,2 + ,51 + ,3 + ,42 + ,49 + ,38 + ,4 + ,2 + ,49 + ,4 + ,33 + ,55 + ,24 + ,3 + ,2 + ,31 + ,4 + ,39 + ,51 + ,25 + ,3 + ,2 + ,43 + ,3 + ,27 + ,46 + ,37 + ,3 + ,2 + ,31 + ,3 + ,35 + ,37 + ,33 + ,3 + ,2 + ,28 + ,4 + ,23 + ,43 + ,30 + ,4 + ,2 + ,43 + ,4 + ,32 + ,41 + ,22 + ,3 + ,2 + ,31 + ,3 + ,22 + ,45 + ,28 + ,2 + ,2 + ,51 + ,3 + ,17 + ,39 + ,24 + ,4 + ,2 + ,58 + ,4 + ,35 + ,38 + ,33 + ,2 + ,2 + ,25 + ,5 + ,34 + ,41 + ,37) + ,dim=c(7 + ,195) + ,dimnames=list(c('Teamwork33' + ,'geslacht' + ,'leeftijd' + ,'opleiding' + ,'Neuroticisme' + ,'Extraversie' + ,'Openheid') + ,1:195)) > y <- array(NA,dim=c(7,195),dimnames=list(c('Teamwork33','geslacht','leeftijd','opleiding','Neuroticisme','Extraversie','Openheid'),1:195)) > 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] "Teamwork33" > x[,par1] [1] 4 4 5 2 3 5 4 4 4 4 5 4 3 4 3 3 4 3 4 2 4 3 4 4 3 3 2 4 5 4 2 5 4 3 4 2 3 [38] 3 3 5 4 4 4 2 4 3 4 1 2 2 4 3 4 3 3 5 3 2 1 2 5 4 4 3 4 4 2 3 4 3 2 4 4 3 [75] 5 1 3 3 5 2 3 3 4 2 4 3 3 3 2 3 2 4 4 2 1 5 4 4 4 3 3 1 5 3 3 2 4 4 3 4 4 [112] 2 3 3 3 4 5 3 3 2 3 1 4 4 4 3 5 2 2 3 3 2 1 3 5 4 4 4 3 5 3 3 3 3 4 2 2 4 [149] 3 3 2 3 3 4 1 1 5 4 3 3 4 3 2 1 1 5 4 3 4 5 4 4 2 3 4 3 4 3 4 1 2 3 3 5 4 [186] 3 3 3 3 4 3 2 4 2 > 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]) 1 2 3 4 5 12 31 67 64 20 > colnames(x) [1] "Teamwork33" "geslacht" "leeftijd" "opleiding" "Neuroticisme" [6] "Extraversie" "Openheid" > colnames(x)[par1] [1] "Teamwork33" > x[,par1] [1] 4 4 5 2 3 5 4 4 4 4 5 4 3 4 3 3 4 3 4 2 4 3 4 4 3 3 2 4 5 4 2 5 4 3 4 2 3 [38] 3 3 5 4 4 4 2 4 3 4 1 2 2 4 3 4 3 3 5 3 2 1 2 5 4 4 3 4 4 2 3 4 3 2 4 4 3 [75] 5 1 3 3 5 2 3 3 4 2 4 3 3 3 2 3 2 4 4 2 1 5 4 4 4 3 3 1 5 3 3 2 4 4 3 4 4 [112] 2 3 3 3 4 5 3 3 2 3 1 4 4 4 3 5 2 2 3 3 2 1 3 5 4 4 4 3 5 3 3 3 3 4 2 2 4 [149] 3 3 2 3 3 4 1 1 5 4 3 3 4 3 2 1 1 5 4 3 4 5 4 4 2 3 4 3 4 3 4 1 2 3 3 5 4 [186] 3 3 3 3 4 3 2 4 2 > 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/1gxsd1293311349.tab") + } + } > m Conditional inference tree with 2 terminal nodes Response: Teamwork33 Inputs: geslacht, leeftijd, opleiding, Neuroticisme, Extraversie, Openheid Number of observations: 194 1) Extraversie <= 40; criterion = 1, statistic = 22.749 2)* weights = 55 1) Extraversie > 40 3)* weights = 139 > postscript(file="/var/www/html/freestat/rcomp/tmp/2gxsd1293311349.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/396sy1293311349.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 4 3.460432 0.5395683 2 4 3.460432 0.5395683 3 5 3.460432 1.5395683 4 2 2.727273 -0.7272727 5 3 2.727273 0.2727273 6 5 3.460432 1.5395683 7 4 3.460432 0.5395683 8 4 3.460432 0.5395683 9 4 3.460432 0.5395683 10 4 3.460432 0.5395683 11 5 3.460432 1.5395683 12 4 3.460432 0.5395683 13 3 3.460432 -0.4604317 14 4 2.727273 1.2727273 15 3 3.460432 -0.4604317 16 3 2.727273 0.2727273 17 4 3.460432 0.5395683 18 3 3.460432 -0.4604317 19 4 3.460432 0.5395683 20 2 2.727273 -0.7272727 21 4 3.460432 0.5395683 22 3 2.727273 0.2727273 23 4 3.460432 0.5395683 24 4 3.460432 0.5395683 25 3 3.460432 -0.4604317 26 3 3.460432 -0.4604317 27 2 3.460432 -1.4604317 28 4 3.460432 0.5395683 29 5 3.460432 1.5395683 30 4 3.460432 0.5395683 31 2 3.460432 -1.4604317 32 5 3.460432 1.5395683 33 4 3.460432 0.5395683 34 3 3.460432 -0.4604317 35 4 3.460432 0.5395683 36 2 2.727273 -0.7272727 37 3 3.460432 -0.4604317 38 3 3.460432 -0.4604317 39 3 2.727273 0.2727273 40 5 3.460432 1.5395683 41 4 3.460432 0.5395683 42 4 2.727273 1.2727273 43 4 3.460432 0.5395683 44 2 2.727273 -0.7272727 45 4 3.460432 0.5395683 46 3 3.460432 -0.4604317 47 4 3.460432 0.5395683 48 1 3.460432 -2.4604317 49 2 3.460432 -1.4604317 50 2 2.727273 -0.7272727 51 4 3.460432 0.5395683 52 3 3.460432 -0.4604317 53 4 2.727273 1.2727273 54 3 3.460432 -0.4604317 55 3 3.460432 -0.4604317 56 5 3.460432 1.5395683 57 3 3.460432 -0.4604317 58 2 2.727273 -0.7272727 59 1 3.460432 -2.4604317 60 2 3.460432 -1.4604317 61 5 3.460432 1.5395683 62 4 3.460432 0.5395683 63 4 3.460432 0.5395683 64 3 3.460432 -0.4604317 65 4 3.460432 0.5395683 66 4 3.460432 0.5395683 67 2 2.727273 -0.7272727 68 3 3.460432 -0.4604317 69 4 2.727273 1.2727273 70 3 3.460432 -0.4604317 71 2 2.727273 -0.7272727 72 4 2.727273 1.2727273 73 4 3.460432 0.5395683 74 3 2.727273 0.2727273 75 5 3.460432 1.5395683 76 1 2.727273 -1.7272727 77 3 3.460432 -0.4604317 78 3 3.460432 -0.4604317 79 5 3.460432 1.5395683 80 2 3.460432 -1.4604317 81 3 2.727273 0.2727273 82 3 3.460432 -0.4604317 83 4 3.460432 0.5395683 84 2 3.460432 -1.4604317 85 4 3.460432 0.5395683 86 3 2.727273 0.2727273 87 3 2.727273 0.2727273 88 3 3.460432 -0.4604317 89 2 3.460432 -1.4604317 90 3 3.460432 -0.4604317 91 2 2.727273 -0.7272727 92 4 2.727273 1.2727273 93 4 3.460432 0.5395683 94 2 3.460432 -1.4604317 95 1 2.727273 -1.7272727 96 5 2.727273 2.2727273 97 4 3.460432 0.5395683 98 4 3.460432 0.5395683 99 4 3.460432 0.5395683 100 3 3.460432 -0.4604317 101 3 3.460432 -0.4604317 102 1 2.727273 -1.7272727 103 5 3.460432 1.5395683 104 3 3.460432 -0.4604317 105 3 3.460432 -0.4604317 106 2 2.727273 -0.7272727 107 4 3.460432 0.5395683 108 4 2.727273 1.2727273 109 3 3.460432 -0.4604317 110 4 2.727273 1.2727273 111 4 3.460432 0.5395683 112 2 2.727273 -0.7272727 113 3 2.727273 0.2727273 114 3 3.460432 -0.4604317 115 3 3.460432 -0.4604317 116 4 3.460432 0.5395683 117 5 2.727273 2.2727273 118 3 3.460432 -0.4604317 119 3 3.460432 -0.4604317 120 2 2.727273 -0.7272727 121 3 3.460432 -0.4604317 122 1 2.727273 -1.7272727 123 4 3.460432 0.5395683 124 4 3.460432 0.5395683 125 4 2.727273 1.2727273 126 3 3.460432 -0.4604317 127 5 3.460432 1.5395683 128 2 3.460432 -1.4604317 129 2 3.460432 -1.4604317 130 3 3.460432 -0.4604317 131 3 2.727273 0.2727273 132 2 3.460432 -1.4604317 133 1 3.460432 -2.4604317 134 3 2.727273 0.2727273 135 5 3.460432 1.5395683 136 4 3.460432 0.5395683 137 4 3.460432 0.5395683 138 4 3.460432 0.5395683 139 3 3.460432 -0.4604317 140 5 2.727273 2.2727273 141 3 3.460432 -0.4604317 142 3 2.727273 0.2727273 143 3 2.727273 0.2727273 144 3 3.460432 -0.4604317 145 4 3.460432 0.5395683 146 2 2.727273 -0.7272727 147 2 2.727273 -0.7272727 148 4 2.727273 1.2727273 149 3 2.727273 0.2727273 150 3 3.460432 -0.4604317 151 2 2.727273 -0.7272727 152 3 3.460432 -0.4604317 153 3 3.460432 -0.4604317 154 4 3.460432 0.5395683 155 1 2.727273 -1.7272727 156 1 3.460432 -2.4604317 157 5 3.460432 1.5395683 158 4 3.460432 0.5395683 159 3 3.460432 -0.4604317 160 3 3.460432 -0.4604317 161 4 3.460432 0.5395683 162 3 3.460432 -0.4604317 163 2 2.727273 -0.7272727 164 1 2.727273 -1.7272727 165 1 2.727273 -1.7272727 166 5 3.460432 1.5395683 167 4 3.460432 0.5395683 168 3 3.460432 -0.4604317 169 4 3.460432 0.5395683 170 5 3.460432 1.5395683 171 4 3.460432 0.5395683 172 4 2.727273 1.2727273 173 2 3.460432 -1.4604317 174 3 3.460432 -0.4604317 175 4 3.460432 0.5395683 176 3 3.460432 -0.4604317 177 4 3.460432 0.5395683 178 3 3.460432 -0.4604317 179 4 3.460432 0.5395683 180 1 2.727273 -1.7272727 181 2 3.460432 -1.4604317 182 3 3.460432 -0.4604317 183 3 3.460432 -0.4604317 184 5 3.460432 1.5395683 185 4 3.460432 0.5395683 186 3 3.460432 -0.4604317 187 3 3.460432 -0.4604317 188 3 2.727273 0.2727273 189 3 3.460432 -0.4604317 190 4 3.460432 0.5395683 191 3 3.460432 -0.4604317 192 2 2.727273 -0.7272727 193 4 2.727273 1.2727273 194 2 3.460432 -1.4604317 > 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/41xrj1293311349.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/54yp71293311349.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/6qy6v1293311349.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/71p5y1293311349.tab") + } > > try(system("convert tmp/2gxsd1293311349.ps tmp/2gxsd1293311349.png",intern=TRUE)) character(0) > try(system("convert tmp/396sy1293311349.ps tmp/396sy1293311349.png",intern=TRUE)) character(0) > try(system("convert tmp/41xrj1293311349.ps tmp/41xrj1293311349.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.972 0.764 5.135