Autor: | Rainer Opgen-Rhein, Ludwig Fahrmeir, Korbinian Strimmer |
---|---|
Rok vydání: | 2005 |
Předmět: |
Genetics
education.field_of_study Markov chain Population Bayesian probability Inference Markov chain Monte Carlo Reversible-jump Markov chain Monte Carlo Biology Coalescent theory Bayes' theorem symbols.namesake symbols Quantitative Biology::Populations and Evolution education Algorithm Ecology Evolution Behavior and Systematics |
Zdroj: | BMC Evolutionary Biology. 5:6 |
ISSN: | 1471-2148 |
DOI: | 10.1186/1471-2148-5-6 |
Popis: | Coalescent theory is a general framework to model genetic variation in a population. Specifically, it allows inference about population parameters from sampled DNA sequences. However, most currently employed variants of coalescent theory only consider very simple demographic scenarios of population size changes, such as exponential growth. Here we develop a coalescent approach that allows Bayesian non-parametric estimation of the demographic history using genealogies reconstructed from sampled DNA sequences. In this framework inference and model selection is done using reversible jump Markov chain Monte Carlo (MCMC). This method is computationally efficient and overcomes the limitations of related non-parametric approaches such as the skyline plot. We validate the approach using simulated data. Subsequently, we reanalyze HIV-1 sequence data from Central Africa and Hepatitis C virus (HCV) data from Egypt. The new method provides a Bayesian procedure for non-parametric estimation of the demographic history. By construction it additionally provides confidence limits and may be used jointly with other MCMC-based coalescent approaches. |
Databáze: | OpenAIRE |
Externí odkaz: |