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