[R-meta] question on forest plot formatring
Hello Wolfgang, Thank you very much for the clear explanation. It is indeed very helpful! Thanks a lot again. Kind regards, Gabriel On Mon, Jun 12, 2023 at 2:29?PM Viechtbauer, Wolfgang (NP) <
wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:
Dear Gabriel,
Maybe something like this?
library(metafor)
### copy BCG vaccine meta-analysis data into 'dat'
dat <- dat.bcg
### calculate log risk ratios and corresponding sampling variances (and use
### the 'slab' argument to store study labels as part of the data frame)
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat,
slab=paste(author, year, sep=", "))
### fit random-effects model
res <- rma(yi, vi, data=dat)
### a little helper function to add Q-test, I^2, and tau^2 estimate info
mlabfun <- function(text, res) {
list(bquote(paste(.(text),
"Q = ", .(formatC(res$QE, digits=2, format="f")),
", df = ", .(res$k - res$p),
", p ", .(metafor:::.pval(res$QEp, digits=2, showeq=TRUE, sep=" ")),
"; ",
I^2, " = ", .(formatC(res$I2, digits=1, format="f")), "%, ",
tau^2, " = ", .(formatC(res$tau2, digits=2, format="f")))))}
### set up forest plot (with 2x2 table counts added; the 'rows' argument is
### used to specify in which rows the outcomes will be plotted)
forest(res, xlim=c(-16, 4.6), at=log(c(0.05, 0.25, 1, 4)), atransf=exp,
ilab=cbind(tpos, tneg, cpos, cneg), ilab.xpos=c(-9.5,-8,-6,-4.5),
cex=0.75, ylim=c(-1, 27), order=alloc, rows=c(3:4,9:15,20:23),
mlab="RE Model for All Studies",
psize=1, header="Author(s) and Year")
text(-16, -1.8, mlabfun("", res)[[1]], pos=4, cex=0.75)
Modify this as needed.
As for the weight column -- if you want to place it somewhere else, I
would suggest to just include the weights (which you can obtain via
weights()) via 'ilab'. Then you have full control (via 'ilab.xpos') where
to place the values.
Best,
Wolfgang
-----Original Message----- From: R-sig-meta-analysis [mailto:
r-sig-meta-analysis-bounces at r-project.org] On
Behalf Of Michael Dewey via R-sig-meta-analysis Sent: Saturday, 10 June, 2023 16:07 To: Gabriel Cotlier Cc: Michael Dewey; R Special Interest Group for Meta-Analysis Subject: Re: [R-meta] question on forest plot formatring Since presumably my two original suggestions did not work I think I will leave this to someone more expert with metafor graphics. Michael On 10/06/2023 07:01, Gabriel Cotlier wrote:
Hello Michael, Thanks a lot for your response. Since I have a lot of lines in my forest plot (with a large number of correlations ) therefore I am using a small font size in order to make the summary of the model's results line more readable, I thought that maybe it is possible to print in the plot the model's summary results in two lines instead of one, just as a possible option...? Since I mentioned I was using this example as base code:
<
think maybe there exist some parameters in the function "mlabfun" could
be a
possible solution?
Regarding spacing the columns of at right hand side of the plot that is "weight %" and CI fisher's z" I could not find a parameter to space them a tab to make the plot more redabel. Thanks a lot. Kind regards, Gabriel