The Expected Sample Allele Frequencies from Populations of Changing Size via Orthogonal Polynomials

04/05/2023
by   Lynette Caitlin Mikula, et al.
0

In this article, discrete and stochastic changes in (effective) population size are incorporated into the spectral representation of a biallelic diffusion process for drift and small mutation rates. A forward algorithm inspired by Hidden-Markov-Model (HMM) literature is used to compute exact sample allele frequency spectra for three demographic scenarios: single changes in (effective) population size, boom-bust dynamics, and stochastic fluctuations in (effective) population size. An approach for fully agnostic demographic inference from these sample allele spectra is explored, and sufficient statistics for step-wise changes in population size are found. Further, convergence behaviours of the polymorphic sample spectra for population size changes on different time scales are examined and discussed within the context of inference of the effective population size. Joint visual assessment of the sample spectra and the temporal coefficients of the spectral decomposition of the forward diffusion process is found to be important in determining departure from equilibrium. Stochastic changes in (effective) population size are shown to shape sample spectra particularly strongly.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
12/17/2021

Forward-backward algorithms with a biallelic mutation-drift model: Orthogonal polynomials, and a coalescent/urn-model based approach

Inference of the marginal likelihood of sample allele configurations usi...
research
12/13/2017

Geometry of the sample frequency spectrum and the perils of demographic inference

The sample frequency spectrum (SFS), which describes the distribution of...
research
01/10/2020

Convergence of Bayesian Estimators for Diffusions in Genetics

A number of discrete time, finite population size models in genetics des...
research
06/15/2020

A demographic microsimulation model with an integrated household alignment method

Many dynamic microsimulation models have shown their ability to reasonab...
research
08/16/2021

Agentization of Two Population-Driven Models of Mathematical Biology

Single species population models and discrete stochastic gene frequency ...
research
08/11/2022

Uncertainty-Aware Blob Detection with an Application to Integrated-Light Stellar Population Recoveries

Context. Blob detection is a common problem in astronomy. One example is...

Please sign up or login with your details

Forgot password? Click here to reset