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

2	88	95	95	81
3	98	92	87	90
4	89	95	91	91
5	81	97	84	87
6	88	92	92	89
7	95	104	84	85
8	79	79	84	81
9	91	78	69	80
10	80	79	84	75
11	91	86	84	84
12	85	92	85	91
13	94	78	89	87
14	86	86	83	90
15	82	71	86	72
16	107	90	91	96
17	97	93	90	88
18	78	80	77	73
19	91	90	88	86
20	74	77	88	67
21	95	86	88	89
22	97	89	102	94
23	112	88	103	88
24	80	86	84	77
25	98	97	87	89
26	102	81	88	104
27	100	88	93	89
28	97	82	101	107
29	92	82	88	96
30	114	100	113	111
31	88	91	95	84
32	96	111	99	83
33	91	102	85	82
34	114	105	105	105
35	91	92	105	100
36	114	100	96	111
37	89	82	88	90
38	103	91	93	87
39	83	71	93	90
40	74	71	69	67
41	88	89	83	84
42	88	82	83	84
43	98	105	93	96
44	74	71	69	67
45	86	91	86	90
46	93	105	90	83
47	85	82	87	78
48	91	86	92	100
49	75	72	75	67
50	97	98	89	94
51	88	80	88	89
52	79	84	87	89
53	106	109	113	111
54	102	80	102	88
55	98	112	101	94
56	98	86	87	111
57	114	112	105	111
58	95	90	101	105
59	98	96	79	97
60	82	84	86	83
61	97	94	95	105
62	94	112	104	69
63	84	81	79	80
64	88	105	86	87
65	83	77	81	67
66	98	105	92	92
67	106	112	95	102
68	94	91	97	90
69	91	92	87	98
70	89	72	94	94
71	102	100	105	101
72	103	103	94	96
73	96	89	88	94
74	89	86	105	105
75	112	112	113	105
76	74	71	69	84
77	74	73	69	72
78	102	99	100	96
79	106	85	113	94
80	91	80	88	107
81	94	88	104	103
82	98	95	94	95
83	85	86	91	73
84	108	92	103	110
85	90	96	95	84
86	74	71	69	67
87	81	88	86	73
88	104	98	108	84
89	90	108	83	79
90	95	103	102	91
91	85	81	94	92
92	114	96	113	107
93	114	112	113	111
94	93	95	93	89
95	104	92	105	110
96	92	103	88	98
97	100	95	90	95
98	92	94	92	84
99	77	79	81	73
100	86	91	86	86
101	87	87	85	98
102	91	82	86	90
103	102	97	97	97
104	100	92	103	106
105	113	112	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')