getsigmas function

Estimate individual-level residual error from the elimination phase