Marjoram Paul, Molitor John, Plagnol Vincent, Tavare Simon
Department of Preventive Medicine, Keck School of Medicine, University of Southern California, Los Angeles, CA 90089, USA.
Proc Natl Acad Sci U S A. 2003 Dec 23;100(26):15324-8. doi: 10.1073/pnas.0306899100. Epub 2003 Dec 8.
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.
许多用于从后验分布生成观测值的随机模拟方法都依赖于知道似然函数。然而,对于许多复杂的概率模型,获得这样的似然函数要么是不可能的,要么在计算上是令人望而却步的。在此,我们提出一种马尔可夫链蒙特卡罗方法,用于在不使用似然函数的情况下从后验分布生成观测值。它也可用于频率论应用,特别是用于最大似然估计。通过群体遗传学中祖先推断的一个例子来说明该方法。讨论中强调了一些未解决的问题。