### load BCG vaccine data
data(dat.bcg)
### calculate log risk ratios and corresponding sampling variances
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
dat.bcg$yi <- dat$yi
dat.bcg$vi <- dat$vi
### random-effects model
res <- rma(yi, vi, data=dat.bcg, method="REML")
ranktest(res)
Run the code above in your browser using DataLab