if (require("marginaleffects")) {
model <- lm(Sepal.Width ~ Species * Petal.Length, data = iris)
get_marginaleffects(model, trend = "Petal.Length", by = "Species")
get_marginaleffects(model, trend = "Petal.Length", by = "Petal.Length")
get_marginaleffects(model, trend = "Petal.Length", by = c("Species", "Petal.Length"))
}
Run the code above in your browser using DataLab