Розуміння виводу з логістичної регресії

Узагальнені лінійні моделі в R

Richard Erickson

Instructor

Як пояснювати результати логістичної регресії?

  • Лінійна регресія проста:

    • Додає перехоплення
    • Множить нахили
  • Пуасонівська регресія:

    • Потрібна експоненціальна трансформація
    • Після трансформації схожа на лінійну
  • Логістична регресія???
Узагальнені лінійні моделі в R

Відношення шансів (odds ratio)

  • Не така пряма, як пуасонівська з експонентою
  • Використовують для порівняння відносних шансів двох подій
Узагальнені лінійні моделі в R

Приклад відношення шансів

  • Нечесна монета:
    • Порівняйте орла з решкою
    • Орел випадає 3 рази на 1 решку
    • Шанси 3 до 1
    • Відношення шансів 3,0
  • Часто в спорті/азартних іграх
  • Медичні дослідження
Узагальнені лінійні моделі в R

Логістичний вивід відношення шансів

Лог-шанси («logit»):

$\phi(x) = \text{ln}(\frac{p(x)}{1 - p(x)}) = \beta_0 + \beta_1 x$

Шанси, беремо експоненту ($e^{x}$):

$\frac{p(x)}{1 - p(x)} = e^{\beta_0 + \beta_1x}$

Узагальнені лінійні моделі в R

Відношення шансів для неперервної змінної

Відношення шансів (OR) для неперервної змінної:

$\text{OR} = \frac{e^{\beta_0 + \beta_1(x +1)}}{e^{\beta_0 + \beta_1x}} = e^{\beta_1}$

Узагальнені лінійні моделі в R

Інтерпретація

Значення OR:

  • OR = 1: коефіцієнт не впливає
  • OR < 1: коефіцієнт зменшує шанси
  • OR > 1: коефіцієнт збільшує шанси
Узагальнені лінійні моделі в R

Приклад про рак

Чоловіки-некурці проти курців (Pesch та ін., 2012)

  • OR: 103,5 (95% ДІ 74,8–143,2)
  • $>$ 100 до 1 шанси захворіти на рак у курців!
  • У медицині часто наводять 95% довірчі інтервали, а не p-значення
  • Загальна тенденція відходити від p-значень
Узагальнені лінійні моделі в R

Отримання з GLM

glm_out <- glm(y ~ x, family = 'binomial')

coef(glm_out)

exp(coef(glm_out))

confint(glm_out)

exp(confint(glm_out))
Узагальнені лінійні моделі в R

Tidyverse

library(broom)

glm_out <- glm(y ~ x, family = 'binomial')

tidy(glm_out, exponentiate = TRUE, conf.int= TRUE)
Узагальнені лінійні моделі в R

Давайте потренуємось!

Узагальнені лінійні моделі в R

Preparing Video For Download...