Error: Package `marginaleffects` required for this function to work.
  Please install it by running `install.packages("marginaleffects")`.

model <- stan_glm(Sepal.Width ~ Species, data = iris, refresh = 0)
contrasts <- modelbased::estimate_contrasts(model, contrast = "Species")
means <- modelbased::estimate_means(model, by = "Species")