Hi, When using glm_fit with an input volume that has many voxels with a 0 value (due to a previous thresholding for quality reasons), and the voxel removes might be distributed randomly across our population, is there a way not to take into account those voxels in the glm_fit model? I saw the --pune flag, but I am not sure if it is the one to use. Thanks Celine