R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-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(1117 + ,698 + ,1120 + ,700 + ,1140 + ,641 + ,1143 + ,857 + ,1152 + ,864 + ,1152 + ,864 + ,1176 + ,735 + ,1176 + ,735 + ,1257 + ,785 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,768 + ,1280 + ,800 + ,1280 + ,768 + ,1280 + ,800 + ,1280 + ,768 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,1024 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,768 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1280 + ,800 + ,1366 + ,768 + ,1366 + ,768 + ,1366 + ,768 + ,1366 + ,768 + ,1366 + ,768 + ,1366 + ,768 + ,1366 + ,768 + ,1366 + ,768 + ,1408 + ,880 + ,1408 + ,880 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1440 + ,900 + ,1503 + ,845 + ,1600 + ,900 + ,1600 + ,900 + ,1600 + ,900 + ,1600 + ,900 + ,1680 + ,1050 + ,1680 + ,1050 + ,1680 + ,1050 + ,1680 + ,1050) + ,dim=c(2 + ,110) + ,dimnames=list(c('Width' + ,'Height') + ,1:110)) > y <- array(NA,dim=c(2,110),dimnames=list(c('Width','Height'),1:110)) > 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' > z <- as.data.frame(t(y)) > z[z[,1]*z[,2]>800000 & z[,1]*z[,2]<1700000,] Width Height 4 1143 857 5 1152 864 6 1152 864 7 1176 735 8 1176 735 9 1257 785 10 1280 1024 11 1280 800 12 1280 1024 13 1280 800 14 1280 800 15 1280 1024 16 1280 800 17 1280 800 18 1280 1024 19 1280 800 20 1280 800 21 1280 768 22 1280 800 23 1280 768 24 1280 800 25 1280 768 26 1280 1024 27 1280 800 28 1280 800 29 1280 800 30 1280 800 31 1280 800 32 1280 800 33 1280 1024 34 1280 800 35 1280 800 36 1280 800 37 1280 800 38 1280 800 39 1280 1024 40 1280 800 41 1280 800 42 1280 1024 43 1280 800 44 1280 800 45 1280 800 46 1280 1024 47 1280 800 48 1280 800 49 1280 800 50 1280 800 51 1280 800 52 1280 800 53 1280 800 54 1280 800 55 1280 800 56 1280 800 57 1280 800 58 1280 800 59 1280 1024 60 1280 800 61 1280 800 62 1280 800 63 1280 800 64 1280 800 65 1280 768 66 1280 800 67 1280 800 68 1280 800 69 1280 800 70 1280 800 71 1366 768 72 1366 768 73 1366 768 74 1366 768 75 1366 768 76 1366 768 77 1366 768 78 1366 768 79 1408 880 80 1408 880 81 1440 900 82 1440 900 83 1440 900 84 1440 900 85 1440 900 86 1440 900 87 1440 900 88 1440 900 89 1440 900 90 1440 900 91 1440 900 92 1440 900 93 1440 900 94 1440 900 95 1440 900 96 1440 900 97 1440 900 98 1440 900 99 1440 900 100 1440 900 101 1440 900 102 1503 845 103 1600 900 104 1600 900 105 1600 900 106 1600 900 > postscript(file="/var/www/rcomp/tmp/1tp6w1317754485.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,] 1117 698 [2,] 1280 800 [3,] 1280 800 [4,] 1440 900 [5,] 1680 1050 $n [1] 110 110 $conf [,1] [,2] [1,] 1255.896 784.9353 [2,] 1304.104 815.0647 $out [1] 641 $group [1] 2 $names [1] "Width" "Height" Warning message: In bxp(list(stats = c(1117, 1280, 1280, 1440, 1680, 698, 800, 800, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/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/www/rcomp/tmp/2s0r51317754485.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/www/rcomp/tmp/3bpjs1317754485.tab") > > try(system("convert tmp/1tp6w1317754485.ps tmp/1tp6w1317754485.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.672 0.096 0.753