Behavior of different numerical schemes for random genetic drift
Date
2019-09-01
Journal Title
Journal ISSN
Volume Title
Repository Usage Stats
views
downloads
Citation Stats
Abstract
In the problem of random genetic drift, the probability density of one gene is governed by a degenerated convection-dominated diffusion equation. Dirac singularities will always be developed at boundary points as time evolves, which is known as the fixation phenomenon in genetic evolution. Three finite volume methods: FVM1-3, one central difference method: FDM1 and three finite element methods: FEM1-3 are considered. These methods lead to different equilibrium states after a long time. It is shown that only schemes FVM3 and FEM3, which are the same, preserve probability, expectation and positiveness and predict the correct probability of fixation. FVM1-2 wrongly predict the probability of fixation due to their intrinsic viscosity, even though they are unconditionally stable. Contrarily, FDM1 and FEM1-2 introduce different anti-diffusion terms, which make them unstable and fail to preserve positiveness.
Type
Department
Description
Provenance
Subjects
Citation
Permalink
Published Version (Please cite this version)
Publication Info
Xu, S, M Chen, C Liu, R Zhang and X Yue (2019). Behavior of different numerical schemes for random genetic drift. BIT Numerical Mathematics, 59(3). pp. 797–821. 10.1007/s10543-019-00749-4 Retrieved from https://hdl.handle.net/10161/27448.
This is constructed from limited available data and may be imprecise. To cite this article, please review & use the official citation provided by the journal.
Collections
Scholars@Duke
Shixin Xu
Shixin Xu is an Assistant Professor of Mathematics whose research spans several dynamic and interconnected fields. His primary interests include machine learning and data-driven models for disease prediction, multiscale modeling of complex fluids, neurovascular coupling, homogenization theory, and numerical analysis. His current projects reflect a diverse and impactful portfolio:
- Developing predictive models based on image data to identify hemorrhagic transformation in acute ischemic stroke.
- Conducting electrodynamics modeling of saltatory conduction along myelinated axons to understand nerve impulse transmission.
- Engaging in electrochemical modeling to explore the interactions between electric fields and chemical processes.
- Investigating fluid-structure interactions with mass transport and reactions, crucial for understanding physiological and engineering systems.
These projects demonstrate his commitment to addressing complex problems through interdisciplinary approaches that bridge mathematics with biological and physical sciences.
Unless otherwise indicated, scholarly articles published by Duke faculty members are made available here with a CC-BY-NC (Creative Commons Attribution Non-Commercial) license, as enabled by the Duke Open Access Policy. If you wish to use the materials in ways not already permitted under CC-BY-NC, please consult the copyright owner. Other materials are made available here through the author’s grant of a non-exclusive license to make their work openly accessible.