use datastats to estimate sigma