### load data
dat <- get(data(dat.bourassa1996))
### calculate log(OR) and corresponding sampling variance with 1/2 correction
dat <- escalc(measure="OR", ai=ai, bi=bi, ci=ci, di=di, data=dat, add=1/2, to="all")
### overall association between handedness and eyedness
res <- rma(yi, vi, data=dat, subset=sex=="combined")
res
predict(res, transf=exp, digits=2)
Run the code above in your browser using DataLab