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

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

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