2012-04-16 13 views
9

Znalazłem pakiet "bivpois" dla R, który ocenia model dla dwóch powiązanych procesów poissona (na przykład liczba bramek w domu i drużynie gości w meczu piłki nożnej). Jednak ten pakiet wydaje się przestać być użyteczny w nowszych wersjach R.Bivariate Poisson Regresja w R?

Czy istnieje rozsądny sposób zmodyfikowania funkcji glm() w celu wykonania podobnego procesu, lub uruchomienia tego starszego pakietu w mojej nowej wersji R? Znalazłem bardzo mało literatury na temat tego rodzaju procesów i znalazłem bardzo niewiele jeśli chodzi o łatwą implementację w innych pakietach statystycznych, takich jak STATA.

Wszelkie sugestie będą mile widziane.

+1

Jeśli jesteś otwarty na używanie OpenBUGS, wygląda na to, że autor pakietu może zawierać tutaj kod (i może omówić podejście na stronie 291-293 swojej książki) http://stat-athens.aueb.gr/~ jbn/winbugs_book/Nawet jeśli nie chcesz używać Bayesian, możesz nadal czytać książkę, ponieważ wydaje mi się, że czasami prezentuje również częsty kod R, chociaż nie jestem tego w 100% pewny. Rozważ także wysłanie e-maila do niego. –

+0

Google włączyło to http://www.ats.ucla.edu/stat/r/dae/poissonreg.htm – Seth

+0

@Seth, nie sądzę, że twój link odpowiada na pytanie (OP chce * bivariate * regres Poisson'a, nie zwykły-wanilia regres Poissona) –

Odpowiedz

6

Podczas gdy CRAN nie obsługuje aktualnego pliku binarnego bivpois, można zbudować pakiet z archiwizowanego kodu źródłowego (patrz http://cran.r-project.org/doc/manuals/R-exts.html#Checking-and-building-packages). Budowanie bpopois 0.50-3.1 ze źródła (dostępne pod http://cran.r-project.org/src/contrib/Archive/bivpois/) działa dla mnie na R 2.15.0 Windows x64. Spakowany binarny plik systemu Windows, który zbudowałem, jest dostępny tutaj: http://commondatastorage.googleapis.com/jthetzel-public/bivpois_0.50-3.1.zip.

+0

jthetzel, Dzięki za informację. Poszedłem do przodu i pobrałem starszą wersję R, a mój skompresowany pakiet działał na tym. Tak więc myślę, że następne pytanie brzmi: czy istnieje jakakolwiek funkcja, taka jak "predprob" (w pakiecie "pscl"), czy też muszę ręcznie tworzyć prognozy z modelu (nie są one tak proste, jak się spodziewałem). Dzięki jeszcze raz. – Millsy11

+0

Można również zajrzeć do pakietu VGLM na CRAN. –

+0

Przepraszam, to powinien być pakiet VGAM! –

Powiązane problemy