R version 2.7.2 (2008-08-25) 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. 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. > y <- c(680.2,668.9,611.4,640.8,549.8,541.8,628.6,686.5,611.8,588.4,566.9,563.7,569.9,635.4,590.8,634.3,576.1,351.6,507.5,586.2,666.4,693.6,650.6,654.8,733.5,648.1,678.1,816.2,591,563.5,742.5,694.4,728.6,749,538.9,568.5,692.8,580.5,506.9,612.8,442.9,523.3,596.7,533.7,523.1,559.2,430.7,538.2,612.4,428,522.4,531.1,425.9,410.3,551,555.6,460.2,288.9,392.3,400.5,399,354.9,337.6,379.2,334.1,321.6,449.8,486.3,421.9,405.6,420,432.4,418.1) > x <- c(12710.3,12120.8,12469.5,12054.6,12112.9,9617.2,12645.8,13581.3,12162.3,10969.7,11880,11887.6,12926.9,12300,12092.8,12380.8,12196.9,9455,13168,13427.9,11980.5,11884.8,11691.7,12233.8,14341.4,13130.7,12421.1,14285.8,12864.6,11160.2,14316.2,14388.7,14013.9,13419,12769.6,13315.5,15332.9,14243,13824.4,14962.9,13202.9,12199,15508.9,14199.8,15169.6,14058,13786.2,14147.9,16541.7,13587.5,15582.4,15802.8,14130.5,12923.2,15612.2,16033.7,16036.6,14037.8,15330.6,15038.3,17401.8,14992.5,16043.7,16929.6,15921.3,14417.2,15961,17851.9,16483.9,14215.5,17429.7,17839.5,17629.2) > par5 = 'Uitvoer Brusselse gewest' > par4 = 'Uitvoer Vlaamse gewest' > par3 = 'TRUE' > par2 = 'brown' > 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) Attaching package: 'Hmisc' The following object(s) are masked from package:base : format.pval, round.POSIXt, trunc.POSIXt, units > z <- data.frame(cbind(x,y)) > names(z) <- list(par4,par5) > postscript(file="/var/www/html/rcomp/tmp/1e06w1224541001.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() null device 1 > > system("convert tmp/1e06w1224541001.ps tmp/1e06w1224541001.png") > > > proc.time() user system elapsed 0.938 0.178 1.023