ga_model function

Fit a generalised additive model (gam)