R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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(207.221 + ,247.182 + ,210.506 + ,251.362 + ,216.222 + ,257.784 + ,224.325 + ,264.661 + ,214.419 + ,258.558 + ,189.422 + ,235.162 + ,188.222 + ,235.042 + ,191.395 + ,238.575 + ,195.978 + ,242.254 + ,201.920 + ,249.192 + ,203.583 + ,250.333 + ,191.987 + ,242.470 + ,186.424 + ,244.894 + ,191.568 + ,248.197 + ,200.140 + ,254.511 + ,211.451 + ,261.678 + ,211.524 + ,253.299 + ,182.510 + ,231.770 + ,182.414 + ,233.271 + ,187.827 + ,237.520 + ,192.684 + ,242.082 + ,198.480 + ,243.287 + ,203.089 + ,246.526 + ,199.834 + ,245.622 + ,195.770 + ,247.447 + ,203.555 + ,252.470 + ,211.868 + ,262.976 + ,229.527 + ,267.145 + ,229.139 + ,260.951 + ,198.563 + ,239.171 + ,195.722 + ,238.454 + ,202.196 + ,245.385 + ,205.816 + ,252.487 + ,212.588 + ,257.375 + ,214.320 + ,261.240 + ,220.376 + ,252.222 + ,204.442 + ,250.627 + ,206.903 + ,255.159 + ,214.126 + ,263.948 + ,226.899 + ,268.145 + ,223.532 + ,263.281 + ,195.309 + ,242.967 + ,186.005 + ,241.326 + ,188.906 + ,248.985 + ,191.563 + ,247.599 + ,189.226 + ,247.323 + ,186.412 + ,246.546 + ,178.037 + ,243.985 + ,166.827 + ,241.478 + ,169.362 + ,242.936 + ,174.330 + ,251.178 + ,187.069 + ,254.016 + ,186.530 + ,249.999 + ,158.114 + ,229.521 + ,151.001 + ,228.793 + ,159.612 + ,239.709 + ,161.914 + ,239.494 + ,164.182 + ,251.846 + ,169.701 + ,253.376 + ,171.297 + ,248.284) + ,dim=c(2 + ,60) + ,dimnames=list(c('Vlaanderen' + ,'Wallonie') + ,1:60)) > y <- array(NA,dim=c(2,60),dimnames=list(c('Vlaanderen','Wallonie'),1:60)) > 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/wessaorg/rcomp/tmp/1s8b81355247626.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] [1,] 158.1140 228.793 [2,] 186.4180 242.168 [3,] 195.7460 247.898 [4,] 208.8635 253.696 [5,] 229.5270 268.145 $n [1] 60 60 $conf [,1] [,2] [1,] 191.1676 245.5466 [2,] 200.3244 250.2494 $out [1] 151.001 $group [1] 1 $names [1] "Vlaanderen" "Wallonie" > 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/249ad1355247627.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/3ce0d1355247627.tab") > > try(system("convert tmp/1s8b81355247626.ps tmp/1s8b81355247626.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.976 0.144 1.098