R version 3.0.1 (2013-05-16) -- "Good Sport" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: i686-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(784437 + ,792318 + ,792318 + ,801087 + ,816849 + ,833499 + ,850038 + ,850038 + ,858030 + ,865911 + ,865911 + ,865911 + ,874569 + ,914862 + ,988455 + ,997113 + ,1046175 + ,1086468 + ,1095237 + ,1135530 + ,1144188 + ,1144188 + ,1144188 + ,1168719 + ,1193250 + ,1193250 + ,1242312 + ,1266843 + ,1348206 + ,1421799 + ,1421799 + ,1429680 + ,1446330 + ,1470861 + ,1470861 + ,1487400 + ,1503273 + ,1519812 + ,1519812 + ,1568874 + ,1683537 + ,1691418 + ,1732599 + ,1757130 + ,1757130 + ,1772892 + ,1781550 + ,1789542 + ,1806081 + ,1813962 + ,1813962 + ,1838493 + ,2051280 + ,2067042 + ,2108223 + ,2108223 + ,2157174 + ,2157174 + ,2181705 + ,2189586 + ,2222886 + ,2247306 + ,2336661 + ,2394381 + ,2565987 + ,2615049 + ,2655342 + ,2664111 + ,2679873 + ,2713062 + ,2721054 + ,2762124 + ,2762124 + ,2811186 + ,2966142 + ,3064266 + ,3235761 + ,3292704 + ,3309354 + ,3350424 + ,3358416 + ,3366297 + ,3382836 + ,3399486 + ,3415359 + ,3448548 + ,3628035 + ,3701517 + ,3930066 + ,3995778 + ,4012317 + ,4036848 + ,4061379 + ,4069260 + ,4077141 + ,4077141 + ,4085910 + ,4158615 + ,4314459 + ,4387941 + ,4706622 + ,4739034 + ,4755684 + ,4771446 + ,4771446 + ,4812627 + ,4853697 + ,4853697 + ,4861578 + ,4886109 + ,5033184 + ,5082246 + ,5408808 + ,5466528 + ,5491059 + ,5548002 + ,5621595 + ,5629476 + ,5638134 + ,5638134) + ,dim=c(1 + ,120) + ,dimnames=list(c('Omzet') + ,1:120)) > y <- array(NA,dim=c(1,120),dimnames=list(c('Omzet'),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' > par1 <- 'grey' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa P., (2012), Notched Boxplots (v1.0.6) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_notchedbox1.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > z <- as.data.frame(t(y)) > postscript(file="/var/fisher/rcomp/tmp/1rcv31377062661.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] [1,] 784437 [2,] 1421799 [3,] 2206236 [4,] 4004048 [5,] 5638134 $n [1] 120 $conf [,1] [1,] 1833789 [2,] 2578683 $out numeric(0) $group numeric(0) $names [1] "Omzet" > dev.off() null device 1 > > #Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/fisher/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/fisher/rcomp/tmp/2jkc01377062661.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/fisher/rcomp/tmp/3gohd1377062661.tab") > > try(system("convert tmp/1rcv31377062661.ps tmp/1rcv31377062661.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.972 0.266 1.214