Computing procedures for genetic evaluation including phenotypic, full pedigree, and genomic information
- 1 September 2009
- journal article
- research article
- Published by American Dairy Science Association in Journal of Dairy Science
- Vol. 92 (9), 4648-4655
- https://doi.org/10.3168/jds.2009-2064
Abstract
Currently, genomic evaluations use multiple-step procedures, which are prone to biases and errors. A single-step procedure may be applicable when genomic predictions can be obtained by modifying the numerator relationship matrix A to H = A + A(Delta), where A(Delta) includes deviations from expected relationships. However, the traditional mixed model equations require H-1, which is usually difficult to obtain for large pedigrees. The computations with H are feasible when the mixed model equations are expressed in an alternate form that also applies for singular H and when those equations are solved by the conjugate gradient techniques. Then the only computations involving H are in the form of Aq or A(Delta q), where q is a vector. The alternative equations have a nonsymmetric left-hand side. Computing A(Delta q) is inexpensive when the number of nonzeros in A(Delta) is small, and the product Aq can be calculated efficiently in linear time using an indirect algorithm. Generalizations to more complicated models are proposed. The data included 10.2 million final scores on 6.2 million Holsteins and were analyzed by a repeatability model. Comparisons involved the regular and the alternative equations. The model for the second case included simulated A(Delta). Solutions were obtained by the preconditioned conjugate gradient algorithm, which works only with symmetric matrices, and by the bi-conjugate gradient stabilized algorithm, which also works with nonsymmetric matrices. The convergence rate associated with the nonsymmetric solvers was slightly better than that with the symmetric solver for the original equations, although the time per round was twice as much for the nonsymmetric solvers. The convergence rate associated with the alternative equations ranged from 2 times lower without A(Delta) to 3 times lower for the largest simulated A(Delta). When the information attributable to genomics can be expressed as modifications to the numerator relationship matrix, the proposed methodology may allow the upgrading of an existing evaluation to incorporate the genomic information.Keywords
Funding Information
- Holstein Association,
- EADGENE network of excellence,
- Maison de Relations Internationales (INRA, France).
This publication has 18 references indexed in Scilit:
- A relationship matrix including full pedigree and genomic informationJournal of Dairy Science, 2009
- Distribution and location of genetic effects for dairy traitsJournal of Dairy Science, 2009
- Invited Review: Reliability of genomic predictions for North American Holstein bullsJournal of Dairy Science, 2009
- Efficient Methods to Compute Genomic PredictionsJournal of Dairy Science, 2008
- Genomic-Assisted Prediction of Genetic Value With Semiparametric ProceduresGenetics, 2006
- An indirect approach to the extensive calculation of relationship coefficientsGenetics Selection Evolution, 2002
- Calculation of prediction error variances using sparse matrix methodsJournal of Animal Breeding and Genetics, 1994
- Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear SystemsSIAM Journal on Scientific and Statistical Computing, 1992
- A Simple Method for Computing the Inverse of a Numerator Relationship Matrix Used in Prediction of Breeding ValuesBiometrics, 1976
- Extension of the Gauss-Markov Theorem to Include the Estimation of Random EffectsThe Annals of Statistics, 1976