Distribución de la relación entre dos variables aleatorias uniformes independientes

17

Supppse XX e YY están distribuidos uniformemente de manera uniforme en [ 0 , 1 ][0,1] , y son independientes, ¿cuál es el PDF de Z = Y / XZ=Y/X ?

La respuesta de algún libro de texto de teoría de la probabilidad es

f Z ( z ) = { 1 / 2 , si  0 z 1 1 / ( 2 z 2 ) , si  z > 1 0 , de lo contrario .

fZ(z)=1/2,1/(2z2),0,if 0z1if z>1otherwise.

Me pregunto, por simetría, no debe f Z ( 1 / 2 ) = f Z ( 2 )fZ(1/2)=fZ(2) ? Este no es el caso según el PDF anterior.

qed
fuente
¿Cuál es el dominio de XX e YY ?
Sobi
2
¿Por qué esperarías que esto sea cierto? La función de densidad se explica cómo apretada la probabilidad está en el entorno de un punto, y es claramente más difícil para ZZ para estar cerca de 22 a 1 / 21/2 (considerar por ejemplo que ZZ puede ser siempre 1 / 21/2 no importa qué XX es, pero Z < 2Z<2 cuando X > 1 / 2X>1/2 ).
dsaxton
3
No creo que sea un duplicado, esa pregunta está buscando el PDF, aquí tengo el PDF, solo estoy cuestionando su corrección (quizás bastante ingenua).
qed

Respuestas:

19

La lógica correcta es que con X , Y U ( 0 , 1 ) independientes , Z = YX,YU(0,1)X y Z-1=XZ=YXY tiene la mismadistribucióny por lo tanto para0<z<1 P { YZ1=XY0<z<1 Xz}= P { XYz}= P { YX1z }FZ(z)= 1 - F Z ( 1z )

P{YXz}FZ(z)=P{XYz}=P{YX1z}=1FZ(1z)
donde la ecuación con CDF utiliza el hecho de que YX es una variable aleatoria continua y, por lo tanto,P{Za}=P{Z>a}=1-FZ(a). Por lo tanto, el pdf deZsatisface fZ(z)=z-2fZ(z-1),YXP{Za}=P{Z>a}=1FZ(a)Z0 < z < 1. Así
fZ(z)=z2fZ(z1),0<z<1.
f Z ( 12 )=4fZ(2), y noffZ(12)=4fZ(2) Z ( 12 )=fZ( 2 ) como pensabas que debería ser.fZ(12)=fZ(2)
Dilip Sarwate
fuente
14

Esta distribución es simétrica, si la miras de la manera correcta.

La simetría que ha observado (correctamente) es que Y / X y X / Y = 1 / ( Y / X ) deben distribuirse de manera idéntica. Cuando trabajas con razones y potencias, realmente estás trabajando dentro del grupo multiplicativo de los números reales positivos. El análogo de la medida invariante de ubicación d λ = d x en los números reales aditivos R es la medida invariante de escala d μ = d x / xY/XX/Y=1/(Y/X)dλ=dxR dμ=dx/x en el grupo multiplicativo RR of positive real numbers. It has these desirable properties:

  1. dμdμ is invariant under the transformation xaxxax for any positive constant aa: dμ(ax)=d(ax)ax=dxx=dμ.

    dμ(ax)=d(ax)ax=dxx=dμ.
  2. dμdμ is covariant under the transformation xxbxxb for nonzero numbers bb: dμ(xb)=d(xb)xb=bxb1dxxb=bdxx=bdμ.

    dμ(xb)=d(xb)xb=bxb1dxxb=bdxx=bdμ.
  3. dμdμ is transformed into dλdλ via the exponential: dμ(ex)=dexex=exdxex=dx=dλ.

    dμ(ex)=dexex=exdxex=dx=dλ.
    Likewise, dλdλ is transformed back to dμdμ via the logarithm.

(3) establishes an isomorphism between the measured groups (R,+,dλ)(R,+,dλ) and (R,,dμ)(R,,dμ). The reflection xxxx on the additive space corresponds to the inversion x1/xx1/x on the multiplicative space, because ex=1/exex=1/ex.

Let's apply these observations by writing the probability element of Z=Y/XZ=Y/X in terms of dμdμ (understanding implicitly that z>0z>0) rather than dλdλ:

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

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

That is, the PDF with respect to the invariant measure dμdμ is gZ(z)gZ(z), proportional to zz when 0<z10<z1 and to 1/z1/z when 1z1z, close to what you had hoped.


This is not a mere one-off trick. Understanding the role of dμdμ makes many formulas look simpler and more natural. For instance, the probability element of the Gamma function with parameter kk, xk1exdxxk1exdx becomes xkexdμxkexdμ. It's easier to work with dμdμ than with dλdλ when transforming xx 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).

whuber
fuente
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 XX-YY plane, curves of constant Z=Y/XZ=Y/X are lines through the origin. (Y/XY/X is the slope.) One can read off the value of ZZ from a line through the origin by finding its intersection with the line X=1X=1. (If you've ever studied projective space: here XX is the homogenizing variable, so looking at values on the slice X=1X=1 is a relatively natural thing to do.)

Consider a small interval of ZZs, (a,b)(a,b). This interval can also be discussed on the line X=1X=1 as the line segment from (1,a)(1,a) to (1,b)(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](X,Y)U=[0,1]×[0,1], which is the region we're actually interested in. If 0a<b10a<b1, then the area of the triangle is 12(10)(ba)12(10)(ba), so keeping the length of the interval constant and sliding it up and down the line X=1X=1 (but not past 00 or 11), the area is the same, so the probability of picking an (X,Y)(X,Y) in the triangle is constant, so the probability of picking a ZZ in the interval is constant.

However, for b>1b>1, the boundary of the region UU turns away from the line X=1X=1 and the triangle is truncated. If 1a<b1a<b, the projections down lines through the origin from (1,a)(1,a) and (1,b)(1,b) to the upper boundary of UU are to the points (1/a,1)(1/a,1) and (1/b,1)(1/b,1). The resulting area of the triangle is 12(1a1b)(10)12(1a1b)(10). From this we see the area is not uniform and as we slide (a,b)(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)fZ(1/2) corresponds to a line that reaches X=1X=1, but fZ(2)fZ(2) does not, so the desired symmetry does not hold.

Eric Towers
fuente
3

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

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

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

and this is indeed the case.

qed
fuente
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
fuente