A Modified Poisson Regression - Zhou et al Approach in R-INLA

308 views
Skip to first unread message

Prince Michael Amegbor

unread,
Feb 15, 2024, 1:59:40 PM2/15/24
to R-inla discussion group
Hi everyone,

I'm exploring the implementation of a modified Poisson regression with robust variance, following the recommendation by Zhou et al. in 2004. This is particularly useful for modeling a binary outcome with a common prevalence ratio (PR).

In R, you can achieve this using the glm model and the sandwich package for robust variance. Here's an example:

library("sandwich")
library("lmtest")

model <- glm(outcome ~ exposure, family = poisson(link = log), data = dataset)
coeftest(model, vcov = sandwich)

This model, expressed as log(λ) = β₀ + β₁ ⋅ exposure, calculates standard errors with robust variance for more reliable results.


I'm curious if there's a way to implement this approach in the r-inla.
 


Helpdesk (Haavard Rue)

unread,
Feb 17, 2024, 6:57:22 AM2/17/24
to Prince Michael Amegbor, R-inla discussion group
not supported
> --
> You received this message because you are subscribed to the Google
> Groups "R-inla discussion group" group.
> To unsubscribe from this group and stop receiving emails from it, send
> an email to r-inla-discussion...@googlegroups.com.
> To view this discussion on the web, visit
> https://groups.google.com/d/msgid/r-inla-discussion-group/97dd2fa4-69a5-41e8-8366-30290c8ef4c0n%40googlegroups.com
> .

--
Håvard Rue
he...@r-inla.org
Reply all
Reply to author
Forward
0 new messages