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(100 + ,100 + ,100 + ,96.62921348 + ,99.99006359 + ,99.88118812 + ,93.25842697 + ,99.79133545 + ,99.55445545 + ,93.25842697 + ,100.1589825 + ,99.83168317 + ,93.25842697 + ,100.6856121 + ,100.5049505 + ,94.38202247 + ,101.2321145 + ,101.1485149 + ,95.50561798 + ,101.2023052 + ,101.3762376 + ,94.38202247 + ,101.391097 + ,101.5247525 + ,96.62921348 + ,101.5699523 + ,101.8316832 + ,95.50561798 + ,101.977345 + ,102.4455446 + ,95.50561798 + ,101.9972178 + ,102.5445545 + ,95.50561798 + ,101.8879173 + ,102.6633663 + ,95.50561798 + ,101.7189984 + ,102.4752475 + ,95.50561798 + ,102.027027 + ,102.4455446 + ,95.50561798 + ,102.1065183 + ,102.4257426 + ,95.50561798 + ,102.1661367 + ,102.4554455 + ,95.50561798 + ,102.6530207 + ,102.9009901 + ,95.50561798 + ,102.5735294 + ,102.8613861 + ,95.50561798 + ,102.9411765 + ,103.3663366 + ,95.50561798 + ,103.2889507 + ,103.7524752 + ,96.62921348 + ,103.2690779 + ,103.7326733 + ,94.38202247 + ,103.5870429 + ,104.0891089 + ,91.01123596 + ,103.7162162 + ,104.2178218 + ,89.88764045 + ,103.6963434 + ,103.9207921 + ,89.88764045 + ,103.6565978 + ,103.7128713 + ,89.88764045 + ,103.9149444 + ,103.970297 + ,89.88764045 + ,104.0143084 + ,104.1089109 + ,88.76404494 + ,104.2527822 + ,104.1584158 + ,87.64044944 + ,104.7893482 + ,104.7227723 + ,87.64044944 + ,104.5608108 + ,104.7326733 + ,88.76404494 + ,104.9085851 + ,105.2079208 + ,91.01123596 + ,104.6701113 + ,105.0792079 + ,89.88764045 + ,104.6104928 + ,105.0693069 + ,85.39325843 + ,105.0278219 + ,105.5148515 + ,82.02247191 + ,104.9980127 + ,105.3861386 + ,78.65168539 + ,105.0377583 + ,105.4851485 + ,76.40449438 + ,105.5147059 + ,106.039604 + ,78.65168539 + ,106.25 + ,107.029703 + ,79.7752809 + ,106.7567568 + ,107.3267327 + ,80.8988764 + ,107.1641494 + ,107.7623762 + ,79.7752809 + ,108.0186804 + ,108.5346535 + ,77.52808989 + ,108.6248013 + ,109.3267327 + ,75.28089888 + ,108.7937202 + ,109.5742574 + ,75.28089888 + ,109.4992051 + ,110.5544554 + ,74.15730337 + ,109.9165342 + ,111.1683168 + ,77.52808989 + ,110.5127186 + ,111.7524752 + ,82.02247191 + ,110.1748808 + ,111.0693069 + ,84.26966292 + ,110.4431638 + ,111.2475248 + ,82.02247191 + ,110.5822734 + ,111.049505 + ,79.7752809 + ,110.3835453 + ,110.3861386 + ,77.52808989 + ,110.5325914 + ,110.1485149 + ,79.7752809 + ,110.741256 + ,110.2574257 + ,84.26966292 + ,111.0393482 + ,110.6336634 + ,86.51685393 + ,110.3636725 + ,110 + ,87.64044944 + ,110.4630366 + ,110.2277228 + ,87.64044944 + ,110.254372 + ,110.1485149 + ,86.51685393 + ,109.7972973 + ,109.9405941 + ,86.51685393 + ,109.7774245 + ,109.8712871 + ,87.64044944 + ,109.9562798 + ,110.2079208 + ,88.76404494 + ,109.7575517 + ,109.9207921) + ,dim=c(3 + ,60) + ,dimnames=list(c('wl' + ,'icp' + ,'g_index') + ,1:60)) > y <- array(NA,dim=c(3,60),dimnames=list(c('wl','icp','g_index'),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' > #'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/1ea2l1258322726.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] [1,] 74.15730 99.79134 99.55446 [2,] 82.02247 102.13633 102.60396 [3,] 88.76404 104.64030 104.90099 [4,] 95.50562 109.76749 109.93069 [5,] 100.00000 111.03935 111.75248 $n [1] 60 60 60 $conf [,1] [,2] [,3] [1,] 86.0138 103.0837 103.4065 [2,] 91.5143 106.1969 106.3955 $out numeric(0) $group numeric(0) $names [1] "wl" "icp" "g_index" > 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/2ue2l1258322726.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/37qww1258322726.tab") > > system("convert tmp/1ea2l1258322726.ps tmp/1ea2l1258322726.png") > > > proc.time() user system elapsed 0.371 0.156 0.604