I am planning to produces consonance intervals for a meta-analysis Project. I need to use Mantel-Haenszel method to pool binary outcome. Here is an issue I ran in to.
library(metafor)
library(concurve)
#dataset from metafor
data(dat.linde2005)
#pool binary outcomes using RR with fixed effects in linear model with rma.uni function
justsomeRR<-rma.uni(ai=ai,bi=(n1i-ai),ci=ci,di=(n2i-ci),data=dat.linde2005,measure="RR",add=0.5,digits=2,method="FE")
#get consonance intervals
metaf<-curve_meta(justsomeRR,steps=100)
#pool binary outcomes using RR based on fixed effects Mantel-Haenszel methods with rma.mh function
justsomeRR.2<-rma.mh(ai=ai,bi=(n1i-ai),ci=ci,di=(n2i-ci),data=dat.linde2005,measure="RR",add=0.5,digits=2)
#a problem with consonance intervals
metaf2<-curve_meta(justsomeRR.2,steps=100)