# copy data into 'dat' and examine data
dat <- dat.hannum2020
dat
if (FALSE) {
# load metafor package
library(metafor)
# compute effect size
dat <- escalc(measure="PR", xi=xi, ni=ni, data=dat)
# split data into objective and subjective datasets
dat_split <- split(dat, dat$objectivity)
dat_objective <- dat_split[["Objective"]]
dat_subjective <- dat_split[["Subjective"]]
# random-effects model all studies
res_all <- rma(yi, vi, data=dat)
print(res_all, digits=2)
# random-effects model objective
res_objective <- rma(yi, vi, data=dat_objective)
print(res_objective, digits=2)
# random-effects model subjective
res_subjective <- rma(yi, vi, data=dat_subjective)
print(res_subjective, digits=2)
}
Run the code above in your browser using DataLab