Lorsque l’on calcule une régression logistique avec glm(), la fonction summary() ne fournit pas directement les odds ratio. Bien que ces derniers soient aisés à calculer (voir http://www.ats.ucla.edu/stat/r/dae/..., il peut être utile d’avoir une fonction les renvoyant directement.
Le constat est presque identique lorsque l’on réalise une régression multinomiale avec multinom() du package nnet. Mais cette fois-ci, les valeurs de p ne sont même pas accessibles via summary() (voir http://www.ats.ucla.edu/stat/r/dae/....
La fonction odds.ratio(), dont le code est visible à https://github.com/juba/questionr/blob/master/R/odds.ratio.r, permet de se faciliter la vie. On peut lui passer soit le résultat d’une régression logistique exécutée avec glm(family=binomial(logit)) soit le résultat d’une régression multinomiale réalisée avec multinom(). Elle affichera les odds ratio, leur intervalles de confiance ainsi que la valeur de p.
On peut préciser la confiance de l’intervalle avec le paramètre level ainsi que le nombre de décimale à afficher pour les odds ratio avec digits.