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

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

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