Run a Generalized Additive Model on all voxels of a NIfTI image within a mask and return parametric and smooth coefficients tables