### copy data into 'dat' and examine data
dat <- dat.hartmannboyce2018
head(dat, 10)
if (FALSE) {
### load metafor package
library(metafor)
### turn treatment into a factor with the desired ordering
dat$treatment <- factor(dat$treatment, levels=unique(dat$treatment))
### meta-analysis per treatment using the M-H method
lapply(split(dat, dat$treatment), function(x)
rma.mh(measure="RR", ai=x.nrt, n1i=n.nrt,
ci=x.ctrl, n2i=n.ctrl, data=x, digits=2))
### all combined
rma.mh(measure="RR", ai=x.nrt, n1i=n.nrt,
ci=x.ctrl, n2i=n.ctrl, data=dat, digits=2)
}
Run the code above in your browser using DataLab