Quelles sont les définitions des prieurs conjugués semi-conjugués et conditionnels?


12

Quelles sont les définitions des prieurs semi-conjugués et des prieurs conjugués conditionnels ? Je les ai trouvés dans l'analyse des données bayésiennes de Gelman , mais je n'ai pas pu trouver leurs définitions.

Réponses:


16

En utilisant la définition de Bayesian Data Analysis (3e éd.) , Si est une classe de distributions d'échantillonnage et est une classe de distributions antérieures pour , alors le la classe est conjuguée à si p ( y | θ ) P θ P FFp(y|θ)PθPF

p(θ|y)P pour tous p(|θ)F et p()P.

Si est une classe de distributions d'échantillonnage , et est une classe de distributions antérieures pour conditionnelle à , alors la classe est un conjugué conditionnel pour si p ( y | θ , ϕ ) P θ ϕ P FFp(y|θ,ϕ)PθϕPF

p(θ|y,ϕ)P pour tous p(|θ,ϕ)F et p(|ϕ)P.

Les antérieurs conjugués conditionnellement sont pratiques pour construire un échantillonneur de Gibbs car le conditionnel complet sera une famille connue.

J'ai recherché une version électronique de Bayesian Data Analysis (3e éd.) Et je n'ai pas pu trouver une référence à semi-conjugué auparavant. Je suppose qu'il est synonyme de conjugué conditionnellement, mais si vous fournissez une référence à son utilisation dans le livre, je devrais être en mesure de fournir une définition.


+1. Quelle est l'URL de la 3e édition de Bayesian Data Analysis?
Patrick Coulombe

1
Merci! Le semi-conjugué apparaît ici (2e éd.) Books.google.com/… . Au fait, comment avez-vous obtenu l'ebook pour la 3e édition?
Tim

1
Je ne sais pas pourquoi il est écrit semi-conjugué avant, car le prieur est entièrement conjugué. Cette déclaration est supprimée dans la 3e édition. L'ebook peut être acheté ici: crcpress.com/product/isbn/9781439840955 .
jaradniemi

@jaradniemi: Dans le lien que j'ai donné, en plus de p84, il est souligné que le prémium conjugué n'est pas un prieur conjugué.
Tim

1
Dans à quoi chacun des ⋅ fait-il référence et fait-il référence à la même chose?
p(θ|y,ϕ)P pour tous p(|θ,ϕ)F et p(|ϕ)P.
Muno

6

Je voudrais utiliser la normale multivariée comme exemple.

Rappelons que la probabilité est donnée par

P(y1,y2,...,yn|μ,Σ)=(2π)-N2det(Σ)-N2exp(12je=1N(Xje-μ)TΣ-1(Xje-μ))

Afin de trouver un avant cette probabilité, nous pouvons choisir

P(μ,Σ)=Ordinaire(μ;μ0,Λ0)InverseWishart(Σ;ν0,S0)

Je vous assure de ne PAS vous soucier de pour l'instant; ce sont simplement des paramètres de la distribution antérieure.μ0,Λ0,ν0,S0

Ce qui est important, cependant, c'est que ce n'est pas conjugué à la probabilité. Pour voir pourquoi, je voudrais citer une référence que j'ai trouvée en ligne.

notons que et Σ apparaissent ensemble de manière non factorisée dans la vraisemblance; par conséquent, ils seront également couplés ensemble dans la partie postérieureμΣ

La référence est "Machine Learning: A Probabilistic Perspective" par Kevin P. Murphy. Voici le lien . Vous pouvez trouver la citation dans la section 4.6 (déduction des paramètres d'un MVN) en haut de la page 135.

Pour continuer la citation,

Le précédent ci-dessus est parfois appelé semi-conjugué ou conjugué conditionnellement , car les deux conditions, et p ( Σ | μ ) , sont individuellement conjuguées. Pour créer un a priori conjugué complet , nous devons utiliser un a priori où μ et Σ dépendent l'un de l'autre. Nous utiliserons une distribution conjointe du formulairep(μ|Σ)p(Σ|μ)μΣ

p(μ,Σ)=p(Σ)p(μ|Σ)

L'idée ici est que la première distribution antérieure

P(μ,Σ)=Ordinaire(μ;μ0,Λ0)InverseWishart(Σ;ν0,S0)

suppose que et Σ sont séparables (ou indépendants dans un sens). Néanmoins, nous observons que dans la fonction de vraisemblance, μ et Σ ne peuvent pas être factorisés séparément, ce qui implique qu'ils ne seront pas séparables dans le postérieur (Rappel, ( Postérieur ) ( Prior ) ( Probabilité ) ). Cela montre que les antérieurs "indissociables" et les antérieurs "séparables" ne sont pas conjugués. D'autre part, en réécrivantμΣμΣ(Postérieur)(Avant)(Probabilité)

p(μ,Σ)=p(Σ)p(μ|Σ)

de sorte que et Σ dépendent l'un de l'autre (à travers p ( μ | Σ ) ), vous obtiendrez un a priori conjugué, qui est nommé a priori semi-conjugué . J'espère que cela répond à votre question.μΣp(μ|Σ)

ps : Une autre référence vraiment utile que j'ai utilisée est "Un premier cours sur les méthodes statistiques bayésiennes" de Peter D. Hoff. Voici un lien vers le livre. Vous pouvez trouver du contenu pertinent dans la section 7 à partir de la page 105, et il a une très bonne explication (et intuition) sur la distribution normale à une seule variable dans la section 5 à partir de la page 67, qui sera encore renforcée dans la section 7 lorsqu'il traitera MVN.


1

Fp(y|θ,ϕ)PθPFp(θ|y,ϕ)Pp(|θ,ϕ)Fp(θ,ϕ)=p(θ)×p(ϕ)p(θ)Pp(ϕ)P

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.