R version 2.7.0 (2008-04-22)
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(31.58,34.32,27.05,15.70,17.57,21.12,27.69,32.49,16.00,27.88,32.78,26.02,13.94,18.42,21.25,29.54,33.66,14.36,27.32,28.38,25.76,12.23,17.96,21.80,29.31,32.44,14.98,28.89,29.18,25.94,14.75,18.39,21.84,29.26,34.38,16.88,28.05,28.62,24.97,14.62,17.63,22.21,28.69,32.36,16.56,28.73,28.20,21.74,15.04,17.95,22.64,26.16,30.73,13.31,32.00,29.33,18.16,15.50,17.79,23.54,27.12,30.31,9.61,34.53,29.72,16.95,16.10,17.73,23.78,29.40,27.26,9.34,33.47,26.29,16.46,15.44,18.99,23.65,30.99,25.05,7.89,34.09,26.82,16.44,15.14,19.83,23.93,32.96,22.33,1.71,35.47,27.64,18.20,15.42,20.23,24.77,32.20,18.26,0.81,34.59,27.10,16.44,15.69,20.24,26.26,31.67,18.30),dim=c(9,12),dimnames=list(c('2000','2001','2002','2003','2004','2005','2006','2007','2008'),1:12))
>  y <- array(NA,dim=c(9,12),dimnames=list(c('2000','2001','2002','2003','2004','2005','2006','2007','2008'),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 = 'Price'
> xlab = 'Time'
> 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/16hpq1229769538.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]   [,9]
[1,] 27.32 26.290 16.440 13.940 17.570 21.120 26.160 18.260  0.810
[2,] 28.39 27.370 16.705 14.685 17.760 21.820 28.190 23.690  8.615
[3,] 31.79 28.500 19.970 15.280 18.175 23.090 29.355 30.520 13.835
[4,] 34.31 29.525 25.850 15.595 19.410 23.855 31.330 32.465 16.280
[5,] 35.47 29.720 27.050 16.100 20.240 26.260 32.960 34.380 16.880

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 29.08985 27.51709 15.79890 14.86494 17.42242 22.16182 27.92282 26.51766
[2,] 34.49015 29.48291 24.14110 15.69506 18.92758 24.01818 30.78718 34.52234
         [,9]
[1,] 10.33894
[2,] 17.33106

$out
[1] 34.32 32.78 12.23 31.58

$group
[1] 2 2 4 9

$names
[1] "2000" "2001" "2002" "2003" "2004" "2005" "2006" "2007" "2008"

Warning message:
In bxp(list(stats = c(27.32, 28.39, 31.79, 34.31, 35.47, 26.29,  :
  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/21hdr1229769538.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/3xhdi1229769538.tab") 
> 
> system("convert tmp/16hpq1229769538.ps tmp/16hpq1229769538.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.146   0.248   1.201