Bridging trees for posterior inference on ancestral recombination graphs
Abstract
We present a new Markov chain Monte Carlo algorithm, implemented in software Arbores, for inferring the history of a sample of DNA sequences. Our principal innovation is a bridging procedure, previously applied only for simple stochastic processes, in which the local computations within a bridge can proceed independently of the rest of the DNA sequence, facilitating large-scale parallelisation.
- Publication:
-
Proceedings of the Royal Society of London Series A
- Pub Date:
- December 2018
- DOI:
- 10.1098/rspa.2018.0568
- arXiv:
- arXiv:1812.01345
- Bibcode:
- 2018RSPSA.47480568H
- Keywords:
-
- Statistics - Computation;
- Quantitative Biology - Populations and Evolution
- E-Print:
- 23 pages, 9 figures, accepted for publication in Proceedings of the Royal Society A