### copy data into 'dat' and examine data
dat <- dat.michael2013
dat
if (FALSE) {
### load metafor package
library(metafor)
### Data prep
# yi and vi are already provided, but here's how you would use escalc() to obtain
# a raw-mean difference and its variance.
# Note the measure parameter is "MD" for 'raw mean difference'
dat <- metafor::escalc(
measure = "MD",
m1i = Brain_m,
m2i = No_brain_m,
sd1i = Brain_s,
sd2i = No_brain_s,
n1i = Brain_n,
n2i = No_brain_n,
data = dat
)
### meta-analysis using a random-effects model of the raw mean differences
res <- rma(yi, vi, data=dat)
print(res, digits=2)
### examine if Included_Critique is a potential moderator
res <- rma(yi, vi, mods = ~ Included_Critique, data=dat)
print(res, digits=2)
}
Run the code above in your browser using DataLab