Une question sur les paramètres de la distribution gamma en économétrie bayésienne


8

L'article de Wikipédia sur la distribution Gamma répertorie deux méthodes de paramétrage différentes, l'une d'entre elles fréquemment utilisée en économétrie bayésienne avec et , est le paramètre de forme, est le paramètre de taux.α>0β>0αβ

XGamma(α,β).

Dans un manuel d'économétrie bayésienne écrit par Gary Koop, le paramètre de précision suit une distribution Gamma, qui est une distribution antérieure1σ2=h

hGamma(s_2,ν_),

où est moyen et est degrés de liberté selon son Annexe. De plus est une erreur standard avec définitions_2ν_s2

s2=(yiβ^xi)ν.

Ainsi pour moi, ces deux définitions de la distribution Gamma sont complètement différentes, puisque la moyenne et les variances seront différentes. Si nous suivons la définition de wikipedia, la moyenne sera , pas .α/βs_2

Je suis très confus ici, quelqu'un m'aiderait-il à éclaircir les pensées ici?


Je pense que vous faites une confusion: est l'écart-type estimé des données, pas l'écart-type de la distribution gamma. Et ce devrait être le postérieur, pas le prieur. s2
Stéphane Laurent

2
Le Gamma, malheureusement, n'a pas de paramétrage standard unique. Parfois, un Gamma (a, b) a un moyen , parfois et parfois avec le paramètre de forme . (Il ne s'agit pas d'une liste exhaustive.) Ils sont tous équivalents, par exemple, le dans le deuxième cas est égal à l'inverse du dans le premier cas. Vous devez donc porter une attention particulière à la façon dont la fonction de densité est écrite pour voir quel paramétrage est utilisé. aba/babbb
jbowman

Réponses:


3

Pour tous ceux qui luttent encore avec la notation terrible de Koops: Le problème est que Koop n'utilise ni l'échelle ni la paramétrisation du taux , mais plutôt une paramétrisation «moyenne, degrés de liberté» (voir Annexe, Déf. B. 22). La distribution de dans une paramétrisation correcte (forme, taux) est donc utilisant Notation Koops pour les paramètres.h

hGamma(shape=ν_/2,rate=νs_2/2)

2

Je pense que l'article Wikipedia fait référence à une forme spécifique de la distribution gamma connue sous le nom de . Le chi carré est et serait la constante par laquelle la variable aléatoire est multipliée par pour obtenir une variable aléatoire avec la distribution d'une estimation de variance. C'est-à-dire et . C'est s qui est l'erreur standard et non . Dans l'article que vous avez mentionné, le est répertorié dans les cas spéciaux (deuxième puce).χ2Gamma(ν,1/2)s2χ2α=νβ=1/2s2χ2


1

Il est d'usage d'imposer (comme auparavant) soit la distribution gamma à h=1σ2 ou la distribution gamma inverse à σ2. Ensuite, le posteior aura une belle apparence. Je pense que vous pouvez attribuer une distribution gamma àσ2, et toujours tous les calculs pour dériver le marginal en intégrant σ2 va passer.

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.