Yes, it would be ideal if we could only remove the local white matter signal.

Here is the current command:

afni_proc.py -subj_id ${SUBJID[$n]}_${MRQID[$n]} \
-out_dir $restdir/${SUBJID[$n]}/${MRQID[$n]}/output \
-script $restdir/${SUBJID[$n]}/${MRQID[$n]}/proc_script.sh \
-blocks regress \
-copy_anat $fsdir/mri/aseg+orig \
-anat_follower_ROI aaseg anat $fsdir/SUMA/aparc.a2009s+aseg_rank.nii \
-anat_follower_ROI aeseg epi $fsdir/SUMA/aparc.a2009s+aseg_rank.nii \
-anat_follower_ROI FSvent epi $fsdir/SUMA/${SUBJID[$n]}_${MRQID[$n]}_vent.nii \
-anat_follower_ROI FSWe epi $fsdir/SUMA/${SUBJID[$n]}_${MRQID[$n]}_WM.nii \
-anat_follower_erode FSWe FSvent \
-dsets $restdir/${SUBJID[$n]}/${MRQID[$n]}/rRest_psn+orig \
-tcat_remove_first_trs 4 \
-regress_anaticor \
-regress_anaticor_label FSWe \
-regress_anaticor_radius ${WM} \
-regress_motion_file $restdir/${SUBJID[$n]}/${MRQID[$n]}/Rest_mt.1D \
-regress_bandpass 0.01 0.1 \
-regress_ROI FSWe \
-regress_ROI_PC FSvent 3