Distribusi


8

Saya sedang mengerjakan masalah berikut:

Biarkan dan menjadi variabel acak independen dengan kerapatan bersama mana . Biarkan U = \ min (X, Y) dan V = \ maks (X, Y) . Cari kepadatan bersama (U, V) dan karenanya menemukan pdf dari U + V .XYf(x)=αβαxα110<x<βα1U=min(X,Y)V=max(X,Y)(U,V)U+V

Sebagai U+V=X+Y , saya hanya bisa menemukan pdf dari X+Y untuk melihat apa pdf of U+V seharusnya.

Saya mendapatkan pdf dari T=X+Y menjadi

(1)fT(t)=f(ty)f(y)dy=α2β2αmax(tβ,0)min(t,β)(y(ty))α1dy10<t<2β

Tidak yakin apakah integral itu dapat disederhanakan.

Kembali ke pertanyaan aktual, pdf gabungan dari (U,V) diberikan oleh

fU,V(u,v)=2f(u)f(v)10<u<v<β=2α2β2α(uv)α110<u<v<β

Saya melakukan perubahan variabel (U,V)(W,Z) dimana W=U+V dan Z=U . Nilai absolut dari jacobian adalah kesatuan. Juga, 0<u<v<β0<z<w2<β . Jadi pdf marginal dari W adalah

(2)fW(w)=2α2β2α0w/2(z(wz))α1dz10<w<2β

Ada kemungkinan bahwa saya telah membuat beberapa kesalahan dalam dukungan yang tepat dari variabel acak. Mungkin juga integral tidak memiliki solusi dalam hal fungsi dasar. Bagaimanapun, saya tidak bisa melanjutkan dengan integral. Jadi saya bahkan tidak bisa memverifikasi bahwa memiliki pdf yang sama seperti . Tampaknya saya mendapatkan distribusi dan . Dan karena penasaran, apakah distribusi memiliki nama (dalam hal ini saya akan mencari konvolusi dua variabel acak tersebut)?W=U+VT=X+YWTX

Edit.

Melanjutkan dengan integral terakhir saya dapatkan dengan tangan

0w/2(z(wz))α1dz=w2α101/2tα1(1t)α1dt=w2α1I1/2(α,α)B(α,α)
di mana adalah fungsi beta lengkap regularized. Menggunakan properti , kita mendapatkan . Jadi akhirnya kita memilikiIxIx(a,b)=1I1x(b,a)I1/2(α,α)=12
0w/2(z(wz))α1dz=12w2α1B(α,α)

Ini menyiratkan bahwa

fW(w)=α2β2αB(α,α)w2α110<w<2β

Bahwa ini bukan kerapatan dalam kisaran diberikan mudah dilihat. Jadi saya merasa telah membuat kesalahan besar di suatu tempat. Saya telah memeriksa perhitungan saya dengan Mathematica dan mereka tampaknya setuju.w


@ Xi'an Dan jumlah varian beta independen tidak memiliki bentuk pdf tertutup mungkin?
StubbornAtom

@ Xi'an Jadi saya merasa tidak ada yang salah jika saya mengakhiri jawaban saya dengan integral terlepas dari apakah itu memiliki bentuk tertutup dalam hal beberapa fungsi khusus atau tidak?
StubbornAtom

1
Sebagai generalisasi dari stats.stackexchange.com/questions/41467 (kasus di mana ), pertanyaan ini kemungkinan dapat diselesaikan dengan menggunakan satu atau lebih dari berbagai teknik yang dijelaskan dalam utas itu. α=1
whuber

Saya keliru mengatakan bahwa , padahal sebenarnya cukup untuk menjadi kepadatan yang valid. Ini kadang-kadang disebut distribusi fungsi daya . Untuk itu adalah kepadatan beta, dan untuk itu adalah kepadatan yang seragam. α>1α>0fβ=1α=1
StubbornAtom

Jawaban:


5

Karena kita miliki ( ) dan dengan perubahan variabel pada integral kedua dari rhs Demikian pula, ketika

max(tβ,0)min(t,β)(y(ty))α1dy10<t<2β={0t(y(ty))α1dywhen 0tβtββ(y(ty))α1dywhen βt2β
t<β
max(tβ,0)min(t,β)(y(ty))α1dy=0t/2(y(ty))α1dy+t/2t(y(ty))α1dy
z=ty
max(tβ,0)min(t,β)(y(ty))α1dy=20t/2(y(ty))α1dy
t>β , lagi dengan perubahan variabel pada integral kedua dari rhs. Namun saya tidak dapat memulihkan ekspresi fungsional yang sama untuk kepadatan dalam kasus kedua ini , yaitu
tββ(y(ty))α1dy=tβt/2(y(ty))α1dy+t/2β(y(ty))α1dy=2t/2β(y(ty))α1dy
z=ty
20w/2(z(wz))α1dz

Sekarang, seperti yang ditunjukkan dalam pertanyaan, oleh perubahan skala, yang akan menyiratkan bahwa distribusi bunga memiliki kepadatan yang mengubahnya menjadi distribusi Beta diubah pada , maka dengan kepadatan

20w/2(z(wz))α1dzw2(α1)+1=w2α1
f(w)w2α110<w<2β
B(2α,1)(0,2β)
f(w)={2β}2αΓ(2α+1)Γ(2α)w2α110<w<2β=2α{2β}2αw2α110<w<2β

Ini muncul sebagai kontradiksi ketika mempertimbangkan jawaban luar biasa terperinci dari W. Huber , karena Seragam adalah Beta . Dan karena jumlah dua Seragam bukan variabel acak Beta , tetapi rv dengan kepadatan "tenda".B(1,1)B(2,1)

Selain itu: Secara umum jumlah dari varian Beta bukanlah varian Beta yang lain, "penjelasan" langsung ketika melihat Betas sebagai dua Gammas dinormalisasi oleh jumlah mereka. Menambahkan dua Betas melihat jumlah yang berbeda dalam penyebut.

Masalahnya adalah dengan penurunan kerapatan : karena perubahan variabel mengarah ke dan batasan indikatornya adalah Oleh karena itu, sebagai kesimpulan, yaitu (1) dan bukan ekspresi yang diusulkan (2).W=U+V

(U,V)2αβ2[uv]α1I0<u<v<β
(Z,W)=(U,U+V)
(Z,W)2αβ2[z(wz)]α1I0<z<wz<β
0<z2z<wz<βz>wβ0<wandw<2β
W2α2β2αmax{0,wβ}min{β,w/2}[z(wz)]α1dzI0<w<2β

1
Itulah yang saya tanyakan apakah itu setuju dengan (1) atau tidak. Anda mungkin harus menambahkan konstanta yang hilang di dan juga. Terima kasih, tidak heran saya mendapatkan semua hasil yang aneh itu. (Z,W)(U,V)
StubbornAtom
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.