如何进一步划分小组以进行荟萃分析中的简单效果分析？

Question

我使用 metafor 在R中进行元分析。我学会了如何进行子组分析并绘制结果from the website，如下所示：

### 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交替; 80s系统，80s随机，80s交替; 90s系统，90s随机，90s交替。

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