我希望将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感兴趣。