R version 2.8.1 (2008-12-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. > x <- array(list(363 + ,14.3 + ,0.8833 + ,15.89 + ,364 + ,14.2 + ,0.87 + ,16.93 + ,363 + ,15.9 + ,0.8758 + ,20.28 + ,358 + ,15.3 + ,0.8858 + ,22.52 + ,357 + ,15.5 + ,0.917 + ,23.51 + ,357 + ,15.1 + ,0.9554 + ,22.59 + ,380 + ,15 + ,0.9922 + ,23.51 + ,378 + ,12.1 + ,0.9778 + ,24.76 + ,376 + ,15.8 + ,0.9808 + ,26.08 + ,380 + ,16.9 + ,0.9811 + ,25.29 + ,379 + ,15.1 + ,1.0014 + ,23.38 + ,384 + ,13.7 + ,1.0183 + ,25.29 + ,392 + ,14.8 + ,1.0622 + ,28.42 + ,394 + ,14.7 + ,1.0773 + ,31.85 + ,392 + ,16 + ,1.0807 + ,30.1 + ,396 + ,15.4 + ,1.0848 + ,25.45 + ,392 + ,15 + ,1.1582 + ,24.95 + ,396 + ,15.5 + ,1.1663 + ,26.84 + ,419 + ,15.1 + ,1.1372 + ,27.52 + ,421 + ,11.7 + ,1.1139 + ,27.94 + ,420 + ,16.3 + ,1.1222 + ,25.23 + ,418 + ,16.7 + ,1.1692 + ,26.53 + ,410 + ,15 + ,1.1702 + ,27.21 + ,418 + ,14.9 + ,1.2286 + ,28.53 + ,426 + ,14.6 + ,1.2613 + ,30.35 + ,428 + ,15.3 + ,1.2646 + ,31.21 + ,430 + ,17.9 + ,1.2262 + ,32.86 + ,424 + ,16.4 + ,1.1985 + ,33.2 + ,423 + ,15.4 + ,1.2007 + ,35.73 + ,427 + ,17.9 + ,1.2138 + ,34.53 + ,441 + ,15.9 + ,1.2266 + ,36.54 + ,449 + ,13.9 + ,1.2176 + ,40.1 + ,452 + ,17.8 + ,1.2218 + ,40.56 + ,462 + ,17.9 + ,1.249 + ,46.14 + ,455 + ,17.4 + ,1.2991 + ,42.85 + ,461 + ,16.7 + ,1.3408 + ,38.22 + ,461 + ,16 + ,1.3119 + ,40.18 + ,463 + ,16.6 + ,1.3014 + ,42.19 + ,462 + ,19.1 + ,1.3201 + ,47.56 + ,456 + ,17.8 + ,1.2938 + ,47.26 + ,455 + ,17.2 + ,1.2694 + ,44.03 + ,456 + ,18.6 + ,1.2165 + ,49.83 + ,472 + ,16.3 + ,1.2037 + ,53.35 + ,472 + ,15.1 + ,1.2292 + ,58.9 + ,471 + ,19.2 + ,1.2256 + ,59.64 + ,465 + ,17.7 + ,1.2015 + ,56.99 + ,459 + ,19.1 + ,1.1786 + ,53.2 + ,465 + ,18 + ,1.1856 + ,53.24 + ,468 + ,17.5 + ,1.2103 + ,57.85 + ,467 + ,17.8 + ,1.1938 + ,55.69 + ,463 + ,21.1 + ,1.202 + ,55.64 + ,460 + ,17.2 + ,1.2271 + ,62.52 + ,462 + ,19.4 + ,1.277 + ,64.4 + ,461 + ,19.8 + ,1.265 + ,64.65 + ,476 + ,17.6 + ,1.2684 + ,67.71 + ,476 + ,16.2 + ,1.2811 + ,67.21 + ,471 + ,19.5 + ,1.2727 + ,59.37 + ,453 + ,19.9 + ,1.2611 + ,53.26 + ,443 + ,20 + ,1.2881 + ,52.42 + ,442 + ,17.3 + ,1.3213 + ,55.03 + ,363 + ,14.3 + ,0.8833 + ,15.89 + ,364 + ,14.2 + ,0.87 + ,16.93 + ,363 + ,15.9 + ,0.8758 + ,20.28 + ,358 + ,15.3 + ,0.8858 + ,22.52 + ,357 + ,15.5 + ,0.917 + ,23.51 + ,357 + ,15.1 + ,0.9554 + ,22.59 + ,380 + ,15 + ,0.9922 + ,23.51 + ,378 + ,12.1 + ,0.9778 + ,24.76 + ,376 + ,15.8 + ,0.9808 + ,26.08 + ,380 + ,16.9 + ,0.9811 + ,25.29 + ,379 + ,15.1 + ,1.0014 + ,23.38 + ,384 + ,13.7 + ,1.0183 + ,25.29 + ,392 + ,14.8 + ,1.0622 + ,28.42 + ,394 + ,14.7 + ,1.0773 + ,31.85 + ,392 + ,16 + ,1.0807 + ,30.1 + ,396 + ,15.4 + ,1.0848 + ,25.45 + ,392 + ,15 + ,1.1582 + ,24.95 + ,396 + ,15.5 + ,1.1663 + ,26.84 + ,419 + ,15.1 + ,1.1372 + ,27.52 + ,421 + ,11.7 + ,1.1139 + ,27.94 + ,420 + ,16.3 + ,1.1222 + ,25.23 + ,418 + ,16.7 + ,1.1692 + ,26.53 + ,410 + ,15 + ,1.1702 + ,27.21 + ,418 + ,14.9 + ,1.2286 + ,28.53 + ,426 + ,14.6 + ,1.2613 + ,30.35 + ,428 + ,15.3 + ,1.2646 + ,31.21 + ,430 + ,17.9 + ,1.2262 + ,32.86 + ,424 + ,16.4 + ,1.1985 + ,33.2 + ,423 + ,15.4 + ,1.2007 + ,35.73 + ,427 + ,17.9 + ,1.2138 + ,34.53 + ,441 + ,15.9 + ,1.2266 + ,36.54 + ,449 + ,13.9 + ,1.2176 + ,40.1 + ,452 + ,17.8 + ,1.2218 + ,40.56 + ,462 + ,17.9 + ,1.249 + ,46.14 + ,455 + ,17.4 + ,1.2991 + ,42.85 + ,461 + ,16.7 + ,1.3408 + ,38.22 + ,461 + ,16 + ,1.3119 + ,40.18 + ,463 + ,16.6 + ,1.3014 + ,42.19 + ,462 + ,19.1 + ,1.3201 + ,47.56 + ,456 + ,17.8 + ,1.2938 + ,47.26 + ,455 + ,17.2 + ,1.2694 + ,44.03 + ,456 + ,18.6 + ,1.2165 + ,49.83 + ,472 + ,16.3 + ,1.2037 + ,53.35 + ,472 + ,15.1 + ,1.2292 + ,58.9 + ,471 + ,19.2 + ,1.2256 + ,59.64 + ,465 + ,17.7 + ,1.2015 + ,56.99 + ,459 + ,19.1 + ,1.1786 + ,53.2 + ,465 + ,18 + ,1.1856 + ,53.24 + ,468 + ,17.5 + ,1.2103 + ,57.85 + ,467 + ,17.8 + ,1.1938 + ,55.69 + ,463 + ,21.1 + ,1.202 + ,55.64 + ,460 + ,17.2 + ,1.2271 + ,62.52 + ,462 + ,19.4 + ,1.277 + ,64.4 + ,461 + ,19.8 + ,1.265 + ,64.65 + ,476 + ,17.6 + ,1.2684 + ,67.71 + ,476 + ,16.2 + ,1.2811 + ,67.21 + ,471 + ,19.5 + ,1.2727 + ,59.37 + ,453 + ,19.9 + ,1.2611 + ,53.26 + ,443 + ,20 + ,1.2881 + ,52.42 + ,442 + ,17.3 + ,1.3213 + ,55.03) + ,dim=c(4 + ,120) + ,dimnames=list(c('wk>25jB' + ,'ExpBE' + ,'Wssk€-$' + ,'RuwOlB/$') + ,1:120)) > y <- array(NA,dim=c(4,120),dimnames=list(c('wk>25jB','ExpBE','Wssk€-$','RuwOlB/$'),1:120)) > 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/rcomp/tmp/1colv1257595894.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] [,4] [1,] 357.0 11.7 0.87000 15.890 [2,] 393.0 15.1 1.08275 26.305 [3,] 435.5 16.3 1.20285 36.135 [4,] 461.5 17.8 1.26295 53.250 [5,] 476.0 21.1 1.34080 67.710 $n [1] 120 120 120 120 $conf [,1] [,2] [,3] [,4] [1,] 425.62 15.91057 1.176859 32.24863 [2,] 445.38 16.68943 1.228841 40.02137 $out numeric(0) $group numeric(0) $names [1] "wk>25jB" "ExpBE" "Wssk\200-$" "RuwOlB/$" > dev.off() null device 1 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/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/rcomp/tmp/2mb5c1257595894.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/rcomp/tmp/35cb51257595894.tab") > > system("convert tmp/1colv1257595894.ps tmp/1colv1257595894.png") > > > proc.time() user system elapsed 0.410 0.310 0.617