Academic Journals Database
Disseminating quality controlled scientific knowledge

Auto-validating von Neumann rejection sampling from small phylogenetic tree spaces

Author(s): Sainudiin Raazesh | York Thomas

Journal: Algorithms for Molecular Biology
ISSN 1748-7188

Volume: 4;
Issue: 1;
Start page: 1;
Date: 2009;
Original page

Abstract Background In phylogenetic inference one is interested in obtaining samples from the posterior distribution over the tree space on the basis of some observed DNA sequence data. One of the simplest sampling methods is the rejection sampler due to von Neumann. Here we introduce an auto-validating version of the rejection sampler, via interval analysis, to rigorously draw samples from posterior distributions over small phylogenetic tree spaces. Results The posterior samples from the auto-validating sampler are used to rigorously (i) estimate posterior probabilities for different rooted topologies based on mitochondrial DNA from human, chimpanzee and gorilla, (ii) conduct a non-parametric test of rate variation between protein-coding and tRNA-coding sites from three primates and (iii) obtain a posterior estimate of the human-neanderthal divergence time. Conclusion This solves the open problem of rigorously drawing independent and identically distributed samples from the posterior distribution over rooted and unrooted small tree spaces (3 or 4 taxa) based on any multiply-aligned sequence data.
Why do you need a reservation system?      Save time & money - Smart Internet Solutions