Evaluation of population structure inferred by principal component analysis or the admixture model
Research output: Contribution to journal › Journal article › Research › peer-review
Documents
- Fulltext
Submitted manuscript, 1.28 MB, PDF document
Principal component analysis (PCA) is commonly used in genetics to infer and visualize population structure and admixture between populations. PCA is often interpreted in a way similar to inferred admixture proportions, where it is assumed that individuals belong to one of several possible populations or are admixed between these populations. We propose a new method to assess the statistical fit of PCA (interpreted as a model spanned by the top principal components) and to show that violations of the PCA assumptions affect the fit. Our method uses the chosen top principal components to predict the genotypes. By assessing the covariance (and the correlation) of the residuals (the differences between observed and predicted genotypes), we are able to detect violation of the model assumptions. Based on simulations and genome-wide human data, we show that our assessment of fit can be used to guide the interpretation of the data and to pinpoint individuals that are not well represented by the chosen principal components. Our method works equally on other similar models, such as the admixture model, where the mean of the data is represented by linear matrix decomposition.
Original language | English |
---|---|
Article number | iyad157 |
Journal | Genetics |
Volume | 225 |
Issue number | 2 |
Number of pages | 14 |
ISSN | 0016-6731 |
DOIs | |
Publication status | Published - 2023 |
Bibliographical note
Publisher Copyright:
© 2023 The Author(s). Published by Oxford University Press on behalf of The Genetics Society of America. All rights reserved.
- ancient DNA, PCA, population structure, residuals, statistical fit
Research areas
Number of downloads are based on statistics from Google Scholar and www.ku.dk
ID: 371925290