# \dontshow{
op <- options(digits = 4)
# }
load_mgcv()
# simulate a small example
df <- data_sim("eg1", seed = 42)
# fit the GAM
m <- gam(y ~ s(x0) + s(x1) + s(x2) + s(x3), data = df, method = "REML")
# extract the estimate of the constant term
model_constant(m)
# same as coef(m)[1L]
coef(m)[1L]
# \dontshow{
options(op)
# }
Run the code above in your browser using DataLab