R version 2.9.0 (2009-04-17) Copyright (C) 2009 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. > x <- array(list(262.262 + ,455 + ,261.631 + ,456 + ,258.953 + ,517 + ,259.966 + ,525 + ,262.850 + ,523 + ,262.204 + ,519 + ,263.418 + ,509 + ,262.752 + ,512 + ,266.433 + ,519 + ,267.722 + ,517 + ,266.003 + ,510 + ,262.971 + ,509 + ,265.521 + ,501 + ,264.676 + ,507 + ,270.223 + ,569 + ,269.508 + ,580 + ,268.457 + ,578 + ,265.814 + ,565 + ,266.680 + ,547 + ,263.018 + ,555 + ,269.285 + ,562 + ,269.829 + ,561 + ,270.911 + ,555 + ,266.844 + ,544 + ,271.244 + ,537 + ,269.907 + ,543 + ,271.296 + ,594 + ,270.157 + ,611 + ,271.322 + ,613 + ,267.179 + ,611 + ,264.101 + ,594 + ,265.518 + ,595 + ,269.419 + ,591 + ,268.714 + ,589 + ,272.482 + ,584 + ,268.351 + ,573 + ,268.175 + ,567 + ,270.674 + ,569 + ,272.764 + ,621 + ,272.599 + ,629 + ,270.333 + ,628 + ,270.846 + ,612 + ,270.491 + ,595 + ,269.160 + ,597 + ,274.027 + ,593 + ,273.784 + ,590 + ,276.663 + ,580 + ,274.525 + ,574 + ,271.344 + ,573 + ,271.115 + ,573 + ,270.798 + ,620 + ,273.911 + ,626 + ,273.985 + ,620 + ,271.917 + ,588 + ,273.338 + ,566 + ,270.601 + ,557 + ,273.547 + ,561 + ,275.363 + ,549 + ,281.229 + ,532 + ,277.793 + ,526 + ,279.913 + ,511 + ,282.500 + ,499 + ,280.041 + ,555 + ,282.166 + ,565 + ,290.304 + ,542 + ,283.519 + ,527 + ,287.816 + ,510 + ,285.226 + ,514 + ,287.595 + ,517 + ,289.741 + ,508 + ,289.148 + ,493 + ,288.301 + ,490 + ,290.155 + ,469 + ,289.648 + ,478 + ,288.225 + ,528 + ,289.351 + ,534 + ,294.735 + ,518 + ,305.333 + ,506 + ,309.030 + ,502 + ,310.215 + ,516 + ,321.935 + ,528 + ,325.734 + ,533 + ,320.846 + ,536 + ,323.023 + ,537 + ,319.753 + ,524 + ,321.753 + ,536 + ,320.757 + ,587 + ,324.479 + ,597 + ,324.641 + ,581 + ,322.767 + ,564) + ,dim=c(2 + ,90) + ,dimnames=list(c('staatsschuld' + ,'werkloosheid') + ,1:90)) > y <- array(NA,dim=c(2,90),dimnames=list(c('staatsschuld','werkloosheid'),1:90)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par1 = 'grey' > ylab = 'value' > xlab = 'variables' > main = 'Notched Boxplots' > #'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!) > z <- as.data.frame(t(y)) > postscript(file="/var/www/html/rcomp/tmp/10cz91262205262.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > (r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=TRUE,col=par1)) $stats [,1] [,2] [1,] 258.9530 455 [2,] 268.3510 517 [3,] 271.6305 552 [4,] 287.8160 584 [5,] 310.2150 629 $n [1] 90 90 $conf [,1] [,2] [1,] 268.3887 540.8414 [2,] 274.8723 563.1586 $out [1] 321.935 325.734 320.846 323.023 319.753 321.753 320.757 324.479 324.641 [10] 322.767 $group [1] 1 1 1 1 1 1 1 1 1 1 $names [1] "staatsschuld" "werkloosheid" > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/overview.htm','Boxplot statistics','Boxplot overview'),6,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/lower_whisker.htm','lower whisker','definition of lower whisker'),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/lower_hinge.htm','lower hinge','definition of lower hinge'),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/central_tendency.htm','median','definitions about measures of central tendency'),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/upper_hinge.htm','upper hinge','definition of upper hinge'),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/upper_whisker.htm','upper whisker','definition of upper whisker'),1,TRUE) > a<-table.row.end(a) > for (i in 1:length(y[,1])) + { + a<-table.row.start(a) + a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE) + for (j in 1:5) + { + a<-table.element(a,r$stats[j,i]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/2as3i1262205262.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Boxplot Notches',4,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',1,TRUE) > a<-table.element(a,'lower bound',1,TRUE) > a<-table.element(a,'median',1,TRUE) > a<-table.element(a,'upper bound',1,TRUE) > a<-table.row.end(a) > for (i in 1:length(y[,1])) + { + a<-table.row.start(a) + a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE) + a<-table.element(a,r$conf[1,i]) + a<-table.element(a,r$stats[3,i]) + a<-table.element(a,r$conf[2,i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/36m521262205262.tab") > > try(system("convert tmp/10cz91262205262.ps tmp/10cz91262205262.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.377 0.156 0.462