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(26.54,19.32,27.46,29.59,40.22,57.42,62.5,90.6,25.85,19.82,30.23,30.7,44.23,63.52,55.04,92.23,27.6,20.36,32.33,30.52,45.85,61.71,58.34,94.09,25.71,24.31,29.87,32.67,53.38,63.01,61.92,102.79,25.38,25.97,24.87,33.19,53.26,68.18,67.65,109.65,28.57,25.61,25.48,37.13,51.8,72.03,67.68,124.05,27.64,24.67,27.28,35.54,55.3,69.75,70.3,132.69,25.36,25.59,28.24,37.75,57.81,74.41,75.26,135.81,25.9,26.09,29.58,41.84,63.96,74.33,71.44,116.07,26.29,28.37,26.95,42.94,63.77,64.24,76.36,101.42,21.74,27.34,29.08,49.14,59.15,60.03,81.71,75.73,19.2,24.46,28.76,44.61,56.12,59.44,92.6,55.48),dim=c(8,12),dimnames=list(c('Jaar1','jaar2','jaar3','jaar4','jaar5','Jaar6','Jaar7','jaar8'),1:12))
>  y <- array(NA,dim=c(8,12),dimnames=list(c('Jaar1','jaar2','jaar3','jaar4','jaar5','Jaar6','Jaar7','jaar8'),1:12))
>  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/17wjd1229174217.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]   [,5]  [,6]  [,7]    [,8]
[1,] 25.360 19.320 24.870 29.590 40.220 57.42 55.04  55.480
[2,] 25.370 22.335 27.115 31.685 48.825 60.87 62.21  91.415
[3,] 25.875 25.130 28.500 36.335 54.340 63.88 68.99 102.105
[4,] 27.070 26.030 29.725 42.390 58.480 70.89 75.81 120.060
[5,] 28.570 28.370 32.330 49.140 63.960 74.41 92.60 135.810

$n
[1] 12 12 12 12 12 12 12 12

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]      [,8]
[1,] 25.09962 23.44469 27.30956 31.45238 49.93629 59.30981 62.78695  89.03982
[2,] 26.65038 26.81531 29.69044 41.21762 58.74371 68.45019 75.19305 115.17018

$out
[1] 21.74 19.20

$group
[1] 1 1

$names
[1] "Jaar1" "jaar2" "jaar3" "jaar4" "jaar5" "Jaar6" "Jaar7" "jaar8"

Warning message:
In bxp(list(stats = c(25.36, 25.37, 25.875, 27.07, 28.57, 19.32,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> 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/27lij1229174217.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/3uayk1229174217.tab") 
> 
> system("convert tmp/17wjd1229174217.ps tmp/17wjd1229174217.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.391   0.156   0.470