Hi Rick,

I generated the script using afni_proc.py but I have made changes to it. For instance, doing bandpassing is one of them. Instead of generating sinusoidal regressors using 1dBport, I do bandpassing directly using 3dTproject:
3dTproject -polort 0 -input $output_dir/pb05.$SBJ.r*.scale+orig.HEAD \
-censor $output_dir/censor_${SBJ}_combined_2.1D -cenmode ZERO \
-dsort $output_dir/Local_FSWe_rall+orig \
-ort $output_dir/X.nocensor.xmat.1D \
-passband 0.01 0.1 \
-prefix $output_dir/errts.$SBJ.fanaticor -overwrite

I don't remember where but I remember reading one of the help files, where it was highly recommended to bandpass using 3dTproject, so I made the change.

As for the white-matter mask erosion, I should have mentioned that -2 erosion was done before the mask was resampled to the epi image. So the erosion was performed at 1x1x1 and then the mask was resampled to 4x4x4 during the alignment step. Perhaps I can use a -3 erosion? But, even right, the white-matter mask seems okay, doesn't it? There does not seem to be any overlap between the mask and the gray matter?

Finally, if ANATICOR is not the answer for 4x4x4 resolution, what would you recommend is my best option for noise removal?

Thanks again for all the help guys!! Really appreciate it.

Kausar