Behavior of different numerical schemes for random genetic drift

Loading...
Thumbnail Image

Date

2019-09-01

Journal Title

Journal ISSN

Volume Title

Repository Usage Stats

96
views
8
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.

Department

Description

Provenance

Citation

Published Version (Please cite this version)

10.1007/s10543-019-00749-4

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.

Scholars@Duke

Xu

Shixin Xu

Assistant Professor of Mathematics at Duke Kunshan University

Shixin Xu is an Assistant Professor of Mathematics.   His research interests are machine learning and data-driven models for diseases,  multiscale modeling of complex fluids, Neurovascular coupling, homogenization theory, and numerical analysis.  The current projects he is working on are

  • image data-based for the prediction of hemorrhagic transformation in acute ischemic stroke,
  • electrodynamics modeling of saltatory conduction along a myelinated axon
  • electrochemical modeling
  • fluid-structure interaction with mass transportation and reaction

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.