Estimación de máxima verosimilitud de un modelo de Poisson para el fútbol (soccer) de las puntuaciones de

Question

Tengo un conjunto de resultados de fútbol y quiero hacer una probabilty modelo de resultados de fútbol, como se describe en Dixon, Coles (1997, http://www.math.ku.dk/~rolf/docencia/tesis/DixonColes.pdf). Ellos estiman los parámetros basados en máxima verosimilitud y el modelo asume que las variables son independientes de la distribución de Poisson.

Ahora tengo 2 preguntas:

• ¿Cómo debo transformar los datos de entrada en el modelo?
• ¿Qué paquetes en R son los mejores para la estimación de máxima verosimilitud de esta naturaleza?

Answer 1

El bivpois paquete, escrito por Karlis y Ntzoufras, utiliza el EM-algoritmo para la estimación de máxima verosimilitud, en este tipo de bivariante de Poisson modelos (y algunas generalizaciones de ellos). No creo que sea en CRAN más, pero se puede encontrar aquí.

Para obtener más información, consulte la descripción del paquete en el Journal of Statistical Software y los relacionados con el papel por Karlis y Ntzoufrasin en El Estadístico (que es una continuación de la obra por Dixon y Coles).

El paquete contiene ejemplos donde se puede ver cómo dar formato a los datos. Han pasado algunos años desde que jugué un rato con él, pero de lo que recuerdo que era bastante fácil de usar.

Answer 2

Usted debe buscar en el VGAM paquete tiene funciones para adaptarse a las Bradley-Terry modelo que se describe en los enlaces de las preguntas en los comentarios.