JDS
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Journal of Dairy Science Vol. 73 No. 1 163-172
© 1990 by American Dairy Science Association ®
This Article
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Misztal, I.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Misztal, I.

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

Ignacy Misztal

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:


Home page
J ANIM SCIHome page
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]


Home page
J DAIRY SCIHome page
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
Copyright © 1990 by the American Dairy Science Association ®.