• QR Code (size=2)
• QR Code (size=3)
• QR Code (size=4)
• QR Code (size=5)
• QR Code (size=6)

1	36	88
2	56	94
2	48	90
2	32	73
1	44	68
2	39	80
2	34	86
3	41	86
3	50	91
1	39	79
3	62	96
2	52	92
3	37	72
2	50	96
1	41	70
2	55	86
2	41	87
3	56	88
2	39	79
1	52	90
2	46	95
2	44	85
2	41	90
3	50	115
3	50	84
2	44	79
1	52	94
2	54	97
2	44	86
3	52	111
2	37	87
3	52	98
3	50	87
1	36	68
1	50	88
3	52	82
3	55	111
2	31	75
1	36	94
1	49	95
1	42	80
2	37	95
2	41	68
1	30	94
1	52	88
3	30	84
1	44	101
2	66	98
3	48	78
2	43	109
2	57	102
1	46	81
3	54	97
3	48	75
2	48	97
1	62	101
3	58	101
2	58	95
2	62	95
2	46	95
1	34	90
2	66	107
3	52	92
2	55	86
1	55	70
3	57	95
1	56	96
2	55	91
3	56	87
1	54	92
3	55	97
2	46	102
1	52	91
2	32	68
1	44	88
2	46	97
2	59	90
3	46	101
3	46	94
3	54	101
3	66	109
2	56	100
2	59	103
2	57	94
3	52	97
1	48	85
1	44	75
2	41	77
1	50	87
3	48	78
2	48	108
2	59	97
2	46	106
2	54	107
2	55	95
3	54	107
2	59	115
2	44	101
3	54	85
3	52	90
3	66	115
2	44	95
2	57	97
1	39	112
3	60	97
2	45	77
2	41	90
2	50	94
2	39	103
2	43	77
1	48	98
2	37	90
2	58	111
1	46	77
1	43	88
2	44	75
3	34	92
1	30	78
3	50	106
1	39	80
2	37	87
2	55	92
1	39	86
3	36	85
2	43	90
3	50	101
2	55	94
2	43	86
3	60	86
2	48	90
3	30	75
2	43	86
1	39	91
2	52	97
1	39	91
1	39	70
1	56	98
1	59	96
2	46	95
2	57	100
2	50	95
1	54	97
3	50	97
3	60	92
3	59	115
2	41	88
1	48	87
2	59	100
3	60	98
2	56	102
1	51	96

par1 <- as.numeric(par1) #colour
par2<- as.logical(par2) # Notches ?
par3<-as.numeric(par3) # % trim
if(par3>45){par3<-45;warning('trim limited to 45%')}
if(par3<0){par3<-0;warning('negative trim makes no sense. Trim is zero.')}
lotrm<-as.integer(length(y[1,])*par3/100)+1
hitrm<-as.integer(length(y[1,])*(100-par3)/100)
y1<-array(dim=c(dim(y)[1], hitrm-lotrm+1), dimnames=list(dimnames(y)[[1]], 1:(hitrm-lotrm+1) ))
for(i in 1:dim(y)[1]){
tmp<-order(y[i,])
y1[i,]<- y[i, tmp[lotrm:hitrm] ]
}
bitmap(file='test2.png')
pairs(t(y))
dev.off()
y<-y1
z <- as.data.frame(t(y))
bitmap(file='test1.png')
(r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=par2,col=par1))
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,hyperlink('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('lower_whisker.htm','lower whisker','definition of lower whisker'),1,TRUE)
a<-table.element(a,hyperlink('lower_hinge.htm','lower hinge','definition of lower hinge'),1,TRUE)
a<-table.element(a,hyperlink('central_tendency.htm','median','definitions about measures of central tendency'),1,TRUE)
a<-table.element(a,hyperlink('upper_hinge.htm','upper hinge','definition of upper hinge'),1,TRUE)
a<-table.element(a,hyperlink('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,round(r$stats[j,i], digits=2))
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
if (par2){
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,round(r$conf[1,i], digits=2))
a<-table.element(a, round(r$stats[3,i], digits=2))
a<-table.element(a,round(r$conf[2,i], digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
}
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Boxplot Means',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,hyperlink('trimmed_mean.htm','trimmed mean','definition of trimmed mean'),1,TRUE)
a<-table.element(a,hyperlink('unbiased1.htm','unbiased SD','definition of unbiased SD'),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,round(mean(unlist(z[i]), trim=par3/100, na.rm=TRUE), digits=2))
a<-table.element(a,round(sd(unlist(z[i]), na.rm=TRUE), digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')