Hi, It is likely to be the masking that is the difference. Try fslstats -M instead of fslstats -m All the best, Mark On 17 Apr 2018, at 07:48, huyang <[log in to unmask]<mailto:[log in to unmask]>> wrote: Dear FSL experts, As for grand mean scaling, I find that there is a difference between using fslmaths -ing option and manual computation. For example: (1) fslmaths func -ing 10000 func_norm (2) Mean=$(fslstats func -m); fslmaths func -mul 10000 -div $Mean func_norm The above two methods would get unequal results and the latter one could scale grand mean to 10000. I wonder how the -ing option is implemented and why? Best, Yang Hu