2012-04-16 7 views
9

Encontré un paquete 'bivpois' para R que evalúa un modelo para dos procesos de poisson relacionados (por ejemplo, el número de goles por el equipo local y el equipo visitante en un partido de fútbol). Sin embargo, este paquete parece no ser más útil en las versiones más recientes de R.Regresión de Poisson bivariado en R?

¿Existe alguna forma razonable de modificar la función glm() para hacer un proceso similar o ejecutar este paquete anterior en mi nueva versión de R? He encontrado muy poca literatura sobre este tipo de procesos y he encontrado muy poco en términos de fácil implementación en otros paquetes estadísticos como STATA.

Cualquier sugerencia sería muy apreciada.

+1

Si está abierto a usar OpenBUGS, parece que el autor del paquete puede incluir código aquí (y puede discutir el enfoque en la página 291-293 de su libro) http://stat-athens.aueb.gr/~ jbn/winbugs_book/Incluso si no quieres usar Bayesian, igual puedes consultar el libro porque a veces él también presenta el código R frecuentista, aunque no estoy 100% seguro de eso. Además, considera enviar un correo electrónico a él. –

+0

Google lo presentó http://www.ats.ucla.edu/stat/r/dae/poissonreg.htm – Seth

+0

@Seth, no creo que su enlace responda la pregunta (el OP quiere * regresión de Poisson bivariante *, no regresión de Poisson simple-vainilla) –

Respuesta

6

Mientras que CRAN no aloja un binario actual de bivpois, puede compilarlo a partir del código fuente archivado (consulte http://cran.r-project.org/doc/manuals/R-exts.html#Checking-and-building-packages). Construir bivpois 0.50-3.1 desde la fuente (disponible en http://cran.r-project.org/src/contrib/Archive/bivpois/) funciona para mí en R 2.15.0 Windows x64. El binario comprimido de Windows que construí está disponible aquí: http://commondatastorage.googleapis.com/jthetzel-public/bivpois_0.50-3.1.zip.

+0

jthetzel, gracias por la información. Seguí adelante y descargué una versión anterior de R y mi paquete comprimido funciona en eso. Así que supongo que una pregunta de seguimiento es si hay alguna función como 'predprob' (en el paquete 'pscl'), o si debo hacer predicciones del modelo manualmente (no son tan sencillas como esperaba). Gracias de nuevo. – Millsy11

+0

También podría consultar el paquete VGLM en CRAN. –

+0

Lo siento, ese debería ser el paquete VGAM. –

Cuestiones relacionadas