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)