Setup and Model Fitting

library(modelbased)
library(rstanarm)
library(ggplot2)
library(see)

theme_set(theme_modern())

model <- stan_glm(Sepal.Width ~ Species, data = iris, refresh = 0)

Pairwise Contrasts

contrasts <- estimate_contrasts(model)
means <- estimate_means(model)

plot(contrasts, means)