Run a Linear Mixed Effects Model on all voxels of a NIfTI image and return parametric coefficients and residuals