# NOT RUN {
# calculate effect sizes for example dataset
d = metafor::escalc(measure="RR", ai=tpos, bi=tneg,
ci=cpos, di=cneg, data=metadat::dat.bcg)
# test H0: the median is -0.3
# using only R = 100 for speed, but should be much larger (e.g., 2000) in practice
pct_pval( yi = d$yi,
sei = sqrt(d$vi),
mu = -0.3,
pct = 0.5,
R = 100 )
# }
Run the code above in your browser using DataLab