x <- array(list(23 ,13 ,14 ,22 ,11 ,23 ,8 ,1 ,6 ,15 ,20 ,12 ,7 ,20 ,22 ,24 ,4 ,2 ,5 ,23 ,26 ,26 ,22 ,25 ,23 ,24 ,7 ,2 ,20 ,26 ,19 ,16 ,12 ,23 ,21 ,21 ,4 ,2 ,12 ,19 ,17 ,18 ,15 ,20 ,19 ,21 ,4 ,2 ,11 ,19 ,17 ,12 ,9 ,22 ,12 ,19 ,5 ,2 ,12 ,16 ,21 ,18 ,20 ,18 ,24 ,12 ,15 ,1 ,11 ,23 ,18 ,20 ,10 ,22 ,21 ,21 ,5 ,1 ,9 ,22 ,16 ,18 ,12 ,23 ,21 ,25 ,7 ,2 ,13 ,19 ,26 ,24 ,23 ,28 ,26 ,27 ,4 ,2 ,9 ,24 ,20 ,17 ,10 ,19 ,18 ,21 ,4 ,1 ,14 ,19 ,14 ,19 ,11 ,26 ,21 ,27 ,7 ,1 ,12 ,25 ,22 ,12 ,20 ,27 ,22 ,20 ,8 ,1 ,18 ,23 ,23 ,25 ,11 ,23 ,26 ,16 ,4 ,2 ,9 ,31 ,25 ,23 ,22 ,27 ,20 ,26 ,8 ,1 ,15 ,29 ,24 ,22 ,19 ,23 ,20 ,24 ,4 ,2 ,12 ,18 ,24 ,23 ,20 ,23 ,26 ,25 ,5 ,2 ,12 ,17 ,16 ,16 ,16 ,19 ,27 ,25 ,16 ,1 ,12 ,22 ,16 ,16 ,12 ,21 ,27 ,27 ,7 ,1 ,15 ,21 ,20 ,15 ,14 ,25 ,16 ,23 ,4 ,2 ,11 ,24 ,20 ,24 ,14 ,22 ,26 ,22 ,6 ,1 ,13 ,22 ,15 ,18 ,9 ,13 ,20 ,10 ,4 ,1 ,10 ,16 ,22 ,23 ,19 ,12 ,25 ,25 ,5 ,2 ,17 ,22 ,20 ,18 ,17 ,20 ,16 ,18 ,4 ,1 ,13 ,21 ,20 ,19 ,14 ,24 ,20 ,21 ,4 ,1 ,17 ,25 ,24 ,17 ,19 ,23 ,20 ,20 ,6 ,1 ,15 ,22 ,27 ,22 ,20 ,25 ,24 ,18 ,4 ,1 ,13 ,24 ,25 ,22 ,20 ,28 ,24 ,25 ,4 ,1 ,17 ,25 ,13 ,8 ,9 ,24 ,22 ,28 ,4 ,1 ,21 ,29 ,15 ,12 ,10 ,18 ,18 ,27 ,8 ,1 ,12 ,19 ,19 ,22 ,6 ,19 ,21 ,20 ,5 ,2 ,12 ,29 ,20 ,16 ,15 ,24 ,17 ,20 ,4 ,1 ,15 ,25 ,11 ,12 ,9 ,22 ,15 ,20 ,10 ,2 ,8 ,19 ,28 ,28 ,24 ,28 ,28 ,27 ,4 ,2 ,15 ,27 ,21 ,15 ,11 ,24 ,23 ,23 ,4 ,1 ,16 ,25 ,25 ,17 ,4 ,28 ,19 ,23 ,4 ,2 ,9 ,23 ,22 ,16 ,12 ,21 ,15 ,22 ,5 ,2 ,13 ,24 ,24 ,24 ,22 ,25 ,26 ,26 ,5 ,1 ,11 ,25 ,21 ,27 ,16 ,23 ,20 ,21 ,4 ,1 ,9 ,23 ,15 ,10 ,14 ,17 ,11 ,17 ,6 ,1 ,15 ,22 ,22 ,20 ,13 ,27 ,17 ,27 ,4 ,2 ,9 ,32 ,18 ,17 ,13 ,18 ,16 ,16 ,4 ,2 ,15 ,22 ,23 ,20 ,10 ,23 ,21 ,26 ,4 ,1 ,14 ,18 ,20 ,16 ,12 ,18 ,18 ,17 ,4 ,1 ,8 ,19 ,23 ,16 ,13 ,28 ,17 ,24 ,4 ,2 ,11 ,23 ,24 ,22 ,16 ,28 ,21 ,23 ,4 ,2 ,14 ,24 ,19 ,19 ,18 ,22 ,18 ,20 ,6 ,1 ,14 ,19 ,16 ,11 ,10 ,23 ,16 ,10 ,4 ,1 ,12 ,16 ,18 ,11 ,12 ,22 ,13 ,21 ,5 ,1 ,15 ,23 ,28 ,28 ,9 ,28 ,28 ,25 ,4 ,1 ,11 ,17 ,18 ,12 ,7 ,23 ,25 ,28 ,4 ,1 ,11 ,17 ,21 ,22 ,16 ,26 ,24 ,25 ,5 ,2 ,9 ,28 ,15 ,15 ,12 ,20 ,15 ,20 ,10 ,2 ,8 ,24 ,18 ,19 ,15 ,20 ,21 ,20 ,10 ,1 ,13 ,21 ,24 ,12 ,15 ,28 ,11 ,27 ,4 ,1 ,12 ,14 ,23 ,18 ,8 ,28 ,27 ,26 ,4 ,1 ,24 ,21 ,20 ,21 ,14 ,22 ,23 ,19 ,4 ,2 ,11 ,20 ,20 ,21 ,13 ,21 ,21 ,26 ,8 ,1 ,11 ,25 ,24 ,15 ,18 ,21 ,16 ,20 ,4 ,2 ,16 ,20 ,17 ,12 ,11 ,19 ,20 ,22 ,14 ,1 ,12 ,17 ,26 ,25 ,12 ,21 ,21 ,19 ,4 ,2 ,18 ,26 ,18 ,12 ,12 ,21 ,10 ,23 ,5 ,2 ,12 ,17 ,26 ,25 ,24 ,28 ,18 ,28 ,4 ,2 ,14 ,17 ,21 ,17 ,11 ,23 ,20 ,22 ,8 ,2 ,16 ,24 ,20 ,26 ,5 ,27 ,21 ,27 ,4 ,2 ,24 ,30 ,25 ,24 ,17 ,23 ,24 ,14 ,4 ,1 ,13 ,25 ,9 ,18 ,9 ,23 ,26 ,25 ,5 ,1 ,11 ,15 ,23 ,20 ,20 ,23 ,23 ,22 ,8 ,1 ,14 ,25 ,20 ,17 ,17 ,26 ,22 ,24 ,7 ,1 ,16 ,18 ,19 ,11 ,14 ,23 ,13 ,23 ,4 ,1 ,12 ,20 ,26 ,27 ,23 ,27 ,27 ,25 ,4 ,1 ,21 ,32 ,13 ,14 ,10 ,20 ,24 ,28 ,9 ,2 ,11 ,14 ,21 ,22 ,19 ,28 ,19 ,28 ,4 ,1 ,6 ,20 ,14 ,19 ,5 ,19 ,17 ,16 ,4 ,2 ,9 ,25 ,26 ,19 ,16 ,24 ,16 ,25 ,5 ,1 ,14 ,25 ,23 ,18 ,19 ,26 ,20 ,21 ,4 ,1 ,16 ,25 ,19 ,9 ,5 ,20 ,8 ,27 ,4 ,1 ,18 ,35 ,25 ,22 ,15 ,25 ,16 ,21 ,6 ,2 ,9 ,29 ,21 ,17 ,18 ,25 ,17 ,22 ,6 ,1 ,13 ,25 ,24 ,23 ,20 ,27 ,23 ,26 ,4 ,2 ,17 ,21 ,20 ,16 ,17 ,22 ,18 ,21 ,6 ,1 ,11 ,21 ,22 ,23 ,19 ,25 ,24 ,24 ,4 ,1 ,16 ,24 ,20 ,13 ,11 ,26 ,17 ,24 ,6 ,1 ,11 ,26 ,23 ,21 ,12 ,21 ,20 ,23 ,4 ,1 ,11 ,24 ,21 ,17 ,13 ,23 ,22 ,26 ,8 ,2 ,11 ,20 ,16 ,15 ,7 ,24 ,22 ,21 ,5 ,1 ,20 ,24 ,20 ,16 ,8 ,24 ,20 ,24 ,8 ,1 ,10 ,18 ,16 ,19 ,15 ,20 ,18 ,23 ,7 ,1 ,12 ,17 ,25 ,19 ,13 ,22 ,21 ,21 ,4 ,2 ,11 ,22 ,18 ,16 ,18 ,25 ,23 ,20 ,6 ,1 ,14 ,22 ,25 ,23 ,19 ,27 ,28 ,22 ,4 ,1 ,12 ,22 ,21 ,19 ,12 ,22 ,19 ,26 ,5 ,1 ,12 ,24 ,18 ,17 ,12 ,20 ,22 ,23 ,6 ,1 ,12 ,32 ,21 ,20 ,17 ,24 ,17 ,23 ,4 ,2 ,10 ,19 ,22 ,25 ,17 ,25 ,25 ,22 ,4 ,2 ,12 ,21 ,22 ,22 ,11 ,28 ,22 ,25 ,4 ,2 ,10 ,23 ,19 ,18 ,11 ,20 ,21 ,21 ,8 ,2 ,10 ,18 ,18 ,16 ,17 ,22 ,15 ,21 ,9 ,1 ,13 ,19 ,24 ,18 ,5 ,17 ,20 ,25 ,4 ,1 ,12 ,22 ,23 ,15 ,8 ,20 ,25 ,26 ,12 ,2 ,13 ,27 ,22 ,19 ,17 ,23 ,21 ,21 ,4 ,1 ,9 ,21 ,19 ,23 ,18 ,22 ,24 ,24 ,8 ,1 ,14 ,20 ,17 ,20 ,17 ,22 ,23 ,21 ,8 ,2 ,14 ,21 ,22 ,24 ,17 ,23 ,22 ,23 ,4 ,1 ,12 ,20 ,24 ,17 ,10 ,25 ,14 ,24 ,4 ,1 ,18 ,29 ,24 ,20 ,8 ,28 ,11 ,24 ,4 ,1 ,17 ,30 ,20 ,11 ,9 ,24 ,22 ,24 ,15 ,1 ,12 ,10 ,19 ,20 ,13 ,25 ,22 ,25 ,3 ,1 ,15 ,23 ,19 ,8 ,14 ,25 ,6 ,28 ,8 ,1 ,8 ,29 ,20 ,22 ,5 ,21 ,15 ,18 ,4 ,2 ,8 ,19 ,22 ,20 ,16 ,25 ,26 ,28 ,5 ,1 ,12 ,26 ,25 ,23 ,22 ,23 ,26 ,22 ,4 ,1 ,10 ,22 ,21 ,11 ,15 ,20 ,20 ,28 ,3 ,1 ,18 ,26 ,21 ,22 ,14 ,26 ,26 ,22 ,11 ,1 ,15 ,27 ,18 ,10 ,8 ,21 ,15 ,24 ,6 ,1 ,16 ,19 ,17 ,19 ,10 ,24 ,25 ,27 ,4 ,2 ,11 ,24 ,25 ,26 ,18 ,24 ,22 ,21 ,5 ,2 ,10 ,26 ,23 ,22 ,18 ,25 ,20 ,26 ,4 ,2 ,7 ,22 ,15 ,12 ,9 ,20 ,18 ,24 ,16 ,1 ,17 ,23 ,22 ,13 ,15 ,25 ,23 ,25 ,8 ,1 ,7 ,25 ,20 ,19 ,9 ,11 ,22 ,20 ,4 ,2 ,14 ,19 ,23 ,19 ,15 ,24 ,23 ,21 ,4 ,1 ,12 ,20 ,26 ,21 ,21 ,23 ,17 ,23 ,4 ,1 ,15 ,25 ,16 ,11 ,9 ,24 ,20 ,23 ,5 ,1 ,13 ,14 ,22 ,21 ,16 ,24 ,21 ,19 ,8 ,2 ,10 ,19 ,22 ,25 ,15 ,26 ,23 ,22 ,4 ,1 ,16 ,27 ,25 ,27 ,10 ,27 ,25 ,15 ,4 ,2 ,11 ,21 ,14 ,21 ,4 ,21 ,25 ,24 ,4 ,2 ,7 ,21 ,18 ,14 ,12 ,20 ,21 ,18 ,8 ,2 ,15 ,14 ,16 ,16 ,14 ,18 ,22 ,18 ,8 ,1 ,18 ,21 ,22 ,16 ,14 ,23 ,18 ,23 ,4 ,1 ,11 ,23 ,17 ,19 ,18 ,20 ,18 ,17 ,18 ,1 ,13 ,18 ,27 ,24 ,19 ,24 ,18 ,19 ,4 ,2 ,11 ,20 ,21 ,18 ,16 ,20 ,21 ,21 ,5 ,2 ,13 ,19 ,15 ,16 ,7 ,21 ,21 ,12 ,4 ,2 ,12 ,15 ,24 ,20 ,12 ,28 ,25 ,25 ,4 ,2 ,11 ,23 ,22 ,19 ,18 ,24 ,24 ,25 ,4 ,1 ,11 ,26 ,16 ,20 ,13 ,25 ,24 ,24 ,7 ,1 ,13 ,21 ,25 ,27 ,21 ,23 ,28 ,24 ,4 ,2 ,8 ,13 ,24 ,24 ,24 ,24 ,24 ,24 ,6 ,2 ,12 ,24 ,23 ,23 ,17 ,22 ,22 ,22 ,4 ,2 ,9 ,17 ,20 ,20 ,12 ,25 ,22 ,22 ,4 ,1 ,14 ,21 ,18 ,20 ,12 ,20 ,20 ,21 ,6 ,1 ,18 ,28 ,22 ,20 ,10 ,24 ,25 ,23 ,5 ,1 ,15 ,22 ,18 ,15 ,14 ,19 ,13 ,21 ,4 ,1 ,9 ,18 ,20 ,17 ,14 ,25 ,21 ,24 ,8 ,1 ,11 ,27 ,22 ,16 ,13 ,25 ,23 ,22 ,6 ,1 ,17 ,25 ,23 ,20 ,17 ,26 ,18 ,25 ,5 ,2 ,12 ,21) ,dim=c(10 ,148) ,dimnames=list(c('I/ToKnow' ,'I/Accomp.' ,'I/Exp.Stimulation' ,'E/Identified' ,'E/Introjected' ,'E/Ext.Regulation' ,'Amotivation' ,'gender' ,'PE' ,'PS') ,1:148)) y <- array(NA,dim=c(10,148),dimnames=list(c('I/ToKnow','I/Accomp.','I/Exp.Stimulation','E/Identified','E/Introjected','E/Ext.Regulation','Amotivation','gender','PE','PS'),1:148)) for (i in 1:dim(x)[1]) { for (j in 1:dim(x)[2]) { y[i,j] <- as.numeric(x[i,j]) } } par4 = 'yes' par3 = '2' par2 = 'quantiles' 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) library(Hmisc) par1 <- as.numeric(par1) par3 <- as.numeric(par3) x <- data.frame(t(y)) is.data.frame(x) x <- x[!is.na(x[,par1]),] k <- length(x[1,]) n <- length(x[,1]) colnames(x)[par1] x[,par1] 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]) colnames(x) colnames(x)[par1] x[,par1] 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/11zij1292003332.tab") } } m postscript(file="/var/www/html/rcomp/tmp/2cqzn1292003332.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) plot(m) dev.off() postscript(file="/var/www/html/rcomp/tmp/3cqzn1292003332.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() 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) } 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/4mhgp1292003332.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() 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/580ed1292003332.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/6b0dj1292003332.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/7xju71292003332.tab") } try(system("convert tmp/2cqzn1292003332.ps tmp/2cqzn1292003332.png",intern=TRUE)) try(system("convert tmp/3cqzn1292003332.ps tmp/3cqzn1292003332.png",intern=TRUE)) try(system("convert tmp/4mhgp1292003332.ps tmp/4mhgp1292003332.png",intern=TRUE))