x <- array(list(1280
,1024
,1024
,768
,1120
,700
,1024
,768
,1280
,800
,1280
,1024
,1280
,800
,1024
,768
,1280
,800
,1280
,1024
,1280
,800
,1280
,800
,1280
,1024
,1688
,949
,1440
,900
,1600
,1200
,1280
,800
,1280
,800
,1280
,768
,1176
,735
,1280
,800
,1503
,845
,1440
,900
,1366
,768
,1280
,768
,1024
,768
,1280
,800
,2560
,1440
,1280
,768
,1024
,768
,1280
,1024
,1280
,800
,1440
,900
,1280
,800
,1440
,900
,1024
,768
,1440
,900
,1143
,857
,1280
,800
,1440
,900
,1280
,800
,1366
,768
,1024
,768
,1408
,880
,1366
,768
,1176
,735
,1920
,1200
,1257
,785
,1280
,800
,1280
,800
,1440
,900
,1680
,1050
,1440
,900
,1024
,768
,1140
,641
,1280
,1024
,1280
,800
,1280
,800
,1280
,800
,1280
,800
,1440
,900
,1280
,800
,1152
,864
,1280
,1024
,1280
,800
,1440
,900
,1280
,800
,1280
,1024
,1440
,900
,1280
,800
,1280
,800
,1440
,900
,1280
,800
,1280
,1024
,1600
,900
,1024
,768
,1366
,768
,1280
,800
,1280
,800
,1440
,900
,1366
,768
,1280
,800
,1024
,768
,1280
,800
,1440
,900
,1280
,800
,1280
,800
,1408
,880
,1280
,800
,1600
,900
,1600
,900
,1680
,1050
,1440
,900
,1440
,900
,917
,550
,1280
,800
,1760
,990
,1280
,800
,1280
,800
,1280
,800
,1024
,768
,1366
,768
,1440
,900
,1280
,800
,1280
,1024
,1920
,1080
,1024
,768
,1024
,768
,1600
,900
,1117
,698
,1440
,900
,983
,737
,1024
,768
,1024
,640
,1280
,800
,1440
,900
,1280
,800
,1280
,800
,1280
,800
,1440
,900
,1280
,800
,1024
,768
,1024
,768
,1152
,864
,1280
,768
,1024
,768
,1366
,768
,1680
,1050
,1680
,1050
,1280
,800
,1366
,768
,1024
,768
,1440
,900
,1024
,768
,1280
,800
,1280
,800
,1280
,800
,1024
,768
,1280
,800)
,dim=c(2
,139)
,dimnames=list(c('Width'
,'Height')
,1:139))
 y <- array(NA,dim=c(2,139),dimnames=list(c('Width','Height'),1:139))
 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'
z <- as.data.frame(t(y))
z <- z[z[,1]*z[,2]>800000 & z[,1]*z[,2]<1700000,]
postscript(file="/var/wessaorg/rcomp/tmp/1kli91324254013.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/2z6mo1324254013.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/wessaorg/rcomp/tmp/3090v1324254013.tab") 

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

