method

change the fitting method

Set the minimization method.

Syntax: method<algorithm>[<# of trials/evaluations>[<critical delta>] [<method-specific options>]]

where <algorithm> is the method in use and the other arguments are control values for the minimization. Their meanings are explained under the individual methods. The migrad and simplex methods are taken from the CERN Minuit2 package, with documentation located at
http://seal.web.cern.ch/seal/MathLibs/Minuit2/html/index.html. If either of these are used, then the error command will use the Minuit2 minos method to find the confidence regions.

leven

method leven [<# of eval>[<crit delta>] [<crit beta>]] [delay |nodelay]

The default XSPEC minimization method using the modified Levenberg-Marquardt algorithm based on the CURFIT routine from Bevington. <# of eval> is the number of trial vectors before the user is prompted to say whether they want to continue fitting. <crit delta> is the convergence criterion, which is the (absolute, not fractional) difference in fit statistic between successive iterations, less than which the fit is determined to have converged.

<crit beta> refers to the |beta|/N value reported during a fit. This is the norm of the vector of derivatives of the statistic with respect to the parameters divided by the number of parameters. At the best fit this should be zero, and so provides another measure of how well the fit is converging. When this is set to a positive value, it will provide another fit stopping criterion in addition to that of the <crit delta> setting.

Including the string delay as an argument turns on delayed gratification. It is turned off by nodelay. Delayed gratification modifies the way the damping parameter is set and has been shown in many cases to speed up convergence. The default is nodelay.

<# of eval>, <crit delta>, <crit beta>, delay, and nodelay may also be set through the fit command.

This method requires an estimate of the second derivative of the statistic with respect to the parameters. By default, XSPEC calculates these using an analytic expression which assumes that partial 2nd derivatives of the model with respect to its parameters may be ignored. This may be changed by setting the USE_CHAIN_RULE flag to false in the user's startup Xspec.init initialization file. XSPEC will then calculate all second derivatives numerically, which can be noticeably slower.

migrad

method migrad [<# of eval>]

The Minuit2 migrad method. <# of eval> is the number of function evaluations to perform before giving up. Migrad uses an internal convergence criterion.

The current version of Minuit2 included is that from ROOT v5.34. Documentation on Minuit2 can be found at http://seal.web.cern.ch/seal/MathLibs/Minuit2/html/.

If migrad is not working well try experimenting with different hard and soft limits on parameters.

simplex

method simplex [<# of evaluations>]

The Minuit2 simplex method. <# of evaluations>is the number of function evaluations to perform before giving up. Simplex uses an internal convergence criterion. This method is included for historical interest and is almost always outperformed by migrad.


HEASARC Home | Observatories | Archive | Calibration | Software | Tools | Students/Teachers/Public

Last modified: Friday, 23-Aug-2024 13:20:40 EDT