2012-04-16 10 views
9

Ich habe ein Paket 'bivpois' für R gefunden, das ein Modell für zwei verwandte Poisson-Prozesse auswertet (zum Beispiel die Anzahl der Tore von Heim- und Auswärtsmannschaft in einem Fußballspiel). Dieses Paket scheint jedoch in neueren Versionen von R nicht mehr verwendbar zu sein.Bivariate Poisson-Regression in R?

Gibt es eine sinnvolle Möglichkeit, die glm() - Funktion zu ändern, um einen ähnlichen Prozess durchzuführen, oder dieses ältere Paket auf meiner neuen Version von R auszuführen? Ich habe sehr wenig Literatur über diese Art von Prozessen gefunden und habe sehr wenig in Bezug auf die einfache Implementierung in anderen statistischen Paketen wie STATA gefunden.

Alle Vorschläge würden sehr geschätzt werden.

+1

Wenn Sie offen für die Verwendung von OpenBUGS sind, sieht es so aus, als könnte der Autor des Pakets hier Code enthalten (und vielleicht den Ansatz auf Seite 291-293 seines Buches diskutieren) http://stat-athens.aueb.gr/~ jbn/winbugs_book/Auch wenn du Bayesian nicht benutzen willst, kannst du das Buch immer noch ausprobieren, weil ich denke, dass er manchmal auch Frequency-R-Code präsentiert, obwohl ich mir darüber nicht 100% ig sicher bin. Ziehen Sie auch in Betracht, ihm eine E-Mail zu senden. –

+0

Google stellte dies her http://www.ats.ucla.edu/stat/r/dae/poissonreg.htm – Seth

+0

@Seth, ich glaube nicht, dass Ihr Link die Frage beantwortet (das OP will * bivariate * Poisson-Regression, nicht plain-vanilla Poisson-Regression) –

Antwort

6

Während CRAN keine aktuelle Binärdatei von Bivpois hostet, können Sie das Paket aus dem archivierten Quellcode erstellen (siehe http://cran.r-project.org/doc/manuals/R-exts.html#Checking-and-building-packages). Building Bivpois 0.50-3.1 aus der Quelle (verfügbar unter http://cran.r-project.org/src/contrib/Archive/bivpois/) funktioniert für mich auf R 2.15.0 Windows x64. Die gezippte Windows-Binärdatei, die ich gebaut habe, ist hier verfügbar: http://commondatastorage.googleapis.com/jthetzel-public/bivpois_0.50-3.1.zip.

+0

jthetzel, Danke für die Informationen. Ich ging voran und lud eine ältere Version von R herunter, und mein gezipptes Paket arbeitet daran. Ich denke also, eine Folgefrage ist, ob es irgendeine Funktion wie 'predprob' gibt (im Paket 'pscl'), oder ob ich manuell Vorhersagen vom Modell machen muss (sie sind nicht so geradlinig, wie ich gehofft hatte). Danke noch einmal. – Millsy11

+0

Sie könnten auch in das VGLM-Paket auf CRAN schauen. –

+0

Entschuldigung, das sollte VGAM-Paket sein! –