Rxivist logo

Rxivist combines preprints from bioRxiv with data from Twitter to help you find the papers being discussed in your field. Currently indexing 55,138 bioRxiv papers from 254,399 authors.

Variational Inference of Population Structure in Large SNP Datasets

By Anil Raj, Matthew Stephens, Jonathan K Pritchard

Posted 02 Dec 2013
bioRxiv DOI: 10.1101/001073 (published DOI: 10.1534/genetics.114.164350)

Tools for estimating population structure from genetic data are now used in a wide variety of applications in population genetics. However, inferring population structure in large modern data sets imposes severe computational challenges. Here, we develop efficient algorithms for approximate inference of the model underlying the STRUCTURE program using a variational Bayesian framework. Variational methods pose the problem of computing relevant posterior distributions as an optimization problem, allowing us to build on recent advances in optimization theory to develop fast inference tools. In addition, we propose useful heuristic scores to identify the number of populations represented in a dataset and a new hierarchical prior to detect weak population structure in the data. We test the variational algorithms on simulated data, and illustrate using genotype data from the CEPH-Human Genome Diversity Panel. The variational algorithms are almost two orders of magnitude faster than STRUCTURE and achieve accuracies comparable to those of ADMIXTURE. Furthermore, our results show that the heuristic scores for choosing model complexity provide a reasonable range of values for the number of populations represented in the data, with minimal bias towards detecting structure when it is very weak. Our algorithm, fastSTRUCTURE, is freely available online at http://pritchardlab.stanford.edu/structure.html.

Download data

  • Downloaded 5,820 times
  • Download rankings, all-time:
    • Site-wide: 250 out of 55,165
    • In genetics: 27 out of 3,190
  • Year to date:
    • Site-wide: 7,967 out of 55,165
  • Since beginning of last month:
    • Site-wide: 10,759 out of 55,165

Altmetric data


Downloads over time

Distribution of downloads per paper, site-wide


Sign up for the Rxivist weekly newsletter! (Click here for more details.)


News