A boundary preserving numerical algorithm for the Wright-Fisher model with mutation

Dangerfield, C.E., Kay, D., MacNamara, S., & Burrage, K. (2012) A boundary preserving numerical algorithm for the Wright-Fisher model with mutation. BIT Numerical Mathematics, 52(2), pp. 283-304.

View at publisher


The Wright-Fisher model is an Itô stochastic differential equation that was originally introduced to model genetic drift within finite populations and has recently been used as an approximation to ion channel dynamics within cardiac and neuronal cells. While analytic solutions to this equation remain within the interval [0,1], current numerical methods are unable to preserve such boundaries in the approximation. We present a new numerical method that guarantees approximations to a form of Wright-Fisher model, which includes mutation, remain within [0,1] for all time with probability one. Strong convergence of the method is proved and numerical experiments suggest that this new scheme converges with strong order 1/2. Extending this method to a multidimensional case, numerical tests suggest that the algorithm still converges strongly with order 1/2. Finally, numerical solutions obtained using this new method are compared to those obtained using the Euler-Maruyama method where the Wiener increment is resampled to ensure solutions remain within [0,1].

Impact and interest:

6 citations in Scopus
5 citations in Web of Science®
Search Google Scholar™

Citation counts are sourced monthly from Scopus and Web of Science® citation databases.

These databases contain citations from different subsets of available publications and different time periods and thus the citation count from each is usually different. Some works are not in either database and no count is displayed. Scopus includes citations from articles published in 1996 onwards, and Web of Science® generally from 1980 onwards.

Citations counts from the Google Scholar™ indexing service can be viewed at the linked Google Scholar™ search.

Full-text downloads:

153 since deposited on 01 Jul 2012
28 in the past twelve months

Full-text downloads displays the total number of times this work’s files (e.g., a PDF) have been downloaded from QUT ePrints as well as the number of downloads in the previous 365 days. The count includes downloads for all files if a work has more than one.

ID Code: 51303
Item Type: Journal Article
Refereed: Yes
Keywords: Boundary preserving numerical algorithm, Hölder condition, Ion channels, Split step, Stochastic differential equations, Strong convergence, Wright-Fisher model
DOI: 10.1007/s10543-011-0351-3
ISSN: 0006-3835
Subjects: Australian and New Zealand Standard Research Classification > MATHEMATICAL SCIENCES (010000)
Divisions: Current > Schools > School of Mathematical Sciences
Current > QUT Faculties and Divisions > Science & Engineering Faculty
Copyright Owner: Copyright 2012 Springer
Copyright Statement: The original publication is available at SpringerLink
Deposited On: 01 Jul 2012 22:55
Last Modified: 04 Jul 2012 05:44

Export: EndNote | Dublin Core | BibTeX

Repository Staff Only: item control page