I often use scipy.optimize.leastsq() for my Ph.D. dissertation, but I have no idea how I can get the Jacobian score from the data returned by leastsq() . I need to know the Jacobian estimate that is used to minimize, to compare with the minimum approximation at the minimum.
Does anyone have a formula on how to get it?
It can be a little complicated when you check how, for example, the covariance matrix is ββcomputed inside leastsq()
source share