首页 > 解决方案 > 什么是简单混合模型代码 (R) 中的“+0”?

问题描述

我正在尝试比较使用包 lme4 在混合模型分析中遇到的差异。也许我的统计背景不够清晰,但我无法弄清楚代码中的“+0”是什么以及由此产生的差异(对于没有 +0 的模型)意味着什么。

这是我的 +0 示例:

lmer(Yield ~ Treatment + 0 + (1|Batch) + (1|Irrigation), data = D)

对比:

lmer(Yield ~ Treatment  + (1|Batch) + (1|Irrigation), data = D)

有人对 +0 是什么以及它对结果的作用有一个聪明的解释吗?

标签: rstatisticsmixed-models

解决方案


具有+ 0通常意味着“没有整体截距”的模型(在固定效应中)。默认情况下,模型包含截距,您也可以使用+ 1.

大多数关于回归建模的讨论都建议包含一个截距,除非有充分的理由相信当预测变量都为零时结果将为 0(某些物理过程可能是真的?)。

相比:

fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
fm2 <- lmer(Reaction ~ Days + 0 + (Days | Subject), sleepstudy)
summary(fm1)
summary(fm2)

注意固定效应


推荐阅读