Este es un adenda a la muy buena respuesta de @Macro que expone exactamente lo que hay que saber para determinar la varianza del el producto de dos variables aleatorias correlacionadas. Como \begin {align} \operatorname {var}(XY) &= E \left [(XY)^2 \right ] - \left (E[XY] \right )^2 \tag {1} \\ &= E[(XY)^2] - \left ( \operatorname {cov}(X,Y)+E[X]E[Y] \right )^2 \\ &= E[X^2Y^2] - \left ( \operatorname {cov}(X,Y)+E[X]E[Y] \right )^2 \tag {2} \\ &= \left ( \operatorname {cov}(X^2,Y^2)+E[X^2]E[Y^2] \right ) - \left ( \operatorname {cov}(X,Y)+E[X]E[Y] \right )^2 \tag {3} \\ \end {align} donde cov(X,Y)cov(X,Y) , E[X]E[X] , E[Y]E[Y] , E[X2]E[X2] y E[Y2]E[Y2] se puede suponer que son cantidades conocidas, tenemos que ser capaces de determinar el valor de E[X2Y2]E[X2Y2] en (2)(2) o cov(X2,Y2)cov(X2,Y2) en (3)(3) . Esto no es fácil de hacer en general, pero, como ya se ha señalado, si XX y YY son independiente variables aleatorias, entonces cov(X,Y)=cov(X2,Y2)=0cov(X,Y)=cov(X2,Y2)=0 . De hecho, dependencia, no la correlación (o la falta de ella) es la cuestión clave. Que sepamos que cov(X,Y)cov(X,Y) es igual a 00 en lugar de algún valor no nulo no lo hace, por sí mismo, ayuda en el menos en nuestros esfuerzos son determinar el valor de E[X2Y2]E[X2Y2] o cov(X2,Y2)cov(X2,Y2) a pesar de que hace simplificar los lados derechos de (2)(2) y (3)(3) un poco.
Cuando XX y YY son dependiente variables aleatorias, entonces en al menos una (bastante común o bastante importante) especial caso, es es posible encontrar el valor de E[X2Y2]E[X2Y2] con relativa facilidad.
Supongamos que XX y YY son conjuntamente normal variables aleatorias con coeficiente de correlación ρρ . Entonces, acondicionado en X=xX=x El condicional densidad de YY es una densidad normal con media E[Y]+ρ√var(Y)var(X)(x−E[X])E[Y]+ρ√var(Y)var(X)(x−E[X]) y la varianza var(Y)(1−ρ2)var(Y)(1−ρ2) . Así, \begin {align}E[X^2Y^2 \mid X] &= X^2E[Y^2 \mid X] \\ &= X^2 \left [ \operatorname {var}(Y)(1- \rho ^2) + \left (E[Y] + \rho\left.\left.\sqrt { \frac { \operatorname {var}(Y)}{ \operatorname {var}(X)}} \right (X-E[X] \right ) \right )^2 \right ] \end {align} que es un Cuartita función de XX , digamos que g(X)g(X) y la Ley de la Iteración Expectativa nos dice que E[X2Y2]=E[E[X2Y2∣X]]=E[g(X)] donde el lado derecho de (4) puede calcularse a partir del conocimiento de los 3er y 4to momento de X -- resultados estándar que se pueden encontrar en muchos textos y libros de referencia (lo que significa que soy demasiado perezoso para buscarlos e incluirlos en esta respuesta).
Más adición: En una respuesta ahora borrada, @Hydrologist da la varianza de XY como Var[xy]=(E[x])2Var[y]+(E[y])2Var[x]+2E[x]Cov[x,y2]+2E[y]Cov[x2,y]+2E[x]E[y]Cov[x,y]+Cov[x2,y2]−(Cov[x,y])2 y afirma que esta fórmula procede de dos artículos publicados hace medio siglo en JASA. Esta fórmula es una transcripción incorrecta de los resultados de los artículos citados por Hydrologist. En concreto, Cov[x2,y2] es un error de transcripción de E[(x−E[x])2(y−E[y])2] en el artículo de la revista, y de forma similar para Cov[x2,y] y Cov[x,y2] .