模型拟合:glm vs glmmPQL

时间:2013-04-02 09:43:40

标签: r glm anova

我正在拟合关于缺席 - 存在数据的模型,我想检查随机因素是否重要。 要做到这一点,我应该将glmm与glm进行比较,如果我理解正确的话,用LR测试检查哪一个是最重要的。

但如果我执行ANOVA(glm,glmm),我会得到Deviance Table的分析,而没有比较模型的输出。

如何获得我想要的输出,从而比较两种模型?

提前致谢, 柯恩

2 个答案:

答案 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()的唯一原因。

请参阅此主题以获得验证和模糊解释:

https://stat.ethz.ch/pipermail/r-help/2002-July/022987.html

答案 1 :(得分:0)

anova不能与glmmPQL一起使用,您需要使用lme4软件包中的glmer才能使用anova