lme4包中提取固定效应标准误的R代码

Question

我希望能够提取使用lmer()进行回归分析的标准误差。例如，

library(lme4)
data(iris)
head(iris)
model <- lmer(Sepal.Length ~ Sepal.Width + (1 | Species), data=iris);       
summary(model)

查看汇总的固定效应，您可以看到标准误差。

Fixed effects:
            Estimate Std. Error t value
(Intercept)   3.4062     0.6683   5.097
Sepal.Width   0.7972     0.1062   7.506

我希望能够自动提取这部分内容，而不是手动输入。

- Blake Shurtz

---

为什么这个问题被认为是离题的？它是关于统计编程的一个问题。 - Blake Shurtz

Answer 1

summary(model)$coef[, 2, drop = FALSE]

或者

sqrt(diag(vcov(model)))

- Mark White

---

@amoeba 我想我会回复并将其迁移到SO。 - Mark White