Exact simulation of the Wright-Fisher diffusion

[thumbnail of WRAP-exact-simulation-Wright-Fisher-diffusion-Jenkins-2017.pdf]
Preview
PDF
WRAP-exact-simulation-Wright-Fisher-diffusion-Jenkins-2017.pdf - Published Version - Requires a PDF viewer.

Download (728kB) | Preview
[thumbnail of WRAP_ea-wf_2016-07-22.pdf] PDF
WRAP_ea-wf_2016-07-22.pdf - Accepted Version
Embargoed item. Restricted access to Repository staff only - Requires a PDF viewer.

Download (756kB)
[thumbnail of SHERPA_RoMEO-Annals-of-Applied-Probability-15_11_2018.pdf] PDF
SHERPA_RoMEO-Annals-of-Applied-Probability-15_11_2018.pdf - Permissions Correspondence
Embargoed item. Restricted access to Repository staff only - Requires a PDF viewer.

Download (133kB)

Request Changes to record.

Abstract

The Wright–Fisher family of diffusion processes is a widely used class of evolutionary models. However, simulation is difficult because there is no known closed-form formula for its transition function. In this article, we demonstrate that it is in fact possible to simulate exactly from a broad class of Wright–Fisher diffusion processes and their bridges. For those diffusions corresponding to reversible, neutral evolution, our key idea is to exploit an eigenfunction expansion of the transition function; this approach even applies to its infinite-dimensional analogue, the Fleming–Viot process. We then develop an exact rejection algorithm for processes with more general drift functions, including those modelling natural selection, using ideas from retrospective simulation. Our approach also yields methods for exact simulation of the moment dual of the Wright–Fisher diffusion, the ancestral process of an infinite-leaf Kingman coalescent tree. We believe our new perspective on diffusion simulation holds promise for other models admitting a transition eigenfunction expansion.

Item Type: Journal Article
Subjects: Q Science > QA Mathematics
Divisions: Faculty of Science, Engineering and Medicine > Science > Statistics
Library of Congress Subject Headings (LCSH): Monte Carlo method, Diffusion processes, Perfect simulation (Statistics) , Population genetics
Journal or Publication Title: Annals of Applied Probability
Publisher: Institute of Mathematical Statistics
ISSN: 1050-5164
Official Date: 19 July 2017
Dates:
Date
Event
19 July 2017
Published
22 July 2016
Accepted
Volume: 27
Number: 3
Page Range: pp. 1478-1509
DOI: 10.1214/16-AAP1236
Status: Peer Reviewed
Publication Status: Published
Access rights to Published version: Restricted or Subscription Access
Date of first compliant deposit: 25 July 2016
Date of first compliant Open Access: 27 November 2017
Funder:
RIOXX Funder/Project Grant:
Project/Grant ID
RIOXX Funder Name
Funder ID
EP/L018497/1
Engineering and Physical Sciences Research Council
Related URLs:
Open Access Version:
URI: https://wrap.warwick.ac.uk/79729/

Export / Share Citation


Request changes or add full text files to a record

Repository staff actions (login required)

View Item View Item