# NOT RUN {
if (require("bayestestR") && require("rstanarm")) {
set.seed(123)
m <- stan_glm(Sepal.Length ~ Petal.Width * Species, data = iris, refresh = 0)
result <- point_estimate(m, centrality = "median")
result
plot(result)
}
# }
Run the code above in your browser using DataLab