Solving Large Test-Day Models by Iteration on Data and Preconditioned Conjugate Gradient

¹ Animal Production Research, Agricultural Research Centre, FIN-31600 Jokioinen, Finland

A preconditioned conjugate gradient method was implemented into an iteration on a program for data estimation of breeding values, and its convergence characteristics were studied. An algorithm was used as a reference in which one fixed effect was solved by Gauss-Seidel method, and other effects were solved by a second-order Jacobi method. Implementation of the preconditioned conjugate gradient required storing four vectors (size equal to number of unknowns in the mixed model equations) in random access memory and reading the data at each round of iteration. The preconditioner comprised diagonal blocks of the coefficient matrix. Comparison of algorithms was based on solutions of mixed model equations obtained by a single-trait animal model and a single-trait, random regression test-day model. Data sets for both models used milk yield records of primiparous Finnish dairy cows. Animal model data comprised 665,629 lactation milk yields and random regression test-day model data of 6,732,765 test-day milk yields. Both models included pedigree information of 1,099,622 animals. The animal model {random regression test-day model} required 122 {305} rounds of iteration to converge with the reference algorithm, but only 88 {149} were required with the preconditioned conjugate gradient. To solve the random regression test-day model with the preconditioned conjugate gradient required 237 megabytes of random access memory and took 14% of the computation time needed by the reference algorithm.

Key Words: iteration on data • preconditioned conjugate gradient • test-day model

Submitted on January 5, 1999
Accepted on June 4, 1999

This article has been cited by other articles:

				J. Bohmanova, F. Miglior, and J. Jamrozik Use of test-day records beyond three hundred five days for estimation of three hundred five-day breeding values for production traits and somatic cell score of Canadian Holsteins J Dairy Sci, October 1, 2009; 92(10): 5314 - 5325. [Abstract] [Full Text] [PDF]

				Z. Liu, J. Jaitner, F. Reinhardt, E. Pasman, S. Rensing, and R. Reents Genetic Evaluation of Fertility Traits of Dairy Cattle Using a Multiple-Trait Animal Model J Dairy Sci, November 1, 2008; 91(11): 4333 - 4343. [Abstract] [Full Text] [PDF]

				J. Bohmanova, F. Miglior, J. Jamrozik, I. Misztal, and P. G. Sullivan Comparison of Random Regression Models with Legendre Polynomials and Linear Splines for Production Traits and Somatic Cell Score of Canadian Holstein Cows J Dairy Sci, September 1, 2008; 91(9): 3627 - 3638. [Abstract] [Full Text] [PDF]

				A. Legarra and I. Misztal Technical Note: Computing Strategies in Genome-Wide Selection J Dairy Sci, January 1, 2008; 91(1): 360 - 366. [Abstract] [Full Text] [PDF]

				K. Meyer and M. Kirkpatrick Up hill, down dale: quantitative genetics of curvaceous traits Phil Trans R Soc B, July 29, 2005; 360(1459): 1443 - 1455. [Abstract] [Full Text] [PDF]

				P. R. C. Nobre, I. Misztal, S. Tsuruta, J. K. Bertrand, L. O. C. Silva, and P. S. Lopes Genetic evaluation of growth in Nellore cattle by multiple-trait and random regression models J Anim Sci, April 1, 2003; 81(4): 927 - 932. [Abstract] [Full Text] [PDF]