R代码，用于提取lme4包中的固定效果标准错误

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

我想将其提取出来，而不必手动输入。

Answer 1

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

或

sqrt(diag(vcov(model)))