Comprendre la sortie d'une régression logistique

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

Richard Erickson

Instructor

Communiquer les résultats d'une régression logistique ?

Régression linéaire : simple
- Ajouter l'ordonnée à l'origine
- Multiplier par les pentes
Régression de Poisson :
- Nécessite une transformation exponentielle
- Semblable à la régression linéaire après transformation
Régression logistique ???

Pièce biaisée :
- Comparer pile à face
- Pile 3 fois pour 1 face
- Cotes de 3 contre 1
- Odds ratio : 3,0
Souvent utilisé en sports/jeux d'argent
Études médicales

Log-cotes (« logit ») :

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

Cotes, prendre l'exponentielle ($e^{x}$) :

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

Odds ratio (OR) pour une variable continue :

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

Valeurs d'OR :

Hommes non-fumeurs vs fumeurs (Pesch et al., 2012)

OR : 103,5 (IC 95 % : 74,8–143,2)
100 contre 1 de probabilité d'avoir un cancer chez les fumeurs !
En médecine, on présente souvent des intervalles de confiance à 95 % plutôt que des valeurs p
Tendance générale à s'éloigner des valeurs p

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

coef(glm_out)

exp(coef(glm_out))

confint(glm_out)

exp(confint(glm_out))

library(broom)

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

tidy(glm_out, exponentiate = TRUE, conf.int= TRUE)

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