# NOT RUN {
d = metafor::escalc(measure="RR", ai=tpos, bi=tneg,
ci=cpos, di=cneg, data=metadat::dat.bcg)
# calculate calibrated estimates
d$calib = calib_ests( yi = d$yi,
sei = sqrt(d$vi) )
# look at 5 studies with the largest calibrated estimates
d = d[ order(d$calib, decreasing = TRUE), ]
d$trial[1:5]
# look at kernel density estimate of calibrated estimates
plot(density(d$calib))
# }
Run the code above in your browser using DataLab