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

79	77	81	73	77
80	78	81	73	77
80	78	83	73	78
81	79	83	75	79
81	79	83	77	80
82	79	83	78	81
82	80	84	79	81
83	80	84	80	82
83	80	84	80	82
84	80	84	81	82
85	81	84	81	82
85	81	84	82	82
85	81	85	83	82
85	82	85	83	83
86	82	85	83	83
86	82	86	84	83
86	82	86	84	83
87	82	86	84	83
88	82	86	84	83
88	84	86	84	84
88	84	86	84	84
88	85	86	84	84
88	86	87	84	84
88	86	87	85	84
88	86	87	86	85
89	86	87	86	86
89	86	87	87	86
89	86	87	87	86
89	86	88	87	86
90	86	88	87	87
90	87	88	88	88
91	88	88	88	88
91	88	88	88	88
91	88	88	89	89
91	88	88	89	89
91	89	88	89	89
91	89	88	89	89
91	89	88	89	89
91	90	89	89	89
91	90	89	89	90
92	90	90	90	90
92	91	90	90	90
92	91	90	90	90
93	91	91	90	90
93	91	91	90	91
94	91	91	90	91
94	92	92	90	91
94	92	92	91	92
94	92	92	91	92
95	92	92	91	92
95	92	93	92	92
95	92	93	92	92
95	92	93	94	92
96	92	93	94	92
96	93	93	94	93
97	94	94	94	93
97	94	94	94	93
97	95	94	94	93
97	95	94	95	93
97	95	95	95	93
98	95	95	96	94
98	95	95	96	94
98	96	95	96	94
98	96	95	96	95
98	96	96	96	95
98	97	97	97	96
98	97	97	97	96
98	97	99	98	96
100	98	100	98	97
100	98	101	98	97
100	99	101	100	97
102	100	101	100	98
102	100	102	101	98
102	100	102	102	98
102	102	102	103	98
102	103	102	104	99
103	103	103	105	99
103	103	103	105	100
104	104	103	105	100
104	105	104	105	100
106	105	104	105	100
106	105	105	106	101
106	105	105	106	102
107	105	105	107	104

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