Dataset Information

Fast hierarchical Bayesian analysis of population structure.

ABSTRACT: We present fastbaps, a fast solution to the genetic clustering problem. Fastbaps rapidly identifies an approximate fit to a Dirichlet process mixture model (DPM) for clustering multilocus genotype data. Our efficient model-based clustering approach is able to cluster datasets 10-100 times larger than the existing model-based methods, which we demonstrate by analyzing an alignment of over 110 000 sequences of HIV-1 pol genes. We also provide a method for rapidly partitioning an existing hierarchy in order to maximize the DPM model marginal likelihood, allowing us to split phylogenetic trees into clades and subclades using a population genomic model. Extensive tests on simulated data as well as a diverse set of real bacterial and viral datasets show that fastbaps provides comparable or improved solutions to previous model-based methods, while being significantly faster. The method is made freely available under an open source MIT licence as an easy to use R package at https://github.com/gtonkinhill/fastbaps.

SUBMITTER: Tonkin-Hill G

PROVIDER: S-EPMC6582336 | biostudies-literature | 2019 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Fast hierarchical Bayesian analysis of population structure.

Tonkin-Hill Gerry G Lees John A JA Bentley Stephen D SD Frost Simon D W SDW Corander Jukka J

Nucleic acids research 20190601 11

We present fastbaps, a fast solution to the genetic clustering problem. Fastbaps rapidly identifies an approximate fit to a Dirichlet process mixture model (DPM) for clustering multilocus genotype data. Our efficient model-based clustering approach is able to cluster datasets 10-100 times larger than the existing model-based methods, which we demonstrate by analyzing an alignment of over 110 000 sequences of HIV-1 pol genes. We also provide a method for rapidly partitioning an existing hierarchy ...[more]

PMID: 31076776

Similar Datasets

Project description:When assessing the human risks due to exposure to environmental chemicals, traditional dose-response analyses are not straightforward when there are numerous high-quality epidemiological studies of priority cancer and non-cancer health outcomes. Given this wealth of information, selecting a single "best" study on which to base dose-response analyses is difficult and would potentially ignore much of the available data. Therefore, systematic approaches are necessary for the analysis of these rich databases. Examples are meta-analysis (and further, meta-regression), which are well established methods that consider and incorporate information from multiple studies into the estimation of risks due to exposure to environmental contaminants. In this paper, we propose a hierarchical, Bayesian meta-analysis approach for the dose-response analysis of multiple epidemiological studies. This paper is the second of two papers detailing this approach; the first covered "pre-analysis" steps necessary to prepare the data for dose-response modeling. This paper focuses on the hierarchical Bayesian approach to dose-response modeling and extrapolation of risk to populations of interest using the association between bladder cancer and oral inorganic arsenic (iAs) exposure as an illustrative case study. In particular, this paper addresses the modeling of both case-control and cohort studies with a flexible, logistic model in a hierarchical Bayesian framework that estimates study-specific slopes, as well as a pooled slope across all studies. This approach is akin to a random effects model in which no assumption is made a priori that there is a single, common slope for all included studies. Further, this paper also details extrapolation of the estimates of logistic slope to extra risk in a target population using a lifetable analysis and basic assumptions about background iAs exposure levels. In this case, the target population was the general United States population and information on all-cause mortality and incidence and mortality from bladder cancer was used to perform the lifetable analysis. The methods herein were developed for general use in investigating the association between any pollutant and observed health-effects in epidemiological studies. In order to demonstrate these methods, inorganic arsenic was chosen as a case study given the large epidemiological database that exists for this contaminant.

Dataset Information

Fast hierarchical Bayesian analysis of population structure.

Publications

Fast hierarchical Bayesian analysis of population structure.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets