Derivative-Free Restricted Maximum Likelihood Estimation in Animal Models with a Sparse Matrix Solver

¹ Agricultural Research Service, USDA, R. L. Hruska US Meat Animal Research Center, University of Nebraska, Lincoln 68583-0908

Estimation of (co)variance components by derivative-free REML requires repeated evaluation of the log-likelihood function of the data. Gaussian elimination of the augmented mixed model coefficient matrix is often used to evaluate the likelihood function, but it can be costly for animal models with large coefficient matrices. This study investigated the use of a direct sparse matrix solver to obtain the log-likelihood function. The sparse matrix package SPARSPAK was used to reorder the mixed model equations once and then repeatedly to solve the equations by Cholesky factorization to generate the terms required to calculate the likelihood. The animal model used for comparison contained 19 fixed levels, 470 maternal permanent environmental effects, and 1586 direct and 1586 maternal genetic effects, resulting in a coefficient matrix of order 3661 with .3% nonzero elements after including numerator relationships. Compared with estimation via Gaussian elimination of the unordered system, utilization of SPARSPAK required 605 and 240 times less central processing unit time on mainframes and personal computers, respectively. The SPARSPAK package also required less memory and provided solutions for all effects in the model.

Key Words: restricted maximum likelihood • derivative-free • animal model • sparse matrices

Submitted on May 3, 1991
Accepted on June 13, 1991

This article has been cited by other articles:

				Z. B. Johnson and R. A. Nugent III Heritability of body length and measures of body density and their relationship to backfat thickness and loin muscle area in swine J Anim Sci, August 1, 2003; 81(8): 1943 - 1949. [Abstract] [Full Text] [PDF]

				J. Carvalheira, E. J. Pollak, R. L. Quaas, and R. W. Blake An Autoregressive Repeatability Animal Model for Test-Day Records in Multiple Lactations J Dairy Sci, August 1, 2002; 85(8): 2040 - 2045. [Abstract] [Full Text] [PDF]