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

58	51	77	77
53	54	68	70
59	59	76	76
166	56	76	77
51	52	69	73
64	64	71	71
52	57	65	64
65	66	70	75
62	62	92	101
61	61	76	75
61	61	119	124
54	59	65	66
50	50	66	70
63	61	101	100
58	60	75	73
39	41	79	76
71	71	64	65
52	52	69	69
68	63	88	86
56	54	65	67
54	53	80	80
63	59	78	80
54	55	85	82
54	56	82	85
75	75	74	73
56	57	102	107
78	75	65	64
60	59	73	74
64	63	75	70
64	62	57	58
52	51	68	69
62	61	71	71
55	54	71	76
56	57	97	98
50	50	80	76
50	55	66	66
63	64	69	70
61	60	69	70
53	52	55	56
60	55	59	61
56	56	62	66
53	53	70	68
57	59	84	86
57	56	69	71
56	56	88	87
56	57	103	101
50	50	63	63
52	52	84	90
55	55	79	79
47	47	67	67
45	45	83	83
62	63	96	94
53	51	75	76
52	51	65	66
57	55	78	77
64	64	69	73
59	55	87	89
55	57	83	84
76	77	90	91
62	62	85	83
68	68	66	68
55	55	88	86
45	45	54	58
68	68	69	68
44	44	56	58
62	61	96	95
56	53	76	75
50	47	61	61
53	53	62	64
64	62	71	68
52	53	66	67
53	55	81	82
54	55	68	68
64	66	80	78
55	55	70	70
55	55	76	75
59	55	88	93
70	67	89	86
57	58	74	71
47	47	83	80
47	45	90	91
48	44	79	81

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,signif(r$stats[j,i], digits=4))
}
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,signif(r$conf[1,i], digits=4))
a<-table.element(a, signif(r$stats[3,i], digits=4))
a<-table.element(a,signif(r$conf[2,i], digits=4))
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,signif(mean(z[i], trim=par3/100, na.rm=TRUE), digits=4))
a<-table.element(a,signif(sd(z[i], na.rm=TRUE), digits=4))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')