## Setting up dummy sampling for fast execution of example
## Please use 4 chains and 20x more warmup & iter in practice
.user_mc_options <- options(RBesT.MC.warmup=50, RBesT.MC.iter=100,
RBesT.MC.chains=2, RBesT.MC.thin=1)
set.seed(546346)
map_crohn <- gMAP(cbind(y, y.se) ~ 1 | study,
family = gaussian,
data = transform(crohn, y.se = 88 / sqrt(n)),
weights = n,
tau.dist = "HalfNormal", tau.prior = 44,
beta.prior = cbind(0, 88)
)
## Recover user set sampling defaults
options(.user_mc_options)
Run the code above in your browser using DataLab