パッケージmetaforでフォレストプロットを結合する方法は？

Question

は、このコードは、次の組み合わせ森林プロットを与える（http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroupsで与えられる）以下のコードを

library(metafor) 

### decrease margins so the full space is used 
par(mar=c(4,4,1,2)) 

### fit random-effects model (use slab argument to define study labels) 
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, measure="RR", 
      slab=paste(author, year, sep=", "), method="REML") 

### set up forest plot (with 2x2 table counts added; rows argument is used 
### to specify exactly in which rows the outcomes will be plotted) 
forest(res, xlim=c(-16, 6), at=log(c(.05, .25, 1, 4)), atransf=exp, 
     ilab=cbind(dat.bcg$tpos, dat.bcg$tneg, dat.bcg$cpos, dat.bcg$cneg), 
     ilab.xpos=c(-9.5,-8,-6,-4.5), cex=.75, ylim=c(-1, 27), 
     order=order(dat.bcg$alloc), rows=c(3:4,9:15,20:23), 
     xlab="Relative Risk", mlab="RE Model for All Studies", psize=1) 

### set font expansion factor (as in forest() above) and use bold italic 
### font and save original settings in object 'op' 
op <- par(cex=.75, font=4) 

### add text for the subgroups 
text(-16, c(24,16,5), pos=4, c("Systematic Allocation", 
           "Random Allocation", 
           "Alternate Allocation")) 

### switch to bold font 
par(font=2) 

### add column headings to the plot 
text(c(-9.5,-8,-6,-4.5), 26, c("TB+", "TB-", "TB+", "TB-")) 
text(c(-8.75,-5.25),  27, c("Vaccinated", "Control")) 
text(-16,    26, "Author(s) and Year",  pos=4) 
text(6,     26, "Relative Risk [95% CI]", pos=2) 

### set par back to the original settings 
par(op) 

### 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=.75, atransf=exp, mlab="RE Model for Subgroup") 
addpoly(res.r, row= 7.5, cex=.75, atransf=exp, mlab="RE Model for Subgroup") 
addpoly(res.a, row= 1.5, cex=.75, atransf=exp, mlab="RE Model for Subgroup") 

仮定します

だから私の質問は、総相対リスクとのそれに対応する結果を削除する方法ですプロットの最下部？言い換えれば、私はちょうどいくつかの森林プロットをプロットし、同じプロットでそれらを表示したいが、それらを組み合わせて結果の結果を表示したくない。

Answer 1

forest(...)の部分には、addfit=FALSEを使用してください。次に、ylimの値を調整する必要があります。この場合、ylim=c(1, 27)が動作するはずです。