用于在 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

我想把它拉出来而不是手动输入。

Answer 1

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

或者

sqrt(diag(vcov(model)))