Функції зв'язку — порівняння probit і logit

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

Richard Erickson

Instructor

Навіщо функції зв'язку?

  • Зрозуміти й змоделювати GLM
  • Probit vs logit як приклад
Узагальнені лінійні моделі в R

Чому probit?

  • Показати функцію зв'язку
  • Використовується в окремих галузях (напр., токсикологія)
  • Декому подобається більше
Узагальнені лінійні моделі в R

Що таке probit?

  • Від "pro"bability un"it"
  • Токсикологія: Честер Бліс, 1934
  • Обчислювально простіший за logit
  • Модель відома як probit-аналіз, probit-регресія або probit-модель
Узагальнені лінійні моделі в R

Рівняння probit

Модель бінарних даних

$Y = \text{Binomial}(p)$

Пов'язана з лінійним рівнянням

$\Phi^{-1} (p) = \beta_0 + \beta_1 x + \epsilon $

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

Функція probit

Ґрунтується на кумулятивному нормальному розподілі

$\Phi(z) = \frac{1}{\sqrt{2 \pi}} \int_{-\infty}^z e^{- \frac{1}{2} z^2} \, dz$

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

Кумулятивна ймовірність для probit і logit

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

Підігнати probit в R

  • Параметр family для glm()
    • Рядок: glm(..., family = "binomial")
    • Функція: glm(..., family = binomial())
  • Типово: binomial(link = "logit")
  • Probit: binomial(link = "probit")
  • Дотримуйтесь інструкцій DataCamp
Узагальнені лінійні моделі в R

Імітація з probit

Перетворіть зі шкали probit на шкалу ймовірностей:

p = pnorm(-0.2)

Використайте ймовірність з біноміальним розподілом

rbinom(n = 10, size = 1, prob = p)
Узагальнені лінійні моделі в R

Імітація з logit

Перетворіть зі шкали logit на шкалу ймовірностей:

p = plogis(-.2)

Використайте ймовірність з біноміальним розподілом

rbinom(n = 10, size = 1, prob = p)
Узагальнені лінійні моделі в R

Коли використовувати probit vs logit?

  • Здебільшого залежить від галузі
  • Товстіші хвости у logit
  • Обидва підходи прийнятні
Узагальнені лінійні моделі в R

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

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

Preparing Video For Download...