Estoy introduciendo esto como una "respuesta" a mi propia pregunta aunque probablemente no sea una gran respuesta. Pero espero que esto ayude a mostrar lo que busco. Lo editaré o lo borraré si me lo piden. Edición: Esto es todavía un trabajo en progreso. Probablemente debería mencionar que Stat-R, cardinal, StasK, y Macro me han guiado en cada paso del camino hasta ... donde me he detenido ahora. Por favor, siéntase libre de editar/corregir este borrador.
Estoy buscando ejemplos de distribuciones conjuntas de los rvs dependientes $X$ y $Y$ tal que (1) la distribución puede ser caracterizada, al menos aproximadamente, (2) lo siguiente puede ser calculado exactamente: $$ E[e^{aXY}] $$
Ejemplo 1. Sea: $$ X = x + \epsilon $$ $$ Y = y + \epsilon $$ $$ \epsilon \sim \mathcal{N}(0,\sigma^{2}) $$ donde $x$ y $y$ son números reales fijos (los medios conocidos de $X$ y $Y$ ).
Observación: $$ e^{\epsilon} \sim \mathcal{LN}(0,\sigma^{2}) $$ $$ E[e^{\epsilon}]=e^{\sigma^{2}/2} $$ Se deduce de una sustitución de $x+\epsilon$ y $y+\epsilon$ en $XY$ eso: $$ E[e^{aXY}] = e^{axy} \times E[e^{a(x+y)\epsilon} e^{a\epsilon^{2}}] $$
Edición: El $\epsilon^{2}$ no es despreciable, como se explica en los comentarios. Así,
$$ E[e^{aXY}] = e^{axy} \times E[e^{a(x+y)\epsilon}] \times E[e^{a\epsilon^{2}}] + e^{2axy} \times COV[e^{a(x+y)\epsilon},e^{a\epsilon^{2}}] $$
La primera expectativa de la derecha: $$ E[e^{a(x+y)\epsilon}] = e^{a^{2}(x+y)^{2}\sigma^{2}/2} $$
La segunda expectativa en el rhs presenta el cuadrado de una Normal, que es un Chi-cuadrado. Editar: Me han mostrado, en los comentarios, cómo calcular la expectativa explotando el hecho de que es una evaluación de la MGF de un chi-cuadrado, ya que $(\epsilon/\sigma)^{2}\sim\chi_{1}^{2}$ . Por lo tanto, $$ E[e^{a\epsilon^{2}}] = (1-2a/\sigma^{2})^{-k/2} = 1/\sqrt{1-2a/\sigma^{2}}, \ a < \sigma^{2}/2 $$ donde $k$ denota los grados de libertad, que en este caso es $k=1$ ya que el chi-cuadrado se obtiene de $1$ Normal al cuadrado. Si $a>\sigma^{2}/2$ entonces la expectativa diverge.
Por último, el término de covarianza. Creo que aquí tengo que calcular esto a partir de la integral, espero que algunos resultados conocidos puedan ser utilizados en el camino. Si alguien sabe cómo terminar esto sería genial, pero estoy atascado.
0 votos
No estoy ligado a la distribución lognormal, es sólo un ejemplo que me vino a la mente... ¿Es la exponencial-exponencial de una normal una lognormal de base cuadrada? (o esta última frase es una total tontería)
0 votos
¿Intentaste usar es.wikipedia.org/wiki/Función generadora de momentos ?
1 votos
Gracias Stat-R. No sé mucho sobre estadísticas. He mirado la página de mgf. Si entiendo, se trata de expresiones de la forma $e^{tX}$ , donde $X$ es un vector aleatorio pero $t$ es determinista, así que en efecto es un exponencial de una suma de rvs, en lugar de un exponencial de productos. Así que no estoy seguro de cómo proceder a partir de ahí...
0 votos
Gracias Macro. Sólo busco ejemplos en los que la distribución de $e^{aXY}$ se puede caracterizar, al menos aproximadamente, y calcular los dos primeros momentos de forma explícita. No tengo ninguna suposición sobre $X$ y $Y$ , excepto que son positivos y continuos (no discretos, se me olvidó insistir en eso). Anoté el primer momento cuando $X$ y $Y$ son LogNormal, pero es complicado. Pero eso no significa que no se pueda resolver: escribir el valor esperado de la Normal es muy lioso también, hasta que se sabe qué cambios de variable aplicar.
0 votos
@PatrickT: El mgf de una lognormal no existe para ningún positivo $t$ .
0 votos
Gracias cardenal, ¡creo que ahora entiendo lo que significa!
0 votos
He borrado este párrafo, porque parece que ahí me he equivocado de camino. Lo copio aquí para que se entiendan los comentarios que se referían a ese párrafo: "Por ejemplo, el producto $Z=XY$ de dos variables lognormales dependientes también es lognormal, si no me equivoco, y presumiblemente multiplicar por un número positivo no cambiaría eso, pero ¿qué pasa con la distribución de la exponencial de una lognormal? Ahora bien, eso sería la exponencial de la exponencial de una normal... ¿qué sabemos de esa distribución?"