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

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

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