[R-meta] Bivariate generalized linear mixed model with {metafor}

[Arthur Albuquerque]

Hi all,

I need some help to figure out how to fit a bivariate generalized linear mixed model using metafor.

In the past year, the Journal of Clinical Epidemiology has posted several articles on a controversy between using risk ratio or odds ratio in meta-analyses. Summary of the controversy here:

George A. Wells , Commentary on Controversy and Debate 4 paper series:
Questionable utility of the relative risk in clinical research, Journal of Clinical Epidemiology (2021), doi: https://doi.org/10.1016/j.jclinepi.2021.09.016

One of the articles (https://doi.org/10.1016/j.jclinepi.2021.08.004) suggested fitting a bivariate generalized linear mixed model (BGLMM),  which "obtains effect estimates conditioning on baseline risks with the estimated model parameters, including the correlation parameter.”

They fitted this model using the PROC NLMIXED command in SAS. I would like to fit this model using metafor, could anyone help me by sending the appropriate code of this model with metafor::rma()?

Kind regards,

Arthur M. Albuquerque

Medical student
Universidade Federal do Rio de Janeiro, Brazil


	[[alternative HTML version deleted]]