Add SE calculation using delta method approximation
Code improvements
fix missing SE computation