Bug fix float overflow in estimate_sigma

Hi I noticed, that in HCP data the sigma was all zeros. This change fixes the problem. GGT!