R version 2.7.0 (2008-04-22) 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(29,34,34,30,21,23,18,27,20,25,24,23,18,22,22,18,19,18,22,22,20,14,15,14,19,23,23,17,21,17,24,18,17,18,16,15,21,17,11,13,11,13,13,14,18,13,16,13,15,20,17,18,17,20,18,18,21,19,17,20,19,24,21,21,16,18,17,16,11,14,11,20,22,18,17,18,15,19,18,15,16,17,12,16,14,18,17,10,18,16,13,13,13,13,5,8,5,4,4,1,6,4,5,4,11) > x <- c(93.9,97.5,98.1,89.6,98.4,102,99.2,101.8,108.3,106.7,108.2,94.2,95.1,98.1,93.2,94,97.2,95,90.5,91.6,90.5,79.9,74.9,74.3,75.9,77.7,86.9,90.7,91,89.5,92.5,94.1,98.5,96.8,91.2,97.1,104.9,110.9,104.8,94.1,95.8,99.3,101.1,104,99,105.4,107.1,110.7,117.1,118.7,126.5,127.5,134.6,131.8,135.9,142.7,141.7,153.4,145,137.7,148.3,152.2,169.4,168.6,161.1,174.1,179,190.6,190,181.6,174.8,180.5,196.8,193.8,197,216.3,221.4,217.9,229.7,227.4,204.2,196.6,198.8,207.5,190.7,201.6,210.5,223.5,223.8,231.2,244,234.7,250.2,265.7,287.6,283.3,295.4,312.3,333.8,347.7,383.2,407.1,413.6,362.7) > par5 = '' > par4 = '' > par3 = 'TRUE' > par2 = 'black' > par1 = 'red' > ylab = '' > main = 'vergelijking reeks 4 met reeks 3' > #'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)) Warning message: In cbind(x, y) : number of rows of result is not a multiple of vector length (arg 1) > names(z) <- list(par4,par5) > postscript(file="/var/www/html/rcomp/tmp/19k5b1224447672.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/19k5b1224447672.ps tmp/19k5b1224447672.png") > > > proc.time() user system elapsed 1.776 0.294 1.842