Quantcast
Channel: Active questions tagged r - Stack Overflow
Viewing all articles
Browse latest Browse all 201894

How to visualize the percentage of each study contributed to the overall effect size summary of a binary outcome meta analysis using R?

$
0
0

I am practicing meta-analysis using the built-in "dat.bcg" dataset of "metafor" package. There I have performed all the primary analysis except the 'percentage of study contribution'. I have not found how to visualize this. I am attaching one image for your reference. enter image description here

I want to visualize this above-mentioned image.

Up until now I have been running the following codes.

enter code herelibrary("metafor") enter code heredata("dat.bcg", package = "metafor") enter code hereprint(dat.bcg, row.names = FALSE) enter code hereres <- rma(yi, vi, data = dat) enter code hereres enter code hereres = rma(ai = tpos, bi = tneg, ci = cpos, di = cneg, data = dat, measure = "OR") enter code hereres enter code hereconfint(res) enter code hereforest(res, slab = paste(dat$author, dat$year, sep = ", "), xlim = c(-16, 6), at = log(c(0.05, 0.25, 1, 4)), atransf = exp, ilab = cbind(dat$tpos, dat$tneg, dat$cpos, dat$cneg), ilab.xpos = c(-9.5, -8, -6, -4.5), cex = 0.75) enter code hereop <- par(cex = 0.75, font = 2) enter code heretext(c(-9.5, -8, -6, -4.5), 15, c("TB+", "TB-", "TB+", "TB-")) enter code heretext(c(-8.75, -5.25), 16, c("Vaccinated", "Control")) enter code heretext(-16, 15, "Author(s) and Year", pos = 4) enter code heretext(6, 15, "Relative Risk [95% CI]", pos = 2) enter code herepar(op)


Viewing all articles
Browse latest Browse all 201894

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>