|
|
||||||||
Department of Animal Sciences, University of Illinois, Urbana 61801
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, the traces were based on a sample of inverse elements, resulting in up to 25 times lower cost per round. The largest model required 8 MWords (64 Mbytes) memory and used 5 h central processing unit time for one inverse or 10 h CPU time for complete run. Estimates were .67, 6.08, 4.62, and 3.47 for variance of herd by sire, permanent environment, animal, and residual, respectively.
This article has been cited by other articles:
![]() |
F. Forabosco, R. Bozzi, P. Boettcher, F. Filippini, P. Bijma, and J. A. M. Van Arendonk Relationship between profitability and type traits and derivation of economic values for reproduction and survival traits in Chianina beef cows J Anim Sci, September 1, 2005; 83(9): 2043 - 2051. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. S. Clay, B. T. McDaniel, and C. H. Brown Variances of and Correlations Among Progeny Tests for Reproductive Traits of Cows Sired by AI Bulls J Dairy Sci, July 1, 2004; 87(7): 2307 - 2313. [Abstract] [Full Text] [PDF] |
||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |