updated lmfit - na.rm, and sd_e calculation