# generate and view data
gam_data2 <- gamSim(eg=6)
head(gam_data2)

# run random intercept model
gamm_intercept <- gam(y ~ s(x0) + s(fac, bs="re"), data=gam_data2)

# examine model output
summary(gamm_intercept)$s.table