We present a new Markov chain Monte Carlo algorithm, implemented in the 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 parallelization.
|Number of pages
|Proceedings of the Royal Society A: Mathematical Physical and Engineering Sciences
|Early online date
|12 Dec 2018
|E-pub ahead of print - 12 Dec 2018