# copy data into 'dat' and examine data
dat <- dat.axfors2021
dat
if (FALSE) {
# load metafor package
library(metafor)
# calculate log odds ratios and corresponding sampling variances
dat <- escalc(measure="OR", ai=hcq_arm_event, n1i=hcq_arm_total,
ci=control_arm_event, n2i=control_arm_total, data=dat)
# meta-analysis Hydroxychloroquine
res_hcq <- rma(yi, vi, subset=(hcq_cq=="hcq"), slab = id, data=dat)
print(res_hcq, digits=2)
# meta-analysis Chloroquine
res_cq <- rma(yi, vi, subset=(hcq_cq=="cq"), slab = id, data=dat)
print(res_cq, digits=2)
}
Run the code above in your browser using DataLab