Run a Generalized Additive Mixed Effects Model on all voxels of a NIfTI image return coefficients and residuals