Casal2 / CASAL2

Integrated Population Dynamics Model (Casal2)
https://casal2.github.io
GNU General Public License v2.0
27 stars 24 forks source link

Add MINUIT2 as a minimiser to casal2 to compute Hessian, matrices of var/cov #277

Closed mkienzle closed 2 years ago

mkienzle commented 5 years ago

MINUIT2 is the popular minimisation library developed at CERN (www.cern.ch and used in High Energy Physics to estimate model parameters. This library, as well as many other mathematical libraries written in C++, are open source, freely available and can be downloaded at www.cern.ch following the links included in the present text. This library will return the hessian matrix. Its usage is fully documented on CERN’s webpages and examples are provides with ROOT (their modular scientific software toolkit) source code, for example: https://root.cern.ch/root/html/tutorials/fit/Ifit.C.html. According to Scott Rasmussen, basic integration of this library into CASAL2 looks straight forward.

If you are interested in how MINUIT has been used in the context of fisheries stock assessment software, you can have a look at this paper and this paper and this GitHub repo. This series of papers describe an implementation of the delay difference that uses MINUIT2. A range of models ( > 100) with various processes and parameters fitted with this software to determine which hypotheses, about the working of this fishery, were most supported by the data. This mathematical library was used to estimate up to 49 parameters. The implementation was tested until the software was shown to be capable of estimating all the parameters of interest for the fishing industry to understand its dynamics.

Craig44 commented 2 years ago

We will move this to a list of future features list, rather than keeping as an open issue.