我正在拟合一个关于缺勤数据的模型,我想检查随机因素是否显着。为此,应该将 glmm 与 glm 进行比较,并使用 LR 测试检查哪个最重要,如果我理解正确的话。
但是,如果我执行 ANOVA(glm,glmm) ,我会得到偏差表的分析,并且没有比较模型的输出。
如何获得我想要的输出,从而比较两个模型?
在此先感谢,科恩
在某个地方,您对此使用有错误的印象anova()
。下面re
适合glmmPQL()
在 MASS 包中使用。fe
适合使用glm()
从基地:
> anova(re,fe)
#Error in anova.glmmPQL(re, fe) : 'anova' is not available for PQL fits
该消息似乎anova.glmmPQL()
是创建的唯一原因。
请参阅此线程以进行验证和模糊解释:
根本anova
不适合glmmPQL
你需要使用glmer
from lme4 package 才能使用anova
。