sigma.est function to estimate sigma^2 for dispersion