A likelihood-free estimator of population structure bridging admixture models and principal components analysis
Abstract We introduce a simple and computationally efficient method for fitting the admixture model of genetic population structure, called<jats:monospace>ALStructure</jats:monospace>. The strategy of<jats:monospace>ALStructure</jats:monospace>is to first estimate the low-dimensional linear subspace of the population admixture components and then search for a model within this subspace that is consistent with the admixture model’s natural probabilistic constraints. Central to this strategy is the observation that all models belonging to this constrained space of solutions are risk-minimizing and have equal likelihood, rendering any additional optimization unnecessary. The low-dimensional linear subspace is estimated through a recently introduced principal components analysis method that is appropriate for genotype data, thereby providing a solution that has both principal components and probabilistic admixture interpretations. Our approach differs fundamentally from other existing methods for estimating admixture, which aim to fit the admixture model directly by searching for parameters that maximize the likelihood function or the posterior probability. We observe that<jats:monospace>ALStructure</jats:monospace>typically outperforms existing methods both in accuracy and computational speed under a wide array of simulated and real human genotype datasets. Throughout this work we emphasize that the admixture model is a special case of a much broader class of models for which algorithms similar to<jats:monospace>ALStructure</jats:monospace>may be successfully employed..
Medienart: |
Preprint |
---|
Erscheinungsjahr: |
2022 |
---|---|
Erschienen: |
2022 |
Enthalten in: |
bioRxiv.org - (2022) vom: 11. Aug. Zur Gesamtaufnahme - year:2022 |
---|
Sprache: |
Englisch |
---|
Beteiligte Personen: |
Cabreros, Irineo [VerfasserIn] |
---|
Links: |
---|
Themen: |
---|
doi: |
10.1101/240812 |
---|
funding: |
|
---|---|
Förderinstitution / Projekttitel: |
|
PPN (Katalog-ID): |
XBI000218243 |
---|
LEADER | 01000caa a22002652 4500 | ||
---|---|---|---|
001 | XBI000218243 | ||
003 | DE-627 | ||
005 | 20230429101624.0 | ||
007 | cr uuu---uuuuu | ||
008 | 200312s2022 xx |||||o 00| ||eng c | ||
024 | 7 | |a 10.1101/240812 |2 doi | |
035 | |a (DE-627)XBI000218243 | ||
035 | |a (biorXiv)10.1101/240812 | ||
040 | |a DE-627 |b ger |c DE-627 |e rakwb | ||
041 | |a eng | ||
100 | 1 | |a Cabreros, Irineo |e verfasserin |4 aut | |
245 | 1 | 0 | |a A likelihood-free estimator of population structure bridging admixture models and principal components analysis |
264 | 1 | |c 2022 | |
336 | |a Text |b txt |2 rdacontent | ||
337 | |a Computermedien |b c |2 rdamedia | ||
338 | |a Online-Ressource |b cr |2 rdacarrier | ||
520 | |a Abstract We introduce a simple and computationally efficient method for fitting the admixture model of genetic population structure, called<jats:monospace>ALStructure</jats:monospace>. The strategy of<jats:monospace>ALStructure</jats:monospace>is to first estimate the low-dimensional linear subspace of the population admixture components and then search for a model within this subspace that is consistent with the admixture model’s natural probabilistic constraints. Central to this strategy is the observation that all models belonging to this constrained space of solutions are risk-minimizing and have equal likelihood, rendering any additional optimization unnecessary. The low-dimensional linear subspace is estimated through a recently introduced principal components analysis method that is appropriate for genotype data, thereby providing a solution that has both principal components and probabilistic admixture interpretations. Our approach differs fundamentally from other existing methods for estimating admixture, which aim to fit the admixture model directly by searching for parameters that maximize the likelihood function or the posterior probability. We observe that<jats:monospace>ALStructure</jats:monospace>typically outperforms existing methods both in accuracy and computational speed under a wide array of simulated and real human genotype datasets. Throughout this work we emphasize that the admixture model is a special case of a much broader class of models for which algorithms similar to<jats:monospace>ALStructure</jats:monospace>may be successfully employed. | ||
650 | 4 | |a Biology |7 (dpeaa)DE-84 | |
650 | 4 | |a 570 |7 (dpeaa)DE-84 | |
700 | 1 | |a Storey, John D. |e verfasserin |4 aut | |
773 | 0 | 8 | |i Enthalten in |t bioRxiv.org |g (2022) vom: 11. Aug. |
773 | 1 | 8 | |g year:2022 |g day:11 |g month:08 |
856 | 4 | 0 | |u https://doi.org/10.1534/genetics.119.302159 |z lizenzpflichtig |3 Volltext |
856 | 4 | 0 | |u http://dx.doi.org/10.1101/240812 |z kostenfrei |3 Volltext |
912 | |a GBV_XBI | ||
912 | |a SSG-OLC-PHA | ||
951 | |a AR | ||
952 | |j 2022 |b 11 |c 08 |