Archivo de la etiqueta: variables aleatorias continuas

Probabilidad I: Transformaciones de V.A.’s Continuas

Por Octavio Daniel Ríos García

Introducción

En la entrada pasada comenzamos a abordar el problema de encontrar la distribución de la transformación de una v.a. conocida. En particular, analizamos un método para el caso en el que la v.a. conocida es discreta, y sólamente para ese caso. No obstante, este método nos servirá como «base» para dar paso al caso continuo. Primero presentaremos un método que consiste en manipular directamente la función de distribución; muy parecido al método que vimos en la entrada pasada. Después, analizaremos un método más especializado que permite encontrar la función de densidad de la transformación de una v.a. sin necesidad de manipular la función de distribución.

Motivación del primer método

Sea X una v.a. y g:RR una función Borel-medible. En la entrada pasada ya describimos el proceso para obtener los eventos de g(X) en términos de eventos que involucran a X. De hecho, vimos que para cada AB(R) se cumple que

(g(X)A)=(Xg1[A]).

¡Atención! En la entrada pasada centramos nuestra atención en las v.a.’s discretas, pero la igualdad anterior es cierta para cualquier variable aleatoria. Por ello, también aplica para las v.a.’s continuas. En particular, para cada yR se cumple que (,y]B(R), por lo que

(g(X)y)=(g(X)(,y])=(Xg1[(,y]]).

Por lo tanto, se tiene que

P(g(X)y)=P(Xg1[(,y]]).

Es decir, si definimos a Y=g(X) y FY:RR es la función de distribución de Y, entonces lo anterior quiere decir que para cada yR,

FY(y)=P(Xg1[(,y]]);

por lo que es posible obtener la distribución de Y en términos de la probabilidad de un evento que involura a X, cuya distribución sí conocemos.

Primer método: manipular la función de distribución

Con la discusión anterior llegamos a que si X es una v.a. (cuya distribución es conocida), g:RR es una función Borel-medible, y Y es la v.a. definida como Y=g(X), entonces la función de distribución de Y, FY:RR, puede obtenerse como

FY(y)=P(Xg1[(,y]])para cada yR.

Por ello, el problema consistirá en encontrar el conjunto g1[(,y]], y así encontrar la probabilidad de (Xg1[(,y]]).

Ejemplo 1. Sea X una v.a. con función de densidad fX:RR dada por

fX(x)=12e|x|,para cada xR.

Una v.a. con esta función de densidad es conocida como una v.a. con distribución Laplace, o distribución doble exponencial. Su función de distribución FX:RR está dada por

FX(x)={12exsi x<0,112exsi x0.

Sea g:RR la función dada por g(x)=|x| para cada xR. De este modo, defínase Y=g(X)=|X|. Para obtener la función de distribución Y, podemos seguir un método similar al que usamos en la entrada anterior. Sea yR. Un primer detalle que podemos observar sobre Y es que no toma valores negativos. Por ello, si y<0, se tiene que (Yy)=, y en consecuencia, P(Yy)=0 para y<0. Por otro lado, para y0 se tiene que que

ω(Yy)Y(ω)y|X(ω)|yyX(ω)yω(yXy),

por lo que para cada y0 se tiene que (Yy)=(yXy). Como esos dos eventos son iguales, se sigue que P(Yy)=P(yXy). , por lo queAdemás, nota que

P(yXy)=P(Xy)P(X<y)()=P(Xy)P(Xy)=FX(y)FX(y),

donde el paso () es válido debido a que X es una v.a. continua. Por ello, podemos concluir que para cada y0,

FY(y)=FX(y)FX(y).

Por lo tanto, la función de distribución de Y queda como sigue:

FY(y)={0si y<0,FX(y)FX(y)si y0.

De aquí podemos obtener una expresión explícita. Para cada y0 se tiene que y0, así que

FX(y)FX(y)=(112ey)12ey=112ey12ey=1ey.

En conclusión, la función de distribución de Y queda así:

FY(y)={0si y<0,1eysi y0.

Probablemente te resulte familiar: ¡Es la función de distribución de una v.a. exponencial! Este ejemplo exhibe que algunas transformaciones de algunas v.a.’s «famosas» resultan en otras v.a.’s «famosas». En este caso, vimos que si X es una v.a. que sigue una distribución Laplace, entonces |X| sigue una distribución exponencial. Más adelante veremos muchas más distribuciones importantes, y veremos cómo se relacionan entre sí mediante transformaciones.


Segundo método: teorema de cambio de variable

Existe un método más especializado para obtener la función de densidad de la transformación de una v.a. continua. La razón por la que decimos que es más especializado es porque funciona para transformaciones que cumplen ciertas condiciones.


Teorema. Sea X:ΩR una v.a. continua con función de densidad fX:RR, y sea g:X[Ω]R una función diferenciable y estrictamente creciente o decreciente. Entonces la función de densidad de Y=g(X) está dada por

fY(y)={fX(g1(y))|ddy[g1(y)]|si y(gX)[Ω],0en otro caso,

donde g1:g[R]R es la inversa de g, y (gX)[Ω] es la imagen directa de Ω bajo gX. Esto es, (gX)[Ω]={yRωΩ:(gX)(ω)=y}, que corresponde al conjunto de valores que toma la v.a. Y=g(X).


Demostración. Demostraremos el caso en el que g es estrictamente creciente. Para ello, sea yR. Primero, recuerda que

(Yy)=(Xg1[(,y]]).

Por un lado, se tiene el caso en el que y(gX)[Ω]; es decir, y es uno de los valores que toma la v.a. Y (pues (gX)[Ω]=Y[Ω]). En este caso, el valor g1(y) está bien definido, ya que g1:(gX)[Ω]R es una función cuyo dominio es la imagen de g. De este modo, para cada ωΩ tendremos que

()Y(ω)yX(ω)g1(y).

Como g es una función estrictamente creciente, su inversa g1:(gX)[Ω]R también es estrictamente creciente, y por lo tanto, la desigualdad en () «no se voltea».

De lo anterior se sigue que P(Yy)=P(Xg1(y)) para cada y(gX)[Ω]. En consecuencia, se tiene que

FY(y)=FX(g1(y)).

Podemos diferenciar ambos lados de la igualdad respecto a y, y por la regla de la cadena obtenemos

fY(y)=fX(g1(y))ddy[g1(y)]=fX(g1(y))|ddy[g1(y)]|,

donde el último paso se obtiene de que g1 es estrictamente creciente, y por lo tanto, su derivada es positiva.

Por otro lado, resta el caso en el que y(gX)[Ω]; es decir, cuando y no es uno de los valores que puede tomar Y. En este caso, simplemente fY vale 0, pues la densidad de una v.a. continua es 0 en aquellos valores que no toma. De este modo, Y tiene densidad fY:RR dada por

fY(y)={fX(g1(y))|ddy[g1(y)]|si y(gX)[Ω],0en otro caso,

que es justamente lo que queríamos demostrar.

El caso para g estrictamente decreciente es casi análogo, por lo que te lo dejamos de tarea moral.

◻

Es importante notar que el teorema anterior no funciona para cualquier g:RR Borel-medible, sólamente para aquellas que cumplen las hipótesis del teorema. Bajo estas hipótesis, el teorema permite obtener la densidad de la transformación de una v.a. de manera más eficiente que los otros métodos que hemos abordado.

Ejemplo 2. Sea Z una v.a. con densidad fZ:RR dada por

fZ(z)=12πez2/2,para cada zR.

Se dice que una v.a. con esa función de densidad sigue una distribución normal estándar. Observa que fZ(z)>0 para todo zR, por lo que Z[Ω]=R. Es decir, Z puede tomar cualquier valor en R.

Sea W=eZ. La función exp:RR+ dada por exp(x)=ex es estrictamente creciente y diferenciable, por lo que podemos usar el teorema anterior para obtener la función de densidad de W. Así, tenemos que

fW(w)=fZ(exp1(w))|ddw[exp1(w)]|,

donde exp1:R+R es la inversa de la función exponencial exp. De hecho, la inversa de exp es la función ln:R+R, el logaritmo natural. Ahora, como Z[Ω]=R, se tiene que (expZ)[Ω]=R+, pues la función exp toma únicamente valores positivos.

En consecuencia, para wR+ se tiene

fW(w)=fZ(ln(w))|ddw[ln(w)]|=fZ(ln(w))|1w|=1w2πexp((ln(w))22)

y así, tenemos que W tiene densidad fW:RR dada por

fW(w)={1w2πexp((ln(w))22)si w>0,0en otro caso.

Como nota adicional, fW es la densidad de una v.a. cuya distribución es conocida como log-normal.


Para concluir, es importante mencionar que la transformación g del teorema sólamente necesita ser diferenciable y estrictamente creciente sobre X[Ω] (por eso es que en el enunciado la pusimos como g:X[Ω]R). Por ejemplo, la función g:R+{0}R dada por g(x)=x2 es una función creciente sobre su dominio. Por ello, si X es una v.a. continua que toma únicamente valores no-negativos, entonces puede aplicarse el teorema para obtener la densidad de g(X). En resumidas cuentas, el teorema puede aplicarse siempre y cuando la transformación g sea diferenciable y estrictamente creciente sobre el conjunto de valores que puede tomar X.

Tarea moral

Los siguientes ejercicios son opcionales. Es decir, no formarán parte de tu calificación. Sin embargo, te recomiendo resolverlos para que desarrolles tu dominio de los conceptos abordados en esta entrada.

  1. Verifica que la función de distribución de la v.a. X del Ejemplo 1 es la función que te dimos. Es decir, obtén la función de distribución de X a partir de su función de densidad.
  2. Demuestra el caso en el que g es estrictamente decreciente del teorema de cambio de variable.
  3. Retoma el segundo Ejemplo 2, pero esta vez comienza con la v.a. W, cuya función de densidad es fW(w)={1w2πexp((ln(w))22)si w>0,0en otro caso,y encuentra la función de densidad de Z=ln(W) usando el teorema.

Más adelante…

El teorema de esta entrada es muy útil para obtener la densidad (y, en consecuencia, la distribución) de muchas transformaciones de v.a.’s continuas. Por ello, nos será de utilidad en el futuro relativamente cercano, cuando veamos las distribuciones de probabilidad más conocidas. Por otro lado, te será de utilidad mucho más adelante en materias posteriores, pues este teorema puede generalizarse al caso en el que la transformación tiene como dominio a Rn y como codominio a R con nN+ y n2 (por ejemplo, g:R2R dada por g(x,y)=x+y).

En la siguiente entrada comenzaremos el estudio de un concepto asociado a las v.a.’s llamado el valor esperado de una variable aleatoria.

Entradas relacionadas

Probabilidad I: Variables Aleatorias Continuas

Por Octavio Daniel Ríos García

Introducción

En la entrada anterior comenzamos el estudio de los dos tipos de v.a.’s que nos interesan, y vimos el caso de las v.a.’s discretas. Sin embargo, puede que te preguntes exactamente por qué es necesaria esta distinción. Por ello, en esta entrada presentaremos las propiedades de las v.a.’s continuas para que compares ambos tipos, y puedas apreciar sus diferencias.

A diferencia de las v.a.’s discretas, las v.a.’s continuas pueden tomar una cantidad infinita no numerable de valores distintos. Es decir, el conjunto de valores que puede tomar una v.a. continua tiene la misma cardinalidad que los números reales.

Definición de una v.a. continua

En la introducción comentamos que las v.a.’s continuas, a grandes rasgos, son aquellas cuyo conjunto de valores posibles es infinito no numerable. Sin embargo, no las definiremos a partir de este hecho, sino a partir de una propiedad que debe cumplir su función de distribución.


Definición. Una variable aleatoria X es absolutamente continua si y sólamente si existe una función f:RR integrable y no-negativa tal que para cada xR se cumple que

FX(x)=xf(t)dt.

Llamamos a f la función de densidad (o simplemente la densidad) de X.


A este tipo de v.a.’s se les llama absolutamente continuas debido a que su función de distribución es una función absolutamente continua. Es común encontrarse con el término «continua» en vez de «absolutamente continua«, aunque en contextos más formales existen diferencias entre ambos términos. De manera equivalente, una v.a. X es continua si existe una función f:RR integrable y no-negativa tal que para cada xR se cumple que

P(Xx)=xf(t)dt,

por la definición de FX. Es decir, X es una v.a. continua si la función de distribución de X puede escribirse como la integral de una función integrable y no-negativa f, que es llamada la densidad de X. Es decir, las probabilidades de los eventos que involucran a una v.a. continua se pueden expresar en términos del área debajo de su función de densidad.

En consecuencia, si X es una v.a. continua con densidad f:RR, como FX(x)1 cuando x, se tiene que

f(t)dt=1.

Además, para cada a, bR tales que a<b se tiene que P(X(a,b])=FX(b)FX(a), por lo que

P(X(a,b])=bf(t)dtaf(t)dt=abf(t)dt.

Es importante notar que en la definición pedimos que la densidad de una v.a. continua X debe ser una función integrable, lo cual garantiza que la función de distribución de X es continua.

Ejemplo 1. Sea λR tal que λ>0 y sea X una v.a. con distribución F:RR dada por

F(x)={1eλxsi x0,0en otro caso.

Gráficamente, para algunos valores de λ, F se ve como sigue:

Figura. Gráfica de la función F(x)=1eλx para distintos valores de λ.

Si tomamos a f:RR como sigue

f(x)={λeλxsi x0,0en otro caso,

se cumplirá que f es la densidad de X. Para verlo, tomamos xR. Primero, observa que si x<0, entonces

xf(t)dt=x0dt=0=F(x),

donde F(x)=0 se cumple por la definición de F. Por otro lado, para el caso en el que x0 se tiene que

xf(t)dt=00dt+0xλeλtdt=0+(eλt)|0x=eλ0eλx=1eλx,

por lo que queda demostrado que para cada xR se cumple que F(x)=xf(t)dt. Por lo tanto, f es la densidad de X.

Figura. Gráfica de la función de densidad de X, f(x)=λeλx, para algunos valores de λ.

Cuando una v.a. tiene la función de distribución (y, en consecuencia, la densidad) del ejemplo anterior, se dice que sigue una distribución exponencial. Esta es una de las muchas distribuciones importantes que veremos más adelante.

¿Las v.a.’s continuas tienen función de masa de probabilidad?

Como seguramente ya notaste, la relación que existe entre la densidad y la función de distribución de una v.a. continua se parece mucho a la relación entre la masa de probabilidad y la distribución de una v.a. discreta. En el caso de las discretas, para obtener el valor de F(x) para xR a partir de la función de masa de probabilidad, lo que se hace es sumar todas las probabilidades de los valores menores o iguales a x. Por otro lado, en el caso de las continuas lo que se hace es integrar la función de densidad desde hasta x.

Sin embargo, hay un detalle muy importante en el que difieren las v.a.’s continuas de las discretas, que desarrollaremos a continuación. Sea X una v.a. continua. ¿Recuerdas la siguiente propiedad? La vimos en la entrada pasada. Para cada aR, se cumple que

P(X=a)=FX(a)FX(a),

donde no olvides que FX(a) es el límite de F(x) cuando x tiende a a por la izquierda. Sin embargo, en el caso de una v.a. continua, sabemos que existe f:RR integrable y no-negativa tal que para cada xR, se cumple que

FX(x)=xf(t)dt.

Debido a que f es una función integrable, la función FX es continua. En particular, es continua por la izquierda, por lo que para cada aR se cumple que

P(X=a)=FX(a)FX(a)=FX(a)FX(a)=0.

Esto seguramente te resulta confuso, ¿la probabilidad de que la v.a. X tome cualquier valor real es 0? ¡Así es! Sin embargo, nosotros habíamos dicho que en la probabilidad, medíamos qué tan «probable» es que pase un evento con una calificación del 0 al 1. En particular, habíamos acordado que 0 representa lo más improbable posible. ¿Esto significa que es imposible que una v.a. continua tome algún valor fijo? ¡No! Ten mucho cuidado, nosotros dijimos que cuando un evento tiene probabilidad 0 esto significa que es lo más improbable posible de acuerdo con la medida de probabilidad que se está utilizando. Esto puede interpretarse como que sí es imposible (como pasa con los puntos de probabilidad 0 en una v.a. discreta, o los puntos donde la densidad de una v.a. continua vale 0), o puede significar que es muy improbable, pero no imposible.

Además, debido a esto último, cuando X es una v.a. continua se cumple que

P(X(a,b])=P(X(a,b))=P(X[a,b))=P(X[a,b]),

y que

P(Xa)=P(X<a),P(Xa)=P(X>a),

por lo que con las v.a.’s continuas no es necesario preguntarse si la desigualdad es estricta o no. Mucho cuidado, con las discretas sí debes de tener cuidado con eso, porque en las discretas hay valores aR para los cuales P(X=a)>0.

Pese a que la probabilidad de los eventos (X=a) es 0 para cada aR cuando X es una v.a. continua, la variable aleatoria sí puede tomar cualquiera de los valores en los que su función de densidad es mayor a 0.

Partiendo de una función de densidad

En la entrada pasada vimos que puede definirse la distribución de una v.a. discreta a partir de una función de masa de probabilidad. De manera muy similar, puede definirse la distribución de una v.a. continua a partir de una función de densidad. Sin embargo, hay que establecer las propiedades que debe de satisfacer una función para poder uitlizarla como función de densidad.

Sea f:RR una función integrable. Si se cumple que

(1)f(x)0para cada xR,(2)f(t)dt=1,

entonces f es la densidad de alguna v.a. continua. Para confirmar este hecho, define F:RR como

F(x)=xf(t)dt,para cada xR,

y demuestra (tarea moral) que F es una función de distribución. Esto es, que F es no-decreciente, continua por la derecha y su límite a es 1 y a es 0.

Ejemplo 2. La función de densidad de una v.a. continua debe de ser integrable. Sin embargo, no tiene por qué ser continua, ya que funciones integrables que tienen discontinuidades. Por ejemplo, considera a f:RR dada por

f(x)={1si x[1.5,1],1si x[1,1.5],0en otro caso.

Gráficamente:

Figura. Gráfica de la función f. Observa que claramente f no es una función continua.

Para ver que f es la función de densidad de alguna v.a. continua, hay que demostrar que f satisface:

f(x)0para cada xRyf(t)dt=1.

Primero, observa que por la definición de f, para cada xR se cumple que f(x)=0 o f(x)=1, por lo que f es no-negativa. Por otro lado, veamos cuánto vale la integral de f sobre R.

f(t)dt=1.5f(t)dt+1.51f(t)dt+11f(t)dt+11.5f(t)dt+1.5f(t)dt=1.50dt+1.511dt+110dt+11.51dt+1.50dt=1.511dt+11.51dt=(t|1.51)+(t|11.5)=(1(1.5))+(1.51)=(1+1.5)+(1.51)=0.5+0.5=1,

por lo que f(t)dt=1, así que f es la función de densidad de alguna v.a. continua.

A partir de f podemos obtener la función de distribución que le corresponde, donde para cada xR, se define F como

F(x)=xf(t)dt.

Sin embargo, como f está definida por pedazos, hay que tener cuidado con cada uno de los casos para x. En este caso, f tiene 5 casos que analizar. Primero, cuando x<1.5, tenemos que

F(x)=xf(t)dt=x0dt=0,

pues f(x)=0 cuando x<1.5. Luego, cuando 1.5x1, tenemos que

F(x)=xf(t)dt=1.50dt+1.5x1dt=0+(t|1.5x)=x(1.5)=x+1.5,

¡observa con cuidado cómo los integrandos difieren debido a que f está definida por pedazos! Continuando con el ejemplo, cuando 1<x<1, se tiene que

F(x)=xf(t)dt=1.50dt+1.511dt+110dt=0+(t|1.51)+0=(1)(1.5)=0.5,

es decir, F permanece constante entre 1 y 1, y toma el valor 0.5. Cuando 1x1.5, la integral queda como sigue

F(x)=xf(t)dt=1.50dt+1.511dt+110dt+1x1dt=0+(t|1.51)+0+(t|1x)=((1)(1.5))+(x1)=0.5+x1=x0.5,

y así obtenemos la función de distribución asociada a f, que es

F(x)={0si x<1.5,x+1.5si 1.5x1,0.5si 1<x<1,x0.5si 1x1.5,1si 1.5x,

que gráficamente se ve como la siguiente figura

Figura. Gráfica de la función de distribución F. Observa que, informalmente, puede apreciarse que satisface las propiedades de una función de distribución (como ya lo esperábamos): es no-decreciente, continua por la derecha y sus límites a e son 0 y 1, respectivamente.

Teniendo la función de distribución, es posible calcular muchas probabilidades. Si U es una v.a. cuya distribución es la función F que obtuvimos, entonces podemos obtener:

P(U43)=43+1.5=43+32=8+96=160.16666,

por lo que P(U43)16.666%. Además, como U es una v.a. continua, se tiene que P(U43)=P(U<43), pues P(U=43)=0. Esto pasa con cualquier valor, como ya mencionamos anteriormente. Por ello, al obtener probabilidades de que una v.a. continua esté dentro de algún intervalo, puedes no preocuparte por los extremos. Otro ejemplo:

P(X(1.2,1.4])=F(1.4)F(1.2)=(1.40.5)((1.2)+1.5)=0.90.3=0.6,

que es igual a P(X(1.2,1.4)), pues P(X=1.4)=0.


Recuperando la densidad a partir de la distribución

Habrá situaciones en las que tendremos la función de distribución de una v.a. continua, y necesitaremos su función de densidad. El siguiente ejemplo exhibe una metodología para obtener la densidad de una v.a. continua a partir de su función de distribución.

Ejemplo 3. Sea Z una v.a. con distribución G:RR dada por

G(z)={0si x<0,z2si 0z<12,13(1z)2si 12z<1,1si z1,

que gráficamente se ve como sigue:

Figura. Gráfica de la función G. Parece que sí es una función de distribución… ¡No se te olvide verificarlo!

Ahora, para obtener la función de densidad de Z, digamos, g. Para hacerlo, hay que hacer el procedimiento inverso al que seguimos en el ejemplo anterior. No obstante, la densidad g será una función definida por pedazos, así como G (y como la densidad del ejemplo anterior).

Además, en lugar de integrar, ahora tenemos que encontrar la derivada G. Es decir, la función g que buscamos debe de satisfacer que

G(z)=g(z),para cada zR.

Sin embargo, observa que la función no es diferenciable sobre todo su dominio. Gráficamente, tiene un piquito en z=12, que es un indicador de que la función no es diferenciable en ese punto. Lo mismo pasa en z=1.

Sin embargo, es posible obtener la derivada de la función G en los pedazos en los que sí es diferenciable. Podemos trabajar de manera similar a como lo hicimos en el ejemplo anterior, analizando cada uno de los casos que contempla la regla de correspondencia de G. De este modo, sea z<0. Queremos encontrar g1, una función no-negativa, tal que

G(z)=zg1(t)dt,

y como z<0, G(z)=0, por lo que buscamos g1 tal que

zg1(t)dt=0,

que debe de ser una función constante para que se cumpla que G(t)=g1(t) para t<0. Además, g1 debe de ser no-negativa, pues buscamos que sea parte de una función de densidad, la cual debe de ser no negativa. De este modo, la única función g1 que satisface este hecho es la función dada por g1(t)=0, pues es constante y su integral de a z es 0, con z<0. De este modo, g, la derivada de G que buscamos, cumple que g(z)=0 para cada z<0.

Después, para z[0,12), tenemos que

G(z)=z2,

y lo que queremos es encontrar una función g2 tal que

G(z)=0g1(t)dt+0zg2(t)dt,

donde g1 es la que obtuvimos en el paso anterior (gracias al paso anterior sabemos que la derivada de G que estamos construyendo vale 0 de a 0). Por ello, buscamos g2 tal que

z2=0g1(t)dt+0zg2(t)dt=0+0zg2(t)dt=0zg2(t)dt,

por lo que g2 es la función dada por g2(t)=2z para cada z[0,12); ya que ddz(z2)=2z.

Ahora, para z[12,1), se tiene que

G(z)=13(1z)2,

y ahora buscamos una función g3 tal que

G(z)=0g1(t)dt+012g2(t)dt+12zg3(t)dt.

Desarrollando el lado derecho de esta última igualdad obtenemos la siguiente expresión

0g1(t)dt+012g2(t)dt+12zg3(t)dt=0+((12)202)+12zg3(t)dt=14+12zg3(t)dt,

y, por otro lado, desarrollando el lado izquierdo, obtenemos que

G(z)=13(1z)2=133z2=2(33z)2=3z12

por lo que la función g3 que buscamos debe de cumplir que

3z12=14+12zg3(t)dt,

que desarrollando un poco, nos dice que g3 debe de satisfascer

3z234=12zg3(t)dt.

Ahora, esto es algo que no se notó mucho en el paso anterior en el que obtuvimos g2, pero la expresión de la izquierda incluye las constantes de integración. Esto es de esperarse, pues se trata del teorema fundamental del cálculo.

Bien, ahora observa ddz(3z234)=32; así que g3(t)=32 es la función que queremos. Más aún, podemos confirmarlo integrando g3 de 12 a z:

12z32dt=3212z1dt=32(z12)=3z234,

por lo que g3(t)=32 es precisamente la función que buscamos.

Finalmente, para z1, tenemos que G(z)=1, por lo que nos queda encontrar g4 tal que

1=0g1(t)dt+012g2(t)dt+121g3(t)dt+1g4(t)dt,=00dt+0122tdtdt+12132dt+1g4(t)dt=0+14+(3234)+1g4(t)dt=14+(634)+1g4(t)dt=14+34+1g4(t)dt=1+1g4(t)dt

por lo que g4 debe de satisfacer que

1=1+1g4(t)dt,

o equivalentemente, que 1g4(t)dt=0. Similarmente a g1, la única función no-negativa que satisface esto es g4(z)=0, para cada z1. Así, colocando cada una de las funciones que hemos obtenido en el caso que le corresponde, obtenemos que la densidad de Z es la función g:RR dada por

g(z)={0si x<0,2zsi 0z<12,32si 12z<1,0si z1,

que gráficamente se ve como sigue:

Figura. Gráfica de g, la densidad de Z. Observa que g no es continua, algo que ya dejamos claro que no es necesario. Sin embargo, g claramente es una función integrable.

Usos de las variables aleatorias continuas

Las v.a.’s continuas tienen una cantidad no-numerable de valores que pueden tomar. Esto se debe a que R y muchos de sus subconjuntos son ejemplos de conjuntos infinitos no-numerables. Por ejemplo, los intervalos [0,1], [0,), (5,3] son ejemplos de conjuntos con cardinalidad infinita no-numerable. Por ello, las v.a.’s continuas se utilizan en fenómenos cuyo resultado amerita usar la precisión de los números reales.

Un primer ejemplo son los fenómenos en donde el resultado es un valor de tiempo. El tiempo que esperas hasta que llega un autobús a la parada en la que lo tomas; el tiempo que tarda en fallar algún aparato electrónico (el ejemplo clásico es el tiempo que tarda un foco en fundirse); el tiempo de vida que le queda a una persona (espeluznante, pero a los actuarios les interesa debido a los seguros de vida); etcétera. Cuando se habla de tiempos, suelen usarse v.a.’s continuas que tengan a [0,) como soporte para asegurar que el modelo contempla únicamente valores no-negativos. Las v.a.’s con distribución exponencial (que mencionamos en el primer ejemplo de esta entrada) son un ejemplo de este tipo de v.a.’s.

Otro ejemplo donde se usan v.a.’s continuas son los fenómenos financieros. El precio de un activo en un momento dado; la paridad cambiaria entre divisas; el valor que tomará la tasa de interés o de rendimiento en un instrumento financiero (como una anualidad) en el futuro; etcétera. La teoría de los procesos estocásticos es la rama de la probabilidad que se encarga del estudio de fenómenos aleatorios a través del tiempo, y resulta fundamental para el análisis de fenómenos financieros como los que aquí mencionamos.

Tarea moral

Los siguientes ejercicios son opcionales. Es decir, no formarán parte de tu calificación. Sin embargo, te recomiendo resolverlos para que desarrolles tu dominio de los conceptos abordados en esta entrada.

  1. Demuestra que si f:RR es una función que satisface f(x)0para cada xRyf(t)dt=1,entonces f es la función de densidad de alguna v.a. continua. Sugerencia: Define F:RR como sigue: para cada xR, F(x)=f(t)dt,y demuestra que F es una función de distribución.
  2. Demuestra que la función g:RR dada porg(x)={x+1si x[1,0),1xsi x[0,1],0en otro casoes la función de densidad de alguna v.a. continua. La gráfica de g puede apreciarse en la siguiente figura:
  3. Verifica que la función G del Ejemplo 3 es una función de distribución.

Más adelante…

Es importante que entiendas las diferencias que existen entre las v.a.’s discretas y las continuas. Repasa esta entrada y la anterior las veces que sea necesario para que no te confundas entre ambas. En la siguiente entrada veremos qué es lo que resulta de aplicarle funciones (transformaciones) a una v.a.

En cuanto a las aplicaciones, verás muchísimas más en materias posteriores, como en las materias de estadística. Por nuestro lado veremos varias distribuciones de probabilidad importantes (como la distribución exponencial) con las que seguramente te encontrarás en ese tipo de materias.

Entradas relacionadas