Re: [R] splitting a factor in an analysis of deviance table (negative binomial model)

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

[R] splitting a factor in an analysis of deviance table (negative binomial model)

Dears useRs, I have 2 factors, (for the sake of explanation - A and B), with 4 levels each. I've already fitted a negative binomial generalized linear model to my data, and now I need to split the factors in two distinct analysis of deviance table:  - A within B1, A within B2, A within B3 and A