Use of direct and iterative solvers for estimation of SNP effects in genome-wide selection

The aim of this study was to compare iterative and direct solvers for estimation of marker effects in genomic selection. One iterative and two direct methods were used: Gauss-Seidel with Residual Update, Cholesky Decomposition and Gentleman-Givens rotations. For resembling different scenarios with r...

Full description

Bibliographic Details
Main Authors: Eduardo da Cruz Gouveia Pimentel, Mehdi Sargolzaei, Henner Simianer, Flávio Schramm Schenkel, Zengting Liu, Luiz Alberto Fries, Sandra Aidar de Queiroz
Format: Article
Language:English
Published: Sociedade Brasileira de Genética 2010-01-01
Series:Genetics and Molecular Biology
Subjects:
Online Access:http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1415-47572010000100033
_version_ 1818511229152395264
author Eduardo da Cruz Gouveia Pimentel
Mehdi Sargolzaei
Henner Simianer
Flávio Schramm Schenkel
Zengting Liu
Luiz Alberto Fries
Sandra Aidar de Queiroz
author_facet Eduardo da Cruz Gouveia Pimentel
Mehdi Sargolzaei
Henner Simianer
Flávio Schramm Schenkel
Zengting Liu
Luiz Alberto Fries
Sandra Aidar de Queiroz
author_sort Eduardo da Cruz Gouveia Pimentel
collection DOAJ
description The aim of this study was to compare iterative and direct solvers for estimation of marker effects in genomic selection. One iterative and two direct methods were used: Gauss-Seidel with Residual Update, Cholesky Decomposition and Gentleman-Givens rotations. For resembling different scenarios with respect to number of markers and of genotyped animals, a simulated data set divided into 25 subsets was used. Number of markers ranged from 1,200 to 5,925 and number of animals ranged from 1,200 to 5,865. Methods were also applied to real data comprising 3081 individuals genotyped for 45181 SNPs. Results from simulated data showed that the iterative solver was substantially faster than direct methods for larger numbers of markers. Use of a direct solver may allow for computing (co)variances of SNP effects. When applied to real data, performance of the iterative method varied substantially, depending on the level of ill-conditioning of the coefficient matrix. From results with real data, Gentleman-Givens rotations would be the method of choice in this particular application as it provided an exact solution within a fairly reasonable time frame (less than two hours). It would indeed be the preferred method whenever computer resources allow its use.
first_indexed 2024-12-10T23:30:32Z
format Article
id doaj.art-5777983032564472a367ed5e16974ac8
institution Directory Open Access Journal
issn 1415-4757
1678-4685
language English
last_indexed 2024-12-10T23:30:32Z
publishDate 2010-01-01
publisher Sociedade Brasileira de Genética
record_format Article
series Genetics and Molecular Biology
spelling doaj.art-5777983032564472a367ed5e16974ac82022-12-22T01:29:25ZengSociedade Brasileira de GenéticaGenetics and Molecular Biology1415-47571678-46852010-01-01331198204Use of direct and iterative solvers for estimation of SNP effects in genome-wide selectionEduardo da Cruz Gouveia PimentelMehdi SargolzaeiHenner SimianerFlávio Schramm SchenkelZengting LiuLuiz Alberto FriesSandra Aidar de QueirozThe aim of this study was to compare iterative and direct solvers for estimation of marker effects in genomic selection. One iterative and two direct methods were used: Gauss-Seidel with Residual Update, Cholesky Decomposition and Gentleman-Givens rotations. For resembling different scenarios with respect to number of markers and of genotyped animals, a simulated data set divided into 25 subsets was used. Number of markers ranged from 1,200 to 5,925 and number of animals ranged from 1,200 to 5,865. Methods were also applied to real data comprising 3081 individuals genotyped for 45181 SNPs. Results from simulated data showed that the iterative solver was substantially faster than direct methods for larger numbers of markers. Use of a direct solver may allow for computing (co)variances of SNP effects. When applied to real data, performance of the iterative method varied substantially, depending on the level of ill-conditioning of the coefficient matrix. From results with real data, Gentleman-Givens rotations would be the method of choice in this particular application as it provided an exact solution within a fairly reasonable time frame (less than two hours). It would indeed be the preferred method whenever computer resources allow its use.http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1415-47572010000100033breeding valuegenomic selectionmixed model equationsnumerical method
spellingShingle Eduardo da Cruz Gouveia Pimentel
Mehdi Sargolzaei
Henner Simianer
Flávio Schramm Schenkel
Zengting Liu
Luiz Alberto Fries
Sandra Aidar de Queiroz
Use of direct and iterative solvers for estimation of SNP effects in genome-wide selection
Genetics and Molecular Biology
breeding value
genomic selection
mixed model equations
numerical method
title Use of direct and iterative solvers for estimation of SNP effects in genome-wide selection
title_full Use of direct and iterative solvers for estimation of SNP effects in genome-wide selection
title_fullStr Use of direct and iterative solvers for estimation of SNP effects in genome-wide selection
title_full_unstemmed Use of direct and iterative solvers for estimation of SNP effects in genome-wide selection
title_short Use of direct and iterative solvers for estimation of SNP effects in genome-wide selection
title_sort use of direct and iterative solvers for estimation of snp effects in genome wide selection
topic breeding value
genomic selection
mixed model equations
numerical method
url http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1415-47572010000100033
work_keys_str_mv AT eduardodacruzgouveiapimentel useofdirectanditerativesolversforestimationofsnpeffectsingenomewideselection
AT mehdisargolzaei useofdirectanditerativesolversforestimationofsnpeffectsingenomewideselection
AT hennersimianer useofdirectanditerativesolversforestimationofsnpeffectsingenomewideselection
AT flavioschrammschenkel useofdirectanditerativesolversforestimationofsnpeffectsingenomewideselection
AT zengtingliu useofdirectanditerativesolversforestimationofsnpeffectsingenomewideselection
AT luizalbertofries useofdirectanditerativesolversforestimationofsnpeffectsingenomewideselection
AT sandraaidardequeiroz useofdirectanditerativesolversforestimationofsnpeffectsingenomewideselection