minor bug fix in parameter estimation, I was dividing theta by two in the correction function proper but not in the second stage of parameter fitting