Fecha/hora: 16/06/2024 10:48:21

Biblioteca(s): | INIA Las Brujas.
Fecha corriente: | 09/11/2017
Fecha de la última actualización: | 25/11/2019
Autoria: | MASUDA, Y; MISZTAL, I.; LEGARRA, A.; TSURUTA, S.; LOURENCO, D.A.L.; FRAGOMENI, B.O.; AGUILAR, I.
Afiliación: | Y. MASUDA, Department of Animal and Dairy Science, University of Georgia; I. MISZTAL, Department of Animal and Dairy Science, University of Georgia; A. LEGARRA, INRA (Institut National de la Recherche Agronomique); S. TSURUTA, Department of Animal and Dairy Science, University of Georgia; D.A.L. LOURENCO, Department of Animal and Dairy Science, University of Georgia; B.O. FRAGOMENI, Department of Animal and Dairy Science, University of Georgia; IGNACIO AGUILAR GARCIA, INIA (Instituto Nacional de Investigación Agropecuaria), Uruguay.
Título: | Technical note: Avoiding the direct inversion of the numerator relationship matrix for genotyped animals in single-step genomic best linear unbiased prediction solved with the preconditioned conjugate gradient.
Año de publicación: | 2017
Fuente/Imprenta: | Journal of Animal Science, 2017, v. 95(1): 49-52.
DOI: | 10.2527/jas.2016.0699
Idioma: | Inglés

Notas:

Article history: Received: July 05, 2016; Accepted: Aug 16, 2016; Published: February 2, 2017.
This research was partially funded by the United States Department of Agriculture?s National Institute of Food and Agriculture (Agriculture and Food Research Initiative competitive grant 2015-67015-22936).

Contenido:

ABSTRACT.
This paper evaluates an efficient implementation to multiply the inverse of a numerator relationship matrix for genotyped animals () by a vector (q). The computation is required for solving mixed model equations in single-step genomic BLUP (ssGBLUP) with the preconditioned conjugate gradient (PCG). The inverse can be decomposed into sparse matrices that are blocks of the sparse inverse of a numerator relationship matrix (A&#8722;1) including genotyped animals and their ancestors. The elements of A&#8722;1 were rapidly calculated with the Henderson?s rule and stored as sparse matrices in memory. Implementation of  was by a series of sparse matrix?vector multiplications. Diagonal elements of , which were required as preconditioners in PCG, were approximated with a Monte Carlo method using 1,000 samples. The efficient implementation of was compared with explicit inversion of A22 with 3 data sets including about 15,000, 81,000, and 570,000 genotyped animals selected from populations with 213,000, 8.2 million, and 10.7 million pedigree animals, respectively. The explicit inversion required 1.8 GB, 49 GB, and 2,415 GB (estimated) of memory, respectively, and 42 s, 56 min, and 13.5 d (estimated), respectively, for the computations. The efficient implementation required <1 MB, 2.9 GB, and 2.3 GB of memory, respectively, and <1 sec, 3 min, and 5 min, respectively, for setting up. Only <1 sec was required for the multiplication in each PCG iteration for any data sets. When the equations in ssGBLUP are solved with the PCG algorithm,  is no longer a limiting factor in the computations.

Copyright © 2016. American Society of Animal Science.

Palabras-Clave: | COMPUTATION; GENOMIC SELECTION; INVERSION; NUMERATOR RELATIONSHIP MATRIX; PRECONDITIONED CONJUGATE GRADIENT; SPARSE MATRIX.
Categoría del tema: | --
Registro original: | INIA Las Brujas (LB)
Ejemplares: | LB - (101491-1) - PXI - AP - DD - PP/JAS/2017 -  -  - DIS