# NOT RUN {
if (require("modelbased") && require("rstanarm")) {
model <- stan_glm(Sepal.Width ~ Species, data = iris, refresh = 0)
contrasts <- estimate_contrasts(model)
means <- estimate_means(model)
plot(contrasts, means)
}
# }
Run the code above in your browser using DataLab