Расчет прогнозируемой вероятности с GEE в R
Я использую пакет "geepack" для запуска регрессии Пуассона (link="log") и получения относительных рисков. В качестве альтернативы, я мог бы запустить биномиальную регрессию (link="identity") и получить различия в риске. Но в любом случае я бы хотел получить абсолютный риск со стандартными ошибками (для расчета доверительного интервала). Многие люди используют логистическую регрессию, а затем используют пакет "поля" для получения прогнозируемой вероятности, однако я не могу использовать glm, поскольку у меня есть кластеры, и я не могу создать смешанную модель, поскольку она не сходится таким образом. Есть идеи?