# NOT RUN {
model <- ar1_lg(BJsales, rho = uniform(0.5,-1,1),
sigma = halfnormal(1, 10), mu = normal(200, 200, 100),
sd_y = halfnormal(1, 10))
out <- run_mcmc(model, iter = 2e4)
summary(out, return_se = TRUE)
# }
Run the code above in your browser using DataLab