Студент т как смесь гауссов

23

Используя t-распределение Стьюдента с k>0 степенями свободы, параметр местоположения и параметр шкалы имеют плотностьсls

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

как показать, что распределение Стьюдента может быть записано как смесь гауссовских распределений, если , , и интегрирование плотности соединения чтобы получить предельную плотность ? Каковы параметры результирующего распределения как функции ?tXN(μ,σ2)τ=1/σ2Γ(α,β)f(x,τ|μ)f(x|μ)tμ,α,β

Я затерялся в исчислении, интегрируя условную плотность сустава с гамма-распределением.

Салих Укан
источник

Ответы:

31

PDF нормального распределения

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

но с точки зрения этоτ=1/σ2

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

PDF-файл гамма-распределения

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

Их произведение, слегка упрощенное простой алгеброй, поэтому

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

Его внутренняя часть, очевидно, имеет вид , что делает ее кратной гамма-функции при интегрировании во всем диапазоне от τ = 0 до τ = . Следовательно, этот интеграл является немедленным (полученным, если знать, что интеграл гамма-распределения равен единице), давая предельное распределениеexp(constant1×τ)×τconstant2dττ=0τ=

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

Пытаясь соответствовать шаблону , отведенный для распределения показывает , есть ошибка в вопросе: PDF для распределения Стьюдента Студенческого фактически пропорционаленt

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

(the power of (xl)/s is 2, not 1). Matching the terms indicates k=2α, l=μ, and s=1/αβ.


Notice that no Calculus was needed for this derivation: everything was a matter of looking up the formulas of the Normal and Gamma PDFs, carrying out some trivial algebraic manipulations involving products and powers, and matching patterns in algebraic expressions (in that order).

whuber
источник
10
Inspired by this answer I made an animation of the t distribution as a mixture of normal distributions. It is available here: sumsar.net/blog/2013/12/t-as-a-mixture-of-normals
Rasmus Bååth
1
@whuber: Technically, for that kind of matching there is always an implicit use of calculus in your recognition that you can integrate out the gamma density using its known integral form. (This is the statistician's equivalent of hiding the broccoli by mixing it in with the meat and potatoes.) A clever way of hiding the calculus!
Reinstate Monica
1

I don't know the steps of the calculation, but I do know the results from some book (cannot remember which one...). I usually keep it in mind directly... :-) The Student t distribution with k degree freedom can be regarded as a Normal distribution with variance mixture Y, where Y follows inverse gamma distribution. More precisely, X~t(k),X=Y*Φ,where Y~IG(k/2,k/2),Φ is standard normal rv. I hope this could help you in some sense.

Jingjings
источник
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
источник