【问题标题】:R code for pulling out fixed effect standard errors in lme4 package用于提取 lme4 包中的固定效果标准错误的 R 代码
【发布时间】:2018-08-22 18:50:48
【问题描述】:

我想提取使用 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

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

【问题讨论】:

  • 为什么这是题外话?这是一个关于统计编程的问题。

标签: r lme4


【解决方案1】:
summary(model)$coef[, 2, drop = FALSE]

或者

sqrt(diag(vcov(model)))

【讨论】:

  • @amoeba 我想我会回应,我们会将其迁移到 SO
猜你喜欢
  • 1970-01-01
  • 1970-01-01
  • 2016-01-14
  • 1970-01-01
  • 2020-08-11
  • 2018-07-31
  • 1970-01-01
  • 1970-01-01
  • 2021-10-16
相关资源
最近更新 更多