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.

23 pages, 9 figures, accepted for publication in Proceedings of the Royal Society A