data(brokentrans)
# For \alpha = 10^{-9}, we obtain significant difference
# between 1-2, 2-3, but not 1-3.
tapply(brokentrans$x, brokentrans$g, mean)
pairwise.t.test(brokentrans$x, brokentrans$g, pool.sd=FALSE)
Run the code above in your browser using DataLab