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(3013.24,2638.53,1995.37,2196.72,2921.44,3494.17,4290.89,4116.68,2987.1,2720.25,1946.81,2350.44,2981.85,3667.03,4443.91,3844.49,2995.55,2745.88,1765.9,2440.25,3080.58,3813.06,4502.64,3720.98,2833.18,2735.7,1635.25,2408.64,3106.22,3917.96,4356.98,3674.4,2848.96,2811.7,1833.42,2472.81,3119.31,3895.51,4591.27,3857.62,2794.83,2799.43,1910.43,2407.6,3061.26,3801.06,4696.96,3801.06,2845.26,2555.28,1959.67,2454.62,3097.31,3570.12,4621.4,3504.37,2915.02,2304.98,1969.6,2448.05,3161.69,3701.61,4562.84,3032.6,2892.63,2214.95,2061.41,2497.84,3257.16,3862.27,4202.52,3047.03,2604.42,2065.81,2093.48,2645.64,3277.01,3970.1,4296.49,2962.34,2641.65,1940.49,2120.88,2756.76,3295.32,4138.52,4435.23,2197.82,2659.81,2042,2174.56,2849.27,3363.99,4199.75,4105.18,2014.45),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/1469x1229173824.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,] 2604.42 1940.490 1635.250 2196.720 2921.440 3494.170 4105.180 2014.450
[2,] 2727.32 2140.380 1871.925 2408.120 3070.920 3684.320 4293.690 2997.470
[3,] 2847.11 2596.905 1964.635 2451.335 3112.765 3837.665 4439.570 3589.385
[4,] 2951.06 2740.790 2077.445 2571.740 3267.085 3944.030 4577.055 3822.775
[5,] 3013.24 2811.700 2174.560 2756.760 3363.990 4199.750 4696.960 4116.680

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2745.061 2323.054 1870.896 2376.707 3023.293 3719.210 4310.325 3212.958
[2,] 2949.159 2870.756 2058.374 2525.963 3202.237 3956.120 4568.815 3965.812

$out
[1] 2849.27

$group
[1] 4

$names
[1] "jaar1" "jaar2" "jaar3" "jaar4" "jaar5" "jaar6" "jaar7" "jaar8"

Warning message:
In bxp(list(stats = c(2604.42, 2727.32, 2847.11, 2951.06, 3013.24,  :
  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/2sp061229173824.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/3h5fj1229173824.tab") 
> 
> system("convert tmp/1469x1229173824.ps tmp/1469x1229173824.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.385   0.171   0.463