R version 3.2.3 (2015-12-10) -- "Wooden Christmas-Tree" Copyright (C) 2015 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-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(13408 + ,7820 + ,9079 + ,8307 + ,7865 + ,10028 + ,9054 + ,7143 + ,8006 + ,7638 + ,7600 + ,2904 + ,13224 + ,8079 + ,9678 + ,7746 + ,9007 + ,8362 + ,7458 + ,7753 + ,7352 + ,7117 + ,6971 + ,3304 + ,11812 + ,6867 + ,8296 + ,6489 + ,7784 + ,7506 + ,6514 + ,6323 + ,6201 + ,7169 + ,6744 + ,2087 + ,10668 + ,6406 + ,7730 + ,7105 + ,7694 + ,7160 + ,6820 + ,6025 + ,5877 + ,7191 + ,5778 + ,2273 + ,11321 + ,6759 + ,7150 + ,6363 + ,6442 + ,6453 + ,6228 + ,5325 + ,6504 + ,6817 + ,5789 + ,1894 + ,11068 + ,7174 + ,8269 + ,7060 + ,6681 + ,8953 + ,7815 + ,5925 + ,6805 + ,7044 + ,7169 + ,2824 + ,10717 + ,5245 + ,6237 + ,5871 + ,5508 + ,15801 + ,1236 + ,2656 + ,3425 + ,3533 + ,4287 + ,1380 + ,8584 + ,5522 + ,6423 + ,5173 + ,5583 + ,5716 + ,4752 + ,4977 + ,4999 + ,5285 + ,5747 + ,1713 + ,9923 + ,6737 + ,7433 + ,6388 + ,6855 + ,7658 + ,6585 + ,6847 + ,6353 + ,7361 + ,6929 + ,1714 + ,11798 + ,8378 + ,8131 + ,7676 + ,7505 + ,8168 + ,6455 + ,6141 + ,6554 + ,6888 + ,5339 + ,1624 + ,9187 + ,5047 + ,5289 + ,4169 + ,3862 + ,4253 + ,3768 + ,3066 + ,4108 + ,3890 + ,3420 + ,1221 + ,5984 + ,4064 + ,5151 + ,4027 + ,3530 + ,4819 + ,3855 + ,3584 + ,4322 + ,4154 + ,4656 + ,1464 + ,7780 + ,5060 + ,6084 + ,4778 + ,4989 + ,4903 + ,4142 + ,4101 + ,4595 + ,5034 + ,5407 + ,1782 + ,8395 + ,5291 + ,6116 + ,4210 + ,4621 + ,5299 + ,4293 + ,4542 + ,3831 + ,4360 + ,4088 + ,1508 + ,6743 + ,4159 + ,5105 + ,4283 + ,4019 + ,4206 + ,3948 + ,3407 + ,3701 + ,4159 + ,4208 + ,2622 + ,6229 + ,4432 + ,4986 + ,4226 + ,4349 + ,4688 + ,4002 + ,3381 + ,4250 + ,4154 + ,4350 + ,2713) + ,dim=c(1 + ,192) + ,dimnames=list(c('V1') + ,1:192)) > y <- array(NA,dim=c(1,192),dimnames=list(c('V1'),1:192)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par2 = 'no' > par1 = 'grey' > ylab = 'value' > xlab = 'variables' > main = 'Notched Boxplots' > par2 <- 'no' > par1 <- 'grey' > #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Mon, 28 Dec 2015 11:07:50 +0000) > #Author: root > #To cite this work: Wessa P., (2015), Notched Boxplots (v1.0.7) 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 > # > if(par2=='yes') { + z <- na.omit(as.data.frame(t(y))) + } else { + z <- as.data.frame(t(y)) + } > postscript(file="/var/wessaorg/rcomp/tmp/19v9h1451838344.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,] 1221.0 [2,] 4209.0 [3,] 5901.0 [4,] 7182.5 [5,] 11321.0 $n [1] 192 $conf [,1] [1,] 5561.942 [2,] 6240.058 $out [1] 13408 13224 11812 15801 11798 $group [1] 1 1 1 1 1 $names [1] "V1" > 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/2h14v1451838344.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/3hhjr1451838344.tab") > > try(system("convert tmp/19v9h1451838344.ps tmp/19v9h1451838344.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.518 0.102 0.623