Switched over to sigma from var (internal code). Newton method gets std errors from Hessian.