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

88	95	95	81	84
88	92	87	90	97
88	95	91	91	88
88	97	84	87	87
88	92	92	89	92
88	104	84	85	91
88	79	84	81	80
88	78	69	80	73
88	79	84	75	73
88	86	84	84	77
88	92	85	91	86
88	78	89	87	82
88	86	83	90	82
88	71	86	72	73
88	90	91	96	99
88	93	90	88	90
88	80	77	73	78
88	90	88	86	100
88	77	88	67	73
88	86	88	89	95
88	89	102	94	93
88	88	103	88	91
88	86	84	77	81
88	97	87	89	98
88	81	88	104	88
88	88	93	89	90
88	82	101	107	85
88	82	88	96	79
88	100	113	111	107
88	91	95	84	82
88	111	99	83	98
88	102	85	82	84
88	105	105	105	99
88	92	105	100	89
88	100	96	111	105
88	82	88	90	94
88	91	93	87	92
88	71	93	90	81
88	71	69	67	NA
88	89	83	84	86
88	82	83	84	90
88	105	93	96	101
88	71	69	67	73
88	91	86	90	88
88	105	90	83	90
88	82	87	78	84
88	86	92	100	89
88	72	75	67	73
88	98	89	94	98
88	80	88	89	82
88	84	87	89	77
88	109	113	111	107
88	80	102	88	84
88	112	101	94	93
88	86	87	111	96
88	112	105	111	107
88	90	101	105	95
88	96	79	97	89
88	84	86	83	82
88	94	95	105	92
88	112	104	69	100
88	81	79	80	75
88	105	86	87	94
88	77	81	67	83
88	105	92	92	92
88	112	95	102	102
88	91	97	90	98
88	92	87	98	94
88	72	94	94	83
88	100	105	101	107
88	103	94	96	107
88	89	88	94	91
88	86	105	105	83
88	112	113	105	100
88	71	69	84	77
88	73	69	72	74
88	99	100	96	92
88	85	113	94	100
88	80	88	107	92
88	88	104	103	93
88	95	94	95	89
88	86	91	73	86
88	92	103	110	97
88	96	95	84	82
88	71	69	67	73
88	88	86	73	93
88	98	108	84	97
88	108	83	79	93
88	103	102	91	89
88	81	94	92	86
88	96	113	107	89
88	112	113	111	107
88	95	93	89	90
88	92	105	110	104
88	103	88	98	96
88	95	90	95	93
88	94	92	84	104
88	79	81	73	83
88	91	86	86	92
88	87	85	98	83
88	82	86	90	83
88	97	97	97	96
88	92	103	106	84
88	112	102	106	107

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')