Markov chain Monte Carlo without likelihoods

Many stochastic simulation approaches for generating observations from a posterior distribution depend on knowing a likelihood function. However, for many complex probability models, such likelihoods are either impossible or computationally prohibitive to obtain. Here we present a Markov chain Monte Carlo method for generating observations from a posterior distribution without the use of likelihoods. It can also be used in frequentist applications, in particular for maximum-likelihood estimation. The approach is illustrated by an example of ancestral inference in population genetics. A number of open problems are highlighted in the discussion.

[1]  M. Feldman,et al.  Population growth of human Y chromosomes: a study of Y chromosome microsatellites. , 1999, Molecular biology and evolution.

[2]  D. Balding,et al.  Handbook of statistical genetics , 2004 .

[3]  A. von Haeseler,et al.  Inference of population history using a likelihood approach. , 1998, Genetics.

[4]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[5]  B. Bainbridge,et al.  Genetics , 1981, Experientia.

[6]  L. L. Cam,et al.  Sufficiency and Approximate Sufficiency , 1964 .

[7]  Bradley P. Carlin,et al.  BAYES AND EMPIRICAL BAYES METHODS FOR DATA ANALYSIS , 1996, Stat. Comput..

[8]  S. Tavaré,et al.  The effects of rate variation on ancestral inference in the coalescent. , 2000, Genetics.

[9]  H. Niederreiter,et al.  Monte Carlo and Quasi-Monte Carlo Methods in Scientific Computing , 1995 .

[10]  D. Balding,et al.  Approximate Bayesian computation in population genetics. , 2002, Genetics.

[11]  J. Felsenstein,et al.  Inching toward reality: An improved likelihood model of sequence evolution , 2004, Journal of Molecular Evolution.

[12]  J. Kingman.  On the genealogy of large populations , 1982 .

[13]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .

[14]  S. Tavaré,et al.  The age of a unique event polymorphism. , 2000, Genetics.

[15]  C. Fischer.  Handbook of statistical genetics: , 2002, Human Genetics.

[16]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[17]  K. Pearson.  Biometrika , 1902, The American Naturalist.

[18]  S. Pääbo,et al.  Extensive mitochondrial diversity within a single Amerindian tribe. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Brian D. Ripley,et al.  Stochastic Simulation , 2005 .

[20]  P. Donnelly,et al.  Inferring coalescence times from DNA sequence data. , 1997, Genetics.

[21]  W. Li,et al.  Estimating the age of the common ancestor of a sample of DNA sequences. , 1997, Molecular biology and evolution.