To solve difficult optimization problems, we have developed a computer
program called GENOUD (GENetic Optimization Using Derivatives) that
combines evolutionary algorithm methods with a derivative-based,
quasi-Newton method. GENOUD can work even when the most often used
optimization methods completely fail. The objective function for a
nonlinear model may not be globally concave, making it difficult for
gradient-based optimization methods to find any optimum at all. Multiple
local optima may exist so there is no guarantee that gradient-based
methods will converge to the global optimum. We discuss the theoretical
basis for expecting GENOUD to have a high probability of finding global
optima. We conduct Monte Carlo experiments using scalar Normal mixture
densities to illustrate this capability. We also use a real-data example
(the four-dimensional Hopf model), which has many parameters and multiple
local optima, to compare the performance of GENOUD to that of the
Gauss-Newton algorithm in SAS's PROC MODEL.