Approximate maximum likelihood estimation for population genetic inference

Johanna Bertl, Gregory Ewing, Carolin Kosiol, Andreas Futschik

Research output: Contribution to journalArticlepeer-review

1 Downloads (Pure)


In many population genetic problems, parameter estimation is obstructed by an intractable likelihood function. Therefore, approximate estimation methods have been developed, and with growing computational power, sampling-based methods became popular. However, these methods such as Approximate Bayesian Computation (ABC) can be inefficient in high-dimensional problems. This led to the development of more sophisticated iterative estimation methods like particle filters. Here, we propose an alternative approach that is based on stochastic approximation. By moving along a simulated gradient or ascent direction, the algorithm produces a sequence of estimates that eventually converges to the maximum likelihood estimate, given a set of observed summary statistics. This strategy does not sample much from low-likelihood regions of the parameter space, and is fast, even when many summary statistics are involved. We put considerable efforts into providing tuning guidelines that improve the robustness and lead to good performance on problems with high-dimensional summary statistics and a low signal-to-noise ratio. We then investigate the performance of our resulting approach and study its properties in simulations. Finally, we re-estimate parameters describing the demographic history of Bornean and Sumatran orang-utans.
Original languageEnglish
Pages (from-to)291-312
JournalStatistical Applications in Genetics and Molecular Biology
Issue number5-6
Publication statusPublished - 2 Nov 2017


  • Approximate inference
  • Isolation-migration model
  • Maximum likelihood estimation
  • Orang-utans
  • Population genetics
  • Stochastic approximation


Dive into the research topics of 'Approximate maximum likelihood estimation for population genetic inference'. Together they form a unique fingerprint.

Cite this