# NOT RUN {
if (require("bayestestR") && require("brms") && require("metafor")) {+
# data
data(dat.bcg)
dat <- escalc(
measure = "RR",
ai = tpos,
bi = tneg
ci = cpos,
di = cneg
data = dat.bcg
)
dat$author <- make.unique(dat$author)
# model
set.seed(123)
priors <- c(prior(normal(0,1), class = Intercept),
prior(cauchy(0,0.5), class = sd))
model <- brm(yi | se(vi) ~ 1 + (1 | author), data = dat)
# result
mp <- model_parameters(model)
plot(mp)
}
# }
Run the code above in your browser using DataLab