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

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

我知道我可以使用以下方法提取summary的随机效果部分:
SleepStudy_RE <- as.data.frame(VarCorr(SleepStudy))

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

谢谢你。

最佳答案

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

09-06 22:22