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

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

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(z[i], trim=par3/100, na.rm=TRUE), digits=2))
a<-table.element(a,round(sd(z[i], na.rm=TRUE), digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')