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