33 Vectores aleatorios discretos

Definition 33.1.

Se dice que un vector aleatorio X=(X1,,Xn)X=(X_{1},\ldots,X_{n}) es discreto si para cada componente i=1,2,,ni=1,2,\ldots,n, XiX_{i} es una variable aleatoria discreta.

Observemos que, en tal caso, existiran conjuntos A1,,AnA_{1},\ldots,A_{n}\subset\mathbb{R} numerables tales que

P(X1A1)=1,P(X2A2)=1,,P(XnAn)=1P(X_{1}\in A_{1})=1,P(X_{2}\in A_{2})=1,\ldots,P(X_{n}\in A_{n})=1

Si consideramos el conjunto A=A1××AnA=A_{1}\times\cdots\times A_{n} se tiene que P(XA)=1P(X\in A)=1.

Por tanto, el soporte de un vector aleatorio discreto (que es un conjunto de n\mathbb{R}^{n}) es numerable.

Vamos a analizar las distribuciones multivariantes discretas comenzando por un ejemplo para el caso bivariante.

Example 33.2.

Los agricultores de cierta region pueden utilizar tanto abonos quimicos como organicos.

Las probabilidades para cada numero de abonos quimicos empleados (X1)(X_{1}) combinados con cada numero de abonos organicos (X2)(X_{2}) se recogen en la siguiente tabla:

0 11 22
0 0.10.1 0.40.4 0.10.1
1 0.20.2 0.20.2 0

Esta tabla recoge todas las probabilidades conjuntas relativas al vector (X1,X2)(X_{1},X_{2}). Por ejemplo, vemos que

P(X1=1,X2=0)=0.4,P(X1=0,X2=1)=0.2P(X_{1}=1,X_{2}=0)=0.4,P(X_{1}=0,X_{2}=1)=0.2
Definition 33.3.

La funcion de masa de probabilidad conjunta de un vector aleatorio bivariante, (X1,X2)(X_{1},X_{2}) es la funcion

fX1,X2:2[0,1]f_{X_{1},X_{2}}\colon\mathbb{R}^{2}\to[0,1]

definida como

fX1,X2(k1,k2)=P(X1=k1,X2=k2)f_{X_{1},X_{2}}(k_{1},k_{2})=P(X_{1}=k_{1},X_{2}=k_{2})
Proposition 33.4.

La funcion de masa de un vector aleatorio bivariante verifica las propiedades

  1. 1.

    fX1,X2(k1,k2)0f_{X_{1},X_{2}}(k_{1},k_{2})\geq 0 para todo k1,k2k_{1},k_{2}\in\mathbb{R},

  2. 2.

    k1SX1k2SX2fX1,X2(k1,k2)=1\sum_{k_{1}\in S_{X_{1}}}\sum_{k_{2}\in S_{X_{2}}}f_{X_{1},X_{2}}(k_{1},k_{2})=1

La distribucion de probabilidades de un vector aleatorio bivariante puede resumirse en una tabla de doble entrada que refleje su funcion de masa conjunta.

Definition 33.5.

La distribucion marginal de X1X_{1} es la distribucion de probabilidades de la primera variable del vector aleatorio considerada individualmente.

En el caso discreto esta distribucion marginal puede describirse mediante la funcion de masa marginal:

f1(k1)\displaystyle f_{1}(k_{1})
=fX1(k1)=P(X1=k1)=k2SX2P(X1=k1,X2=k2)\displaystyle{}=f_{X_{1}}(k_{1})=P(X_{1}=k_{1})=\sum{k_{2}\in S_{X_{2}}}P(X_{1% }=k_{1},X_{2}=k_{2})
=k2SX2fX1,X2(k1,k2)\displaystyle{}=\sum_{k_{2}\in S_{X_{2}}}f_{X_{1},X_{2}}(k_{1},k_{2})

Analogamente, la funcion de masa marginal de X2X_{2} viene dada por

f2(k2)=k1SX1fX1,X2(k1,k2)f_{2}(k_{2})=\sum_{k_{1}\in S_{X_{1}}}f_{X_{1},X_{2}}(k_{1},k_{2})

Las funciones de masa de probabilidad marginales se calculan facilmente en los margenes de la tabla sumando las probabiliaddes por filas o por columnas.

Example 33.6.

La distribucion marginal de la cantidad de abonos quimicos empleados en la region del ejemplo anterior es

fX1(0)\displaystyle f_{X_{1}}(0)
=P(X1=0)=P(X1=0,X2=0)+P(X1=0,X2=1)=0.1+0.2=0.3\displaystyle{}=P(X_{1}=0)=P(X_{1}=0,X_{2}=0)+P(X_{1}=0,X_{2}=1)=0.1+0.2=0.3
fX1(1)\displaystyle f_{X_{1}}(1)
=P(X1=1)=0.4+0.2=0.6\displaystyle{}=P(X_{1}=1)=0.4+0.2=0.6
fX1(2)\displaystyle f_{X_{1}}(2)
=P(X1=2)=0.1+0=0.1\displaystyle{}=P(X_{1}=2)=0.1+0=0.1

y la distribucion marginal del numero de abonos organicos utilizados,

fX2(0)\displaystyle f_{X_{2}}(0)
=P(X2=0)=0.1+0.4+0.1=0.6\displaystyle{}=P(X_{2}=0)=0.1+0.4+0.1=0.6
fX2(1)=P(X2=1)=0.2+0.2+0=0.4\displaystyle f_{X_{2}}(1)=P(X_{2}=1)=0.2+0.2+0=0.4

Los calculos resultan mas sencillos sumando a los margenes de la tabla.

Usando estas distribuciones marginales podemos calcular, por ejemplo, la esperanza y la varianza de cada una de estas variables aleatorias:

E(X1)\displaystyle E(X_{1})
=0.8,E(X12)\displaystyle{}=0.8,E(X^{2}_{1})
=1,V(X1)\displaystyle=1,V(X_{1})
=0.36\displaystyle{}=0.36
E(X2)\displaystyle E(X_{2})
=0.4,E(X22)\displaystyle{}=0.4,E(X^{2}_{2})
=0.4,V(X2)\displaystyle=0.4,V(X_{2})
=0.24\displaystyle{}=0.24

Ademas de su distribucion marginal, se pueden considerar otras distribuciones de probabilidad para X1X_{1}: las distribuciones condicionadas por un valor particular de X2X_{2}, ya que, en muchos casos, el hecho de saber que la variable X2X_{2} toma un determinado valor modifica la distribucion de probabilidades sobre X1X_{1}.

Definition 33.7.

La funcion de masa de probabilidad de X1X_{1} condicionada a X2=k2X_{2}=k_{2}, denotada como X1X2=k2X_{1}\mid X_{2}=k_{2}, es

fX1(k1X2=k2)\displaystyle f_{X_{1}}(k_{1}\mid X_{2}=k_{2})
=P(X1=k1X2=k2)=P(X1=k1,X2=k2)P(X2=k2)\displaystyle{}=P(X_{1}=k_{1}\mid X_{2}=k_{2})=\frac{P(X_{1}=k_{1},X_{2}=k_{2}% )}{P(X_{2}=k_{2})}
=fX1,X2(k1,k2)fX2(k2)\displaystyle{}=\frac{f_{X_{1},X_{2}}(k_{1},k_{2})}{f_{X_{2}}(k_{2})}
Proposition 33.8.

La funcion de masa condicionada proporcionada otra distribucion de probabilidades para X1X_{1}, y como tal verfica las propiedades

  1. 1.

    fX1(k1X2=k2)0f_{X_{1}}(k_{1}\mid X_{2}=k_{2})\geq 0 para todo k1SX1,k2SX2k_{1}\in S_{X_{1}},k_{2}\in S_{X_{2}}.

  2. 2.

    k1SX1fX1(k1X2=k2)=1\sum_{k_{1}\in S_{X_{1}}}f_{X_{1}}(k_{1}\mid X_{2}=k_{2})=1.

A partir de esta funcion de probabilidad condicionada podemos calcular, por ejemplo, esperanzas y varianzas condicionadas:

E(X1X2=k2)\displaystyle E(X_{1}\mid X_{2}=k_{2})
=k1SX1k1f1(k1k2)\displaystyle{}=\sum_{k_{1}\in S_{X_{1}}}k_{1}\cdot f_{1}(k_{1}\mid k_{2})
V(X1X2=k2)\displaystyle V(X_{1}\mid X_{2}=k_{2})
=E(X12X2=k2)E2(X1X2=k2)\displaystyle{}=E(X^{2}_{1}\mid X_{2}=k_{2})-E^{2}(X_{1}\mid X_{2}=k_{2})

De manera similar, la funcion de masa de probabilidad condicionada a X1=k1X_{1}=k_{1} viene dada por

fX2(k2X1=k1)=fX1,X2(k1,k2)fX1(k1)f_{X_{2}}(k_{2}\mid X_{1}=k_{1})=\frac{f_{X_{1},X_{2}}(k_{1},k_{2})}{f_{X_{1}}% (k_{1})}
Example 33.9.

Para el caso de los abonos de la region del ejemplo anterior,

  1. 1.

    Determinar la distribucion de la cantidad de abonos quimicos, X1X_{1}, sabiendo que el numero de abonos organicos empleados es X2=0X_{2}=0.

    La funcion de probabilidad de X1X_{1} condicionada por X2=0X_{2}=0 es

    fX1(0X2=0)\displaystyle f_{X_{1}}(0\mid X_{2}=0)
    =P(X1=0X2=0)=fX1,X2(0,0)fX2(0)=0.10.6=16\displaystyle{}=P(X_{1}=0\mid X_{2}=0)=\frac{f_{X_{1},X_{2}}(0,0)}{f_{X_{2}}(0% )}=\frac{0.1}{0.6}=\frac{1}{6}
    fX1(1X2=0)\displaystyle f_{X_{1}}(1\mid X_{2}=0)
    =P(X1=1X2=0)=fX1,X2(1,0)fX2(0)=0.40.6=46\displaystyle{}=P(X_{1}=1\mid X_{2}=0)=\frac{f_{X_{1},X_{2}}(1,0)}{f_{X_{2}}(0% )}=\frac{0.4}{0.6}=\frac{4}{6}
    fX1(2X2=0)\displaystyle f_{X_{1}}(2\mid X_{2}=0)
    =P(X1=2X2=0)=fX1,X2(2,0)fX2(0)=0.10.6=16\displaystyle{}=P(X_{1}=2\mid X_{2}=0)=\frac{f_{X_{1},X_{2}}(2,0)}{f_{X_{2}}(0% )}=\frac{0.1}{0.6}=\frac{1}{6}

    o expresado en forma matricial

    X(012164616)X\equiv\begin{pmatrix}0&1&2\\ \frac{1}{6}&\frac{4}{6}&\frac{1}{6}\\ \end{pmatrix}
  2. 2.

    Calcular E(X1X2=0)E(X_{1}\mid X_{2}=0) y V(X1X2=0)V(X_{1}\mid X_{2}=0)

    Como hemos dicho, tambien se modifican la esperanza y la varianza de X1X_{1} (con respecto a los valores marginales)

    E(X1X2=0)\displaystyle E(X_{1}\mid X_{2}=0)
    =016+46+216=1\displaystyle{}=0\cdot\frac{1}{6}+\frac{4}{6}+2\cdot\frac{1}{6}=1
    E(X12X2=0)\displaystyle E(X^{2}_{1}\mid X_{2}=0)
    =016+146+416=43\displaystyle{}=0\cdot\frac{1}{6}+1\cdot\frac{4}{6}+4\cdot\frac{1}{6}=\frac{4}% {3}
  3. 3.

    Establecer la distribucion de la cantidad de abonos organicos, X2X_{2}, cuando se sabe que el numero de abonos quimicos utilizados es X1=1X_{1}=1.

    Del mismo modo, la distribucion de probabilidades del numero de abonos organicos tambien varia cuando sabemos cuantos abonos quimicos se utilizan.

    Asi, la funcion de masa de X2X_{2} condicionada por X1=1X_{1}=1 es

    fX2(0X1=1)\displaystyle f_{X_{2}}(0\mid X_{1}=1)
    =P(X2=0X1=1)=fX1,X2(1,0)fX1(1)=0.40.6=23\displaystyle{}=P(X_{2}=0\mid X_{1}=1)=\frac{f_{X_{1},X_{2}}(1,0)}{f_{X_{1}}(1% )}=\frac{0.4}{0.6}=\frac{2}{3}
    fX2(1X1=1)\displaystyle f_{X_{2}}(1\mid X_{1}=1)
    =P(X2=1X1=1)=fX1,X2(1,1)fX1(1)=0.20.6=13\displaystyle{}=P(X_{2}=1\mid X_{1}=1)=\frac{f_{X_{1},X_{2}}(1,1)}{f_{X_{1}}(1% )}=\frac{0.2}{0.6}=\frac{1}{3}

    o, resumidamente,

    X(012313)X\equiv\begin{pmatrix}0&1\\ \frac{2}{3}&\frac{1}{3}\\ \end{pmatrix}
  4. 4.

    Calcular E(X2X1=1)E(X_{2}\mid X_{1}=1) y V(X2X1=1)V(X_{2}\mid X_{1}=1).

    Del apartado anterior se deduce que

    E(X2X1=1)\displaystyle E(X_{2}\mid X_{1}=1)
    =023+113=23\displaystyle{}=0\cdot\frac{2}{3}+1\cdot\frac{1}{3}=\frac{2}{3}
    E(X22X1=1)\displaystyle E(X^{2}_{2}\mid X_{1}=1)
    =023+113=13\displaystyle{}=0\cdot\frac{2}{3}+1\cdot\frac{1}{3}=\frac{1}{3}
    V(X22X1=1)\displaystyle V(X^{2}_{2}\mid X_{1}=1)
    =13(13)2=29\displaystyle{}=\frac{1}{3}-\left(\frac{1}{3}\right)^{2}=\frac{2}{9}
Example 33.10.

Para el experimento aleatorio consistente en lanzar al aire una moneda equilibrada tres veces, consideremos las variables aleatorias:

X\displaystyle X
numero de caras\displaystyle{}\equiv\text{numero de caras}
Y\displaystyle Y
numero de cruces que preceden a la primera cara\displaystyle{}\equiv\text{numero de cruces que preceden a la primera cara}
  1. 1.

    Determinar la distribucion conjunta del vector aleatorio (X,Y)(X,Y).

  2. 2.

    Calcular E(X)E(X), E(Y)E(Y), V(X)V(X) y V(Y)V(Y).

  3. 3.

    Calcular E(XY=1)E(X\mid Y=1) y V(XY=1)V(X\mid Y=1).

  4. 4.

    Calcular E(YX=2)E(Y\mid X=2) y V(YX=2)V(Y\mid X=2)

  1. 1.
    X Y
    CCC 3 0
    CCY 2 0
    CXC 2 0
    XCC 2 1
    CXX 1 0
    XCX 1 1
    XXC 1 2
    XXX 0 3
    0 1 2 3
    0 0 18\frac{1}{8} 28\frac{2}{8} 18\frac{1}{8} 48\frac{4}{8}
    1 0 18\frac{1}{8} 18\frac{1}{8} 0 28\frac{2}{8}
    2 0 18\frac{1}{8} 0 0 18\frac{1}{8}
    3 18\frac{1}{8} 0 0 0 18\frac{1}{8}
    18\frac{1}{8} 38\frac{3}{8} 38\frac{3}{8} 18\frac{1}{8} 11
    Table 1: caption
  2. 2.

    E(X)=1.5E(X)=1.5, E(Y)=0.875E(Y)=0.875, V(X)=0.75V(X)=0.75, V(Y)=1.1V(Y)=1.1.

  3. 3.

    E(XY=1)=112+212=1.5E(X\mid Y=1)=1\cdot\frac{1}{2}+2\cdot\frac{1}{2}=1.5, E(X2Y=1)=12+2=2.5E(X^{2}\mid Y=1)=\frac{1}{2}+2=2.5, V(XY=1)=2.51.52V(X\mid Y=1)=2.5-1.5^{2}

  4. 4.

    YX=2(0123231300)=(012313)Y\mid X=2\equiv\begin{pmatrix}0&1&2&3\\ \frac{2}{3}&\frac{1}{3}&0&0\\ \end{pmatrix}=\begin{pmatrix}0&1\\ \frac{2}{3}&\frac{1}{3}\\ \end{pmatrix}. E(YX=2)=130.3333E(Y\mid X=2)=\frac{1}{3}\simeq 0.3333, V(YX=2)=290.2222V(Y\mid X=2)=\frac{2}{9}\simeq 0.2222