Извлекать коэффициенты из glm в R

Я выполнил логистическую регрессию со следующим результатом:

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

Мне нужно подобрать коэффициент для age, и в настоящее время я использую следующий код:

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

Это работает, но здесь мне не нужен критический код, могу ли я использовать имя коэффициента (т.е. (Intercept) или age)

Ответы

Ответ 1

Существует функция выделения, называемая coef, чтобы получить коэффициенты от моделей:

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

Ответ 2

Я нашел его, здесь

Посмотрите на структуру данных, созданную 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