Likelihoods from summary statistics: Recent divergence between species

Scotland C. Leman, Yuguo Chen, Jason E. Stajich, Mohamed A.F. Noor, Marcy K. Uyenoyama

Research output: Contribution to journalArticlepeer-review


We describe an importance-sampling method for approximating likelihoods of population parameters based on multiple summary statistics. In this first application, we address the demographic history of closely related members of the Drosophila pseudoobscura group. We base the maximum-likelihood estimation of the time since speciation and the effective population sizes of the extant and ancestral populations on the pattern of nucleotide variation at DPS2002, a noncoding region tightly linked to a paracentric inversion that strongly contributes to reproductive isolation. Consideration of summary statistics rather than entire nucleotide sequences permits a compact description of the genealogy of the sample. We use importance sampling first to propose a genealogical and mutational history consistent with the observed array of summary statistics and then to correct the likelihood with the exact probability of the history determined from a system of recursions. Analysis of a subset of the data, for which recursive computation of the exact likelihood was feasible, indicated close agreement between the approximate and exact likelihoods. Our results for the complete data set also compare well with those obtained through Metropolis-Hastings sampling of fully resolved genealogies of entire nucleotide sequences.

Original languageEnglish (US)
Pages (from-to)1419-1436
Number of pages18
Issue number3
StatePublished - Nov 2005

ASJC Scopus subject areas

  • General Medicine


Dive into the research topics of 'Likelihoods from summary statistics: Recent divergence between species'. Together they form a unique fingerprint.

Cite this