# NOT RUN {
library(rstanarm)
library(bayestestR)
library(see)
library(ggplot2)
model <- stan_glm(
Sepal.Length ~ Petal.Width + Species + Sepal.Width,
data = iris,
chains = 2, iter = 200
)
result <- hdi(model, ci = c(0.5, 0.75, 0.9, 0.95))
data <- data_plot(result, data = model)
p <- data %>%
ggplot(aes(x = x, y = y, height = height, group = y, fill = fill)) +
ggridges::geom_ridgeline_gradient()
p
p + add_plot_attributes(data)
# }
Run the code above in your browser using DataLab