x <- array(list(41
,34
,38
,30
,39
,33
,32
,28
,30
,29
,35
,31
,31
,34
,33
,35
,34
,32
,37
,35
,35
,35
,29
,37
,39
,41
,31
,39
,34
,27
,36
,31
,36
,40
,35
,38
,37
,40
,38
,37
,38
,36
,31
,37
,36
,40
,34
,35
,38
,43
,35
,37
,39
,40
,38
,42
,33
,33
,37
,28
,32
,37
,33
,37
,36
,32
,32
,36
,38
,26
,38
,37
,39
,36
,38
,33
,32
,39
,32
,40
,32
,38
,33
,30
,31
,34
,31
,36
,39
,35
,38
,33
,37
,41
,39
,40
,39
,42
,32
,37
,41
,36
,32
,37
,36
,39
,35
,39
,33
,33
,37
,35
,33
,33
,33
,36
,34
,36
,33
,34
,31
,37
,28
,36
,27
,36
,32
,32
,37
,34
,31
,33
,34
,32
,37
,27
,34
,35
,30
,37
,32
,39
,33
,32
,29
,30
,31
,31
,36
,25
,33
,31
,29
,29
,31
,32
,35
,39
,33
,37
,37
,31
,32
,25
,34
,26
,33
,30
,38
,28
,32
,37
,35
,40
,33
,37
,38
,32
,28
,40
,37
,35
,35
,35
,38
,32
,39
,35
,33
,41
,34
,43
,36
,34
,38
,32
,38
,36
,32
,42
,32
,38
,38
,35
,32
,34
,30
,27
,32
,32
,33
,30
,34
,34
,38
,31
,32
,32
,32
,36
,37
,40
,32
,36
,39
,43
,34
,41
,29
,35
,34
,34
,37
,45
,36
,36
,35
,36
,34
,33
,30
,39
,28
,35
,38
,31
,34
,28
,34
,36
,35
,33
,31
,36
,35
,38
,34
,37
,31
,37
,35
,40
,37
,39
,36
,35
,35
,34
,30
,36
,27
,32
,39
,32
,40
,36
,35
,36
,37
,36
,38
,37
,36
,35
,31
,42
,38
,33
,34
,37
,39
,42
,38
,36
,41
,36
,34
,36
,27
,33
,39
,33
,30
,36
,37
,37
,37
,32
,34
,35
,31
,35
,28
,37
,31
,38
,37
,28
,27
,33
,33
,33
,36
,32
,37
,45
,38
,38
,35
,38
,37
,39
,37
,43
,33
,39
,32
,37
,34
,39
,33
,36
,31
,30
,38
,40
,39
,38
,33
,39
,34
,38
,29
,43
,32
,42
,33
,32
,33
,41
,31
,37
,36
,31
,36
,34
,32
,39
,35
,44
,41
,40
,32
,35
,28
,31
,29
,34
,30
,34
,39
,37
,36
,23
,37
,40
,35
,28
,35
,36
,31
,36
,37
,44
,34
,41
,32
,35
,36
,29
,38
,34
,36
,31
,37
,40
,35
,33
,36
,34
,37
,35
,32
,39
,28
,35
,33
,36
,39
,34
,40
,40
,32
,40
,38
,37
,35
,34
,41
,35
,39
,36
,36
,45
,35
,35
,43
,39
,42
,39)
,dim=c(4
,110)
,dimnames=list(c('Y2010C'
,'UC'
,'Y2010S'
,'US')
,1:110))
 y <- array(NA,dim=c(4,110),dimnames=list(c('Y2010C','UC','Y2010S','US'),1:110))
 for (i in 1:dim(x)[1])
 {
 	for (j in 1:dim(x)[2])
 	{
 		y[i,j] <- as.numeric(x[i,j])
 	}
 }
par20 = ''
par19 = ''
par18 = ''
par17 = ''
par16 = ''
par15 = ''
par14 = ''
par13 = ''
par12 = ''
par11 = ''
par10 = ''
par9 = ''
par8 = ''
par7 = ''
par6 = ''
par5 = ''
par4 = ''
par3 = ''
par2 = ''
par1 = 'grey'
ylab = 'value'
xlab = 'variables'
main = 'Notched Boxplots'
z <- as.data.frame(t(y))
postscript(file="/var/fisher/rcomp/tmp/1e6oj1349775497.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
(r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=TRUE,col=par1))
dev.off()

#Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/fisher/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/fisher/rcomp/tmp/2a3p71349775497.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/fisher/rcomp/tmp/38eid1349775498.tab") 

try(system("convert tmp/1e6oj1349775497.ps tmp/1e6oj1349775497.png",intern=TRUE))