# NOT RUN {
# load data (Source: ?mgcv::gam)
library(mgcv)
n <- 200
sig <- 2
dat <- gamSim(1, n = n, scale = sig, verbose = FALSE)
dat$x4 <- as.factor(c(
rep("A", 50), rep("B", 50), rep("C", 50),
rep("D", 50)
))
fit_gam <- gam(y ~ s(x0) + s(I(x1^2)) + s(x2) +
offset(x3) + x4, data = dat) # fit model
tmp <- plot(fit_gam, pages = 1) # plot output
tmp <- no_plot(fit_gam) # no plot output
# }
Run the code above in your browser using DataLab