GAMs with GCV/AIC/REML smoothness estimation and GAMMs by PQL
Routines for GAMs and other generalized ridge regression
with multiple smoothing parameter selection by GCV, REML or
UBRE/AIC. Also GAMMs by REML or PQL. Includes a gam() function.