EWF : simulating exact paths of the Wright--Fisher diffusion

DOI10.1093/BIOINFORMATICS/BTAD017MaRDI QIDQ6423296zbMATH OpenFDO

Authors Paul A. Jenkins, Jere Koskela, Dario Spanò

Publication date 13 January 2023

Full work available at URL https://doi.org/10.1093/BIOINFORMATICS/BTAD017

Diffusion processes (60J60) Applications of Brownian motions and diffusion theory (population genetics, absorption problems, etc.) (60J70) Probabilistic methods, stochastic differential equations (65C99) Population dynamics (general) (92D25)

Abstract: The Wright--Fisher diffusion is important in population genetics in modelling the evolution of allele frequencies over time subject to the influence of biological phenomena such as selection, mutation, and genetic drift. Simulating paths of the process is challenging due to the form of the transition density. We present EWF, a robust and efficient sampler which returns exact draws for the diffusion and diffusion bridge processes, accounting for general models of selection including those with frequency-dependence. Given a configuration of selection, mutation, and endpoints, EWF returns draws at the requested sampling times from the law of the corresponding Wright--Fisher process. Output was validated by comparison to approximations of the transition density via the Kolmogorov--Smirnov test and QQ plots. All software is available at https://github.com/JaroSant/EWF

This page was built for publication: EWF : simulating exact paths of the Wright--Fisher diffusion

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6423296)