Displaying 1 result from an estimated 1 matches for "blevel4".
Did you mean:
level4
2010 Jan 19
1
splitting a factor in an analysis of deviance table (negative binomial model)
...re(MASS)
my.data <- data.frame("A"=rep(c(rep("Alevel1",4),rep("Alevel2",4),rep("Alevel3",4),rep("Alevel4",4)),4),
"B"=c(rep("Blevel1",16),rep("Blevel2",16),rep("Blevel3",16),rep("Blevel4",16)),
"value"=rnegbin(64, 10, 10))
# fitting the model with interaction (a + b + a:b)
model <- glm.nb(value ~ A*B, data=my.data)
anova(model, test="F")
Df Deviance Resid. Df Resid. Dev F Pr(>F)
NULL 63 80.639...