Importance sampling for Lambda-coalescents in the infinitely many sites model

Abstract

We present and discuss new importance sampling schemes for the approximate computation of the sample probability of observed genetic types in the infinitely many sites model from population genetics. More specifically, we extend the 'classical framework', where genealogies are assumed to be governed by Kingman's coalescent, to the more general class of Lambda-coalescents and develop further Hobolth et. al.'s (2008) idea of deriving importance sampling schemes based on 'compressed genetrees'. The resulting schemes extend earlier work by Griffiths and Tavaré (1994), Stephens and Donnelly (2000), Birkner and Blath (2008) and Hobolth et. al. (2008). We conclude with a performance comparison of classical and new schemes for Beta- and Kingman coalescents.

(38 pages, 40 figures)

Cited by in corpus (8)

• Robust model selection between population growth and multiple merger coalescents
• Multi-locus data distinguishes between population growth and multiple merger coalescents
• Computational inference beyond Kingman's coalescent
• Inference and rare event simulation for stopped Markov processes via reverse-time sequential Monte Carlo
• Bayesian non-parametric inference for $\Lambda$-coalescents: consistency and a parametric method
• Asymptotic behaviour of sampling and transition probabilities in coalescent models under selection and parent dependent mutations
• Weak convergence of the scaled jump chain and number of mutations of the Kingman coalescent
• Large-sample analysis of cost functionals for inference under the coalescent