Distribusi konvolusi variabel normal kuadrat dan chi-kuadrat?

14

masalah berikut muncul baru-baru ini saat menganalisis data. Jika variabel acak X mengikuti distribusi normal dan Y mengikuti distribusi χn2 (dengan n dof), bagaimana Z=X2+Y2 didistribusikan? Hingga sekarang saya datang dengan pdf dari : ψ 2 n ( x )Y2

ψn2(x)=F(x)x=(0xtn/21et/22n/2Γ(n/2)dt)x=12n/2Γ(n/2)(x)n/21ex/2(x)x=12n/21Γ(n/2)xn/41ex/2

serta beberapa penyederhanaan untuk integral konvolusi ( memiliki pdf dengan m dof):χ 2 mX2χm2

Kmn(t):=(χm2ψn2)(t)=0tχm2(x)ψn2(tx)dx=(2(n+m)2+1Γ(m2)Γ(n2))10t(tx)n41xm21exp((tx+x)/2)dx

Apakah seseorang melihat cara yang baik untuk menghitung integral ini untuk t nyata atau apakah itu harus dihitung secara numerik? Atau apakah saya melewatkan solusi yang lebih sederhana?

Leo Szilard
sumber
1
Jika tidak dikuadratkan, saya akan memiliki beberapa saran khusus. Saya tidak berpikir yang satu ini akan traktable (juga tidak perlu mencerahkan bahkan jika itu terbukti traktat). Saya akan tergoda untuk melihat pendekatan komputasi, seperti konvolusi numerik atau simulasi, tergantung pada apa yang ingin Anda lakukan dengan hasilnya. Y
Glen_b -Reinstate Monica
Sangat tidak mungkin menurut saya bahwa integral dapat dilakukan.
Dave31415
1
@ Dave31415 Untuk dan m genap, integral dapat secara eksplisit dihitung untuk nilai integral positif dari n dan m . Ini akan sama dengan kombinasi linear dari fungsi eksponensial dan kesalahan dengan koefisien polinomial dalam nmnm . Evaluasi dapat dilakukan melalui substitusix=t-u2. Misalnya, dengann=2,m=4kita memperoleh 1tx=tu2n=2,m=4. 14e18(2t+1)2(et2(2π(4t+3)(erfi(2t122)+erfi(122))+4e18)4et2+18(2t+1))
whuber
Bagus. Untuk angka ganjil, Anda mungkin dapat memperkirakannya dengan rata-rata hasil untuk membatasi angka genap? Atau mungkin tidak.
Dave31415
1
Terima kasih atas balasan Anda! Untuk beberapa kasus genap, saya mendapatkan hasil yang serupa yang melibatkan fungsi Dawson, tetapi sepertinya saya harus melakukan lebih banyak pekerjaan untuk solusi umum ...
Leo Szilard

Jawaban:

2

Dalam kasus ini membantu, variabel adalah variabel acak gamma umum (lihat misalnya, Stacy 1962). Pertanyaan Anda menanyakan distribusi jumlah variabel acak chi-kuadrat dan variabel acak gamma umum. Sepengetahuan saya, kerapatan variabel yang dihasilkan tidak memiliki ekspresi bentuk tertutup. Oleh karena itu, konvolusi yang Anda peroleh adalah integral dengan solusi tanpa bentuk tertutup. Saya pikir Anda akan terjebak dengan solusi numerik untuk yang satu ini.Y2


Stacy, EW (1962). Generalisasi Distribusi Gamma. Annals of Statistics Matematika 33 (3) , hlm. 1187-1192.

Pasang kembali Monica
sumber
0

Ini hanya petunjuk. Pearson tipe III dapat dikuadratkan. Terkadang konvolusi dapat ditemukan dengan menggabungkan sesuatu dengan dirinya sendiri. Saya berhasil melakukan ini untuk menggabungkan ND dan GD , yang mana saya menggabungkan Pearson III dengan dirinya sendiri. Bagaimana ini bekerja dengan ND 2 dan Chi-Squared, saya tidak yakin. Tapi, Anda meminta petunjuk, dan ini adalah petunjuk umum. Saya harap itu cukup untuk memulai.2

Carl
sumber
2
Bisakah Anda menjelaskan bagaimana ini menjawab pertanyaan? Tampaknya tidak berhubungan langsung.
Whuber
Konvolusi Pearson tipe III dengan dirinya sendiri dapat dilakukan. Untuk beberapa alasan, menggabungkan satu hal dengan dirinya sendiri lebih mudah untuk diselesaikan daripada menggabungkan satu hal dengan yang lain. Sebagai contoh, saya memecahkan konvolusi Pearson tipe III dan mendapatkan konvolusi ND dengan GD, masalah terkait.
Carl
Tampaknya tidak membantu, akan segera dihapus.
Carl