[R-sig-ME] glht hypothesis test, full cross, 2 treatments.

[Colin Wahl]

I would like to compare between all levels of each treatment cross. I
have two treatments, one with 4 levels and the other with two levels.
I am using glmer with binomial errors.
My code is as follows:
glht(modelEPT, c1)

glht(modelEPT, linfct=mcp(wsh="Tukey"))

This gives me contrasts for the main effects (4 levels)

c1 <- rbind ("C: F vs. NF Rip Effect" = c(0,0,0,0,1,0,0,0),
     "D: F vs. NF Rip Effect" = c(0,0,0,0,1,1,0,0),

     "F: F vs. NF Rip Effect" = c(0,0,0,0,1,0,1,0),

     "G: F vs. NF Rip Effect" = c(0,0,0,0,1,0,0,1))

summary(glht(modelEPT, c1))





--
CW