Jacobian from lsqnonlin with multistart
显示 更早的评论
Dear all,
For my research, I'm running multiple lsqnonlin-fits using Matlab's MultiStart function. My objective function does not supply a jacobian.
If I run a single lsqnonlin-fit (from one start point; no MultiStart), I'm able to retrieve the (approximated) jacobian using: [x,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(...) I need the jacobian to calculate a confidence interval for my estimated parameters. The run command (to run MultiStart), however, does not return a jacobian.
Is there a way to get the jacobian of the best fit of a MultiStart-run (again using lsqnonlin as the local solver)?
Thanks in advance!
Bart
回答（4 个）
Cody
2015-10-24
0 个投票
yes I'd also like to know how this can be accomplished
Carlos Henriques
2017-4-22
0 个投票
I'm doing it in this way:
you can extract the best X0 from multistart and then use it to run the lsqnonlin again.
For instance:
[params,fval,eflag,output,manymins] = run(ms,problem,20); % manymins have all the X0 used by multistart, you can take the better
[params,resnorm,residual,exitflag,fit_out,lambda,jacobian] = lsqnonlin(fun,manymins(1).X0{1},lb,ub,options);
cheers
类别
在 帮助中心 和 File Exchange 中查找有关 Global or Multiple Starting Point Search 的更多信息
2017-4-22
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
