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