getsigma function

Compute overall residual variability from elimination phase