Displaying 1 result from an estimated 1 matches for "blevel1".
Did you mean:
level1
2010 Jan 19
1
splitting a factor in an analysis of deviance table (negative binomial model)
...B within A4
Here is a code that illustrates my problem:
# inputing my data
require(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 D...