2018/12/04 by Kari Heine, Alex Beskos, Ajay Jasra +2 · 1 citation
Mathematics · Biochemistry, Genetics and Molecular Biology · #stat.CO #q-bio.PE
paper · pdf · doi:10.1098/rspa.2018.0568
23 pages, 9 figures, accepted for publication in Proceedings of the Royal Society A
arxiv created 2018/12/04 · arxiv updated 2019/03/06
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.