2012-04-16 32 views
9

Ho trovato un pacchetto "bivpois" per R che valuta un modello per due processi di poisson correlati (ad esempio, il numero di goal della casa e la squadra in trasferta in una partita di calcio). Tuttavia, questo pacchetto sembra non essere più utilizzabile nelle versioni più recenti di R.Regressione di Poisson bivariata in R?

Esiste un modo ragionevole per modificare la funzione glm() per eseguire un processo simile o eseguire questo pacchetto precedente sulla mia nuova versione di R? Ho trovato pochissima letteratura su questo tipo di processi e ho trovato molto poco in termini di facile implementazione in altri pacchetti statistici come STATA.

Qualsiasi suggerimento sarebbe molto apprezzato.

+1

Se si è aperti all'utilizzo di OpenBUGS sembra che l'autore del pacchetto possa includere il codice qui (e potrebbe discutere l'approccio a pagina 291-293 del suo libro) http://stat-athens.aueb.gr/~ jbn/winbugs_book/Anche se non vuoi usare Bayesian potresti ancora controllare il libro perché penso che a volte presenti anche il codice R frequentista, anche se non ne sono sicuro al 100%. Inoltre, considera di inviarlo via email. –

+0

Google ha rivelato questo http://www.ats.ucla.edu/stat/r/dae/poissonreg.htm – Seth

+0

@Seth, non penso che il tuo link risponda alla domanda (l'OP vuole * bivariata * regressione di Poisson, non semplice-regressione di Poisson alla vaniglia) –

risposta

6

Mentre CRAN non ospita un file binario corrente di bivpo, è possibile creare il pacchetto dal codice sorgente archiviato (vedere http://cran.r-project.org/doc/manuals/R-exts.html#Checking-and-building-packages). Building bivpois 0.50-3.1 from source (disponibile a http://cran.r-project.org/src/contrib/Archive/bivpois/) funziona per me su R 2.15.0 Windows x64. Il file binario di Windows con zip che ho creato è disponibile qui: http://commondatastorage.googleapis.com/jthetzel-public/bivpois_0.50-3.1.zip.

+0

jthetzel, Grazie per l'informazione. Sono andato avanti e ho scaricato una versione precedente di R e il mio pacchetto zippato funziona su questo. Quindi suppongo che una domanda successiva sia se esiste una funzione come "predprob" (nel pacchetto "pscl"), o se devo fare manualmente delle previsioni dal modello (non sono così semplici come speravo). Grazie ancora. – Millsy11

+0

Si potrebbe anche esaminare il pacchetto VGLM su CRAN. –

+0

Siamo spiacenti, questo dovrebbe essere il pacchetto VGAM! –

Problemi correlati