Fonctions de liaison : probit comparé au logit

Modèles linéaires généralisés (GLM) en R

Richard Erickson

Instructor

Pourquoi des fonctions de liaison ?

  • Comprendre et simuler les GLM
  • Probit vs logit comme exemple
Modèles linéaires généralisés (GLM) en R

Pourquoi le probit ?

  • Illustrer la fonction de liaison
  • Utilisée dans certains domaines (p. ex., toxicologie)
  • Préférée par certaines personnes
Modèles linéaires généralisés (GLM) en R

Qu'est-ce qu'un probit ?

  • « Pro »bability « it » (probit)
  • Toxicologie par Chester Bliss en 1934
  • Plus simple à calculer que le logit
  • Modèle connu comme analyse probit, régression probit ou modèle probit
Modèles linéaires généralisés (GLM) en R

Équation du probit

Modèle pour données binaires

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

Lié à une équation linéaire

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

Modèles linéaires généralisés (GLM) en R

Fonction probit

Basée sur la loi normale cumulée

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

Modèles linéaires généralisés (GLM) en R

Probabilité cumulée : probit et logit

Modèles linéaires généralisés (GLM) en R

Ajuster un probit en R

  • Option family pour glm()
    • Caractère : glm(..., family = "binomial")
    • Fonction : glm(..., family = binomial())
  • Par défaut : binomial(link = "logit")
  • Probit : binomial(link = "probit")
  • Suivre les consignes pour DataCamp
Modèles linéaires généralisés (GLM) en R

Simuler avec probit

Convertir de l'échelle probit à l'échelle de probabilité :

p = pnorm(-0.2)

Utiliser la probabilité avec une loi binomiale

rbinom(n = 10, size = 1, prob = p)
Modèles linéaires généralisés (GLM) en R

Simuler avec logit

Convertir de l'échelle logit à l'échelle de probabilité :

p = plogis(-.2)

Utiliser la probabilité avec une loi binomiale

rbinom(n = 10, size = 1, prob = p)
Modèles linéaires généralisés (GLM) en R

Quand utiliser probit vs logit ?

  • Largement propre au domaine
  • Queues plus épaisses pour le logit
  • Les deux sont défendables
Modèles linéaires généralisés (GLM) en R

Passons à la pratique !

Modèles linéaires généralisés (GLM) en R

Preparing Video For Download...