Combien de distributions sont dans le GLM?


11

J'ai identifié plusieurs endroits dans les manuels où le GLM est décrit avec 5 distributions (à savoir, Gamma, gaussienne, binomiale, gaussienne inverse et Poisson). Ceci est également illustré dans la fonction familiale dans R.

Parfois, je rencontre des références au GLM où des distributions supplémentaires sont incluses ( exemple ). Quelqu'un peut-il expliquer pourquoi ces 5 sont spéciaux ou sont toujours dans le GLM, mais parfois d'autres le sont?

D'après ce que j'ai appris jusqu'à présent, les distributions GLM dans la famille exponentielle tiennent toutes sous la forme: où est le paramètre de dispersion et est le paramètre canonique.ϕθ

f(y;θ,ϕ)=exp{yθb(θ)ϕ+c(y,ϕ)}
ϕθ

Aucune distribution ne peut-elle être transformée pour s'intégrer dans le GLM?


5
De toute évidence, la distribution uniforme n'appartient pas à la famille exponentielle.
Zhanxiong

Bonne question. Par exemple, qu'en est-il du lognormal?
Michael M

1
@Zhanxiong, l'uniforme n'est-il pas un cas particulier de distribution bêta, et la distribution bêta est dans la famille exponentielle?
shf8888

@ shf8888 AFAIK c'est seulement une distribution de famille exponentielle dans la limite, quand elle converge vers la distribution gamma.
shadowtalker

@Zhanxiong, merci d'avoir clarifié! Toutes mes excuses, vous avez raison, avec des limites inconnues n'est pas une distribution familiale exponentielle.
shf8888

Réponses:


4

Comme vous l'indiquez, la qualification pour utiliser une distribution dans un GLM est qu'elle soit de la famille exponentielle (note: ce n'est pas la même chose que la distribution exponentielle! Bien que la distribution exponentielle, en tant que distribution gamma, fasse elle-même partie de la famille exponentielle). Les cinq distributions que vous listez sont toutes de cette famille, et plus important encore, sont des distributions TRÈS courantes, elles sont donc utilisées comme exemples et explications.

Comme le note Zhanxiong, la distribution uniforme (avec des bornes inconnues) est un exemple classique de distribution familiale non exponentielle. shf8888 confond la distribution uniforme générale, à n'importe quel intervalle, avec un uniforme (0, 1). La distribution uniforme (0,1) est un cas particulier de la distribution bêta, qui est une famille exponentielle. Les autres distributions non exponentielles des familles sont les modèles de mélange et la distribution t.

Vous avez la définition de la famille exponentielle correcte, et le paramètre canonique est très important pour utiliser GLM. Pourtant, j'ai toujours trouvé un peu plus facile de comprendre la famille exponentielle en l'écrivant comme:

f(x;θ)=a(θ)g(x)exp[b(θ)R(x)]

Il existe une manière plus générale d'écrire ceci, avec un vecteur au lieu d'un scalaire ; mais le cas unidimensionnel explique beaucoup de choses. Plus précisément, vous devez être capable de factoriser la partie non exponentiée de votre densité en deux fonctions, une de paramètre inconnu mais pas de données observées et une de et non ; et de même pour la partie exponentiée. Il peut être difficile de voir comment, par exemple, la distribution binomiale peut être écrite de cette façon; mais avec quelques jonglages algébriques, cela devient clair finalement.θθθxxθ

Nous utilisons la famille exponentielle car elle facilite beaucoup de choses: par exemple, trouver des statistiques suffisantes et tester des hypothèses. Dans GLM, le paramètre canonique est souvent utilisé pour trouver une fonction de lien. Enfin, une illustration connexe de la raison pour laquelle les statisticiens préfèrent utiliser la famille exponentielle dans presque tous les cas tente de faire une inférence statistique classique sur, disons, une distribution uniforme ( , ) où les deux et sont inconnus . Ce n'est pas impossible, mais c'est beaucoup plus compliqué et impliqué que de faire de même pour les distributions familiales exponentielles.θ1θ2θ1θ2


La distribution bêta avec les deux paramètres inconnus est toujours une famille exponentielle (mais une famille exponentielle à 2 paramètres). Qu'est-ce qui vous fait penser que ce n'est pas le cas? www2.stat.duke.edu/courses/Spring11/sta114/lec/… ou wikipedia
DavidR

Merci d'avoir signalé cela, j'ai changé mon commentaire ... vous avez raison! Je ne sais vraiment pas ce que je voulais dire
Henry
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.