La distribución de Skellam es la distribución de probabilidad discreta de la diferencia de dos variables aleatorias estadísticamente independientes y cada una de ellas distribuida según Poisson con sus respectivos valores esperados y . Es útil para describir las estadísticas de la diferencia de dos imágenes con ruido fotónico simple , así como para describir la distribución de puntos en deportes donde todos los puntos anotados son iguales, como el béisbol , el hockey y el fútbol .
La distribución también es aplicable a un caso especial de la diferencia de variables aleatorias de Poisson dependientes, pero sólo al caso obvio donde las dos variables tienen una contribución aleatoria aditiva común que se cancela mediante la diferenciación: véase Karlis y Ntzoufras (2003) para obtener detalles y una aplicación.
La función de masa de probabilidad para la distribución de Skellam para una diferencia entre dos variables aleatorias independientes distribuidas por Poisson con medias y está dada por:
donde I k ( z ) es la función de Bessel modificada de primera especie. Como k es un entero tenemos que I k ( z )= I |k| ( z ).
Derivación
La función de masa de probabilidad de una variable aleatoria distribuida por Poisson con media μ está dada por
para (y cero en caso contrario). La función de masa de probabilidad de Skellam para la diferencia de dos recuentos independientes es la convolución de dos distribuciones de Poisson: ( Skellam , 1946)
Dado que la distribución de Poisson es cero para valores negativos del recuento , la segunda suma solo se toma para aquellos términos donde y . Se puede demostrar que la suma anterior implica que
de modo que:
donde I k (z) es la función de Bessel modificada de primera especie. El caso especial para está dado por Irwin (1937):
Utilizando los valores límite de la función de Bessel modificada para argumentos pequeños, podemos recuperar la distribución de Poisson como un caso especial de la distribución de Skellam para .
Propiedades
Como es una función de probabilidad discreta, la función de masa de probabilidad de Skellam está normalizada:
Sabemos que la función generadora de probabilidad (pgf) para una distribución de Poisson es:
De ello se deduce que la pgf, , para una función de masa de probabilidad de Skellam será:
Obsérvese que la forma de la función generadora de probabilidad implica que la distribución de las sumas o las diferencias de cualquier número de variables independientes distribuidas según Skellam también se distribuyen según Skellam. A veces se afirma que cualquier combinación lineal de dos variables distribuidas según Skellam también se distribuye según Skellam, pero esto claramente no es cierto, ya que cualquier multiplicador distinto de este cambiaría el soporte de la distribución y alteraría el patrón de momentos de una manera que ninguna distribución de Skellam puede satisfacer.
La función generadora de momentos viene dada por:
que da como resultado los momentos brutos m k . Definir:
Entonces los momentos crudos m k son
Los momentos centrales M k son
La media , la varianza , la asimetría y el exceso de curtosis son respectivamente:
La función generadora de cumulantes viene dada por:
Lo que produce los cumulantes :
Para el caso especial cuando μ 1 = μ 2 , una expansión asintótica de la función de Bessel modificada del primer tipo produce para μ grande:
(Abramowitz y Stegun 1972, p. 377). Además, para este caso especial, cuando k también es grande y del orden de la raíz cuadrada de 2μ, la distribución tiende a una distribución normal :
Estos resultados especiales pueden extenderse fácilmente al caso más general de medias diferentes.
Límites de peso por encima de cero
Si , con , entonces
Los detalles se pueden encontrar en la distribución de Poisson#Razas de Poisson
Referencias
- Abramowitz, Milton; Stegun, Irene A., eds. (junio de 1965). Handbook of mathematics functions with formulas, graphs, and mathematics tables (Republicación íntegra y sin modificaciones [der Ausg.] 1964, 5.ª edición de Dover). Dover Publications. págs. 374–378. ISBN 0486612724. Recuperado el 27 de septiembre de 2012 .
- Irwin, JO (1937) "La distribución de frecuencia de la diferencia entre dos variables independientes que siguen la misma distribución de Poisson". Journal of the Royal Statistical Society : Series A , 100 (3), 415–416. JSTOR 2980526
- Karlis, D. y Ntzoufras, I. (2003) "Análisis de datos deportivos utilizando modelos de Poisson bivariados". Journal of the Royal Statistical Society, Serie D , 52 (3), 381–393. doi :10.1111/1467-9884.00366
- Karlis D. y Ntzoufras I. (2006). Análisis bayesiano de las diferencias de datos de recuento. Statistics in Medicine , 25, 1885–1905. [1]
- Skellam, JG (1946) "La distribución de frecuencias de la diferencia entre dos variables de Poisson pertenecientes a diferentes poblaciones". Journal of the Royal Statistical Society, Serie A , 109 (3), 296. JSTOR 2981372
Véase también