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(0.527 + ,0.7461 + ,0.0479 + ,0.472 + ,0.7775 + ,0.0595 + ,0.000 + ,0.779 + ,0.0546 + ,0.052 + ,0.7744 + ,0.0575 + ,0.313 + ,0.7905 + ,0.0515 + ,0.364 + ,0.7719 + ,0.0496 + ,0.363 + ,0.7811 + ,0.0528 + ,-0.155 + ,0.7557 + ,0.0573 + ,0.052 + ,0.7637 + ,0.0575 + ,0.568 + ,0.7595 + ,0.0588 + ,0.668 + ,0.7471 + ,0.063 + ,1.378 + ,0.7615 + ,0.0674 + ,0.252 + ,0.7487 + ,0.0675 + ,-0.402 + ,0.7389 + ,0.0734 + ,-0.05 + ,0.7337 + ,0.0664 + ,0.555 + ,0.751 + ,0.0662 + ,0.05 + ,0.7382 + ,0.0632 + ,0.15 + ,0.7159 + ,0.0532 + ,0.45 + ,0.7542 + ,0.0568 + ,0.299 + ,0.7636 + ,0.0618 + ,0.199 + ,0.7433 + ,0.0502 + ,0.496 + ,0.7658 + ,0.021 + ,0.444 + ,0.7627 + ,0.04 + ,-0.393 + ,0.748 + ,0.03 + ,-0.444 + ,0.7692 + ,0.0473 + ,0.198 + ,0.785 + ,0.0514 + ,0.494 + ,0.7913 + ,0.0581 + ,0.133 + ,0.772 + ,0.0624 + ,0.388 + ,0.788 + ,0.0449 + ,0.484 + ,0.807 + ,0.0422 + ,0.278 + ,0.8268 + ,0.0488 + ,0.369 + ,0.8244 + ,0.0518 + ,0.165 + ,0.8487 + ,0.0519 + ,0.155 + ,0.8572 + ,0.0506 + ,0.087 + ,0.8214 + ,0.0465 + ,0.414 + ,0.8827 + ,0.0483 + ,0.36 + ,0.9216 + ,0.046 + ,0.975 + ,0.8865 + ,0.0472 + ,0.27 + ,0.8816 + ,0.0433 + ,0.359 + ,0.8884 + ,0.0497 + ,0.169 + ,0.9466 + ,0.0537 + ,0.381 + ,0.918 + ,0.0419 + ,0.154 + ,0.9337 + ,0.0454 + ,0.486 + ,0.9559 + ,0.0582 + ,0.925 + ,0.9626 + ,0.0549 + ,0.728 + ,0.9434 + ,0.0328 + ,-0.014 + ,0.8639 + ,0.0511 + ,0.046 + ,0.7996 + ,0.0624 + ,-0.819 + ,0.668 + ,0.0641 + ,-1.674 + ,0.6572 + ,0.0643 + ,-0.788 + ,0.6928 + ,0.0842 + ,0.279 + ,0.6438 + ,0.0823 + ,0.396 + ,0.6454 + ,0.0317 + ,-0.141 + ,0.6873 + ,0.0272 + ,-0.019 + ,0.7265 + ,0.03 + ,0.099 + ,0.7912 + ,0.0347 + ,0.742 + ,0.8114 + ,0.0388 + ,0.005 + ,0.8281 + ,0.0343 + ,0.448 + ,0.8393 + ,0.0406) + ,dim=c(3 + ,59) + ,dimnames=list(c('Amerikaanse_inflatie' + ,'USdollar' + ,'Rentevoet_auto') + ,1:59)) > y <- array(NA,dim=c(3,59),dimnames=list(c('Amerikaanse_inflatie','USdollar','Rentevoet_auto'),1:59)) > 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/1mwj31260543800.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] [,3] [1,] -0.444 0.64380 0.02720 [2,] 0.051 0.74835 0.04515 [3,] 0.278 0.77750 0.05150 [4,] 0.449 0.83370 0.05915 [5,] 0.975 0.95590 0.07340 $n [1] 59 59 59 $conf [,1] [,2] [,3] [1,] 0.196132 0.7599436 0.04862022 [2,] 0.359868 0.7950564 0.05437978 $out [1] 1.3780 -0.8190 -1.6740 -0.7880 0.9626 0.0210 0.0842 0.0823 $group [1] 1 1 1 1 2 3 3 3 $names [1] "Amerikaanse_inflatie" "USdollar" "Rentevoet_auto" > 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/2dcer1260543800.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/33m971260543800.tab") > > system("convert tmp/1mwj31260543800.ps tmp/1mwj31260543800.png") > > > proc.time() user system elapsed 0.375 0.152 0.647