Student t as mixture of gaussian

23

Utilisation de la distribution t de Student avec k>0 degrés de liberté, paramètre de localisation et paramètre d'échelle ayant une densitésls

Γ(k+12)Γ(k2kπs2){1+k1(xls)}(k+1)/2,

comment montrer que la distribution Student peut être écrite comme un mélange de distributions gaussiennes en laissant , , et intégrant la densité conjointe pour obtenir la densité marginale ? Quels sont les paramètres de la distribution résultante, en tant que fonctions de ?tXN(μ,σ2)τ=1/σ2Γ(α,β)f(x,τ|μ)f(x|μ)tμ,α,β

Je me suis perdu dans le calcul en intégrant la densité conditionnelle conjointe à la distribution gamma.

Salih Ucan
la source

Réponses:

31

Le PDF d'une distribution normale est

fμ,σ(x)=12πσe(xμ)22σ2dx

mais en termes de il estτ=1/σ2

gμ,τ(x)=τ2πeτ(xμ)22dx.

Le PDF d'une distribution Gamma est

hα,β(τ)=1Γ(α)eτβτ1+αβαdτ.

Leur produit, légèrement simplifié avec l'algèbre facile, est donc

fμ,α,β(x,τ)=1βαΓ(α)2πeτ((xμ)22+1β)τ1/2+αdτdx.

Sa partie intérieure a évidemment la forme , ce qui en fait un multiple d'une fonction Gamma lorsqu'elle est intégrée sur toute la plage τ = 0 à τ = . Cette intégrale est donc immédiate (obtenue en sachant que l'intégrale d'une distribution gamma est l'unité), donnant la distribution marginaleexp(constant1×τ)×τconstant2dττ=0τ=

fμ,α,β(x)=βΓ(α+12)2πΓ(α)1(β2(xμ)2+1)α+12.

Essayer de faire correspondre le modèle fourni pour la distribution montre qu'il y a une erreur dans la question: le PDF pour la distribution t de Student est en fait proportionnel àt

1ks(11+k1(xls)2)k+12

(la puissance de est 2 , pas 1 ). La correspondance des termes indique k = 2 α , l = μ et s = 1 / (xl)/s21k=2αl=μ .s=1/αβ


Notez qu'aucun calcul n'était nécessaire pour cette dérivation: tout consistait à rechercher les formules des PDF normaux et gamma, à effectuer des manipulations algébriques triviales impliquant des produits et des puissances, et à associer des modèles dans des expressions algébriques (dans cet ordre).

whuber
la source
10
Inspiré par cette réponse, j'ai fait une animation de la distribution t comme un mélange de distributions normales. Il est disponible ici: sumsar.net/blog/2013/12/t-as-a-mixture-of-normals
Rasmus Bååth
1
@whuber: Techniquement, pour ce type de correspondance, il y a toujours une utilisation implicite du calcul dans votre reconnaissance que vous pouvez intégrer la densité gamma en utilisant sa forme intégrale connue. (C'est l'équivalent pour le statisticien de cacher le brocoli en le mélangeant avec la viande et les pommes de terre.) Une façon intelligente de cacher le calcul!
Rétablir Monica le
1

tkYYXt(k)XYΦ,where Y~IG(k/2,k/2),Φ is standard normal rv. I hope this could help you in some sense.

Jingjings
la source
0

To simplify we assume mean 0. Using representation, we show the result for integer degrees of freedom.

1/τX=Y
is equivalent to a Gaussian mixture with that prior: conditioned on τ, Y is Gaussian with precision τ, and the prior τ is as desired. Then it remains to show that 1/τX is a t-distribution. We can write
τΓ(α,β)β2Γ(α,2)β2χ2(2α)
using a well-known result about gammas and Chi-squares (decompose a gamma as a sum of exponentials and combine the exponentials to normals to Chi squares) This in turn implies that
YX1(β/2)χ2(2α)
=Xαβχ2α2/(2α)
which is a scaled t with k=2α and s=1/αβ by variance of t. We can recenter our representation at μ and l would follow.

Romil Sirohi
la source