Siswa t sebagai campuran dari gaussian

23

Menggunakan distribusi-t siswa dengan k>0 derajat kebebasan, parameter lokasi l dan parameter skala s memiliki kepadatan

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

bagaimana menunjukkan bahwa distribusi- Student tdapat ditulis sebagai campuran dari distribusi Gaussian dengan membiarkan XN(μ,σ2) , τ=1/σ2Γ(α,β) , dan mengintegrasikan kepadatan bersama f(x,τ|μ) untuk mendapatkan kepadatan marginal f(x|μ) ? Apa parameter dari t yang dihasilkant-distribusi, sebagai fungsi ?μ,α,β

Saya tersesat dalam kalkulus dengan mengintegrasikan kepadatan bersyarat bersama dengan distribusi Gamma.

Salih Ucan
sumber

Jawaban:

31

PDF dari distribusi Normal adalah

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

tetapi dalam hal ituτ=1/σ2

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

PDF distribusi Gamma adalah

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

Produk mereka, sedikit disederhanakan dengan aljabar mudah, karenanya

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

Bagian dalamnya ternyata memiliki bentuk , membuatnya menjadi beberapa dari fungsi Gamma ketika terintegrasi atas lengkap τ = 0 untuk τ = . Karena itu integral itu langsung (diperoleh dengan mengetahui integral dari distribusi Gamma adalah kesatuan), memberikan distribusi marjinalexp(constant1×τ)×τconstant2dττ=0τ=

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

Mencoba untuk mencocokkan pola disediakan untuk menunjukkan distribusi ada kesalahan dalam pertanyaan: yang PDF untuk distribusi t Student sebenarnya sebanding dengant

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

(kekuatan adalah 2 , bukan 1 ). Mencocokkan istilah menunjukkan k = 2 α , l = μ , dan s = 1 / (xl)/s21k=2αl=μ .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
sumber
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: Secara teknis, untuk pencocokan semacam itu selalu ada penggunaan kalkulus secara implisit dalam pengakuan Anda bahwa Anda dapat mengintegrasikan kepadatan gamma menggunakan bentuk integral yang dikenal. (Ini setara dengan ahli statistik dalam menyembunyikan brokoli dengan mencampurkannya dengan daging dan kentang.) Cara pintar menyembunyikan kalkulus!
Pasang kembali Monica
1

Saya tidak tahu langkah-langkah perhitungannya, tetapi saya tahu hasil dari beberapa buku (tidak ingat yang mana ...). Saya biasanya mengingatnya secara langsung ... :-) Mahasiswat distribusi dengan k derajat kebebasan dapat dianggap sebagai distribusi normal dengan campuran varians Ydimana Ymengikuti distribusi gamma terbalik. Lebih tepatnya,X~t(k),X=Y*Φ,dimana Y~sayaG(k/2,k/2),Φadalah standar normal rv. Saya harap ini bisa membantu Anda dalam beberapa hal.

Jingjings
sumber
0

Untuk menyederhanakan kita anggap berarti 0. Dengan menggunakan representasi, kami menunjukkan hasilnya untuk derajat kebebasan integer.

1/τX=Y
setara dengan campuran Gaussian dengan yang sebelumnya: dikondisikan pada τ, Y adalah Gaussian dengan presisi τ, dan sebelumnya τseperti yang diinginkan. Maka tetap menunjukkan itu1/τXadalah distribusi-t. Kita bisa menulis
τΓ(α,β)β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
sumber