如何从lme4中提取固定效果的一部分摘要？

Question

我希望将summary()的固定效果部分提取为data.frame。我正在使用lme4来运行以下模型：

SleepStudy <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)
summary(SleepStudy)

我知道我可以使用以下内容提取summary的随机效果部分：

SleepStudy_RE <- as.data.frame(VarCorr(SleepStudy))

是否有类似的固定效果代码，包括估算，标准误差，自由度和精确的p值？

谢谢。

Answer 1

coef(summary(fitted_model))应该这样做。

library(lme4)
SleepStudy <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)
coef(summary(SleepStudy))
##              Estimate Std. Error  t value
## (Intercept) 251.40510  9.7467163 25.79383
## Days         10.46729  0.8042214 13.01543

如果你想要p值，你需要lmerTest（你需要重新拟合模型）：

library(lmerTest)
SleepStudy <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)
coef(summary(SleepStudy))
##              Estimate Std. Error       df  t value Pr(>|t|)
## (Intercept) 251.40510  9.7467163  22.8102 25.79383        0
## Days         10.46729  0.8042214 161.0036 13.01543        0

在这种情况下，我不知道为什么p值正好为零;也许需要了解lmerTest维护者。

您可能也对broom package感兴趣。