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