Метаанализ: как убрать суммарный многоугольник на лесном участке?

Question

Я meta-analysing некоторые исследования и рисую лесной участок для моих результатов. Однако я бы хотел скрыть / удалить / удалить итоговый многоугольник внизу графика. Я уже удалил текст "RE Model", используя mlab=''. Кроме того, как я могу визуализировать p-value для взаимодействия в анализе подгруппы?

Большое спасибо!

### fit random-effects model (use slab argument to define study labels)
res <- rma(n1i=n1, m1i=m1, sd1i=sd1, n2i=n2, m2i=m2, sd2i=sd2, data=madata,  measure="MD", method = "REML",
           slab =paste(Author, Year, sep = "   "))



#FOREST PLOT WITH THE RE MODEL PER GROUP 
#png(height=400, width=800, pointsize = 12,  filename="plot_w_cairo.png", type="cairo")
forest(res, xlim=c(-40, 55), 
       ilab=cbind(madata$n1, madata$m1, madata$sd1, madata$n2, madata$m2, madata$sd2),
       ilab.xpos=c(-25,-22,-19,-15,-12,-9), cex=.75, ylim=c(-1, 15),
       order=order(madata$Cage), rows=c(11:8, 6:2))
op <- par(cex=.75, font=4)
text(-40, c(6.5,11.5), c("Static", "Expandable"), pos=4)
text(-35, -1, pos=4, cex=0.75, bquote(paste("Model for All Studies (Q = ",
     .(formatC(res$QE, digits=2, format="f")), ", df = ", .(res$k - res$p),
     ", p = ", .(formatC(res$QEp, digits=2, format="f")), "; ", I^2, " = ",
     .(formatC(res$I2, digits=1, format="f")), "%)")))
par(font=2)
text(c(-25,-22,-19,-15,-12,-9), 13.5, c("n1", "m1", "sd1", "n2", "m2", "sd2"))
text(c(-22,-11.5),     14, c("Preoperative", "Postopoperative"))
text(-40,            13.5, "Author and Year",  pos=4)
text(53,             14, "MD [95% CI]", pos=2)
res.e <- rma(n1i=n1, m1i=m1, sd1i=sd1, n2i=n2, m2i=m2, sd2i=sd2, data=madata, measure="MD",
            subset=(Cage=="Expandable"), method="REML")
res.s <- rma(n1i=n1, m1i=m1, sd1i=sd1, n2i=n2, m2i=m2, sd2i=sd2, data=madata, measure="MD",
            subset=(Cage=="Static"), method="REML")
text(-40, 7.5, pos=4, cex=0.75, bquote(paste("RE Model for Expandable Group (Q = ",
     .(formatC(res.e$QE, digits=2, format="f")), ", df = ", .(res.e$k - res.e$p),
     ", p = ", .(formatC(res.e$QEp, digits=2, format="f")), "; ", I^2, " = ",
     .(formatC(res.e$I2, digits=1, format="f")), "%)")))
text(-40, 1, pos=4, cex=0.75, bquote(paste("RE Model for Static Group (Q = ",
     .(formatC(res.s$QE, digits=2, format="f")), ", df = ", .(res.s$k - res.s$p),
     ", p = ", .(formatC(res.s$QEp, digits=2, format="f")), "; ", I^2, " = ",
     .(formatC(res.s$I2, digits=1, format="f")), "%)")))
addpoly(res.e, row=7.5, cex=0.75,  mlab="", font=6)
addpoly(res.s, row= 1, cex=0.75, mlab="", font=6)
title("Forest Plot ODI in MIS", font=2)