# 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