# NOT RUN {
### copy data into 'dat'
dat <- dat.yusuf1985
### to select a table for the analysis
tab <- "6" ### either: 6, 9, 10, 11, 12a, 12b
### to double-check total counts as reported in article
apply(dat[dat$table==tab,4:7], 2, sum, na.rm=TRUE)
### meta-analysis using Peto's one-step method
res <- rma.peto(ai=ai, n1i=n1i, ci=ci, n2i=n2i, data=dat, subset=(table==tab))
res
predict(res, transf=exp, digits=2)
# }
Run the code above in your browser using DataLab