title: Genome-wide regression and prediction with the BGLR statistical package.
abstract: Many modern genomic data analyses require implementing regressions where the number of parameters (p, e.g., the number of marker effects) exceeds sample size (n). Implementing these large-p-with-small-n regressions poses several statistical and computational challenges, some of which can be confronted using Bayesian methods. This approach allows integrating various parametric and nonparametric shrinkage and variable selection procedures in a unified and consistent manner. The BGLR R-package implements a large collection of Bayesian regression models, including parametric variable selection and shrinkage methods and semiparametric procedures (Bayesian reproducing kernel Hilbert spaces regressions, RKHS). The software was originally developed for genomic applications; however, the methods implemented are useful for many nongenomic applications as well. The response can be continuous (censored or not) or categorical (either binary or ordinal). The algorithm is based on a Gibbs sampler with scalar updates and the implementation takes advantage of efficient compiled C and Fortran routines. In this article we describe the methods implemented in BGLR, present examples of the use of the package, and discuss practical issues emerging in real-data analysis.
authors.affiliation: Socio Economía Estadística e Informática, Colegio de Postgraduados 56230, México perpdgo@colpos.mx.
authors.forename: Paulino
authors.initials: P
authors.lastname: Pérez
authors.affiliation: Department of Biostatistics, Section on Statistical Genetics, University of Alabama, Birmingham, Alabama 35294.
authors.forename: Gustavo
authors.initials: G
authors.lastname: de los Campos
mesh_terms: D000465:Algorithms; D000818:Animals; D001499:Bayes Theorem; D003627:Data Interpretation, Statistical; D016678:Genome; D051379:Mice; D040641:Quantitative Trait Loci; D012044:Regression Analysis; D012984:Software; D014908:Triticum / Q000235:genetics
keywords: Bayesian methods; GenPred; R; genome-wide regression; regression; reproducing kernel Hilbert spaces regressions, RKHS; semiparametric regression; shared data resource; shrinkage; variable selection; whole-genome prediction; whole-genome regression