He oído que las razones o inversas de variables aleatorias a menudo son problemáticas, al no tener expectativas. ¿Porqué es eso?


Respuestas:


25

Me gustaría ofrecer una explicación muy simple e intuitiva. Equivale a mirar una imagen: el resto de esta publicación explica la imagen y saca conclusiones de ella.

Esto es a lo que se reduce: cuando hay una "masa de probabilidad" concentrada cerca de , habrá demasiada probabilidad cerca de , lo que hace que su expectativa sea indefinida.X=01/X±


En lugar de ser completamente general, concentrémonos en variables aleatorias que tienen densidades continuas en una vecindad de . Supongamos que . Visualmente, estas condiciones significan que la gráfica de encuentra sobre el eje alrededor de :XfX0fX(0)0f0

Figura que muestra el gráfico de una densidad y el área debajo de ella.

La continuidad de alrededor de implica que para cualquier altura positiva menor que y suficientemente pequeña , podemos tallar un rectángulo debajo de este gráfico que se centra alrededor de , tiene un ancho y una altura , como se muestra. Esto corresponde a expresar la distribución original como una mezcla de una distribución uniforme (con peso ) y lo que queda. 0 p f X ( 0 ) ϵ x = 0 2 ϵ p p × 2 ϵ = 2 p ϵfX0pfX(0)ϵx=02ϵpp×2ϵ=2pϵ

Figura que muestra el gráfico como una mezcla.

En otras palabras, podemos pensar que surge de la siguiente manera:X

  1. Con probabilidad , dibuje un valor de una distribución Uniforme .( - ϵ , ϵ )2pϵ(ϵ,ϵ)

  2. De lo contrario, extraiga un valor de la distribución cuya densidad sea proporcional a . (Esta es la función dibujada en amarillo a la derecha).fXpI(ϵ,ϵ)

( es la función del indicador).I

Paso muestra que para cualquier , la probabilidad de que está entre y excede . De manera equivalente, esta es la posibilidad de que exceda de . Para decirlo de otra manera: escribir para la función de supervivencia de0 < u < ϵ X 0 u p u / 2 1 / X 1 / u S 1 / X(1)0<u<ϵX0upu/21/X1/uS1/X

S(x)=Pr(1/X>x),

la imagen muestra para todos los .S(x)>p/(2x)x>1/ϵ

Ya hemos terminado, porque este hecho sobre implica que la expectativa no está definida. S Compare las integrales involucradas en el cálculo de la expectativa de la parte positiva de , :( 1 / X ) + = max ( 0 , 1 / X )1/X(1/X)+=max(0,1/X)

E[(1/X)+]=0S(x)dx>1/ϵxS(x)dx>1/ϵxp2xdx=p2log(xϵ).

(Este es un argumento puramente geométrico: cada integral representa una región bidimensional identificable y todas las desigualdades surgen de inclusiones estrictas dentro de esas regiones. De hecho, ni siquiera necesitamos saber que la integral final es un logaritmo: hay formas geométricas simples Los argumentos que muestran esta integral divergen.)

Dado que el lado derecho diverge como , diverge. La situación con la parte negativa de es la misma (porque el rectángulo está centrado alrededor de ), y el mismo argumento muestra que la parte negativa de diverge. En consecuencia, la expectativa de sí no está definida.E [ ( 1 / X ) + ] 1 / X 0 1 / X 1 / XxE[(1/X)+]1/X0 01/ /X1/ /X

Por cierto, el mismo argumento muestra que cuando tiene la probabilidad concentrada en un lado de , como cualquier distribución exponencial o gamma (con un parámetro de forma menor que ), entonces la expectativa positiva diverge, pero la expectativa negativa es cero. En este caso la expectativa está definida, pero es infinita.0 1X0 01


3
¿Tengo razón al sospechar que la suposición es crucial para el resultado? Quiero decir, tenemos casos donde tiene momentos al menos para algún rango de parámetros involucrados, y parece que es en casos donde , como Gamma / Inverse-Gamma1 / X f X ( 0 ) = 0FX(0 0)0 01/ /XFX(0 0)=0 0
Alecos Papadopoulos

3
@ Alecos No, esa suposición no es crucial. Eso y la continuidad de en hacen que el argumento sea simple, pero ninguno de los dos es esencial. Considere una con densidad proporcional a para y . Esto es continuo en pero no tiene expectativas. 0 X f X - 1 / log ( x ) 0 < x < 1 / e f X ( 0 ) = 0 0 1 / XF0 0XFX-1/ /Iniciar sesión(X)0 0<X<1/ /miFX(0 0)=0 00 01/ /X
whuber

15

Las razones e inversas son principalmente significativas con variables aleatorias no negativas, por lo que asumiré casi con seguridad. Entonces, si es una variable discreta que toma el valor cero con probabilidad positiva, estaremos dividiendo con cero con una probabilidad positiva, lo que explica por qué la expectativa de no existirá.X 1 / XX0 0X1/ /X

Ahora observe el caso de distribución continua, con una variable aleatoria con función de densidad . Asumiremos que y que es continua (al menos en cero). Luego hay un tal que para . El valor esperado de viene dado por Ahora cambiemos la variable de integración a , tenemos , obteniendo f ( x ) f ( 0 ) > 0 f ϵ > 0 f ( x ) > ϵ 0 x < ϵ 1 / X E 1X0 0F(X)f(0)>0fϵ>0f(x)>ϵ0x<ϵ1/Xu = 1 / x d u = - 1

E1X=01xf(x)dx
u=1/xE 1du=1x2dxf ( u ) > ϵ [ 0 , ϵ ) f ( 1
E1X=0uf(1u)(1u)2du=01uf(1u)du
Ahora, suponiendo que en entonces en , usando esto tenemos mostrando que la expectativa no existe. Un ejemplo que cumple este supuesto es la distribución exponencial con tasa 1.f(u)>ϵ[0,ϵ)(1/ϵ,)E1f(1u)>1/ϵ(1/ϵ,)
E1X>ϵ1/ϵ1udu=

Hemos dado una respuesta para inversas, ¿qué pasa con las razones? Sea la razón de dos variables aleatorias no negativas. Si son independientes, podemos escribir así que esto se reduce al primer caso y no hay mucho nuevo que decir . ¿Qué pasa si son dependientes, con factorización de densidad conjunta como Entonces obtenemos (usando la misma sustitución que arriba) y podemos razonar como anteriormente en la integral interna. El resultado será que si la densidad condicional (dadaE Z = E YZ=Y/X f(x,y)=f(xy)g(y)EY

EZ=EYX=EYE1x
f(x,y)=f(xy)g(y)
y y 1 / X Y / X
EYX=0y01xf(xy)dxg(y)dy=0y01uf(1uy)dug(y)dy
y) es positivo y continuo en cero, para un conjunto de 's con probabilidad marginal positiva, la expectativa será infinita. Supongo que no será fácil encontrar ejemplos en los que la expectativa marginal de sea ​​infinita, pero la expectativa de la relación sea ​​finita, a menos que haya una correlación perfecta. ¡Sería bueno ver algunos ejemplos!y1/XY/X
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.