gam_data = gamSim(eg=5)
head(gam_data)