dat.yusuf1985
character
table number
id character
trial id number
trial character
trial name or first author
ai numeric
number of deaths/reinfarctions in treatment group
n1i numeric
number of patients in treatment group
ci numeric
number of deaths/reinfarctions in control group
n2i numeric
number of patients in control group
}rma.peto
function.### load data
dat <- get(data(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