我正在拟合关于缺席 - 存在数据的模型,我想检查随机因素是否重要。 要做到这一点,我应该将glmm与glm进行比较,如果我理解正确的话,用LR测试检查哪一个是最重要的。
但如果我执行ANOVA(glm,glmm),我会得到Deviance Table的分析,而没有比较模型的输出。
如何获得我想要的输出,从而比较两种模型?
提前致谢, 柯恩
答案 0 :(得分:0)
在某个地方,您对使用anova()
给出了错误的印象。在MASS包中使用re
符合glmmPQL()
以下。 fe
适合使用glm()
来自基地:
> anova(re,fe)
#Error in anova.glmmPQL(re, fe) : 'anova' is not available for PQL fits
该消息似乎是创建anova.glmmPQL()
的唯一原因。
请参阅此主题以获得验证和模糊解释:
答案 1 :(得分:0)
anova
不能与glmmPQL
一起使用,您需要使用lme4软件包中的glmer
才能使用anova
。