R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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(116222
+ ,344744
+ ,31899
+ ,492865
+ ,110924
+ ,338653
+ ,31384
+ ,480961
+ ,103753
+ ,327532
+ ,30650
+ ,461935
+ ,99983
+ ,326225
+ ,30400
+ ,456608
+ ,93302
+ ,318672
+ ,30003
+ ,441977
+ ,91496
+ ,317756
+ ,29896
+ ,439148
+ ,119321
+ ,337302
+ ,31557
+ ,488180
+ ,139261
+ ,349420
+ ,31883
+ ,520564
+ ,133739
+ ,336923
+ ,30830
+ ,501492
+ ,123913
+ ,330758
+ ,30354
+ ,485025
+ ,113438
+ ,321002
+ ,29756
+ ,464196
+ ,109416
+ ,320820
+ ,29934
+ ,460170
+ ,109406
+ ,327032
+ ,30599
+ ,467037
+ ,105645
+ ,324047
+ ,30378
+ ,460070
+ ,101328
+ ,316735
+ ,29925
+ ,447988
+ ,97686
+ ,315710
+ ,29471
+ ,442867
+ ,93093
+ ,313427
+ ,29567
+ ,436087
+ ,91382
+ ,310527
+ ,29419
+ ,431328
+ ,122257
+ ,330962
+ ,30796
+ ,484015
+ ,139183
+ ,339015
+ ,31475
+ ,509673
+ ,139887
+ ,341332
+ ,31708
+ ,512927
+ ,131822
+ ,339092
+ ,31917
+ ,502831
+ ,116805
+ ,323308
+ ,30871
+ ,470984
+ ,113706
+ ,325849
+ ,31512
+ ,471067
+ ,113012
+ ,330675
+ ,32362
+ ,476049
+ ,110452
+ ,332225
+ ,31928
+ ,474605
+ ,107005
+ ,331735
+ ,31699
+ ,470439
+ ,102841
+ ,328047
+ ,30363
+ ,461251
+ ,98173
+ ,326165
+ ,30386
+ ,454724
+ ,98181
+ ,327081
+ ,30364
+ ,455626
+ ,137277
+ ,346764
+ ,32806
+ ,516847
+ ,147579
+ ,344190
+ ,33423
+ ,525192
+ ,146571
+ ,343333
+ ,33071
+ ,522975
+ ,138920
+ ,345777
+ ,33888
+ ,518585
+ ,130340
+ ,344094
+ ,34805
+ ,509239
+ ,128140
+ ,348609
+ ,35489
+ ,512238
+ ,127059
+ ,354846
+ ,37259
+ ,519164
+ ,122860
+ ,356427
+ ,37722
+ ,517009
+ ,117702
+ ,353467
+ ,38764
+ ,509933
+ ,113537
+ ,355996
+ ,39594
+ ,509127
+ ,108366
+ ,352487
+ ,40004
+ ,500857
+ ,111078
+ ,355178
+ ,40715
+ ,506971
+ ,150739
+ ,374556
+ ,44028
+ ,569323
+ ,159129
+ ,375021
+ ,45564
+ ,579714
+ ,157928
+ ,375787
+ ,44277
+ ,577992
+ ,147768
+ ,372720
+ ,44976
+ ,565464
+ ,137507
+ ,364431
+ ,45406
+ ,547344
+ ,136919
+ ,370490
+ ,47379
+ ,554788
+ ,136151
+ ,376974
+ ,49200
+ ,562325
+ ,133001
+ ,377632
+ ,50221
+ ,560854
+ ,125554
+ ,378205
+ ,51573
+ ,555332
+ ,119647
+ ,370861
+ ,53091
+ ,543599
+ ,114158
+ ,369167
+ ,53337
+ ,536662
+ ,116193
+ ,371551
+ ,54978
+ ,542722
+ ,152803
+ ,382842
+ ,57885
+ ,593530
+ ,161761
+ ,381903
+ ,67099
+ ,610763
+ ,160942
+ ,384502
+ ,67169
+ ,612613
+ ,149470
+ ,392058
+ ,69796
+ ,611324
+ ,139208
+ ,384359
+ ,70600
+ ,594167
+ ,134588
+ ,388884
+ ,71982
+ ,595454
+ ,130322
+ ,386586
+ ,73957
+ ,590865
+ ,126611
+ ,387495
+ ,75273
+ ,589379
+ ,122401
+ ,385705
+ ,76322
+ ,584428
+ ,117352
+ ,378670
+ ,77078
+ ,573100
+ ,112135
+ ,377367
+ ,77954
+ ,567456
+ ,112879
+ ,376911
+ ,79238
+ ,569028
+ ,148729
+ ,389827
+ ,82179
+ ,620735
+ ,157230
+ ,387820
+ ,83834
+ ,628884
+ ,157221
+ ,387267
+ ,83744
+ ,628232
+ ,146681
+ ,380575
+ ,84861
+ ,612117
+ ,136524
+ ,372402
+ ,86478
+ ,595404
+ ,132111
+ ,376740
+ ,88290
+ ,597141)
+ ,dim=c(4
+ ,72)
+ ,dimnames=list(c('min25'
+ ,'25-50'
+ ,'50+'
+ ,'totaal')
+ ,1:72))
>  y <- array(NA,dim=c(4,72),dimnames=list(c('min25','25-50','50+','totaal'),1:72))
>  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))
> postscript(file="/var/www/html/freestat/rcomp/tmp/1gv1o1257701767.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,]  91382.0 310527.0 29419.0 431328.0
[2,] 111001.0 330716.5 30850.5 472836.0
[3,] 123386.5 350953.5 36374.0 516928.0
[4,] 139051.5 376942.5 56431.5 571211.5
[5,] 161761.0 392058.0 88290.0 628884.0

$n
[1] 72 72 72 72

$conf
         [,1]   [,2]     [,3]   [,4]
[1,] 118163.4 342346 31610.70 498610
[2,] 128609.6 359561 41137.30 535246

$out
numeric(0)

$group
numeric(0)

$names
[1] "min25"  "25-50"  "50+"    "totaal"

> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/2bxiw1257701767.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/freestat/rcomp/tmp/34tzh1257701767.tab") 
> 
> system("convert tmp/1gv1o1257701767.ps tmp/1gv1o1257701767.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.568   0.222   0.639