32 Distribucion de transformaciones de variables aleatorias

Dada una variable aleatoria XX y una funcion g:g\colon\mathbb{R}\to\mathbb{R} en ocasiones puede interesarnos saber como se distribuye g(X)g(X). Observemos que g(X):Ωg(X)\colon\Omega\to\mathbb{R} es una nueva variable aleatoria. Si llamamos Y=g(X)Y=g(X) nos interesa saber como se distribuye la variable aleatoria YY.

Si la variable aleatoria X:ΩX\colon\Omega\to\mathbb{R} es discreta, entonces para cualquier funcion gg la variable aleatoria Y=g(X)Y=g(X) tambien sera discreta.

Su funcion de masa de probabilidad es facil de calcular:

fY(k)=P(Y=k)=P[g(X)=k]=P[X=g1(k)]f_{Y}(k)=P(Y=k)=P[g(X)=k]=P[X=g^{-1}(k)]

aunque al realizar los calculos debe ponerse especial atencion en los casos en los que la funcion gg no sea inyectiva.

Example 32.1.

Dada la variable aleatoria HH con funcion de distribucion

FH(t)={0 si t<3,0.2 si 3t1,0.5 si 1t<1,0.6 si 1t<2,0.9 si 2t<3,1 si t3F_{H}(t)=\begin{cases}0\quad&\text{ si }t<-3,\\ 0.2&\text{ si }-3\leq t\leq-1,\\ 0.5&\text{ si }-1\leq t<1,\\ 0.6&\text{ si }1\leq t<2,\\ 0.9&\text{ si }2\leq t<3,\\ 1&\text{ si }t\geq 3\end{cases}
  1. 1.

    Determinar la funcion de masa de probabilidad de J=H4+5J=H^{4}+5

    Escrita abreviadamente en forma matricial, la distribucion de HH es

    H(311230.20.30.10.30.1)H\equiv\begin{pmatrix}-3&-1&1&2&3\\ 0.2&0.3&0.1&0.3&0.1\\ \end{pmatrix}

    y por tanto la distribucion de J=H4+5J=H^{4}+5 es

    J(621860.40.30.3)J\equiv\begin{pmatrix}6&21&86\\ 0.4&0.3&0.3\\ \end{pmatrix}
  2. 2.

    Calcular P(J>20)P(J>20) y E(J)E(J).

    P(J>20)\displaystyle P(J>20)
    =P(J=21)+P(J=86)=0.3+0.3=0.6\displaystyle{}=P(J=21)+P(J=86)=0.3+0.3=0.6
    E(J)\displaystyle E(J)
    =0.46+0.321+0.386=34.5\displaystyle{}=0.4\cdot 6+0.3\cdot 21+0.3\cdot 86=34.5

Cuando la variable aleatoria X:ΩX\colon\Omega\to\mathbb{R} es continua, entonces dependiendo de como sea gg la variable aleatoria Y=g(X)Y=g(X) puede ser discreta o continua. Ademas, si la transformacion gg presenta discontinuidades YY tambien tendra discontinuidades. Si la transformacion gg es continua YY sera continua.

Vamos a centrarnos en los casos en los que la variable aleatoria XX es continua y gg es una funcion continua. En estos casos, la forma mas inmediata de determinar la distribucion de probabilidad de Y=g(X)Y=g(X) es calcular su funcion de distribucion

FY(t)=P(Yt)=P[g(X)t]F_{Y}(t)=P(Y\leq t)=P[g(X)\leq t]
Example 32.2.

Sea MM una variable aleatoria con distribucion uniforme en el intervalo (0,1)(0,1), MU(0,1)M\sim U(0,1).

  1. 1.

    Determinar la funcion de densidad de la variable aleatoria X=eMX=e^{M}.

    Tenemos que encontrar la funcion de densidad de XX. Para ello, primero determinamos su funcion de distribucion

    FX(t)=P(Xt)=P(eMt)F_{X}(t)=P(X\leq t)=P(e^{M}\leq t)

    El primer paso es determinar cual es el soporte de XX. Para ello lo mas sencillo es representar la funcion

    g(m)=emg(m)=e^{m}

    a lo largo del soporte de MM, que en este caso es el intervalo [0,1][0,1]. El soporte de XX es SX=[1,e]S_{X}=[1,e].

    Observemos ahora que, dado que MU(0,1)M\sim U(0,1), la funcion de distribucion de MM, FMF_{M}, es

    FM(z)=P[Mz]={0 si z<0z si 0z11 si z>1F_{M}(z)=P[M\leq z]=\begin{cases}0\quad&\text{ si }z<0\\ z&\text{ si }0\leq z\leq 1\\ 1&\text{ si }z>1\end{cases}

    Nuestro objetivo es calcular FX(t)=P(Xt)=P(eMt)F_{X}(t)=P(X\leq t)=P(e^{M}\leq t). En primer lugar podemos observar que, puesto que SX=[1,e]S_{X}=[1,e], FX(t)=0t<1F_{X}(t)=0\;\forall t<1 y FX(t)=1t>eF_{X}(t)=1\;\forall t>e.

    Por otra parte, para t[1,e]t\in[1,e], FX(t)=P(Xt)=P(eMt)=P(Mlnt)=ln(t)F_{X}(t)=P(X\leq t)=P(e^{M}\leq t)=P(M\leq\ln t)=\ln(t). Luego la funcion de distribucion de XX es

    FX(t)={0 si t<1ln(t) si 1t<e1 si teF_{X}(t)=\begin{cases}0\quad&\text{ si }t<1\\ \ln(t)&\text{ si }1\leq t<e\\ 1&\text{ si }t\geq e\end{cases}

    y por tanto su funcion de densidad es

    fX(t)={1t si 1te0 si t[1,e]f_{X}(t)=\begin{cases}\frac{1}{t}\quad&\text{ si }1\leq t\leq e\\ 0\text{ si }t\notin[1,e]\end{cases}
  2. 2.

    Calcular P(1.5<X<2)P(1.5<X<2).

    P(1.5<X<2)=FX(2)FX(1.5)=ln(2)ln(1.5)=0.2877P(1.5<X<2)=F_{X}(2)-F_{X}(1.5)=\ln(2)-\ln(1.5)=0.2877
Example 32.3.

Sea GG una variable aleatoria con funcion de densidad

fG(x)={2x si x(0,1)0 si x(0,1)f_{G}(x)=\begin{cases}2x\quad&\text{ si }x\in(0,1)\\ 0&\text{ si }x\notin(0,1)\end{cases}

Determinar la funcion de densidad de las siguientes variables aleatorias:

  1. 1.

    A=3G+1A=3G+1.

  2. 2.

    B=eGB=e^{-G}

  3. 3.

    C=G2C=G^{2}

La funcion de distribucion de GG es

FG(x)=P(Gx)={0 si x0x2 si x(0,1)1 si x1F_{G}(x)=P(G\leq x)=\begin{cases}0\quad&\text{ si }x\leq 0\\ x^{2}&\text{ si }x\in(0,1)\\ 1&\text{ si }x\geq 1\end{cases}
  1. 1.

    Queremos determinar la funcion de densidad de A=3G+1A=3G+1. Para ello comenzaremos por encontrar su funcion de distribucion. El soporte de AA sera el intervalo (1,4)(1,4).

    Por tanto sabemos que FA(t)=0F_{A}(t)=0 para todo t1t\leq 1 y FA(t)=1F_{A}(t)=1 para todo t4t\geq 4. Para t(1,4)t\in(1,4) tenemos que

    FA(t)=P(At)=P(3G+1t)=P(Gt13)=(t13)2F_{A}(t)=P(A\leq t)=P(3G+1\leq t)=P\left(G\leq\frac{t-1}{3}\right)=\left(\frac% {t-1}{3}\right)^{2}

    Juntando los tres tramos tenemos que la expresion completa de la funcion de distribucion de BB es

    FA(t)={0 si t1(t13)2 si 1<t<41 si t4F_{A}(t)=\begin{cases}0\quad&\text{ si }t\leq 1\\ \left(\frac{t-1}{3}\right)^{2}&\text{ si }1<t<4\\ 1&\text{ si }t\geq 4\end{cases}

    y derivando FA(t)F_{A}(t) obtenemos que la funcion de densidad de AA es

    fA(t)={2(t1)9 si t(1,4)0 si t(1,4)f_{A}(t)=\begin{cases}\frac{2(t-1)}{9}\quad&\text{ si }t\in(1,4)\\ 0&\text{ si }t\notin(1,4)\end{cases}
  2. 2.

    Ahora queremos determinar la funcion de densidad de B=eGB=e^{-G}. El soporte de BB sera el intervalo (e1,1)(e^{-1},1). Por tanto sabemos que FB(t)=0F_{B}(t)=0 para todo te1t\leq e^{-1} y FB(t)=1F_{B}(t)=1 para todo t1t\geq 1.

    Para t(e1,1)t\in(e^{-1},1) tenemos que

    FB(t)\displaystyle F_{B}(t)
    =P(Bt)=P(eGt)=P(Gln(t))=P(Gln(t))\displaystyle{}=P(B\leq t)=P(e^{-G}\leq t)=P(-G\leq\ln(t))=P(G\geq-\ln(t))
    =1FG(ln(t))=1ln(t)2\displaystyle{}=1-F_{G}(-\ln(t))=1-\ln(t)^{2}

    Juntando los tres tramos se tiene la funcion de distribucion de BB. Y derivando FB(t)F_{B}(t) obtenemos que la funcion de densidad de BB es

    fB(t)={2ln(t)t si t(e1,1)0 si t(e1,1)f_{B}(t)=\begin{cases}\frac{-2\ln(t)}{t}\quad&\text{ si }t\in(e^{-1},1)\\ 0&\text{ si }t\notin(e^{-1},1)\end{cases}
  3. 3.

    Por ultimo, tenemos que encontrar la distribucion de C=G2C=G^{2}. Como se puede comprobar, el soporte de C=X2C=X^{2} es tambien el intervalo (0,1)(0,1).

    Para t(0,1)t\in(0,1) se tiene

    FC(t)=P(Ct)=P(G2t)=P(G3)=FG(t)=t2=tF_{C}(t)=P(C\leq t)=P(G^{2}\leq t)=P(G\leq\sqrt{3})=F_{G}(\sqrt{t})=\sqrt{t}^{% 2}=t

    La expresion completa de la funcion de distribucion de CC es por tanto

    FC(t)={0 si t0t si 0<t<11 si t1F_{C}(t)=\begin{cases}0\quad&\text{ si }t\leq 0\\ t&\text{ si }0<t<1\\ 1\text{ si }t\geq 1\end{cases}

    y derivando FB(t)F_{B}(t) obtenemos que la funcion de densidad de BB es

    fB(t)={1 si t(0,1)0 si t(0,1)f_{B}(t)=\begin{cases}1\quad&\text{ si }t\in(0,1)\\ 0&\text{ si }t\notin(0,1)\end{cases}

    Notese que CC sigue una distribucion uniforme en el intervalo (0,1)(0,1), CU(0,1)C\sim U(0,1).

Example 32.4.

Sea KK una variable aleatoria con funcion de densidad

fK(x)={x218 si x(3,3)0 si x(3,3)f_{K}(x)=\begin{cases}\frac{x^{2}}{18}\quad&\text{ si }x\in(-3,3)\\ 0&\text{ si }x\notin(-3,3)\end{cases}

Determinar la funcion de densidad de la variable aleatoria Y=|K|Y=\left|K\right|.

Integrando fK(x)f_{K}(x) se obtiene que la funcion de distribucion de KK es

FK(t)=P(Kt)={0 si t<3t3+2754 si 3t<31 si t3F_{K}(t)=P(K\leq t)=\begin{cases}0\quad&\text{ si }t<-3\\ \frac{t^{3}+27}{54}&\text{ si }-3\leq t<3\\ 1&\text{ si }t\geq 3\end{cases}

En este caso debemos encontrar la funcion de densidad de Y=|K|Y=\left|K\right|.

Para determinar el soporte de YY y la relacion entre KK e YY representamos la funcion g(x)=|x|g(x)=\left|x\right| a lo largo del soporte de KK, que es el intervalo (3,3)(-3,3).

El grafico de la pagina siguiente se aprecia claramente que el soporte de YY va a ser el intervalo SY=[0,3)S_{Y}=[0,3).

Por tanto sabemos que

FY(t)=0 para todo t<0F_{Y}(t)=0\text{ para todo }t<0

y que

FY(t)=1 para todo t3F_{Y}(t)=1\text{ para todo }t\geq 3

Tambien se observa que la transformacion no es monotona. Para t[0,3)t\in[0,3) tenemos que

FY(t)\displaystyle F_{Y}(t)
=P(Tt)=P(|K|t)=P(tKt)=FK(t)FK(t)\displaystyle{}=P(T\leq t)=P(\left|K\right|\leq t)=P(-t\leq K\leq t)=F_{K}(t)-% F_{K}(-t)
=t3+2754(t)3+2754=2t354=t327\displaystyle{}=\frac{t^{3}+27}{54}-\frac{(-t)^{3}+27}{54}=\frac{2t^{3}}{54}=% \frac{t^{3}}{27}

Juntando los tres tramos se obtiene la funcion de distribucion de YY. Derivando FY(t)F_{Y}(t) obtenemos la funcion de densidad de YY:

fY(t)={t29 si t[0,3)0 si t[0,3)f_{Y}(t)=\begin{cases}\frac{t^{2}}{9}\quad&\text{ si }t\in[0,3)\\ 0&\text{ si }t\notin[0,3)\end{cases}