### copy data into 'dat' and examine data
dat <- dat.debruin2009
dat
if (FALSE) {
### load metafor package
library(metafor)
### calculate proportions and corresponding sampling variances
dat <- escalc(measure="PR", xi=xi, ni=ni, data=dat)
dat
### random-effects model
res <- rma(yi, vi, data=dat)
print(res, digits=2)
### mixed-effects meta-regression model with all predictors/covariates
res <- rma(yi, vi, mods = ~ scq + ethnicity + patients + select + sens, data=dat)
print(res, digits=3)
### mixed-effects meta-regression model with scq and ethnicity as predictors/covariates
res <- rma(yi, vi, mods = ~ scq + ethnicity, data=dat)
print(res, digits=3)
}
Run the code above in your browser using DataLab