Distribusi jika Beta dan chi-kuadrat dengan derajat

8

Seandainya X memiliki distribusi beta, Beta(1,K1) dan Y mengikuti chi-squared dengan 2Kderajat. Selain itu, kami menganggap ituX dan Y independen.

Apa distribusi produk .Z=XY

Perbarui
Upaya saya:

fZ=y=y=+1|y|fY(y)fX(zy)dy=0+1B(1,K1)2KΓ(K)1yyK1ey/2(1z/y)K2dy=1B(1,K1)2KΓ(K)0+ey/2(yz)K2dy=1B(1,K1)2KΓ(K)[2K1ez/2Γ(K1,yz2)]0=2K1B(1,K1)2KΓ(K)ez/2Γ(K1,z/2)

Apakah itu benar? jika ya, bagaimana kita menyebut distribusi ini?

tam
sumber
2
Jika ini adalah pekerjaan rumah atau belajar mandiri, silakan tambahkan tag yang sesuai. Kami tidak (biasanya) menyelesaikan masalah seperti itu untuk Anda, melainkan membantu memandu Anda untuk solusi sendiri, yang secara umum akan memberi Anda pemahaman yang lebih baik tentang bagaimana menyelesaikan masalah tersebut di masa depan.
Jbowman
Tidak yakin tapi mungkin ini bisa membantu: en.wikipedia.org/wiki/Noncentral_beta_distribution
Sudahkah Anda mencoba membuat variabel kedua? MengatakanW=X+Y? Maka Anda bisa mendapatkan distribusi bersamaW,Z dan mengintegrasikan W untuk mendapatkan distribusi Z.
1
Saya tidak melihat di mana Anda menggunakan fakta bahwa fungsi kerapatan Beta adalah nol pada komplemen interval [0,1].
whuber
@whuber saya pikir saya menemukan kesalahan. Apakah Anda ingin memberikan jawaban lengkap atau saya melakukannya sendiri?
tam

Jawaban:

9

Setelah beberapa komentar yang berharga, saya dapat menemukan solusinya:

Kita punya fX(x)=1B(1,K1)(1x)K2 dan fY(y)=12KΓ(K)yK1ey/2.

Kami juga punya 0x1. Jadi, jikax=zy, kita mendapatkan 0zy1 yang menyiratkan itu zy.

Karenanya:

fZ=y=y=+1|y|fY(y)fX(zy)dy=z+1B(1,K1)2KΓ(K)1yyK1ey/2(1z/y)K2dy=1B(1,K1)2KΓ(K)z+ey/2(yz)K2dy=1B(1,K1)2KΓ(K)[2K1ez/2Γ(K1,yz2)]z=2K1B(1,K1)2KΓ(K)ez/2Γ(K1)=12ez/2
di mana kesetaraan terakhir berlaku sejak itu B(1,K1)=Γ(1)Γ(K1)Γ(K).

Begitu Z mengikuti distribusi parameter eksponensial 12; atau yang setara,Zχ22.

tam
sumber
8

Ada solusi statistik yang menyenangkan dan alami untuk masalah ini untuk nilai integral K, menunjukkan bahwa produk tersebut memiliki χ2(2)distribusi. Itu hanya bergantung pada hubungan-hubungan yang dikenal dan mudah diketahui di antara fungsi-fungsi variabel normal standar.

Kapan K adalah bagian integral, sebuah Beta(1,K1) distribusi muncul sebagai rasio

XX+Z
dimana X dan Z independen, X mempunyai sebuah χ2(2) distribusi, dan Z mempunyai sebuah χ2(2K2)distribusi. (Lihat artikel Wikipedia tentang distribusi Beta misalnya.)

Apa saja χ2(n) distribusi adalah jumlah dari kuadrat dari nindependen standar Normal variates. Karena itu,X+Z didistribusikan sebagai panjang kuadrat dari a 2+2K2=2K vektor dengan distribusi multinormal standar di R2K dan X/(X+Z) adalah panjang kuadrat dari dua komponen pertama ketika vektor diproyeksikan secara radial ke unit sphere S2K1.

Proyeksi multinormal standar n-Vektor ke unit sphere memiliki distribusi seragam karena distribusi multinormal simetris berbentuk bola. (Yaitu, itu tidak tetap di bawah kelompok ortogonal, hasil yang mengikuti segera dari dua fakta sederhana: (a), kelompok ortogonal memperbaiki asal dan menurut definisi tidak mengubah kovarian; dan (b) rerata dan kovarian sepenuhnya menentukan distribusi normal multivariat. Saya menggambarkan ini untuk kasus inin=3di https://stats.stackexchange.com/a/7984 ). Bahkan, simetri bola segera menunjukkan distribusi ini bersyarat pada panjang vektor asli. RasioX/(X+Z)oleh karena itu tidak tergantung pada panjangnya.

Yang disiratkan oleh semua ini adalah penggandaan itu X/(X+Z) oleh yang independen χ2(2K) variabel Y membuat variabel dengan distribusi yang sama dengan X/(X+Z) dikalikan dengan X+Z; kecerdasan, distribusiX, yang memiliki a χ2(2) distribusi.

whuber
sumber
Analogi yang sangat bagus! Saya merasa sedikit tidak yakin tentang paragraf terakhir karena penyederhanaan hanya terjadi karenaX+Zada di kedua sisi perkalian, yang tidak dapat bekerja untuk independen χ2(2K).
Xi'an
1
Tetapi setelah beberapa renungan lebih lanjut di métro Paris, saya menyadari itu karena X/(X+Z) dan (X+Z) independen, menggunakan (X+Z)×X/(X+Z) atau menggunakan Y×X/(X+Z)mengarah pada distribusi yang sama. Selamat!
Xi'an
1
tambahan: alasannya berlaku untuk non-integer K juga, jika seseorang mendefinisikan aχq2 sebagai Gamma Ga(q/2,1/2).
Xi'an
1
@ Xi'an Terima kasih atas komentar yang mengungkapkan itu. Memang, salah satu cara memanfaatkan pengakuan ituX/(X+Z) dan X+Z independen adalah untuk mengejar implikasi bahwa fungsi kepadatannya akan dapat dipisahkan: dan gagasan itu berlaku tanpa modifikasi pada kasus umum non-integral K. Bahkan bagi mereka yang lebih suka menghitung konvolusiXYsecara langsung, wawasan statistik ini menyarankan cara yang sederhana dan efektif untuk melanjutkan integrasi melalui perubahan variabel yang sesuai.
whuber
3

Saya sangat mencela taktik yang biasa digunakan untuk menemukan kepadatan Z=g(X,Y) oleh komputasi pertama menghitung kepadatan bersama Z dan X (atau Y) karena "mudah" untuk menggunakan Jacobian, dan kemudian mendapatkan fZsebagai kepadatan marginal (lih. jawaban Rusty Statistician). Adalah jauh lebih mudah untuk menemukan CDF dariZlangsung dan kemudian bedakan untuk menemukan pdf tersebut. Ini adalah pendekatan yang digunakan di bawah ini.

X dan Y adalah variabel acak independen dengan kepadatan fX(x)=(K1)(1x)K21(0,1)(x) dan fY(y)=12K(K1)!yK1ey/21(0,)(y). Lalu, denganZ=XY, kami punya untuk z>0,

P{Z>z}=P{XY>z}=y=z12K(K1)!yK1ey/2[x=zy1(K1)(1x)K2dx]dy=y=z12K(K1)!yK1ey/2(1zy)K1dy=y=z12K(K1)!(yz)K1ey/2dy=ez/2012K(K1)!tK1et/2dy   on setting yz=t=ez/2on noting that the integral is thatof a Gamma pdf

Sudah terkenal kalau VExponential(λ), kemudian P{V>v}=eλv. Karena ituZ=XY memiliki kerapatan eksponensial dengan parameter λ=12, yang juga merupakan χ2(2) distribusi.

Dilip Sarwate
sumber