Jacobian from lsqnonlin with multistart
5 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaires
Réponses (4)
Carlos Henriques
le 22 Avr 2017
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 commentaires
Voir également
Catégories
En savoir plus sur Global or Multiple Starting Point Search dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!