search for: letters_mod

...# combined variable for pairwise comparisons treat.seas <- (as.numeric(dat$treat) * 10) + as.numeric(dat$seas) dat$treat.seas <- factor(treat.seas) mod.2 <- glm(response ~ 1 + treat.seas + covariate, data = dat) mod.2_glht <- glht(mod.2, linfct = mcp(treat.seas = "Tukey")) letters_mod.2 <- cld(mod.2_glht) Thank you in advance Melanie -- View this message in context: http://r.789695.n4.nabble.com/glht-multiple-comparisons-for-glm-with-2-factors-tp4633365.html Sent from the R help mailing list archive at Nabble.com.