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

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

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