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

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

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