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(108.5
+ ,108.2
+ ,106.2
+ ,110.8
+ ,112.3
+ ,108.8
+ ,107.1
+ ,119.3
+ ,116.6
+ ,110.2
+ ,108.5
+ ,128.1
+ ,115.5
+ ,109.5
+ ,106.5
+ ,127.6
+ ,120.1
+ ,109.5
+ ,108.3
+ ,137.9
+ ,132.9
+ ,116
+ ,125.6
+ ,151.4
+ ,128.1
+ ,111.2
+ ,124
+ ,143.6
+ ,129.3
+ ,112.1
+ ,127.2
+ ,143.4
+ ,132.5
+ ,114
+ ,136.9
+ ,141.9
+ ,131
+ ,119.1
+ ,135.8
+ ,135.2
+ ,124.9
+ ,114.1
+ ,124.3
+ ,133.1
+ ,120.8
+ ,115.1
+ ,115.4
+ ,129.6
+ ,122
+ ,115.4
+ ,113.6
+ ,134.1
+ ,122.1
+ ,110.8
+ ,114.4
+ ,136.8
+ ,127.4
+ ,116
+ ,118.4
+ ,143.5
+ ,135.2
+ ,119.2
+ ,117
+ ,162.5
+ ,137.3
+ ,126.5
+ ,116.5
+ ,163.1
+ ,135
+ ,127.8
+ ,115.4
+ ,157.2
+ ,136
+ ,131.3
+ ,113.6
+ ,158.8
+ ,138.4
+ ,140.3
+ ,117.4
+ ,155.4
+ ,134.7
+ ,137.3
+ ,116.9
+ ,148.5
+ ,138.4
+ ,143
+ ,116.4
+ ,154.2
+ ,133.9
+ ,134.5
+ ,111.1
+ ,153.3
+ ,133.6
+ ,139.9
+ ,110.2
+ ,149.4
+ ,141.2
+ ,159.3
+ ,118.9
+ ,147.9
+ ,151.8
+ ,170.4
+ ,131.8
+ ,156
+ ,155.4
+ ,175
+ ,130.6
+ ,163
+ ,156.6
+ ,175.8
+ ,138.3
+ ,159.1
+ ,161.6
+ ,180.9
+ ,148.4
+ ,159.5
+ ,160.7
+ ,180.3
+ ,148.7
+ ,157.3
+ ,156
+ ,169.6
+ ,144.3
+ ,156.4
+ ,159.5
+ ,172.3
+ ,152.5
+ ,156.6
+ ,168.7
+ ,184.8
+ ,162.9
+ ,162.4
+ ,169.9
+ ,177.7
+ ,167.2
+ ,166.8
+ ,169.9
+ ,184.6
+ ,166.5
+ ,162.6
+ ,185.9
+ ,211.4
+ ,185.6
+ ,168.1
+ ,190.8
+ ,215.3
+ ,193.2
+ ,171.3
+ ,195.8
+ ,215.9
+ ,207.8
+ ,171
+ ,211.9
+ ,244.7
+ ,223.4
+ ,178.7
+ ,227.1
+ ,259.3
+ ,246.4
+ ,187.5
+ ,251.3
+ ,289
+ ,266.3
+ ,211.4
+ ,256.7
+ ,310.9
+ ,264.3
+ ,211.6
+ ,251.9
+ ,321
+ ,255.8
+ ,199.4
+ ,251.2
+ ,315.1
+ ,259.5
+ ,198.7
+ ,270.3
+ ,333.2
+ ,289.1
+ ,209.3
+ ,267.2
+ ,314.1
+ ,288.5
+ ,215.5
+ ,243
+ ,284.7
+ ,243.8
+ ,212.8
+ ,229.9
+ ,273.9
+ ,224.8
+ ,203.2
+ ,187.2
+ ,216
+ ,177
+ ,175.7
+ ,178.2
+ ,196.4
+ ,171.6
+ ,171
+ ,175.2
+ ,190.9
+ ,165.4
+ ,172.5
+ ,192.4
+ ,206.4
+ ,191.1
+ ,183.5
+ ,187
+ ,196.3
+ ,181.1
+ ,185.4
+ ,184
+ ,199.5
+ ,174.9
+ ,180.9
+ ,194.1
+ ,198.9
+ ,198.1
+ ,187.3
+ ,212.7
+ ,214.4
+ ,223.5
+ ,202
+ ,217.5
+ ,214.2
+ ,236.4
+ ,203.3
+ ,200.5
+ ,187.6
+ ,205.1
+ ,205.8
+ ,205.9
+ ,180.6
+ ,208.5
+ ,221.5
+ ,196.5
+ ,172.2
+ ,181.9
+ ,226.6
+ ,206.3
+ ,187.2
+ ,186.9
+ ,236.7)
+ ,dim=c(4
+ ,61)
+ ,dimnames=list(c('Totaal'
+ ,'Graan'
+ ,'Olie'
+ ,'Suiker')
+ ,1:61))
>  y <- array(NA,dim=c(4,61),dimnames=list(c('Totaal','Graan','Olie','Suiker'),1:61))
>  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: Wessa P., (2007), Notched Boxplots (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_notchedbox1.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description
> z <- as.data.frame(t(y))
> postscript(file="/var/www/html/rcomp/tmp/150ap1260375765.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,] 108.5 108.2 106.2 110.8
[2,] 133.6 119.2 117.0 148.5
[3,] 160.7 175.8 148.7 162.5
[4,] 196.5 211.4 198.1 187.3
[5,] 270.3 333.2 289.1 236.7

$n
[1] 61 61 61 61

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 147.9754 157.1481 132.2936 154.6508
[2,] 173.4246 194.4519 165.1064 170.3492

$out
numeric(0)

$group
numeric(0)

$names
[1] "Totaal" "Graan"  "Olie"   "Suiker"

> 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/2wpa51260375765.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/31mj71260375765.tab") 
> 
> system("convert tmp/150ap1260375765.ps tmp/150ap1260375765.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.396   0.168   0.453