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 Response:

[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