R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(588264,462824,431664,336150,641250,862205,1260562,1523124,1680746,1523124,577918,399492,325803,263595,493873,735641,1071790,1344597,1576009,1344597,567562,347630,378790,294755,452367,651596,998102,1376780,1680746,1376780,546859,347630,357986,305101,409838,714939,1082137,1397583,1733732,1397583,756344,546859,536412,577918,694135,987746,1282389,1534503,1796963,1534503,745987,567562,493873,577918,714939,1071790,1460814,1712929,1880998,1712929,588264,409838,325803,420195,661953,1051088,1418286,1586355,1796963,1586355,483527,231412,200263,399492,714939,1092483,1418286,1649698,1848826,1649698,493873,325803,315447,462824,704481,1082137,1439089,1596712,1785585,1596712,493873,325803,347630,431664,620447,977399,1366423,1565653,1775239,1565653,504229,399492,378790,515709,714939,1155825,1555307,1807421,2037699,1807421,526055,442021,420195,620447,819676,1198354,1555307,1754435,2059525,1754435),dim=c(10,12),dimnames=list(c('1','2','3','4','5','6','7','8','9','10'),1:12)) > y <- array(NA,dim=c(10,12),dimnames=list(c('1','2','3','4','5','6','7','8','9','10'),1:12)) > 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 = 'Block Index' > main = 'Notched Boxplots Omzetcijfers Product X' > #'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/wessaorg/rcomp/tmp/146xk1313154188.ps",horizontal=F,onefile=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] [,5] [,6] [,7] [,8] [,9] [1,] 483527.0 231412.0 200263.0 263595.0 409838 651596 998102 1344597 1576009 [2,] 499051.0 336716.5 325803.0 320625.5 557160 798923 1171350 1460354 1707239 [3,] 557210.5 399492.0 368388.0 425929.5 678044 1019417 1392354 1576004 1791274 [4,] 588264.0 452422.5 425929.5 546813.5 714939 1087310 1449952 1681314 1864912 [5,] 588264.0 567562.0 536412.0 620447.0 819676 1198354 1555307 1807421 2059525 [,10] [1,] 1344597 [2,] 1460354 [3,] 1576004 [4,] 1681314 [5,] 1807421 $n [1] 12 12 12 12 12 12 12 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 516519.9 346717.7 322719.6 322763.6 606079.9 887881.8 1265282 1475223 [2,] 597901.1 452266.3 414056.4 529095.4 750008.1 1150952.2 1519427 1676785 [,9] [,10] [1,] 1719358 1475223 [2,] 1863190 1676785 $out [1] 756344 745987 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" Warning message: In bxp(list(stats = c(483527, 499051, 557210.5, 588264, 588264, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/2rtrc1313154188.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/wessaorg/rcomp/tmp/3szmv1313154188.tab") > > try(system("convert tmp/146xk1313154188.ps tmp/146xk1313154188.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.483 0.080 0.558