5 Transformaciones de variables aleatorias. Univariantes y multivariantes.

Supongamos que a partir de la v.a. XX definimos una nueva v.a. Y=h(X)Y=h(X). Queremos saber cual es el efecto de la transformacion sbore la esperanza, la varianza y la distribucion de probabilidad.

Esto ya lo vimos con una transformacion lineal de XX, Y=a+bXY=a+bX, pero, que pasa si la transformacion no es lineal?

Proposición III.4.

Si XX es una v.a. continua y la transformacion hh tiene inversa, entonces se verifica que la densidad de Y=h(X)Y=h(X) se puede obtener como

fY(y)=fX(h1(y))|dh1(y)/dy|,f_{Y}(y)=f_{X}(h^{-1}(y))\left|dh^{-1}(y)/dy\right|,

donde |dh1(y)/dy|\left|dh^{-1}(y)/dy\right| es el modulo de la derivada de la inversa de la transformacion Y=h(X)Y=h(X).

Ejemplo III.5.

Supongamos que XX tiene densidad exe^{-x} si x>0x>0, consideremos la nueva variable Y=bXY=bX, cual es la distribucion de YY?

Para calcular la densidad de YY usando el resultado anterior, tenemos que la transformacion es Y=h(X)=b(X)Y=h(X)=b(X), entonces la transformacion inversa es

X=h1(Y)=YbX=h^{-1}(Y)=\frac{Y}{b}

La derivada con respecto a YY es |dh1(y)/dy|=1/b\left|dh^{-1}(y)/dy\right|=1/b, con lo que la densidad de la nueva variable es

fY(y)=fX(h1(y))|dh1(y)/dy|=1bey/b.f_{Y}(y)=f_{X}(h^{-1}(y))\left|dh^{-1}(y)/dy\right|=\frac{1}{b}e^{-y/b}% \mathchar 46\relax

Supongamos ahora que (X,Y)(X,Y) es una variable aleatoria bidimensional, con una funcion de densidad fX,Y(x,y)f_{X,Y}(x,y).

Consideremos una nueva variable aleatoria bidimensional

U=g1(X,Y) y V=g2(X,Y)U=g_{1}(X,Y)\text{ y }V=g_{2}(X,Y)

con g1g_{1} y g2g_{2} funciones especificas bivariantes. Supongamos ademas que existe la inversa de estas funciones, es decir, que podemos calcular

X=h1(U,V) e Y=h2(U,V)X=h_{1}(U,V)\text{ e }Y=h_{2}(U,V)

Estamos interesados en calcular la funcion de densidad bidimensional del vector (U,V)(U,V).

Existe un resultado similar al del cambio de variable unidimensional, donde se utiliza el Jacobiano de la transformacion:

J=|xuxvyuyv|=|xuyvxvyu|J=\begin{vmatrix}\frac{\partial x}{\partial u}&\frac{\partial x}{\partial v}\\ \frac{\partial y}{\partial u}&\frac{\partial y}{\partial v}\end{vmatrix}=\left% |\frac{\partial x}{\partial u}\frac{\partial y}{\partial v}-\frac{\partial x}{% \partial v}\frac{\partial y}{\partial u}\right|

con

xu=h1(u,v)u,xv=h1(u,v)v,yu=h2(u,v)u,yv=h2(u,v)v\frac{\partial x}{\partial u}=\frac{\partial h_{1}(u,v)}{\partial u},\;\frac{% \partial x}{\partial v}=\frac{\partial h_{1}(u,v)}{\partial v},\;\frac{% \partial y}{\partial u}=\frac{\partial h_{2}(u,v)}{\partial u},\;\frac{% \partial y}{\partial v}=\frac{\partial h_{2}(u,v)}{\partial v}

La densidad del nuevo vector aleatorio (U,V)(U,V) viene dado por

fU,V(u,v)=fX,Y(h1(u,v),h2(u,v))|J|f_{U,V}(u,v)=f_{X,Y}(h_{1}(u,v),h_{2}(u,v))\left|J\right|
Ejemplo III.6.

Supongamos que (X,Y)(X,Y) siguen distribuciones exponenciales independientes con el mismo parametro λ\lambda, es decir

f(x,y)={λ2eλxeλyx>0,y>00en el restof(x,y)=\begin{cases}\lambda^{2}e^{-\lambda x}e^{-\lambda y}&x>0,y>0\\ 0&\text{en el resto}\end{cases}

Supongamos que nos interesa calcular como se distribuye la suma de v.a. exponenciales, es decir, U=X+YU=X+Y, como involucra a las dos variables, tenemos que tener en cuenta un cambio bivariante, y despues calcular la marginal, integrando con respecto a la otra variable. Definamos por ejemplo V=XV=X.

Obten la distribucion marginal de UU.