sigma_hat function

Statistical estimator for sigma