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(7612 + ,7381 + ,6978 + ,6819 + ,6688 + ,6454 + ,6679 + ,6921 + ,7807 + ,7898 + ,7832 + ,7384 + ,7620 + ,7281 + ,6929 + ,6587 + ,6071 + ,5928 + ,5964 + ,6374 + ,7160 + ,7213 + ,6890 + ,6525 + ,6739 + ,6580 + ,6391 + ,6254 + ,6114 + ,5978 + ,6315 + ,6427 + ,7132 + ,7292 + ,7708 + ,7525 + ,7450 + ,7526 + ,7263 + ,7070 + ,6893 + ,6781 + ,7188 + ,7015 + ,8273 + ,8470 + ,8230 + ,8137 + ,8122 + ,8367 + ,8141 + ,7750 + ,7504 + ,7330 + ,7608 + ,7647 + ,8942 + ,8865 + ,8320 + ,8207 + ,8105 + ,8290 + ,8162 + ,8051 + ,7699 + ,7440 + ,7656 + ,7549 + ,9086 + ,8942 + ,8764 + ,8500 + ,8239 + ,8443 + ,8349 + ,8288 + ,7970 + ,7496 + ,7745 + ,7543 + ,9036 + ,9075 + ,8859 + ,8605 + ,8419 + ,8495 + ,8284 + ,7582 + ,7691 + ,7046 + ,7442 + ,7596 + ,8597 + ,8436 + ,7881 + ,7477 + ,7508 + ,7361 + ,7299 + ,6914 + ,6768 + ,6746 + ,7052 + ,7139 + ,7714 + ,7750 + ,7622 + ,7424 + ,7444 + ,7208 + ,7128 + ,7022 + ,6688 + ,6199 + ,6400 + ,6474 + ,7182 + ,7330 + ,7410 + ,7442 + ,7753 + ,7762 + ,7814 + ,7838 + ,7298 + ,7155 + ,7076 + ,7450 + ,8216 + ,8246 + ,8335 + ,8171 + ,8485 + ,8435 + ,8369 + ,8210 + ,7888 + ,8061 + ,8139 + ,7837 + ,8943 + ,8523 + ,8104 + ,7969 + ,7921 + ,7930 + ,7706 + ,7552 + ,7379 + ,6946 + ,7128 + ,7393 + ,8092 + ,8004 + ,7903 + ,7710 + ,7867 + ,7860 + ,7723 + ,7477 + ,7126 + ,7161 + ,7162 + ,7406 + ,7944 + ,8084 + ,8088 + ,7972 + ,8184 + ,7914 + ,7845 + ,7610 + ,7278 + ,6883 + ,7123 + ,7182 + ,7912 + ,7893 + ,7671 + ,7403 + ,7663 + ,7589 + ,7450 + ,7069 + ,6670 + ,6285 + ,6506 + ,6539 + ,7291 + ,7391 + ,7126 + ,6752 + ,6835 + ,6664 + ,6562 + ,6174 + ,5741 + ,5398 + ,5203 + ,5673 + ,6379 + ,6418 + ,6272 + ,6059) + ,dim=c(1 + ,204) + ,dimnames=list(c('NWWZ') + ,1:204)) > y <- array(NA,dim=c(1,204),dimnames=list(c('NWWZ'),1:204)) > 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 = 'Aantal mensen' > main = 'NWWZ, -25jr, Man, Brussel' > 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/1m0051456586148.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,] 5673.0 [2,] 7018.5 [3,] 7500.0 [4,] 7971.0 [5,] 9086.0 $n [1] 204 $conf [,1] [1,] 7394.632 [2,] 7605.368 $out [1] 5398 5203 $group [1] 1 1 $names [1] "NWWZ" > 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/2h9gf1456586148.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/3mgx61456586148.tab") > > try(system("convert tmp/1m0051456586148.ps tmp/1m0051456586148.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.539 0.112 0.642