### copy data into 'dat' and examine data
dat <- dat.baskerville2012
dat
if (FALSE) {
### load metafor package
library(metafor)
### random-effects model
res <- rma(smd, sei=se, data=dat, method="DL")
print(res, digits=2)
### funnel plot
funnel(res, xlab="Standardized Mean Difference", ylim=c(0,0.6))
### rank and regression tests for funnel plot asymmetry
ranktest(res)
regtest(res)
### meta-regression analyses examining various potential moderators
rma(smd, sei=se, mods = ~ score, data=dat, method="DL")
rma(smd, sei=se, mods = ~ alloconc, data=dat, method="DL")
rma(smd, sei=se, mods = ~ blind, data=dat, method="DL")
rma(smd, sei=se, mods = ~ itt, data=dat, method="DL")
rma(smd, sei=se, mods = ~ duration, data=dat, method="DL")
rma(smd, sei=se, mods = ~ tailor, data=dat, method="DL")
rma(smd, sei=se, mods = ~ pperf, data=dat, method="DL")
rma(smd, sei=se, mods = ~ I(meetings * hours), data=dat, method="DL")
}
Run the code above in your browser using DataLab