Message-ID: <14a0912cbca3412db7015f6cd6cb8f18@UM-MAIL3214.unimaas.nl>
Date: 2019-10-22T15:02:58Z
From: Wolfgang Viechtbauer
Subject: [R-meta] metafor - averaging over rma result
In-Reply-To: <OFEFDD8EEE.DC2566D6-ONC125849B.004A7DC8-C125849B.004B6C05@lotus.uzh.ch>
Dear Stephanie,
My suggestion would be to fit a sensible model that uses 'disorder' as a predictor and then plot the model coefficients using the forest() function. Here is an example, not using VRs, but the idea is the same:
library(metafor)
dat <- dat.bourassa1996
# calculate log(OR) and corresponding sampling variances
dat <- escalc(measure="OR", ai=lh.le, bi=lh.re, ci=rh.le, di=rh.re, data=dat)
dat
# fit model
res <- rma(yi, vi, mods = ~ eye_assess - 1, data=dat, subset=sex=="combined")
res
# by removing the intercept, the coefficients are the estimated log(OR)s for each level of 'eye_assess'
# pass coefficients and variance of the coefficients to forest() function
forest(coef(res), diag(vcov(res)), slab=names(coef(res)), atrans=exp,
at = c(log(c(1, 2, 4, 8, 16, 32))), xlab="Odds Ratio (log scale)")
Best,
Wolfgang
-----Original Message-----
From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-project.org] On Behalf Of stephanie.winkelbeiner at bli.uzh.ch
Sent: Tuesday, 22 October, 2019 15:44
To: r-sig-meta-analysis at r-project.org
Subject: [R-meta] metafor - averaging over rma result
Dear experts,
for a meta-analysis, we calculated the variability ratio (VR) of all included studies
using the rma function of the metafor package.
We would like to plot the results in a forest plot grouped by disorder for a less
busy presentation (compared to plotting the VR for every study).
Yet, averaging over the VRs, especially over the confidence intervals seems
worrying.
Do you have any suggestions?
Thank you,
Stephanie