y <- c(2.00,2.00,2.00,2.00,2.00,2.00,2.00,2.00,2.00,2.00,2.00,2.00,2.00,2.21,2.25,2.25,2.45,2.50,2.50,2.64,2.75,2.93,3.00,3.17,3.25,3.39,3.50,3.50,3.65,3.75,3.75,3.90,4.00,4.00,4.00,4.00,4.00,4.00,4.00,4.00,4.00,4.00,4.00,4.00,4.18,4.25,4.25,3.97,3.42,2.75,2.31,2.00,1.66,1.31,1.09,1.00,1.00,1.00,1.00,1.00) x <- c(4.6,4.3,4.3,4.6,5.1,4.8,4.5,4.9,5.1,5.1,5.2,4.5,4.6,4.9,4.6,4.4,3.7,4,4.2,3.9,3.6,3.6,3.2,3.2,3.5,3.6,3.7,3.8,3.8,3.8,3.3,3.3,3.4,3.1,3.5,4.2,4.9,5.1,5.5,5.6,6.4,6.2,7.2,7.8,7.9,7.4,7.5,6.7,5.1,4.6,4.3,3.9,2.6,2.6,1.6,0.9,0.3,1.2,0.8,1) par5 = 'rente' par4 = 'index' par3 = 'TRUE' par2 = 'blue' par1 = 'red' ylab = '' main = 'Back to Back Histogram' #'GNU S' R Code compiled by R2WASP v. 1.0.44 () #Author: Prof. Dr. P. Wessa #To cite this work: AUTHOR(S), (YEAR), 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: Office for Research, Development, and Education #Technical description: Write here your technical program description (don't use hard returns!) if (par3 == 'TRUE') par3 <- TRUE if (par3 == 'FALSE') par3 <- FALSE library(Hmisc) z <- data.frame(cbind(x,y)) names(z) <- list(par4,par5) postscript(file="/var/www/html/rcomp/tmp/1c65d1258369830.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 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() system("convert tmp/1c65d1258369830.ps tmp/1c65d1258369830.png")