R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(8357.00
+ ,7454.00
+ ,8076.00
+ ,7248.00
+ ,7339.00
+ ,7292.00
+ ,7359.00
+ ,7537.00
+ ,7441.00
+ ,8057.00
+ ,8037.00
+ ,8257.00
+ ,8692.00
+ ,8119.00
+ ,8236.00
+ ,7432.00
+ ,7669.00
+ ,7453.00
+ ,7566.00
+ ,7731.00
+ ,7657.00
+ ,8130.00
+ ,8401.00
+ ,8737.00
+ ,9009.00
+ ,7919.00
+ ,8228.00
+ ,7903.00
+ ,7912.00
+ ,7857.00
+ ,7965.00
+ ,8091.00
+ ,8024.00
+ ,8772.00
+ ,8656.00
+ ,8953.00
+ ,9014.00
+ ,8103.00
+ ,8876.00
+ ,8231.00
+ ,8173.00
+ ,8087.00
+ ,8296.00
+ ,8007.00
+ ,8382.00
+ ,9168.00
+ ,9137.00
+ ,9321.00
+ ,9234.00
+ ,8451.00
+ ,9101.00
+ ,8279.00
+ ,8284.00
+ ,8225.00
+ ,8597.00
+ ,8305.00
+ ,8620.00
+ ,9102.00
+ ,9258.00
+ ,9652.00
+ ,9522.00
+ ,8874.00
+ ,9415.00
+ ,8525.00
+ ,8862.00
+ ,8421.00
+ ,8626.00
+ ,8750.00
+ ,8852.00
+ ,9412.00
+ ,9570.00
+ ,9513.00
+ ,9986.00
+ ,8907.00
+ ,9663.00
+ ,8799.00
+ ,8931.00
+ ,8732.00
+ ,8936.00
+ ,9127.00
+ ,9070.00
+ ,9773.00
+ ,9670.00
+ ,9929.00
+ ,10095.00
+ ,9025.00
+ ,9659.00
+ ,8954.00
+ ,9022.00
+ ,8855.00
+ ,9034.00
+ ,9196.00
+ ,9038.00
+ ,9650.00
+ ,9715.00
+ ,10052.00
+ ,10436.00
+ ,9314.00
+ ,9717.00
+ ,8997.00
+ ,9062.00
+ ,8885.00
+ ,9058.00
+ ,9095.00
+ ,9149.00
+ ,9857.00
+ ,9848.00
+ ,10269.00
+ ,10341.00
+ ,9690.00
+ ,10125.00
+ ,9349.00
+ ,9224.00
+ ,9224.00
+ ,9454.00
+ ,9347.00
+ ,9430.00
+ ,9933.00
+ ,10148.00
+ ,10677.00
+ ,10735.00
+ ,9760.00
+ ,10567.00
+ ,9333.00
+ ,9409.00
+ ,9502.00
+ ,9348.00
+ ,9319.00
+ ,9594.00
+ ,10160.00
+ ,10182.00
+ ,10810.00
+ ,11105.00
+ ,9874.00
+ ,10958.00
+ ,9311.00
+ ,9610.00
+ ,9398.00
+ ,9784.00
+ ,9425.00
+ ,9557.00
+ ,10166.00
+ ,10337.00
+ ,10770.00
+ ,11265.00
+ ,10183.00
+ ,10941.00
+ ,9628.00
+ ,9709.00
+ ,9637.00
+ ,9579.00
+ ,9741.00
+ ,9754.00
+ ,10508.00
+ ,10749.00
+ ,11079.00
+ ,11608.00
+ ,10668.00
+ ,10933.00
+ ,9703.00
+ ,9799.00
+ ,9656.00
+ ,9648.00
+ ,9712.00
+ ,9766.00
+ ,10540.00
+ ,10564.00
+ ,10911.00
+ ,11218.00
+ ,10230.00
+ ,10410.00
+ ,9227.00
+ ,9378.00
+ ,9105.00
+ ,9128.00)
+ ,dim=c(1
+ ,175)
+ ,dimnames=list(c('verbruik')
+ ,1:175))
>  y <- array(NA,dim=c(1,175),dimnames=list(c('verbruik'),1:175))
>  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/1ajl01263389234.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]
[1,]  7248
[2,]  8561
[3,]  9227
[4,]  9757
[5,] 11265

$n
[1] 175

$conf
         [,1]
[1,] 9084.154
[2,] 9369.846

$out
[1] 11608

$group
[1] 1

$names
[1] "verbruik"

> 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/2szmx1263389234.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/30y001263389234.tab") 
> 
> try(system("convert tmp/1ajl01263389234.ps tmp/1ajl01263389234.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.374   0.154   0.462