# NOT RUN {
library(modelbased)
model <- lm(Sepal.Length ~ Species + Petal.Width, data = iris)
# By default, 'levels' is set to "Species"
model_emmeans(model)
# One can estimate marginal means at several values of a 'modulate' variable
model_emmeans(model, modulate = "Petal.Width", length = 3)
# }
Run the code above in your browser using DataLab