# Générez des données
gam_data2 <- gamSim(eg=6)
head(gam_data2)

# Faites rouler un modèle avec intercepte aléatoire
gamm_intercept <- gam(y ~ s(x0) + s(fac, bs="re"), data=gam_data2)
summary(gamm_intercept)