Jacobian from lsqnonlin with multistart
조회 수: 3 (최근 30일)
이전 댓글 표시
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
댓글 수: 0
답변 (4개)
Carlos Henriques
2017년 4월 22일
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
댓글 수: 0
참고 항목
카테고리
Help Center 및 File Exchange에서 Global or Multiple Starting Point Search에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!