library(mgcv)
set.seed(10)
data <- gamSim(4)
model <- gam(y ~ fac + s(x2) + s(x2, by = fac) + s(x0), data = data)
get_smooths_difference(model, x2, list(fac = c("1", "2")))
# For details, see vignette
if (FALSE) {
vignette("plot-smooths", package = "tidymv")
}
Run the code above in your browser using DataLab