如何从 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感兴趣。

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

How to extract fixed effects part of summary from lme4?

r

lme4