Moreover, I have one other question : I read the file "fslmeants.cc" in order to understand how the eigenvariate is computed, and I would like to know why at the line "Corr << dat * dat.t() / dat.Ncols()" the correlation matrix is divided by the number of voxels. Thank you in advance, Justine