data(data_mod_cat_2w)
lm_out <- lm(y ~ x*gp + x*site + c1 + c2, data_mod_cat_2w)
out <- cond_effects(
wlevels = c("site", "gp"),
x = "x",
fit = lm_out
)
out
plot(out, facet_grid_cols = "site", graph_type = "tumble")
Run the code above in your browser using DataLab