Apunte

1. Hallar el estimador máximo verosímil de υ (tasa promedio) de una población que distribuye Poisson. Considere el procedimiento descrito en los contenidos de la semana. Se valorará el proceso en detalle. Interprete su resultado.

La Distribución de Poisson es una distribución de variable discreta, que se aplica principalmente en situaciones que buscan medir el número de hechos de cierto tipo que se pueden producir en un intervalo de tiempo o de espacio, bajo supuestos de aleatoriedad y ciertas restricciones, especializándose en la probabilidad de ocurrencia de sucesos con posibilidades muy pequeñas .

La función de probabilidad se representa como:

P(X=x)=(e^(-v)*v^x)/x!

Donde υ es la tasa de incidencia media del evento, o tasa promedio, es decir un parámetro que representa el número de veces que ocurra un fenómeno en un determinado período de tiempo. Por su parte x es el número de ocurrencias del evento, por lo que esta función nos entrega la probabilidad de que un acontecimiento suceda precisamente x veces.

El ejercicio planteado pide encontrar el Estimador de Máxima Verosimilitud de la tasa promedio v, por lo que se debe utilizar el Método de Máxima Verosimilitud. En este sentido cabe recordar que la Función de Verosimilitud para una muestra de tamaño n de variables aleatorias x, que definen a una población con parámetro θ, a través de una función de probabilidad f(x;θ), se estructura como:

L(x_1,x_1,⋯,x_n;θ)=∏_(i=1)^n▒f(x_i;θ)

Por lo que para encontrar el Estimador de Máxima Verosimilitud, lo primero es tomar esta Función de Verosimilitud y transformarla en una función monótona a través de logaritmo natural, para luego derivar y luego igualando a 0. Por lo tanto se tiene:

L(x_1,x_1,⋯,x_n;υ)=∏_(i=1)^n▒(e^(-v)*v^(x_i ))/(x_i !)

L(x_1,x_1,⋯,x_n;υ)=(e^(-v)*v^(x_1 ))/(x_1 !)*(e^(-v)*v^(x_2 ))/(x_2 !)*⋯*(e^(-v)*v^(x_n ))/(x_n !)

L(x_1,x_1,⋯,x_n;υ)=(e^(-nv)*v^∑▒x_i )/(∏▒〖x_i !〗)

Donde ∏▒〖x_i !〗es una constante que se puede denotar por k, por lo que al aplicar logaritmo natural a la función queda de la siguiente forma:

ln(L(x;υ) )=ln((e^(-nv)*v^∑▒x_i )/(∏▒〖x_i !〗))

Dado que por propiedades de los logaritmos se tiene que ln(a*b)=ln⁡(a)+ln(b), y en paralelo que ln(a/b)=ln⁡(a)-ln(b), y además que ln(a^n )=n*ln(a)y ln(e)=1, la expresión queda:

ln⁡(L(x;υ) )=-nv+∑_(i=1)^n▒x_i ln(v)-ln(k)

Que es la fórmula que hay que maximizar a través de la derivada del parámetro, lo que resulta en lo siguiente:

∂ln⁡(L(x;υ) )/∂υ=-n+(∑▒x_i

...