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(15.8368
+ ,15.7995
+ ,8.9
+ ,17.82
+ ,17.5704
+ ,16.6856
+ ,8.6
+ ,16.52
+ ,18.2521
+ ,18.3578
+ ,8.3
+ ,20.54
+ ,16.1967
+ ,15.9527
+ ,8.3
+ ,18.36
+ ,16.643
+ ,16.5409
+ ,8.3
+ ,18.82
+ ,17.729
+ ,18.1595
+ ,8.4
+ ,20.62
+ ,16.4461
+ ,16.8472
+ ,8.5
+ ,18.64
+ ,15.9938
+ ,17.2825
+ ,8.4
+ ,18.2
+ ,16.3735
+ ,18.0423
+ ,8.6
+ ,18.86
+ ,17.8422
+ ,19.6149
+ ,8.5
+ ,19.88
+ ,22.3215
+ ,23.3357
+ ,8.5
+ ,23.14
+ ,22.7867
+ ,23.6504
+ ,8.4
+ ,23.36
+ ,18.2741
+ ,20.3283
+ ,8.5
+ ,19.96
+ ,22.3929
+ ,22.7516
+ ,8.5
+ ,19.2
+ ,23.8993
+ ,24.3477
+ ,8.5
+ ,23.18
+ ,21.3435
+ ,21.4796
+ ,8.5
+ ,21.82
+ ,22.9523
+ ,23.5338
+ ,8.5
+ ,23.46
+ ,21.3744
+ ,22.1762
+ ,8.5
+ ,21.96
+ ,21.1641
+ ,20.8668
+ ,8.5
+ ,22.56
+ ,20.9065
+ ,22.1687
+ ,8.5
+ ,22.14
+ ,17.8774
+ ,19.4846
+ ,8.6
+ ,20
+ ,20.6643
+ ,21.1315
+ ,8.4
+ ,22.66
+ ,22.16
+ ,21.9689
+ ,8.1
+ ,24.48
+ ,19.8136
+ ,19.2208
+ ,8
+ ,22.5
+ ,17.7354
+ ,17.7896
+ ,8
+ ,20.84
+ ,19.6402
+ ,19.2275
+ ,8
+ ,18.5
+ ,20.8444
+ ,19.8992
+ ,8
+ ,23.44
+ ,19.8231
+ ,19.563
+ ,7.9
+ ,21.86
+ ,18.5946
+ ,18.0752
+ ,7.8
+ ,21.22
+ ,21.3506
+ ,20.4386
+ ,7.8
+ ,23.76
+ ,18.5741
+ ,18.0767
+ ,7.9
+ ,21.06
+ ,18.9242
+ ,18.6373
+ ,8.1
+ ,21.2
+ ,17.3434
+ ,17.5794
+ ,8
+ ,20.4
+ ,19.9612
+ ,19.4221
+ ,7.6
+ ,22.58
+ ,19.9321
+ ,19.27
+ ,7.3
+ ,23.3
+ ,19.4646
+ ,18.553
+ ,7
+ ,22.96
+ ,16.1654
+ ,16.8956
+ ,6.8
+ ,20.1
+ ,17.5749
+ ,16.3989
+ ,7
+ ,17.08
+ ,19.7954
+ ,18.9531
+ ,7.1
+ ,22.92
+ ,19.4395
+ ,18.69
+ ,7.2
+ ,21.98
+ ,17.17
+ ,17.0653
+ ,7.1
+ ,20.14
+ ,21.0724
+ ,21.4899
+ ,6.9
+ ,23.1
+ ,17.7518
+ ,17.7761
+ ,6.7
+ ,20.14
+ ,17.5155
+ ,18.1019
+ ,6.7
+ ,19.8
+ ,18.0403
+ ,19.021
+ ,6.6
+ ,20.46
+ ,19.0901
+ ,18.883
+ ,6.9
+ ,21.76
+ ,17.7465
+ ,16.7687
+ ,7.3
+ ,21.18
+ ,19.2021
+ ,18.1
+ ,7.5
+ ,22.64
+ ,15.1416
+ ,15.6998
+ ,7.3
+ ,19.14
+ ,16.2581
+ ,15.0521
+ ,7.1
+ ,16.18
+ ,18.5865
+ ,18.3209
+ ,6.9
+ ,22.78
+ ,17.2094
+ ,16.505
+ ,7.1
+ ,19.62
+ ,17.8387
+ ,16.889
+ ,7.5
+ ,20.56
+ ,19.1235
+ ,18.4477
+ ,7.7
+ ,20.94
+ ,16.5836
+ ,16.3599
+ ,7.8
+ ,19.18
+ ,15.9912
+ ,15.6639
+ ,7.8
+ ,18.92
+ ,16.7044
+ ,16.8692
+ ,7.7
+ ,20.32
+ ,17.4204
+ ,16.6148
+ ,7.8
+ ,20.78
+ ,17.872
+ ,16.6707
+ ,7.8
+ ,22.06
+ ,17.8232
+ ,16.6296
+ ,7.9
+ ,22.82)
+ ,dim=c(4
+ ,60)
+ ,dimnames=list(c('UV'
+ ,'IV'
+ ,'WL'
+ ,'IP')
+ ,1:60))
>  y <- array(NA,dim=c(4,60),dimnames=list(c('UV','IV','WL','IP'),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 = 'waarden'
> xlab = 'variabelen'
> 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/10hof1258017305.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,] 15.14160 15.05210 6.60 16.18
[2,] 17.38190 16.85820 7.30 19.71
[3,] 18.26310 18.33935 7.90 20.89
[4,] 19.94665 19.75705 8.45 22.61
[5,] 22.95230 23.65040 8.90 24.48

$n
[1] 60 60 60 60

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 17.73995 17.74805 7.665426 20.29847
[2,] 18.78625 18.93065 8.134574 21.48153

$out
[1] 23.8993 24.3477

$group
[1] 1 2

$names
[1] "UV" "IV" "WL" "IP"

> 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/2mid11258017305.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/3qcvt1258017305.tab") 
> 
> system("convert tmp/10hof1258017305.ps tmp/10hof1258017305.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.560   0.218   0.633