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