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

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

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