data(detergent)
N <- rep(2, 5)
# BIBD: prepare the contrast matrix = all-pair comparisons for
# the 5 levels of detergent
C <- contrMat(N, type="Tukey")
# adjusted p-values: either using the contrast matrix C
summary(simtest(plates ~ block+detergent, data=detergent, whichf="detergent",
ttype="logical", cmatrix=C))
# or, more easily, type="Tukey"
summary(simtest(plates ~ detergent + block, data=detergent,
whichf="detergent", type="Tukey", ttype="logical"))
Run the code above in your browser using DataLab