A boundary preserving numerical algorithm for the Wright-Fisher model with mutation
Data(s) |
2012
|
---|---|
Resumo |
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]. |
Formato |
application/pdf |
Identificador | |
Publicador |
Springer |
Relação |
http://eprints.qut.edu.au/51303/1/51303.pdf DOI:10.1007/s10543-011-0351-3 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. |
Direitos |
Copyright 2012 Springer The original publication is available at SpringerLink http://www.springerlink.com |
Fonte |
School of Mathematical Sciences; Science & Engineering Faculty |
Palavras-Chave | #010000 MATHEMATICAL SCIENCES #Boundary preserving numerical algorithm #Hölder condition #Ion channels #Split step #Stochastic differential equations #Strong convergence #Wright-Fisher model |
Tipo |
Journal Article |