Estimate generalized additive mixed models via a version of function gamm() from 'mgcv', using 'lme4' for estimation.