从 R 中的 glm 中提取系数答案

【问题标题】：extract coefficients from glm in R从 R 中的 glm 中提取系数
【发布时间】：2011-01-24 08:57:35
【问题描述】：

我进行了逻辑回归，结果如下：

ssi.logit.single.age["coefficients"]
# $coefficients
#  (Intercept)          age 
# -3.425062382  0.009916508

我需要为age取系数，目前我使用以下代码：

ssi.logit.single.age["coefficients"][[1]][2]

有效，但我不喜欢这里的神秘代码，我可以使用系数的名称（即(Intercept)或age）

【问题讨论】：

标签： list r

【解决方案1】：

有一个名为coef的提取函数可以从模型中获取系数：

coef(ssi.logit.single.age)["age"]

【讨论】：

那么如何提取支架错误？我需要计算 95% CI
查看 coef(summary(ssi.logit.single.age)) 以获取包含标准错误的表格。您可以使用 confint.default(ssi.logit.single.age) 获得基于 +/- 1.96 SE 的置信区间。或者更好的是，MASS::confint(ssi.logit.single.age) 会给你配置置信区间，更准确（虽然有点慢）。

【解决方案2】：

我找到了，来自here

看一下summary()产生的数据结构

> names(summary(lm.D9))
  [1] "call"          "terms"         "residuals"     "coefficients"
  [5] "aliased"       "sigma"         "df"            "r.squared"
  [9] "adj.r.squared" "fstatistic"    "cov.unscaled"

现在看一下摘要中系数的数据结构：

> summary(lm.D9)$coefficients
             Estimate Std. Error   t value     Pr(>|t|)
(Intercept)    5.032  0.2202177 22.850117 9.547128e-15
groupTrt      -0.371  0.3114349 -1.191260 2.490232e-01

> class(summary(lm.D9)$coefficients)
[1] "matrix"

> summary(lm.D9)$coefficients[,3]
(Intercept)    groupTrt
   22.850117   -1.191260

【讨论】：

通过str()函数显示R对象的内部结构也很有用，例如：str(summary(lm.D9))在你的例子中
不要开始钻研这样的对象！使用提取器功能；在这种情况下coef()。例如，查看$residuals 组件。您可以通过直接访问$residuals 来获取模型残差，但这些是什么？它们当然不是您可能立即想到的模型残差。在模型上使用resid()，可以选择残差的类型。
?? > all.equal( residuals(lm.D9) , lm.D9$residuals) [1] TRUE（glm 残差是“工作”残差。）
@DWin @lokheart 实际上有一个 GLM，这就是我所提到的，尽管答案是指链接到的 R-Help 帖子中的一个 LM。那里（他的 GLM）mod$residuals 与 resid(mod) 不同，前者不太可能是人们经常需要的东西。 R-help 帖子中的建议和@lokheart 的回答中重复的建议没有帮助，也不推荐。
好的。同意，也不准确。辛普森的建议特别重要，因为系数列表的第三列是 p 值而不是任何系数。残差（）的 glm 默认值不是“工作”，这是存储在模型拟合中的。