Bayesian Inference of Selection in the Wright-Fisher Diffusion Model

by   Jeffrey J. Gory, et al.

The increasing availability of population-level allele frequency data across one or more related populations necessitates the development of methods that can efficiently estimate population genetics parameters, such as the strength of selection acting on the population(s), from such data. Existing methods for this problem in the setting of the Wright-Fisher diffusion model are primarily likelihood-based, and rely on numerical approximation for likelihood computation and on bootstrapping for assessment of variability in the resulting estimates, requiring extensive computation. Recent work (Jenkins and Spano, 2015) has provided a method for obtaining exact samples from general Wright-Fisher diffusion processes, enabling the development of methods for Bayesian estimation in this setting. We develop and implement a Bayesian method for estimating the strength of selection based on the Wright-Fisher diffusion for data sampled at a single time point. The method utilizes the work of Jenkins and Spano (2015) to develop a Markov chain Monte Carlo algorithm to draw samples from the joint posterior distribution of the selection coefficient and the allele frequencies. We demonstrate that when assumptions about the initial allele frequencies are accurate the method performs well for both simulated data and for an empirical data set on hypoxia in flies (Zhou et al. 2011), where we find evidence for strong positive selection in a region of chromosome 2L previously identified by Ronen et al. (2013). We discuss possible extensions of our method to the more general settings commonly encountered in practice, highlighting the advantages of Bayesian approaches to inference in this setting.


page 1

page 2

page 3

page 4


Practical bounds on the error of Bayesian posterior approximations: A nonasymptotic approach

Bayesian inference typically requires the computation of an approximatio...

Statistic Selection and MCMC for Differentially Private Bayesian Estimation

This paper concerns differentially private Bayesian estimation of the pa...

Self-contained Beta-with-Spikes Approximation for Inference Under a Wright-Fisher Model

We construct a reliable estimation of evolutionary parameters within the...

Rare event ABC-SMC^2

Approximate Bayesian computation (ABC) is a well-established family of M...

EWF : simulating exact paths of the Wright–Fisher diffusion

The Wright–Fisher diffusion is important in population genetics in model...

Convergence of Bayesian Estimators for Diffusions in Genetics

A number of discrete time, finite population size models in genetics des...

Distributed Computation for Marginal Likelihood based Model Choice

We propose a general method for distributed Bayesian model choice, using...

Please sign up or login with your details

Forgot password? Click here to reset