  Robust and efficient optimization methods for variance component estimationusing Restricted Maximum Likelihood (REML) models for genetic mapping ofquantitative traits are considered. We show that the standard Newton-AI schememay fail when the optimum is located at one of the constraint boundaries, andwe introduce different approaches to remedy this by taking the constraints intoaccount. We approximate the Hessian of the objective function using the averageinformation matrix and also by using an inverse BFGS formula. The robustnessand efficiency is evaluated for problems derived from two experimental datafrom the same animal populations.
We approximate the Hessian of the objective function using the average information matrix and also by using an inverse BFGS formula.