Розподіл відношення між двома незалежними рівномірними випадковими змінними


17

Supppse X і Y стандартно рівномірно розподілені в [0,1] , і вони незалежні, що таке PDF у форматі Z=Y/X ?

Відповідь з деякого підручника з теорії ймовірностей є

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

Я цікаво, в силу симетрії, не повинні fZ(1/2)=fZ(2) ? Це не так, як зазначено у PDF вище.


Що таке область і Y ? XY
Собі


2
Чому ви очікували, що це правда? Функція щільності говорить вам , наскільки щільно упаковані ймовірність в околиці точки, і це, очевидно , більш важким для , щоб бути поруч з 2 , ніж 1 / 2 (вважають, наприклад , що Z завжди може бути 1 / 2 незалежно від того , що X є, але Z < 2 , коли Х > 1 / 2 ). Z21/2Z1/2XZ<2X>1/2
dsaxton


3
Я не думаю, що це дублікат, це питання шукає PDF, тут у мене є PDF, я просто сумніваюся у його правильності (можливо, досить наївно).
qed

Відповіді:


19

Правильна логіка полягає в тому, що при незалежних , Z = YX,YU(0,1) і Z-1=XZ=YX однаковийрозподілі так для0<z<1 P { YZ1=XY0<z<1 де рівняння з CDF використовує той факт, щоY

P{YXz}=P{XYz}=P{YX1z}FZ(z)=1FZ(1z)
- неперервна випадкова величина і томуP{Za}=P{Z>a}=1-FZ(a). Отже, pdf відZзадовольняє fZ(z)=z-2fZ(z-1),YXP{Za}=P{Z>a}=1FZ(a)Z Отже, f Z ( 1
fZ(z)=z2fZ(z1),0<z<1.
, а неfZ(1fZ(12)=4fZ(2) як ти думав, що має бути.fZ(12)=fZ(2)

14

Цей розподіл є симетричним - якщо ви подивіться на це правильний шлях.

Симетрія, яку ви (правильно) спостерігали, полягає в тому, що і X / Y = 1 / ( Y / X ) повинні бути розподілені однаково. Працюючи з співвідношеннями та потужностями, ви дійсно працюєте в межах мультиплікативної групи позитивних дійсних чисел. Аналогом міри інваріантного розташування d λ = d x на додаткові дійсні числа R є міра інваріантної міри d μ = d x / xY/XX/Y=1/(Y/X)dλ=dxR dμ=dx/x для мультиплікативної групи додатних дійсних чисел. Він має такі бажані властивості:R

  1. є інваріантним при перетворенні x a x для будь-якої позитивної константи a : d μ ( a x ) = d ( a x )dμxaxa

    dμ(ax)=d(ax)ax=dxx=dμ.
  2. є коваріантним при перетворенні x x b для ненульових чисел b : d μ ( x b ) = d ( x b )dμxxbb

    dμ(xb)=d(xb)xb=bxb1dxxb=bdxx=bdμ.
  3. перетворюється на d λ за допомогою експоненціалу: d μ ( e x ) = d e xdμdλ Так самоdλперетворюється назад вdμчерез логарифм.

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

(3) встановлює ізоморфізм між вимірюваними групами та ( R , , d µ ) . Відображення x - x на адитивному просторі відповідає інверсії x 1 / x на мультиплікативному просторі, тому що e - x = 1 / e x .(R,+,dλ)(R,,dμ)xxx1/xex=1/ex

Застосуємо ці спостереження, записавши елемент ймовірності через d μ (розуміючи неявно, що z > 0 ), а не d λ :Z=Y/Xdμz>0dλ

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

dμgZ(z)z0<z11/z1z, close to what you had hoped.


This is not a mere one-off trick. Understanding the role of dμ makes many formulas look simpler and more natural. For instance, the probability element of the Gamma function with parameter k, xk1exdx becomes xkexdμ. It's easier to work with dμ than with dλ when transforming x by rescaling, taking powers, or exponentiating.

The idea of an invariant measure on a group is far more general, too, and has applications in that area of statistics where problems exhibit some invariance under groups of transformations (such as changes of units of measure, rotations in higher dimensions, and so on).


3
Looks like a very insightful answer. It's a pity I don't understand it at the moment. I will check back later.
qed

4

If you think geometrically...

In the X-Y plane, curves of constant Z=Y/X are lines through the origin. (Y/X is the slope.) One can read off the value of Z from a line through the origin by finding its intersection with the line X=1. (If you've ever studied projective space: here X is the homogenizing variable, so looking at values on the slice X=1 is a relatively natural thing to do.)

Consider a small interval of Zs, (a,b). This interval can also be discussed on the line X=1 as the line segment from (1,a) to (1,b). The set of lines through the origin passing through this interval forms a solid triangle in the square (X,Y)U=[0,1]×[0,1], which is the region we're actually interested in. If 0a<b1, then the area of the triangle is 12(10)(ba), so keeping the length of the interval constant and sliding it up and down the line X=1 (but not past 0 or 1), the area is the same, so the probability of picking an (X,Y) in the triangle is constant, so the probability of picking a Z in the interval is constant.

However, for b>1, the boundary of the region U turns away from the line X=1 and the triangle is truncated. If 1a<b, the projections down lines through the origin from (1,a) and (1,b) to the upper boundary of U are to the points (1/a,1) and (1/b,1). The resulting area of the triangle is 12(1a1b)(10). From this we see the area is not uniform and as we slide (a,b) further and further to the right, the probability of selecting a point in the triangle decreases to zero.

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.


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.


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.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.