### copy data into 'dat' and examine data
dat <- dat.hine1989
dat
if (FALSE) {
### load metafor package
library(metafor)
### calculate risk differences and corresponding sampling variances
dat <- escalc(measure="RD", n1i=n1i, n2i=n2i, ai=ai, ci=ci, data=dat)
dat
### meta-analysis of risk differences using a random-effects model
res <- rma(yi, vi, data=dat)
res
}
Run the code above in your browser using DataLab