afni.nimh.nih.gov

It's no surprise that you didn't see any differences because you didn't incorporate the voxel-wise variable into the model. With the assumption that the VBM values do not vary within subject, you probably want to try

-bsVars "group+VBM" \

or

-bsVars "group*VBM" \

Gang