I guess your model specification is incoherent. It's seems to me that you incorporated some kind of lmer() style with "(0 + Micro/Trt/Year)" term. Consider the toy data and code below. Nesting model is m2. Models m1 and m3 are the same so showing that terms in brackets don't are used.
da <- expand.grid(A=factor(1:3), B=factor(1:3))[rep(1:9, 3),] da$y <- rnbinom(da$A, size=10, prob=0.5) m1 <- glm.nb(y~A*B, data=da) anova(m1) m2 <- glm.nb(y~A/B, data=da) anova(m2) m3 <- glm.nb(y~A*B+(A/B), data=da) anova(m3) Bests. Walmes Zeviani. ----- ..ooo0 ................................................................................................... ..(....)... 0ooo... Walmes Zeviani ...\..(.....(.....)... Master in Statistics and Agricultural Experimentation ....\_)..... )../.... [email protected], Lavras - MG, Brasil ............ (_/............................................................................................ -- View this message in context: http://n4.nabble.com/Correct-nested-design-for-GLM-tp1576038p1576483.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
