Importance sampling for the infinite sites model

Research output: Contribution to journalJournal articleResearchpeer-review

Importance sampling or Markov Chain Monte Carlo sampling is required for state-of-the-art statistical analysis of population genetics data. The applicability of these sampling-based inference techniques depends crucially on the proposal distribution. In this paper, we discuss importance sampling for the infinite sites model. The infinite sites assumption is attractive because it constraints the number of possible genealogies, thereby allowing for the analysis of larger data sets. We recall the Griffiths-Tavaré and Stephens-Donnelly proposals and emphasize the relation between the latter proposal and exact sampling from the infinite alleles model. We also introduce a new proposal that takes knowledge of the ancestral state into account. The new proposal is derived from a new result on exact sampling from a single site. The methods are illustrated on simulated data sets and the data considered in Griffiths and Tavaré (1994).

Original languageEnglish
Article number32
JournalStatistical Applications in Genetics and Molecular Biology
Volume7
Issue number1
ISSN1544-6115
DOIs
Publication statusPublished - 1 Jan 2008
Externally publishedYes

    Research areas

  • Ancestral inference, Coalescent, Importance sampling, Infinite sites

ID: 203905609