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

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

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

Будем признательны за любые предложения.