R version 3.3.1 (2016-06-21) -- "Bug in Your Hair" Copyright (C) 2016 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(0.0299 + ,0.0538 + ,0.0098 + ,0.0446 + ,0.0526 + ,0.0499 + ,-0.0202 + ,-0.0202 + ,-0.0662 + ,0.0054 + ,0.0890 + ,0.0098 + ,-0.0244 + ,0.0472 + ,0.0768 + ,0.0098 + ,-0.0784 + ,0.0296 + ,-0.0635 + ,-0.0560 + ,-0.0824 + ,-0.0617 + ,-0.0728 + ,-0.0212 + ,-0.0894 + ,-0.0562 + ,-0.1046 + ,-0.0129 + ,-0.1210 + ,-0.0882 + ,-0.1046 + ,0.0110 + ,-0.0651 + ,-0.1014 + ,-0.0845 + ,-0.0560 + ,-0.0582 + ,-0.0872 + ,-0.0687 + ,-0.0470 + ,-0.1056 + ,-0.4204 + ,-0.0992 + ,-0.0212 + ,-0.1116 + ,-0.0472 + ,-0.0672 + ,-0.1021 + ,-0.1205 + ,-0.0601 + ,-0.0845 + ,-0.0728 + ,-0.0732 + ,-0.0734 + ,-0.0388 + ,-0.0894 + ,-0.0466 + ,-0.0587 + ,-0.0845 + ,-0.0692 + ,-0.1004 + ,0.0078 + ,-0.1093 + ,-0.1168 + ,-0.0425 + ,0.0410 + ,-0.0845 + ,-0.0698 + ,-0.1298 + ,-0.0497 + ,-0.0845 + ,-0.1046 + ,-0.0463 + ,-0.0200 + ,-0.0535 + ,-0.1093 + ,-0.0694 + ,-0.0283 + ,-0.1009 + ,-0.1009 + ,-0.0317 + ,0.0516 + ,-0.1093 + ,-0.0845 + ,-0.1143 + ,NA + ,-0.1009 + ,-0.0461 + ,-0.0877 + ,NA + ,-0.1179 + ,-0.0107 + ,-0.1256 + ,NA + ,-0.1179 + ,-0.1009 + ,-0.0695 + ,NA + ,-0.1179 + ,-0.0845 + ,-0.1009 + ,NA + ,-0.1009 + ,-0.0845 + ,-0.0920 + ,NA + ,-0.0829 + ,-0.0535 + ,-0.1092 + ,NA + ,-0.1009 + ,0.0026 + ,-0.0943 + ,NA + ,-0.0245 + ,-0.0926 + ,-0.0641 + ,NA + ,-0.1030 + ,-0.0687 + ,-0.0845 + ,NA + ,-0.0872 + ,-0.0687 + ,-0.0877 + ,NA + ,-0.0430 + ,-0.1194 + ,-0.0577 + ,NA + ,-0.0872 + ,-0.1194 + ,-0.0834 + ,NA + ,-0.0719 + ,-0.1030 + ,-0.0961 + ,NA + ,-0.0719 + ,-0.0795 + ,-0.0613 + ,NA + ,-0.1030 + ,-0.0795 + ,-0.0804 + ,NA + ,-0.0719 + ,-0.0719 + ,-0.0971 + ,NA + ,-0.0572 + ,-0.0872 + ,-0.0790 + ,NA + ,-0.0872 + ,-0.1030 + ,-0.0707 + ,NA + ,-0.0719 + ,-0.0950 + ,-0.0253 + ,NA + ,-0.0719 + ,-0.0719 + ,-0.1009 + ,NA + ,NA + ,-0.1278 + ,-0.1383 + ,NA + ,NA + ,-0.0872 + ,-0.0712 + ,NA + ,NA + ,-0.0719 + ,-0.1439 + ,NA + ,NA + ,-0.0430 + ,-0.1039 + ,NA + ,NA + ,NA + ,-0.0699 + ,NA + ,NA + ,NA + ,-0.0964 + ,NA + ,NA + ,NA + ,-0.1095 + ,NA + ,NA + ,NA + ,-0.0529 + ,NA + ,NA + ,NA + ,-0.1210 + ,NA + ,NA + ,NA + ,-0.0251 + ,NA + ,NA + ,NA + ,-0.0897 + ,NA + ,NA + ,NA + ,-0.0768 + ,NA + ,NA + ,NA + ,-0.0590 + ,NA + ,NA + ,NA + ,-0.0884 + ,NA + ,NA + ,NA + ,-0.0953 + ,NA + ,NA + ,NA + ,-0.0875 + ,NA + ,NA + ,NA + ,-0.0852 + ,NA + ,NA + ,NA + ,-0.1012 + ,NA + ,NA + ,NA + ,-0.0937 + ,NA + ,NA + ,NA + ,-0.0650 + ,NA + ,NA + ,NA + ,-0.0246 + ,NA + ,NA + ,NA + ,-0.0528 + ,NA + ,NA + ,NA + ,-0.0488 + ,NA + ,NA + ,NA + ,-0.0200 + ,NA + ,NA + ,NA + ,0.0686 + ,NA + ,NA + ,NA + ,-0.0127 + ,NA + ,NA + ,NA) + ,dim=c(4 + ,68) + ,dimnames=list(c('Vaihingen' + ,'Herxheim' + ,'Cuiry-les-Chaudardes' + ,'Eilsleben ') + ,1:68)) > y <- array(NA,dim=c(4,68),dimnames=list(c('Vaihingen','Herxheim','Cuiry-les-Chaudardes','Eilsleben '),1:68)) > 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/1uwfx1477401344.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] [1,] -0.14390 -0.1014 -0.1179 -0.1278 [2,] -0.09875 -0.0617 -0.1009 -0.0950 [3,] -0.08140 -0.0472 -0.0845 -0.0719 [4,] -0.05795 0.0296 -0.0672 -0.0461 [5,] -0.01270 0.0538 -0.0202 0.0110 $n [1] 68 21 41 45 $conf [,1] [,2] [,3] [,4] [1,] -0.08921741 -0.0786788 -0.09281563 -0.08341754 [2,] -0.07358259 -0.0157212 -0.07618437 -0.06038246 $out [1] 0.0299 0.0526 0.0686 -0.4204 0.0098 0.0890 0.0768 0.0446 $group [1] 1 1 1 2 3 3 3 4 $names [1] "Vaihingen" "Herxheim" "Cuiry-les-Chaudardes" [4] "Eilsleben\r\r\r" Warning messages: 1: In bxp(list(stats = c(-0.1439, -0.09875, -0.0814, -0.05795, -0.0127, : some notches went outside hinges ('box'): maybe set notch=FALSE 2: In axis(side = 1, at = 1:4, labels = c("Vaihingen", "Herxheim", : font width unknown for character 0xd 3: In axis(side = 1, at = 1:4, labels = c("Vaihingen", "Herxheim", : font width unknown for character 0xd 4: In axis(side = 1, at = 1:4, labels = c("Vaihingen", "Herxheim", : font width unknown for character 0xd 5: In axis(side = 1, at = 1:4, labels = c("Vaihingen", "Herxheim", : font width unknown for character 0xd 6: In axis(side = 1, at = 1:4, labels = c("Vaihingen", "Herxheim", : font width unknown for character 0xd 7: In axis(side = 1, at = 1:4, labels = c("Vaihingen", "Herxheim", : font width unknown for character 0xd > 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/20rkc1477401344.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/36yv41477401344.tab") > > try(system("convert tmp/1uwfx1477401344.ps tmp/1uwfx1477401344.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.582 0.039 0.627