1

我正在使用metafor在 R 中进行元分析。我学会了如何进行亚组分析并从网站上绘制结果，如下所示：

### fit random-effects model in the three subgroups
res.s <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, measure="RR",
             subset=(alloc=="systematic"), method="REML")
res.r <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, measure="RR",
             subset=(alloc=="random"), method="REML")
res.a <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, measure="RR",
             subset=(alloc=="alternate"), method="REML")

### add summary polygons for the three subgroups
addpoly(res.s, row=18.5, cex=0.75, atransf=exp, mlab="")
addpoly(res.r, row= 7.5, cex=0.75, atransf=exp, mlab="")
addpoly(res.a, row= 1.5, cex=0.75, atransf=exp, mlab="")

我的问题是，我可以进一步划分这三个子组吗？例如，以年划分以及分配：70s系统，70s随机，70s交替；80年代系统，80年代随机，80年代交替；90 年代系统，90 年代随机，90 年代交替。

而且，我希望这是在元回归中处理交互的方法。

4

0 回答 0