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(11956.6
+ ,10519.2
+ ,149
+ ,462
+ ,10723.8
+ ,10414.9
+ ,139
+ ,455
+ ,13938.9
+ ,12476.8
+ ,135
+ ,461
+ ,13979.8
+ ,12384.6
+ ,130
+ ,461
+ ,13807.4
+ ,12266.7
+ ,127
+ ,463
+ ,12973.9
+ ,12919.9
+ ,122
+ ,462
+ ,12509.8
+ ,11497.3
+ ,117
+ ,456
+ ,12934.1
+ ,12142
+ ,112
+ ,455
+ ,14908.3
+ ,13919.4
+ ,113
+ ,456
+ ,13772.1
+ ,12656.8
+ ,149
+ ,472
+ ,13012.6
+ ,12034.1
+ ,157
+ ,472
+ ,14049.9
+ ,13199.7
+ ,157
+ ,471
+ ,11816.5
+ ,10881.3
+ ,147
+ ,465
+ ,11593.2
+ ,11301.2
+ ,137
+ ,459
+ ,14466.2
+ ,13643.9
+ ,132
+ ,465
+ ,13615.9
+ ,12517
+ ,125
+ ,468
+ ,14733.9
+ ,13981.1
+ ,123
+ ,467
+ ,13880.7
+ ,14275.7
+ ,117
+ ,463
+ ,13527.5
+ ,13435
+ ,114
+ ,460
+ ,13584
+ ,13565.7
+ ,111
+ ,462
+ ,16170.2
+ ,16216.3
+ ,112
+ ,461
+ ,13260.6
+ ,12970
+ ,144
+ ,476
+ ,14741.9
+ ,14079.9
+ ,150
+ ,476
+ ,15486.5
+ ,14235
+ ,149
+ ,471
+ ,13154.5
+ ,12213.4
+ ,134
+ ,453
+ ,12621.2
+ ,12581
+ ,123
+ ,443
+ ,15031.6
+ ,14130.4
+ ,116
+ ,442
+ ,15452.4
+ ,14210.8
+ ,117
+ ,444
+ ,15428
+ ,14378.5
+ ,111
+ ,438
+ ,13105.9
+ ,13142.8
+ ,105
+ ,427
+ ,14716.8
+ ,13714.7
+ ,102
+ ,424
+ ,14180
+ ,13621.9
+ ,95
+ ,416
+ ,16202.2
+ ,15379.8
+ ,93
+ ,406
+ ,14392.4
+ ,13306.3
+ ,124
+ ,431
+ ,15140.6
+ ,14391.2
+ ,130
+ ,434
+ ,15960.1
+ ,14909.9
+ ,124
+ ,418
+ ,14351.3
+ ,14025.4
+ ,115
+ ,412
+ ,13230.2
+ ,12951.2
+ ,106
+ ,404
+ ,15202.1
+ ,14344.3
+ ,105
+ ,409
+ ,17056
+ ,16093.4
+ ,105
+ ,412
+ ,16077.7
+ ,15413.6
+ ,101
+ ,406
+ ,13348.2
+ ,14705.7
+ ,95
+ ,398
+ ,16402.4
+ ,15972.8
+ ,93
+ ,397
+ ,16559.1
+ ,16241.4
+ ,84
+ ,385
+ ,16579
+ ,16626.4
+ ,87
+ ,390
+ ,17561.2
+ ,17136.2
+ ,116
+ ,413
+ ,16129.6
+ ,15622.9
+ ,120
+ ,413
+ ,18484.3
+ ,18003.9
+ ,117
+ ,401
+ ,16402.6
+ ,16136.1
+ ,109
+ ,397
+ ,14032.3
+ ,14423.7
+ ,105
+ ,397
+ ,17109.1
+ ,16789.4
+ ,107
+ ,409
+ ,17157.2
+ ,16782.2
+ ,109
+ ,419
+ ,13879.8
+ ,14133.8
+ ,109
+ ,424
+ ,12362.4
+ ,12607
+ ,108
+ ,428
+ ,12683.5
+ ,12004.5
+ ,107
+ ,430
+ ,12608.8
+ ,12175.4
+ ,99
+ ,424
+ ,13583.7
+ ,13268
+ ,103
+ ,433
+ ,12846.3
+ ,12299.3
+ ,131
+ ,456
+ ,12347.1
+ ,11800.6
+ ,137
+ ,459
+ ,13967
+ ,13873.3
+ ,135
+ ,446
+ ,13214.7
+ ,12315
+ ,124
+ ,441)
+ ,dim=c(4
+ ,61)
+ ,dimnames=list(c('UE'
+ ,'IE'
+ ,'WJ'
+ ,'WO')
+ ,1:61))
>  y <- array(NA,dim=c(4,61),dimnames=list(c('UE','IE','WJ','WO'),1:61))
>  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/1foob1258043642.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]
[1,] 10723.8 10414.9   84  385
[2,] 13154.5 12476.8  107  413
[3,] 13979.8 13643.9  117  442
[4,] 15452.4 14423.7  131  461
[5,] 18484.3 17136.2  157  476

$n
[1] 61 61 61 61

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 13514.94 13250.05 112.1448 432.2897
[2,] 14444.66 14037.75 121.8552 451.7103

$out
[1] 18003.9

$group
[1] 2

$names
[1] "UE" "IE" "WJ" "WO"

> 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/2x6681258043642.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/38ndg1258043642.tab") 
> 
> system("convert tmp/1foob1258043642.ps tmp/1foob1258043642.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.385   0.159   0.467