# NOT RUN {
data(iris)
sampler <- create_sampler(Sepal.Length ~
reg(~ Petal.Length + Species, name="beta"), data=iris)
sim <- MCMCsim(sampler, burnin=100, n.iter=400)
(summary(sim))
par_names(sim)
# }
Run the code above in your browser using DataLab