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(100
+ ,100
+ ,100
+ ,100
+ ,100
+ ,100.2738004
+ ,104.011209
+ ,98.53703704
+ ,100
+ ,97.74002424
+ ,107.3498799
+ ,95.08333333
+ ,100
+ ,93.72278828
+ ,105.0360288
+ ,94.34259259
+ ,100
+ ,89.72126218
+ ,104.1953563
+ ,91.44444444
+ ,100
+ ,87.43435522
+ ,105.692554
+ ,94.11111111
+ ,100
+ ,107.2955474
+ ,103.5868695
+ ,91.72222222
+ ,100
+ ,100.0869653
+ ,101.6333066
+ ,92.62037037
+ ,100
+ ,96.55617398
+ ,97.39791833
+ ,97.81481481
+ ,100
+ ,104.2826204
+ ,96.37309848
+ ,94.63888889
+ ,100
+ ,91.2187486
+ ,98.41473179
+ ,93.07407407
+ ,100
+ ,84.9544414
+ ,98.1265012
+ ,97.33333333
+ ,100
+ ,107.7365456
+ ,96.19695757
+ ,97.89814815
+ ,100
+ ,99.56966201
+ ,94.36349079
+ ,100.1944444
+ ,110.5
+ ,107.1081512
+ ,94.92393915
+ ,94.13888889
+ ,112.5
+ ,101.2180753
+ ,96.90152122
+ ,100.9351852
+ ,112.5
+ ,98.27359846
+ ,95.58046437
+ ,100.8888889
+ ,122.5
+ ,99.59939854
+ ,96.23698959
+ ,101.2407407
+ ,125
+ ,118.2301719
+ ,98.24659728
+ ,102.3981481
+ ,125
+ ,96.33511378
+ ,102.2417934
+ ,97.34259259
+ ,132
+ ,109.0685174
+ ,101.2810248
+ ,95.09259259
+ ,137.5
+ ,111.065353
+ ,101.5532426
+ ,96.28703704
+ ,146.5
+ ,98.60687194
+ ,102.570056
+ ,99.44444444
+ ,150
+ ,90.69863997
+ ,101.897518
+ ,95.02777778
+ ,158.5
+ ,109.2093451
+ ,100.968775
+ ,91.50925926
+ ,162.5
+ ,111.8323309
+ ,103.1305044
+ ,94.48148148
+ ,169.5
+ ,111.9956012
+ ,105.7886309
+ ,94.4537037
+ ,175
+ ,97.30800305
+ ,104.0752602
+ ,96.44444444
+ ,175
+ ,106.1773419
+ ,104.6757406
+ ,97.75
+ ,182.5
+ ,104.2130482
+ ,106.0208167
+ ,94.12962963
+ ,187.5
+ ,119.7910588
+ ,108.2145717
+ ,93.2037037
+ ,187.5
+ ,104.3280668
+ ,108.1745396
+ ,93.67592593
+ ,195
+ ,111.2712644
+ ,107.4379504
+ ,92.34259259
+ ,200
+ ,116.9509403
+ ,109.8158527
+ ,95.92592593
+ ,200
+ ,110.1945779
+ ,109.0632506
+ ,93.82407407
+ ,200
+ ,99.50738363
+ ,111.2570056
+ ,98.76851852
+ ,200
+ ,111.1674671
+ ,113.9071257
+ ,100.7685185
+ ,200
+ ,124.3323309
+ ,117.5660528
+ ,99.94444444
+ ,200
+ ,115.9404596
+ ,116.6533227
+ ,103.8148148
+ ,200
+ ,100.304098
+ ,117.8382706
+ ,104.0925926
+ ,200
+ ,117.2993626
+ ,118.0784628
+ ,106.3240741
+ ,200
+ ,118.7446699
+ ,124.3154524
+ ,102.8333333
+ ,200
+ ,119.9245927
+ ,126.1008807
+ ,102.0277778
+ ,200
+ ,128.7776606
+ ,124.5556445
+ ,99.69444444
+ ,200
+ ,119.7512231
+ ,124.5236189
+ ,97.93518519
+ ,209
+ ,134.0909601
+ ,126.2610088
+ ,100.5462963
+ ,212.5
+ ,125.6390547
+ ,119.8959167
+ ,97.91666667
+ ,212.5
+ ,102.5298487
+ ,115.0440352
+ ,98.88888889
+ ,198.5
+ ,127.8485345
+ ,106.6613291
+ ,101.1666667
+ ,171
+ ,125.2384533
+ ,101.93755
+ ,109.0833333
+ ,137.5
+ ,100.1066026
+ ,107.6781425
+ ,104.0833333
+ ,115.5
+ ,91.86621931
+ ,105.9967974
+ ,114.3240741
+ ,100
+ ,89.39921002
+ ,102.3618895
+ ,117.5833333
+ ,83
+ ,92.27355357
+ ,104.4835869
+ ,113.0925926
+ ,65.5
+ ,99.47147538
+ ,105.6044836
+ ,111.7222222
+ ,54.5
+ ,93.37829346
+ ,109.2874299
+ ,111.8981481
+ ,50
+ ,90.87425378
+ ,112.2177742
+ ,110.1296296
+ ,50
+ ,102.4064141
+ ,112.7942354
+ ,109.1944444
+ ,50
+ ,98.58162395
+ ,114.2353883
+ ,144.1388889
+ ,50
+ ,88.85497554
+ ,116.5892714
+ ,149.3981481)
+ ,dim=c(4
+ ,60)
+ ,dimnames=list(c('Rente'
+ ,'Wisselkoers'
+ ,'Reserves'
+ ,'Uitvoer')
+ ,1:60))
>  y <- array(NA,dim=c(4,60),dimnames=list(c('Rente','Wisselkoers','Reserves','Uitvoer'),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/freestat/rcomp/tmp/1yyt81257958629.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,]  50.00  84.95444  94.3635  91.44444
[2,] 100.00  97.52401 101.5933  94.56019
[3,] 134.75 102.46813 105.6485  98.65278
[4,] 199.25 111.91397 112.5060 102.61574
[5,] 212.50 128.77766 126.2610 114.32407

$n
[1] 60 60 60 60

$conf
         [,1]      [,2]     [,3]      [,4]
[1,] 114.5053  99.53291 103.4226  97.00963
[2,] 154.9947 105.40335 107.8745 100.29593

$out
[1] 134.0910 117.5833 144.1389 149.3981

$group
[1] 2 4 4 4

$names
[1] "Rente"       "Wisselkoers" "Reserves"    "Uitvoer"    

> 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/2l7ci1257958629.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/3pn5i1257958629.tab") 
> 
> system("convert tmp/1yyt81257958629.ps tmp/1yyt81257958629.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.564   0.226   0.639