R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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(26,21,21,23,17,23,4,20,16,15,24,17,20,4,19,19,18,22,18,20,6,19,18,11,20,21,21,8,20,16,8,24,20,24,8,25,23,19,27,28,22,4,25,17,4,28,19,23,4,22,12,20,27,22,20,8,26,19,16,24,16,25,5,22,16,14,23,18,23,4,17,19,10,24,25,27,4,22,20,13,27,17,27,4,19,13,14,27,14,22,4,24,20,8,28,11,24,4,26,27,23,27,27,25,4,21,17,11,23,20,22,8,13,8,9,24,22,28,4,26,25,24,28,22,28,4,20,26,5,27,21,27,4,22,13,15,25,23,25,8,14,19,5,19,17,16,4,21,15,19,24,24,28,7,7,5,6,20,14,21,4,23,16,13,28,17,24,4,17,14,11,26,23,27,5,25,24,17,23,24,14,4,25,24,17,23,24,14,4,19,9,5,20,8,27,4,20,19,9,11,22,20,4,23,19,15,24,23,21,4,22,25,17,25,25,22,4,22,19,17,23,21,21,4,21,18,20,18,24,12,15,15,15,12,20,15,20,10,20,12,7,20,22,24,4,22,21,16,24,21,19,8,18,12,7,23,25,28,4),dim=c(7,37),dimnames=list(c('I1','I2','I3','E1','E2','E3','A'),1:37)) > y <- array(NA,dim=c(7,37),dimnames=list(c('I1','I2','I3','E1','E2','E3','A'),1:37)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = '' > par2 = '' > par1 = 'grey' > ylab = '' > xlab = '' > main = '' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > z <- as.data.frame(t(y)) > postscript(file="/var/www/html/rcomp/tmp/14szm1287474556.ps",horizontal=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] [,5] [,6] [,7] [1,] 13 8 4 18 8 14 4 [2,] 19 15 9 23 17 20 4 [3,] 21 18 14 24 21 23 4 [4,] 23 20 17 27 23 25 6 [5,] 26 27 24 28 28 28 8 $n [1] 37 37 37 37 37 37 37 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 19.961 16.70125 11.92200 22.961 19.44150 21.70125 3.480499 [2,] 22.039 19.29875 16.07800 25.039 22.55850 24.29875 4.519501 $out [1] 7 5 11 12 15 10 $group [1] 1 2 4 6 7 7 $names [1] "I1" "I2" "I3" "E1" "E2" "E3" "A" Warning message: In bxp(list(stats = c(13, 19, 21, 23, 26, 8, 15, 18, 20, 27, 4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/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/html/rcomp/tmp/20kfd1287474556.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/html/rcomp/tmp/3tteg1287474556.tab") > try(system("convert tmp/14szm1287474556.ps tmp/14szm1287474556.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.368 0.155 0.572