Dunnett

# Data input

> a<-c(6,8,8,7,7,11,8,10,10,8)

> b<-c(7,6,6,5,6,8,7,7,7,7)

> c<-c(5,5,6,4,5,6,6,6,4,4,5,5,5,5)

> d<-c(17,15,15,18,14)

> e<-c(10,9,10,11,14,17,14,12,15,17,18,19,18,20,16,15,17,16)

# boxplot

> boxplot(a,b,c,d,e, names=c("a","b","c","d","e"))

# Data input

> data<-c(6,8,8,7,7,11,8,10,10,8,7,6,6,5,6,8,7,7,7,7,5,5,6,4,5,6,6,6,4,4,5,5,5,5,17,15,15,18,14,10,9,10,11,14,17,14,12,15,17,18,19,18,20,16,15,17,16)

> group <- rep(1:5, c(10,10,14,5,18))

# Data output

> data

[1] 6 8 8 7 7 11 8 10 10 8 7 6 6 5 6 8 7 7 7 7 5 5 6 4 5 6 6 6 4 4 5 5 5

[34] 5 17 15 15 18 14 10 9 10 11 14 17 14 12 15 17 18 19 18 20 16 15 17 16

> group

[1] 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 5 5 5 5 5 5 5 5 5 5 5

[51] 5 5 5 5 5 5 5

# confirm the data, by comparing this plot with the boxplot you draw above.

> plot(group,data)

# open package

# *please install the "mvtnorm" package before use.

> library(mvtnorm)

# dunnet

> dunnett(data, group)

t p

1:2 1.799640 2.272777e-01

1:3 3.691650 1.993688e-03

1:4 6.482647 8.268556e-08

1:5 7.908992 9.221490e-11