data(litter)
# define the contrast matrix
C <- matrix(c(0,0,0,0,0,0,1,1,1,.75,.384,.887,-1,0,0,.25,.37,.113,
0,-1,0,-.25,.246,-.339,0,0,-1,-.75,-1,-.661,0,0,0,
0,0,0,0,0,0,0,0,0), ncol=7)
# numerate the contrasts
rownames(C) <- paste("C", 1:nrow(C), sep="")
# simultaneous confidence intervals for above comparisons
# in an ANCOVA
summary(simint(weight ~ dose + gesttime + number, data=litter,
alternative="greater", cmatrix=C))
Run the code above in your browser using DataLab