if (par3 == 'TRUE') par3 <- TRUE if (par3 == 'FALSE') par3 <- FALSE library(Hmisc) load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Summary statistics (variable X)',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('',RC.texteval('summary(x)'),' ',sep='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Summary statistics (variable Y)',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('',RC.texteval('summary(y)'),' ',sep='')) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab') summary(x) summary(y) z <- data.frame(cbind(x,y)) names(z) <- list(par4,par5) bitmap(file='plot.png') out <- histbackback(z, probability=par3, main = main, ylab = ylab) barplot(-out$left, col=par1, horiz=TRUE, space=0, add=TRUE, axes=FALSE) barplot(out$right, col=par2, horiz=TRUE, space=0, add=TRUE, axes=FALSE) dev.off()
|