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

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

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