R version 2.8.1 (2008-12-22)
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(9356
+ ,7291
+ ,14271
+ ,4071
+ ,9337
+ ,6820
+ ,14013
+ ,4351
+ ,10149
+ ,8031
+ ,15912
+ ,4871
+ ,9788
+ ,7862
+ ,14290
+ ,4649
+ ,9770
+ ,7357
+ ,14744
+ ,4922
+ ,9911
+ ,7213
+ ,14721
+ ,4879
+ ,9429
+ ,7079
+ ,13918
+ ,4853
+ ,8775
+ ,7012
+ ,13263
+ ,4545
+ ,10189
+ ,7319
+ ,15660
+ ,4733
+ ,10529
+ ,8148
+ ,15629
+ ,5191
+ ,9914
+ ,7599
+ ,15113
+ ,4983
+ ,9790
+ ,6908
+ ,14526
+ ,4593
+ ,9625
+ ,7878
+ ,15132
+ ,4656
+ ,9729
+ ,7407
+ ,14908
+ ,4513
+ ,10589
+ ,7911
+ ,16167
+ ,4857
+ ,9611
+ ,7323
+ ,14122
+ ,4681
+ ,9388
+ ,7179
+ ,13985
+ ,4897
+ ,9510
+ ,6758
+ ,14236
+ ,4547
+ ,9690
+ ,6934
+ ,13921
+ ,4692
+ ,8434
+ ,6696
+ ,12394
+ ,4390
+ ,9844
+ ,7688
+ ,15454
+ ,5341
+ ,10601
+ ,8296
+ ,16146
+ ,5415
+ ,9942
+ ,7697
+ ,15107
+ ,4890
+ ,10229
+ ,7907
+ ,14593
+ ,5120
+ ,9381
+ ,7592
+ ,14695
+ ,4422
+ ,9635
+ ,7710
+ ,14513
+ ,4797
+ ,11228
+ ,9011
+ ,17071
+ ,5689
+ ,9999
+ ,8225
+ ,15179
+ ,5171
+ ,10089
+ ,7733
+ ,15460
+ ,4265
+ ,11622
+ ,8062
+ ,17173
+ ,5215
+ ,10533
+ ,7859
+ ,15938
+ ,4874
+ ,9965
+ ,8221
+ ,15003
+ ,4590
+ ,11567
+ ,8330
+ ,18216
+ ,4994
+ ,11321
+ ,8868
+ ,17847
+ ,4988
+ ,11686
+ ,9053
+ ,18029
+ ,5110
+ ,11747
+ ,8811
+ ,17281
+ ,5141
+ ,10595
+ ,8120
+ ,16706
+ ,4395
+ ,10751
+ ,7953
+ ,16750
+ ,4523
+ ,12199
+ ,8878
+ ,18912
+ ,5306
+ ,11690
+ ,8601
+ ,17763
+ ,5365
+ ,10978
+ ,8361
+ ,16736
+ ,5496
+ ,11753
+ ,9116
+ ,18061
+ ,5647
+ ,10839
+ ,9310
+ ,16713
+ ,5443
+ ,10518
+ ,9891
+ ,16769
+ ,5546
+ ,12183
+ ,10147
+ ,19514
+ ,5912
+ ,11967
+ ,10317
+ ,19251
+ ,5665
+ ,12363
+ ,10682
+ ,19951
+ ,5963
+ ,12359
+ ,10276
+ ,19052
+ ,5861
+ ,12162
+ ,10614
+ ,19555
+ ,5366
+ ,12096
+ ,9413
+ ,19083
+ ,5619
+ ,14325
+ ,11068
+ ,22534
+ ,6721
+ ,12670
+ ,9772
+ ,18854
+ ,6054
+ ,13865
+ ,10350
+ ,19801
+ ,6619
+ ,13563
+ ,10541
+ ,20346
+ ,6856
+ ,12734
+ ,10049
+ ,18169
+ ,6193
+ ,12464
+ ,10714
+ ,19087
+ ,6317
+ ,13389
+ ,10759
+ ,20842
+ ,6618
+ ,13961
+ ,11684
+ ,21602
+ ,6585
+ ,14088
+ ,11462
+ ,22360
+ ,6852
+ ,13143
+ ,10485
+ ,20334
+ ,6586
+ ,13413
+ ,11056
+ ,21215
+ ,6154
+ ,13579
+ ,10184
+ ,20530
+ ,6193
+ ,15388
+ ,11082
+ ,23152
+ ,7606
+ ,13708
+ ,10554
+ ,20134
+ ,6588
+ ,14689
+ ,11315
+ ,21193
+ ,7143
+ ,14883
+ ,10847
+ ,21628
+ ,7629
+ ,13991
+ ,11104
+ ,20823
+ ,7041
+ ,13854
+ ,11026
+ ,20493
+ ,7146
+ ,14364
+ ,11073
+ ,22106
+ ,7200
+ ,15672
+ ,12073
+ ,24178
+ ,7739
+ ,15904
+ ,12328
+ ,24958
+ ,7953
+ ,14016
+ ,11172
+ ,21620
+ ,7082)
+ ,dim=c(4
+ ,72)
+ ,dimnames=list(c('InvEU'
+ ,'InvnietEU'
+ ,'UitvEU'
+ ,'Uivniet-EU')
+ ,1:72))
>  y <- array(NA,dim=c(4,72),dimnames=list(c('InvEU','InvnietEU','UitvEU','Uivniet-EU'),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'
> #'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/rcomp/tmp/1myqe1257280656.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,]  8434.0  6696.0 12394.0 4071.0
[2,]  9912.5  7703.5 15055.0 4825.0
[3,] 11274.5  8706.0 17122.0 5323.5
[4,] 13266.0 10547.5 20042.5 6255.0
[5,] 15904.0 12328.0 24958.0 7953.0

$n
[1] 72 72 72 72

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 10650.06 8176.434 16193.30 5057.227
[2,] 11898.94 9235.566 18050.70 5589.773

$out
numeric(0)

$group
numeric(0)

$names
[1] "InvEU"      "InvnietEU"  "UitvEU"     "Uivniet-EU"

> 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/20oqh1257280656.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/31qju1257280656.tab") 
> 
> system("convert tmp/1myqe1257280656.ps tmp/1myqe1257280656.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.410   0.350   0.698