Autor Tema: Hallar probabilidad a partir de vectores con ayuda de MATLAB

0 Usuarios y 1 Visitante están viendo este tema.

02 Marzo, 2010, 12:03 am
Leído 1812 veces

newone

  • $$\Large \color{#6a84c0}\pi$$
  • Mensajes: 14
  • Karma: +0/-0
Me dan la siguiente función de densidad para resolver ciertos ejercicios en MATLAB:

\( f_x(x)=\begin{Bmatrix} 8*exp(-8*(x-1)) & \mbox{ si }& 1\leq{x\leq{\infty}}\\0 & \mbox{,}& otros\end{matrix} \)

Primero me ponen a encontrar la función inversa de la función de distirbución, luego me ponen a generar un vector de 1000 elementos de números aleatorios con distribución uniforme para después generar un nuevo vector (llamemos a este último vector el vector \( v \)) de números aleatorios con función de probabilidad \( f_x \) substituyendo el vector anterior (el de los números aleatorios con distribución uniforme) en la función inversa hallada al principio. Hasta ahí todo perfecto y sencillo; pero luego me dicen:

"Para el vector \( v \) aplique la siguiente transformación:

\( Y = 8*ln(X)+8 \)

Luego halle la probabilidad \( P(9\leq{y\leq{10}}) \)"

Por lo que sé, la primera parte la puedo hacer remplazando el vector \( v \) en esa transformación, lo que me generaría un nuevo vector de 1000 elementos (llamémoslo el vector \( y \)). Lo que no sé es cómo hallar la probabilidad solicitada con ayuda del vector \( y \) (de hecho ni siquiera sabía que se podía hallar una probabilidad de esta manera). ¿Alguien tiene alguna idea? Cualquier ayuda me sería útil.

PD: me disculparán por el uso pobre de LaTeX, la verdad es algo que no uso frecuentemente.

EDIT: Por cierto, alguien me había dicho que eso se podía hacer usando operadores relaciones, la función sum() y la función length() pero la verdad yo no veo para qué se usarían semejantes cosas y por lo general no hago las cosas a menos de que entienda la lógica detrás de ellas...

02 Marzo, 2010, 05:03 am
Respuesta #1

topo23

  • $$\Large \color{#c88359}\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 937
  • Karma: +0/-0
Si lo quieres hacer todo con lógica deberías chequear que definición te han dado de función de densidad.

Por ejemplo esta defincion http://e-stadistica.bio.ucm.es/glosario/def_fun_densidad.html, se puede aplicar.

En concreto si h es la funcion de densidad de la v.a. Z, entonces \( \mathrm{Pr}(Z \le z) = \int\limits_{-\infty}^z h(t)dt \)
.

02 Marzo, 2010, 05:47 pm
Respuesta #2

newone

  • $$\Large \color{#6a84c0}\pi$$
  • Mensajes: 14
  • Karma: +0/-0
Sí, yo sé hacer ese cálculo (de hecho ya lo hice). Lo que no logro ver es cómo calcular una probabilidad aproximada con tan solo usar ese arreglo \( y \) de 1000 elementos. Se me han ocurrido ciertas ideas pero el margen de error es ciertamente alto a comparación con la integral  :-\.

03 Marzo, 2010, 01:33 am
Respuesta #3

topo23

  • $$\Large \color{#c88359}\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 937
  • Karma: +0/-0
\( Y = 8*ln(X)+8 \)

Luego halle la probabilidad \( P(9\leq{y\leq{10}}) \)"

La verdad no esta claro que es lo que te piden en el enunciado, en todo caso \( P(9\le Y \le 10)=P(9\le 8\ln(X)+8\le 10)=\dots=P(\exp(1/8)\le X\le\exp(2/8)) \)

También podría ser que te piden calcular una aproximación a la probabilidad usando Monte Carlo, para eso podrías usar el vector.

O tal vez una combinación de ambos...
.