Двумерная регрессия Пуассона в R?

Я нашел пакет bivpois для R, который оценивает модель для двух связанных пуассоновских процессов (например, количества голов домашней и гостевой команды в футбольном матче). Однако этот пакет, похоже, больше не может использоваться в более новых версиях R.

Есть ли разумный способ изменить функцию glm () для выполнения аналогичного процесса или запустить этот старый пакет в моей новой версии R? Я нашел очень мало литературы по такого рода процессам и очень мало с точки зрения простой реализации в других статистических пакетах, таких как STATA.

Любые предложения будут высоко ценится.


r poisson
person Millsy11    schedule 16.04.2012    source источник
comment
Если вы готовы использовать OpenBUGS, похоже, что автор пакета может включить сюда код (и может обсудить подход на страницах 291-293 своей книги) stat-athens.aueb.gr/~jbn/winbugs_book Даже если вы не хотите использовать байесовский метод, вы все равно можете почитать книгу, потому что я иногда думаю он также представляет частотный код R, хотя я не уверен на 100% в этом. Также подумайте о том, чтобы написать ему по электронной почте.   -  person Mark Miller    schedule 17.04.2012
comment
Google сделал это ats.ucla.edu/stat/r/dae/ poissonreg.htm   -  person Seth    schedule 17.04.2012
comment
@Seth, я не думаю, что ваша ссылка отвечает на вопрос (OP хочет двумерной регрессии Пуассона, а не простой ванильной регрессии Пуассона)   -  person Ben Bolker    schedule 17.04.2012
comment
Это не совсем bivpois решение, но вы можете проверить multicmp на наличие двумерного распределения COM-Пуассона github.com/diagdavenport/multicmp   -  person radek    schedule 25.02.2019

Ответы (2)


arrow_upward
7
arrow_downward

Хотя CRAN не содержит текущий двоичный файл bivpois, вы можете собрать пакет из заархивированного исходного кода (см. http://cran.r-project.org/doc/manuals/R-exts.html#Пакетыпроверкиисборки). Создание двойных точек 0.50-3.1 из исходного кода (доступно по адресу http://cran.r-project.org/src/contrib/Archive/bivpois/) у меня работает на R 2.15.0 Windows x64. Созданный мной заархивированный двоичный файл Windows доступен здесь: http://commondatastorage.googleapis.com/jthetzel-public/bivpois_0.50-3.1.zip.

person jthetzel    schedule 17.04.2012
comment
jthetzel, спасибо за информацию. Я пошел дальше и загрузил старую версию R, и мой заархивированный пакет работает с этим. Поэтому я предполагаю, что последующий вопрос заключается в том, есть ли какая-либо функция, подобная 'predprob' (в пакете 'pscl'), или мне нужно вручную делать прогнозы из модели (они не так просты, как я надеялся). Спасибо еще раз. - person Millsy11; 30.04.2012
comment
Вы также можете изучить пакет VGLM на CRAN. - person kjetil b halvorsen; 09.02.2017
comment
Извините, это должен быть пакет VGAM! - person kjetil b halvorsen; 09.02.2017

arrow_upward
0
arrow_downward

Вы можете смело ссылаться на моделирование шансов и тестирование неэффективности букмекеров, так как я изменил соответствующие функции внутри bivpois пакета.

person RYO ENG Lian Hu    schedule 05.12.2017
comment
вы не предоставляете никакого кода, я прав (в вашем github)? - person Fierce82; 01.10.2018
comment
Вы можете свободно просматривать исходный код с помощью пакета Rmodel. - person RYO ENG Lian Hu; 01.10.2018