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")
# the additional 10 columns of are for the 10 blocks
C <- cbind( matrix(0, ncol=10, nrow=10), C )
# numerate the contrasts
colnames(C) <- NULL
rownames(C) <- paste("C", 1:nrow(C), sep="")
# adjusted p-values
summary(simtest(plates ~ block+detergent, data=detergent,
ttype="logical", cmatrix=C))
Run the code above in your browser using DataLab