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

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

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