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

82	79	88	75	77
80	78	86	80	77
91	81	84	80	78
84	86	84	84	79
91	84	87	84	80
79	80	88	89	81
92	82	84	96	81
80	79	84	81	82
83	86	93	77	82
82	84	86	90	82
88	91	95	83	82
90	96	95	84	82
94	78	89	84	82
88	80	88	87	83
86	86	83	89	83
83	77	81	90	83
91	79	81	73	83
89	82	86	90	83
87	87	94	94	83
89	86	85	98	84
85	82	105	105	84
88	95	87	78	84
91	102	95	81	84
102	80	85	82	84
100	92	102	88	85
97	82	103	106	86
85	86	101	73	86
88	89	91	84	86
85	92	83	91	86
85	81	85	92	87
81	97	94	87	88
86	91	84	90	88
89	95	86	91	88
102	81	91	104	89
95	103	88	91	89
98	95	102	95	89
98	96	94	97	89
91	86	92	100	89
91	92	105	100	89
93	96	90	83	90
88	105	83	84	90
97	82	90	88	90
100	93	93	89	90
93	88	93	89	90
95	95	84	85	91
96	104	103	88	91
86	88	88	94	91
103	89	86	86	92
88	91	93	87	92
98	91	92	89	92
102	92	92	92	92
97	105	100	96	92
91	99	95	105	92
81	94	88	107	92
90	80	86	73	93
98	88	83	79	93
97	89	101	94	93
100	95	102	94	93
94	88	90	95	93
88	105	104	103	93
89	82	86	87	94
91	92	88	90	94
95	86	87	98	94
95	90	88	89	95
102	97	101	105	95
92	103	97	97	96
98	86	88	98	96
104	98	87	84	96
98	92	87	90	97
96	92	103	83	97
98	97	99	89	97
94	91	87	90	98
97	98	97	94	98
107	90	89	96	98
94	105	91	105	98
91	90	104	86	99
106	85	88	94	99
98	105	93	105	100
106	94	95	96	100
92	92	92	102	100
104	100	105	84	100
103	103	96	96	101
102	100	94	101	102
106	100	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(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')