Sorry by the delay. You could do: > my.data <- expand.grid(A=factor(1:4), B=factor(1:4), rep=1:4) > my.data$y <- rbinom(my.data$A, 10, 0.5) > > model <- glm(cbind(y, 10-y)~A*B, family=binomial, data=my.data) > anova(model, test="Chisq" ) Analysis of Deviance Table
Model: binomial, link: logit Response: cbind(y, 10 - y) Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev P(>|Chi|) NULL 63 53.528 A 3 1.6792 60 51.849 0.6416 B 3 2.9380 57 48.911 0.4013 A:B 9 13.8380 48 35.073 0.1282 > > X <- model.matrix(~A/B, my.data) > > A <- X[,2:4] > B.A1 <- X[,grep("A1:", colnames(X))] > B.A2 <- X[,grep("A2:", colnames(X))] > B.A3 <- X[,grep("A3:", colnames(X))] > B.A4 <- X[,grep("A4:", colnames(X))] > > model2 <- glm(cbind(y, 10-y)~A+B.A1+B.A2+B.A3+B.A4, family=binomial, > data=my.data) > anova(model2, test="Chisq") Analysis of Deviance Table Model: binomial, link: logit Response: cbind(y, 10 - y) Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev P(>|Chi|) NULL 63 53.528 A 3 1.6792 60 51.849 0.64156 B.A1 3 1.2924 57 50.556 0.73093 B.A2 3 4.3126 54 46.244 0.22962 B.A3 3 1.8865 51 44.357 0.59629 B.A4 3 9.2844 48 35.073 0.02574 * --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > Walmes. ----- ..ooo0 ................................................................................................... ..(....)... 0ooo... Walmes Zeviani ...\..(.....(.....)... Master in Statistics and Agricultural Experimentation ....\_)..... )../.... walmeszevi...@hotmail.com, Lavras - MG, Brasil ............ (_/............................................................................................ -- View this message in context: http://r.789695.n4.nabble.com/splitting-a-factor-in-an-analysis-of-deviance-table-negative-binomial-model-tp1017744p2235345.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org 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.