Restricted Maximum Likelihood Estimation of Variance Components in Animal Model Using Sparse Matrix Inversion and a Supercomputer

The estimation of variance components by EM-type REML algorithms requires repeated inversion of the coefficient matrix of the mixed model equations. In the animal model, the rank of the coefficient matrix is usually larger than the number of records. If inversion is done by dense matrix algorithms, equations are limited to fewer than 2000. This study investigated the inversion of large coefficient matrices by direct sparse matrix solvers in analyzing conformation final scores of selected Holstein herds. The largest model contained 36,771 animal effects, 30,739 permanent environmental effects, 6102 herd X sire interactions, 532 herd-time classes, and 55 unknown parent groups, resulting in a coefficient matrix of order 74,199. Inversion was by sparse matrix package SMPAK, partly vectorized for the Cray-2 supercomputer. The convergence rate was accelerated by modified Aitken's extrapolation. For most rounds,