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

89	81	109
90	91	88
82	96	87
82	85	78
71	74	72
80	78	85
62	71	73
77	72	77
80	87	90
87	72	78
96	90	96
80	87	84
85	89	87
89	110	105
89	95	89
73	72	78
88	70	96
86	79	79
83	92	90
91	93	96
98	97	95
80	75	79
84	88	86
90	105	99
92	82	111
86	89	105
117	114	123
88	81	102
86	73	82
78	82	86
104	110	117
89	93	94
99	114	110
84	91	86
89	97	109
95	97	118
58	68	64
77	72	78
77	81	86
98	102	98
58	62	64
81	102	84
92	84	87
94	75	88
87	102	106
75	62	78
80	100	79
84	82	95
87	78	98
114	106	101
97	88	104
102	85	91
80	97	78
112	114	123
104	105	109
75	100	77
82	72	77
92	106	102
102	78	88
73	80	81
78	88	83
74	62	72
92	93	95
97	87	109
92	109	91
91	91	87
102	110	112
89	98	95
98	99	103
117	95	108
92	97	116
70	88	71
63	72	70
91	102	105
106	95	112
81	93	86
98	90	109
88	80	112
95	75	97
105	109	122
90	75	97
58	72	64
86	77	84
112	93	110
71	72	70
98	91	100
87	95	102
98	111	108
117	105	123
98	90	78
98	105	109
90	90	106
84	92	95
86	95	108
68	75	74
86	92	86
77	100	82
83	82	88
90	102	96
107	104	114
102	98	111
90	95	88
113	114	123
86	93	89
105	110	114
110	104	121
82	91	100
89	79	83
99	97	92
113	102	108
91	77	87
92	105	84
86	95	95
91	104	90
82	75	82
86	82	84
86	91	96
102	105	109
84	71	84
80	89	83
96	79	93
117	113	120
88	82	83
99	102	103
89	82	103
86	92	102
98	82	90
58	86	69
98	108	114
58	74	64
81	87	87
90	97	116
112	102	121
84	75	90
84	82	92
75	72	81
78	88	84
94	97	106
75	85	75
82	62	79
87	95	80
59	95	74
68	72	66
80	92	81
94	91	98
87	97	99
70	78	68
113	95	100
99	91	97
90	95	99
102	88	104
80	92	82
98	97	92
98	104	99
95	113	109
104	113	106
86	82	89
84	90	87
85	91	98
90	95	95
86	102	92
97	92	90

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