J'ai trouvé une formule pour le pseudo dans le livre Extending the Linear Model with R, Julian J. Faraway (p. 59).
Est-ce une formule courante pour le pseudo pour les GLM?
la source
J'ai trouvé une formule pour le pseudo dans le livre Extending the Linear Model with R, Julian J. Faraway (p. 59).
Est-ce une formule courante pour le pseudo pour les GLM?
Il existe un grand nombre de pseudo- pour les GLiM. L'excellent site d'aide des statistiques UCLA en a un aperçu complet ici . Celui que vous listez s'appelle le pseudo- McFadden . Par rapport à la typologie UCLA, elle est comme dans le sens où elle indexe l'amélioration du modèle ajusté sur le modèle nul. Certains logiciels statistiques, notamment SPSS, si je me souviens bien, impriment par défaut le pseudo- de McFadden avec les résultats de certaines analyses comme la régression logistique, donc je soupçonne que c'est assez courant, bien que les pseudo- Cox & Snell et Nagelkerke pseudo- s peut être encore plus. Cependant, le pseudo- McFadden n'a pas toutes les propriétés deR 2 R 2 R 2 R 2 R 2 R 2 R 2 R 2 R 2 R 2 R 2 (pas de pseudo- ). Si quelqu'un souhaite utiliser un pseudo- pour comprendre un modèle, je recommande fortement de lire cet excellent fil de discussion: Quelle pseudo- mesure est celle à signaler pour la régression logistique (Cox & Snell ou Nagelkerke)? (Pour ce que ça vaut, est lui - même que les gens se rendent compte glissante, une grande démonstration qui peut être vu dans la réponse de @ whuber ici: Is ? Utile ou dangereux )
Colin Cameron, A., & Windmeijer, F. A. (1997). An R-squared measure of goodness of fit for some common nonlinear regression models. Journal of Econometrics, 77(2), 329-342.
R donne une déviance nulle et résiduelle dans la sortie pour
glm
que vous puissiez faire exactement ce genre de comparaison (voir les deux dernières lignes ci-dessous).Vous pouvez également extraire ces valeurs de l'objet avec
model$null.deviance
etmodel$deviance
la source
summary.glm
La formule que vous avez proposée a été proposée par Maddala (1983) et Magee (1990) pour estimer R au carré sur le modèle logistique. Par conséquent, je ne pense pas que cela s'applique à tous les modèles glm (voir le livre Modern Regression Methods de Thomas P. Ryan à la page 266).
Si vous faites un faux ensemble de données, vous verrez qu'il sous-estime le R au carré ... pour la glm gaussienne par exemple.
Je pense que pour un gla gaussien, vous pouvez utiliser la formule de base (lm) R au carré ...
Et pour la logistique (ou la famille binomiale en r) j'utiliserais la formule que vous avez proposée ...
Jusqu'à présent, pour poisson glm, j'ai utilisé l'équation de ce post.
/programming/23067475/how-do-i-obtain-pseudo-r2-measures-in-stata-when-using-glm-regression
Il y a aussi un super article sur la pseudo R2 disponible sur les portails de recherche ... voici le lien:
https://www.researchgate.net/publication/222802021_Pseudo_R-squared_measures_for_Poisson_regression_models_with_over-_or_underdispersion
J'espère que cette aide.
la source
1-summary(GLM)$deviance/summary(GLM)$null.deviance
et vous verrez que le R2 correspond à la valeur R2 d'une régression OLS régulière, donc la réponse ci-dessus est correcte! Voir aussi mon article ici - stats.stackexchange.com/questions/412580/…Le paquet de R
modEvA
calcule D-Squared comme1 - (mod$deviance/mod$null.deviance)
mentionné par David J. HarrisLe D-Squared ou la déviance expliquée du modèle est présenté dans (Guisan & Zimmermann 2000) https://doi.org/10.1016/S0304-3800(00)00354-9
la source