Distribusi rasio antara dua variabel acak seragam independen

17

Misalkan dan adalah standar yang terdistribusi secara seragam dalam , dan mereka independen, berapakah PDF dari ?XY[0,1]Z=Y/X

Jawaban dari beberapa buku teks teori probabilitas adalah

fZ(z)={1/2,if 0z11/(2z2),if z>10,otherwise.

Saya bertanya-tanya, apakah simetri, bukankah seharusnya ? Ini bukan kasusnya menurut PDF di atas.fZ(1/2)=fZ(2)

qed
sumber
Apa domain X dan Y ?
Sobi
2
Mengapa Anda berharap ini benar? Fungsi kepadatan memberitahu Anda bagaimana padat probabilitas adalah di lingkungan titik, dan itu jelas lebih sulit bagi Z untuk berada di dekat 2 dari 1/2 (mempertimbangkan misalnya bahwa Z selalu dapat 1/2 tidak peduli apa X , tapi Z<2 ketika X>1/2 ).
dsaxton
2
Kemungkinan rangkap dari Distribusi rasio seragam: Apa yang salah?
Xi'an
3
Saya tidak berpikir itu duplikat, pertanyaan itu mencari PDF, di sini saya punya PDF, saya hanya mempertanyakan kebenarannya (mungkin agak naif).
qed

Jawaban:

19

Logika yang benar adalah bahwa dengan independen , Y U ( 0 , 1 )X,YU(0,1) , Z=YX dan Z1=XY memilikidistribusi yangsamadan untuk0<z<1

P{YXz}=P{XYz}=P{YX1z}FZ(z)=1FZ(1z)
mana persamaan dengan CDF menggunakan fakta bahwaYX adalah variabel acak kontinu sehinggaP{Za}=P{Z>a}=1FZ(a). Karenanya pdf dariZmemenuhi Jadi f Z ( 1
fZ(z)=z2fZ(z1),0<z<1.
, dan bukanfZ(1fZ(12)=4fZ(2) seperti yang Anda pikir seharusnya.fZ(12)=fZ(2)
Dilip Sarwate
sumber
14

Distribusi ini simetris - jika Anda melihatnya dengan cara yang benar.

Simetri yang telah Anda amati (dengan benar) adalah bahwa dan X / Y = 1 / ( Y / X ) harus didistribusikan secara identik. Ketika bekerja dengan rasio dan kekuatan, Anda benar-benar bekerja dalam kelompok multiplikasi angka nyata positif. Analog dari ukuran invarian lokasi d λ = d x pada bilangan real aditif R adalah ukuran invarian skala d μ = d x / xY/XX/Y=1/(Y/X)dλ=dxR dμ=dx/xpada kelompok multiplikasi dari bilangan real positif. Ini memiliki sifat-sifat yang diinginkan:R

  1. invarian di bawah transformasi x a x untuk konstanta positif a : d μ ( a x ) = d ( a x )dμxaxa

    dμ(ax)=d(ax)ax=dxx=dμ.
  2. adalah kovarian di bawah transformasi x x b untuk bilangan nol b : d μ ( x b ) = d ( x b )dμxxbb

    dμ(xb)=d(xb)xb=bxb1dxxb=bdxx=bdμ.
  3. ditransformasikan menjadi d λ melalui eksponensial: d μ ( e x ) = d e xdμdλ Demikian juga,dλditransformasikan kembali kedμmelalui logaritma.

    dμ(ex)=dexex=exdxex=dx=dλ.
    dλdμ

(3) membentuk isomorfisme antara kelompok yang diukur dan ( R , , d μ ) . Refleksi x - x pada ruang aditif sesuai dengan inversi x 1 / x pada ruang multiplikasi, karena e - x = 1 / e x .(R,+,dλ)(R,,dμ)xxx1/xex=1/ex

Mari kita terapkan pengamatan ini dengan menulis elemen probabilitas dalam hal d μ (memahami secara implisit bahwa z > 0 ) daripada d λ :Z=Y/Xdμz>0dλ

fZ(z)dz=gZ(z)dμ=12{1dz=zdμ,if 0z11z2dz=1zdμ,if z>1.

Yaitu, PDF sehubungan dengan ukuran invarian dμ adalah , sebanding dengan z ketika 0 < z 1 dan 1 / z ketika 1 z , dekat dengan apa yang Anda harapkan.gZ(z)z0<z11/z1z


Ini bukan hanya trik satu kali. Memahami peran membuat banyak formula terlihat lebih sederhana dan lebih alami. Sebagai contoh, elemen probabilitas fungsi Gamma dengan parameter k , x k - 1 e xdμk menjadi x k e x d μ . Lebih mudah untuk bekerja dengan d μ daripada dengan d λ ketika mentransformasikan x denganmengubahskala, mengambil kekuatan, atau secara eksponensial.xk1exdxxkexdμdμdλx

Gagasan tentang ukuran invarian pada suatu kelompok juga jauh lebih umum, dan memiliki aplikasi dalam bidang statistik tersebut di mana masalah menunjukkan beberapa invarian di bawah kelompok transformasi (seperti perubahan satuan ukuran, rotasi dalam dimensi yang lebih tinggi, dan sebagainya. ).

whuber
sumber
3
Sepertinya jawaban yang sangat berwawasan. Sayang sekali aku tidak memahaminya saat ini. Saya akan periksa kembali nanti.
qed
4

Jika Anda berpikir secara geometris ...

Dalam bidang - Y , kurva konstanta Z = Y / X adalah garis melalui titik asal. ( Y / X adalah kemiringan.) Seseorang dapat membaca nilai Z dari garis melalui titik asal dengan menemukan persimpangannya dengan garis X = 1 . (Jika Anda pernah mempelajari ruang proyektif: di sini X adalah variabel homogenisasi, jadi melihat nilai pada slice X = 1 adalah hal yang relatif wajar untuk dilakukan.)XYZ=Y/XY/XZX=1XX=1

Pertimbangkan interval kecil s, ( a , b ) . Interval ini juga dapat didiskusikan pada garis X = 1 sebagai segmen garis dari ( 1 , a ) hingga ( 1 , b ) . Himpunan garis melalui titik asal melewati interval ini membentuk segitiga padat di bujur sangkar ( X , Y ) U = [ 0 , 1 ] × [ 0 , 1 ]Z(a,b)X=1(1,a)(1,b)(X,Y)U=[0,1]×[0,1], yang merupakan wilayah tempat kami tertarik. Jika , maka luas segitiga adalah 10a<b1, sehingga menjaga panjang interval konstan dan menggesernya ke atas dan ke bawah garisX=1(tetapi tidak melewati0atau1), area tersebut sama, sehingga probabilitas memilih(X,Y)dalam segitiga adalah konstan, sehingga probabilitas memilihZdalam interval adalah konstan.12(10)(ba)X=101(X,Y)Z

b>1UX=11a<b(1,a)(1,b)U(1/a,1)(1/b,1)12(1a1b)(10)(a,b)

Then the same algebra demonstrated in other answers finishes the problem. In particular, returning to the OP's last question, fZ(1/2) corresponds to a line that reaches X=1, but fZ(2) does not, so the desired symmetry does not hold.

Eric Towers
sumber
3

Just for the record, my intuition was totally wrong. We are talking about density, not probability. The right logic is to check that

1kfZ(z)dz=1/k1fZ(z)=12(11k)
,

and this is indeed the case.

qed
sumber
1

Yea the link Distribution of a ratio of uniforms: What is wrong? provides CDF of Z=Y/X. The PDF here is just derivative of the CDF. So the formula is correct. I think your problem lies in the assumption that you think Z is "symmetric" around 1. However this is not true. Intuitively Z should be a skewed distribution, for example it is useful to think when Y is a fixed number between (0,1) and X is a number close to 0, thus the ratio would be going to infinity. So the symmetry of distribution is not true. I hope this help a bit.

lzstat
sumber