Affordable Access

deepdyve-link
Publisher Website

Using recursion to compute the inverse of the genomic relationship matrix

Authors
  • Misztal, I.
  • Legarra, A.
  • Aguilar, I.1, 2, 3
  • 1 Department of Animal and Dairy Science
  • 2 INRA
  • 3 Instituto Nacional de Investigación Agropecuaria
Type
Published Article
Journal
Journal of Dairy Science
Publisher
American Dairy Science Association
Publication Date
Jan 01, 2014
Accepted Date
Feb 10, 2014
Identifiers
DOI: 10.3168/jds.2013-7752
Source
Elsevier
Keywords
License
Unknown

Abstract

Computing the inverse of the genomic relationship matrix using recursion was investigated. A traditional algorithm to invert the numerator relationship matrix is based on the observation that the conditional expectation for an additive effect of 1 animal given the effects of all other animals depends on the effects of its sire and dam only, each with a coefficient of 0.5. With genomic relationships, such an expectation depends on all other genotyped animals, and the coefficients do not have any set value. For each animal, the coefficients plus the conditional variance can be called a genomic recursion. If such recursions are known, the mixed model equations can be solved without explicitly creating the inverse of the genomic relationship matrix. Several algorithms were developed to create genomic recursions. In an algorithm with sequential updates, genomic recursions are created animal by animal. That algorithm can also be used to update a known inverse of a genomic relationship matrix for additional genotypes. In an algorithm with forward updates, a newly computed recursion is immediately applied to update recursions for remaining animals. The computing costs for both algorithms depend on the sparsity pattern of the genomic recursions, but are lower or equal than for regular inversion. An algorithm for proven and young animals assumes that the genomic recursions for young animals contain coefficients only for proven animals. Such an algorithm generates exact genomic EBV in genomic BLUP and is an approximation in single-step genomic BLUP. That algorithm has a cubic cost for the number of proven animals and a linear cost for the number of young animals. The genomic recursions can provide new insight into genomic evaluation and possibly reduce costs of genetic predictions with extremely large numbers of genotypes.

Report this publication

Statistics

Seen <100 times