# show (some) data
head(dat.crisafulli2020)
if (FALSE) {
# compute logarithmic proportions and associated standard errors
library(metafor)
logp <- escalc(measure="PLN",
xi=cases, ni=total, slab=study,
data=dat.crisafulli2020)
# perform meta-analysis
rma01 <- rma.uni(logp)
# show results
rma01
# illustrate in a forest plot
forest(rma01, xlim=c(-12,-5))
}
Run the code above in your browser using DataLab