Autor Tema: Chebyschev

0 Usuarios y 1 Visitante están viendo este tema.

24 Abril, 2017, 08:59 am
Respuesta #50

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 56,005
  • País: es
  • Karma: +0/-0
Hola

La primera cota cómo es este paso

\( E[E[1_{X+y<0}|Y=y]]\leq E\left[\dfrac{\sigma_X^2}{\sigma_X^2+(\mu_X+Y)^2}\right] \). Puede ser

\( P(X+y<0|Y=y)\leq{}\dfrac{\sigma_X^2}{\sigma_X^2+(\mu_X+y)^2} \)?

Pero no tiene que ir dividido por la \( P(Y=y) \)?

Más o menos. Ten en cuenta que por ejemplo para Y discreta tendríamos:

\(
P(X+y<0|Y=y)=\dfrac{P(X+y<0\cap Y=y)}{P(Y=y)} \)

por la independencia de \( X,Y \) en el numerador:

\(
P(X+y<0\cap Y=y)=P(X+Y<0)P(Y=y) \)

y \( P(Y=y)  \)se cancela con el denominador.

Digo más o menos porque escribirlo así, usando \( P(Y=y) \) requiere que esa probabilidad sea no nula (imposible por ejemplo en una variable continua). Entonces de manera más general, por la independencia de \( X \) e \( Y \):

\( E[1_{X+y<0}|Y=y]=E[1_{X+y<0}] \)

y después:

\( E[1_{X+y<0}]=P(X+y<0)\leq \dfrac{\sigma_X^2}{\sigma_X^2+(\mu_X+y)^2} \)

Citar
En la prueba de la segunda, no es que el soporte de \( Y \) estaba acotado superiormente, de forma que \( g(y) \) fuera cóncava, debe haber una errata ahí.


Es que me quedé en ese paso. En lo que escribí simplemente he usado que el soporte de \( Y \) está en \( [-\mu_X,+\infty) \).

Continuando. Si llamamos:

\( g(y)=\dfrac{\sigma_X^2}{\sigma_X^2+(\mu_X+y)^2} \)

Comprueba que su segunda derivada es negativa en el intervalo \( [-\mu_X,-\mu_X+\color{red}\sigma_x/\sqrt{3}\color{black}] \)

y por tanto cóncava.

Entonces si \( Y \) tiene el soporte en \( [-\mu_X,-\mu_X+\color{red}\sigma_x/\sqrt{3}\color{black}] \), por la desigualdad de Jensen:

\( E[g(Y)]\leq g(E[Y])=g(\mu_Y) \)

De ahí directamente:

\( E\left[\dfrac{\sigma_X^2}{\sigma_X^2+(\mu_X+Y)^2}\right]\leq \dfrac{\sigma_X^2}{\sigma_X^2+(\mu_X+\mu_Y)^2} \)

Citar
Y si puedes ser un poco más específico en la prueba de la cota 3.

Continuará...

Saludos.

CORREGIDO

24 Abril, 2017, 02:15 pm
Respuesta #51

Quema

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,974
  • País: uy
  • Karma: +0/-0
  • Sexo: Masculino
Hola

Ahora si me queda claro. Una pregunta, cuando se exige que el soporte de \( Y \) dependa de la media de \( X \), de alguna forma eso no invalida que \( X,Y \) sean independientes?

Saludos

24 Abril, 2017, 04:05 pm
Respuesta #52

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 56,005
  • País: es
  • Karma: +0/-0
Hola

Ahora si me queda claro. Una pregunta, cuando se exige que el soporte de \( Y \) dependa de la media de \( X \), de alguna forma eso no invalida que \( X,Y \) sean independientes?

No. No invalida. Simplemente restringe las condiciones bajo las cuales podemos aplicar nuestra cota.

Pero que sean independientes o no, no tiene nada que ver con que sus medias y varianzas estén más o menos próximas o cumplan tal o cual relación; tiene que ver con si la distribución marginal de una variable varia o no cuando conocemos datos sobre la otra.

Saludos.

24 Abril, 2017, 06:15 pm
Respuesta #53

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 56,005
  • País: es
  • Karma: +0/-0
Hola

 Para la última cota. Con la hipótesis de que \( sop(Y)\subset [-\mu_X,+\infty) \) partimos de que:

\( P(X+Y\leq 0)\leq E[g(Y)] \)

 con:

\( g(y)=\dfrac{\sigma_X^2}{\sigma_X^2+(\mu_X+y)^2} \)

 Notamos que esa función es siempre positiva, tiene un máximo en \( y=-\mu_X \) igual a \( 1 \) y luego decrece. Buscamos una parábola cóncava (una función de grado dos para que luego la esperanza de \( h(Y) \) quede en  función de la media y la varianza) \( h(y) \) tal que:

\( g(y)\leq h(y) \) en el soporte de \( Y \)

 y luego usar que \( E[g(Y)]\leq E[h(Y)] \)

 Esto ya nos obliga a que el soporte sea acotado, ya que una parábola \( h(y) \) cóncava toma valores negativos para \( y \) alto, mientras que \( g(y) \) siempre es positivo.

 Para que \( h(y) \) tenga también el máximo en \( y=-\mu_X \) tomamos:

\(  h(y)=1-a(y+\mu_X)^2 \)

 Operando observamos que:

\( h(y)-g(y)=(\mu_x+y)^2\left(\dfrac{1}{\sigma_X^2-(\mu_X+y)^2}-a\right) \) (**)

 Si el soporte de \( Y \) es \( [-\mu_X,-\mu_X+b] \) para que (**) sea positivo en ese soporte necesitamos que \( h(b)-g(b)=0 \). Despejando obtenemos:

\( a=\dfrac{1}{\sigma_X^2+b^2} \)

 Por tanto:

\(  h(y)=1-a(\mu_x+y)^2=1-\dfrac{1}{\sigma_X^2+b^2}(\mu_X^2+2y\mu_X+y^2) \)

 y

 \( E[h(Y)]=1-\dfrac{1}{\sigma_X^2+b^2}(\mu_X^2+2\mu_Y\mu_X+E[Y^2]) \)

 Recordemos que \( E[Y^2]=\sigma_Y^2+\mu_Y^2 \). Queda entonces:

 \( E[h(Y)]=1-\dfrac{1}{\sigma_X^2+b^2}(\mu_X^2+2\mu_Y\mu_X+\sigma_Y^2+\mu_Y^2) \)

 \( E[h(Y)]=1-\dfrac{1}{\sigma_X^2+b^2}((\mu_X+\mu_Y)^2+\sigma_Y^2) \)

 De ahí la cota 3 (me sobraba una hipótesis que he tachado).

Saludos.


24 Abril, 2017, 07:36 pm
Respuesta #54

Quema

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,974
  • País: uy
  • Karma: +0/-0
  • Sexo: Masculino
Hola

Lo voy a mirar detenidamente. El \( b\geq{}0 \)?

Saludos

24 Abril, 2017, 07:38 pm
Respuesta #55

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 56,005
  • País: es
  • Karma: +0/-0
Hola

Lo voy a mirar detenidamente. El \( b\geq{}0 \)?

Si.

Saludos.

24 Abril, 2017, 07:53 pm
Respuesta #56

Quema

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,974
  • País: uy
  • Karma: +0/-0
  • Sexo: Masculino
Hola,

Como \( h(y) \) es cóncava y \( g(y) \) es casi siempre convexa, parece que queda bastante espacio libre entre las dos. No se tendría que aproximar \( g(y) \) con una función convexa o al menos con tramos convexos para que se aproxime más?

Y si expandimos \( g(y) \) por Taylor? Y utilizamos el polinomio cuadrático de esta expansión? Viendo si cumple que es mayor o igual a \( g(y) \).

Saludos

24 Abril, 2017, 08:08 pm
Respuesta #57

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 56,005
  • País: es
  • Karma: +0/-0
Hola

Citar
Como \( h(y) \) es cóncava y \( g(y) \) es casi siempre convexa, parece que queda bastante espacio libre entre las dos. No se tendría que aproximar \( g(y) \) con una función convexa o al menos con tramos convexos para que se aproxime más?

Lo de casi siempre convexa es relativo; tiene forma de campana y es convexa para valores de \( y \) grandes (o muy pequeños), pero cóncava en su parte central, la de valores más altos y por tanto más significativos.

Ten en cuenta además que tenemos que acotar por una función de grado dos; sino el método no funciona. Ya que al final \( E[h(y)]  \)tiene que quedar en función de \( E[Y] \) e \( E[Y^2] \).

Tendría sentido acotar por una parábola convexa si \( Y \) tuviese su soporte en la zona donde la función \( g(Y) \) a la que nos referimos es convexa, es decir, mucho más a la derecha de \( -\mu_X. \) Se podría explorar esa otra posibilidad.

Saludos.

24 Abril, 2017, 08:13 pm
Respuesta #58

Quema

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,974
  • País: uy
  • Karma: +0/-0
  • Sexo: Masculino
Hola

Si claro, si pudiéramos encontrar una cota óptima, aún tirando hacia la derecha el soporte de \( Y \) sería bueno. Creo que se puede utilizar la desigualdad de Williassen, que es para funciones crecientes y cóncavas sabiendo la media y varianza a la parte convexa de \( g(y) \) y cómo ésta cota es óptima obtendríamos una cota óptima. No?

Saludos

26 Abril, 2017, 12:18 pm
Respuesta #59

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 56,005
  • País: es
  • Karma: +0/-0
Hola

Si claro, si pudiéramos encontrar una cota óptima, aún tirando hacia la derecha el soporte de \( Y \) sería bueno. Creo que se puede utilizar la desigualdad de Williassen, que es para funciones crecientes y cóncavas sabiendo la media y varianza a la parte convexa de \( g(y) \)

La desigualdad de Williassen para una función convexa decreciente \( g(y) \) con soporte en \( [a,+\infty) \) si no me equivoco queda así:

\( E[g(y)]\leq \dfrac{(\mu_Y-a)^2}{\sigma_Y^2+(\mu_Y-a)^2}g\left(\mu_x+\dfrac{\sigma_Y^2}{\mu_Y-a}\right)+
\dfrac{\sigma_Y^2}{\sigma_Y^2+(\mu_Y-a)^2}g(a) \)

En nuestro caso habría que aplicarlo para la función:

\( g(y)=\dfrac{\sigma_X^2}{\sigma_X^2+(\mu_X+y)^2} \)

y \( a=-\mu_X+\sigma_X/\sqrt{3} \), bajo el supuesto de que el soporte de \( Y \) está en \( [a,+\infty). \)

Citar
y cómo ésta cota es óptima obtendríamos una cota óptima. No?

No. O al menos no necesariamente. Eso es lo que traté de explicar aquí:

P.D. En cuanto a la posiblidad de encontrar cotas óptimas conocida la varianza de \( Y \), habría que cambiar de estrategia. En todo lo que hemos hecho usamos de manera indpendiente la desigualdad de Cantelli para la variable \( X \). Está desigualdad es óptima pero la dsitribuciónd e \( X \) que alcanza el óptimo para \( P(X+a\leq 0) \) depende de \( a \). Si la usamos para \( P(X+Y\leq 0) \) solo podemos aspirar a cota óptima cuando \( Y \) puede tomar un sólo valor; si exigimos varianza de \( Y \) no nula esto ya no se da.

Estamos combinando dos acotaciones: una para la variable \( X \), la de Cantelli; otra para la variable \( Y \) (mi primer intento o la de Williassen ahora). Las dos por separado son óptimas. El problema es que la distribución que da la cota óptima para Cantelli depende de el valor de \( y \). A su vez la distribución óptima de Williassen toma dos valores de \( Y \) distintos, entonces al menos para uno de ellos no será el valor de \( y \) que hemos usado para la optimización de la cota de Cantelli.

Saludos.