我想在对数线性模型中验证交互项是否在我的模型中起到相关作用。为了说明,我使用源自Field(2013)的下面的代码:P。846-848。我不明白为什么“更新”不按预期工作。该公式按原样采用:
数据:
catdogconttable <- xtabs(~ Animal + Training + Dance , data=Dataset)
catdogconttable
完全饱和的模型:
caturated <- loglm(~ Animal * Training * Dance, data = catdogconttable, fit = TRUE)
summary(caturated)
省略3倍交互项
threeway <- loglm(~ Animal + Training + Dance + Animal : Training +
Animal : Dance + Dance:Training , data = catdogconttable, fit = TRUE)
summary(threeway)
使用update
threeway1<-update(caturated, .~., -Animal * Training * Dance)
summary(threeway1)
总结anova的差异表明得到的结果不是 相同的
anova(caturated, threeway)
anova(caturated, threeway1)
可用数据
问题是为什么threeway1没有产生相同的输出 三通?