fit<-metabin(dat$e.e,dat$e.tot,dat$c.e,dat$c.tot, data=dat,
sm="RR", method="MH", comb.fixed=T,
comb.random=F,studlab=dat$lab)
Error in (function (yi, vi, sei, weights, ai, bi, ci, di, n1i, n2i, x1i, :
Fisher scoring algorithm did not converge. See 'help(rma)' for possible remedies.
数据无法估计，有两种方法解决
第一种：采用stata软件 重新估计，不会报错
第二种：换用metafor的rma包，但是rma包绘制的森林图不如meta包绘制的精美，需要手手动调整，同时，若估计固定效应，rma.mh估计结果与stata一致
dt.1<-rma.mh(measure="RR", ai=e.e, bi=e.tot-e.e,
ci=c.e, di=c.tot-c.e,n1i=e.tot,n2i = c.tot,
data=dat, method="FE", #固定效应
slab = dat$lab
)
森林图绘制
forest(dt.1,
atransf=exp, #转换成RR
slab=dat$lab,
ilab.xpos = c(-1.1,-0.8,-0.5,-0.2), #位置参数
ilab = cbind(e.e,e.tot,c.e,c.tot), #event total controlevent control total
showweights=T,
header = F,
xlab = "",
mlab = "Fixed-effect Model for All Studies",
psize = 1
)
text(c(-1.6,-1.1,-0.8,-0.5,-0.2,1.5,2), #X位置
dt.1$k+2, #Y轴位置
c("Study", "Events", "Total", "Events", "Total",
"Weight","RR 95%-CI"),
cex=1.2,font = 2)
text(c(-0.9,-0.4),
dt.1$k+3,
c("Experimental", "Control"),
cex=1.2, font=2)
