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