接受-拒绝抽样下的蒙特卡洛积分

Monte Carlo Integration With Acceptance-Rejection

Journal of Computational and Graphical Statistics · 2006
被引 17
ABS 3

中文导读

研究了在拒绝抽样或Metropolis-Hastings抽样下进行蒙特卡洛积分的方法,提出似然估计量并证明其渐近方差更小,通过子采样降低计算成本,并用三个例子展示有效性。

Abstract

This article considers Monte Carlo integration under rejection sampling or Metropolis-Hastings sampling. Each algorithm involves accepting or rejecting observations from proposal distributions other than a target distribution. While taking a likelihood approach, we basically treat the sampling scheme as a random design, and define a stratified estimator of the baseline measure. We establish that the likelihood estimator has no greater asymptotic variance than the crude Monte Carlo estimator under rejection sampling or independence Metropolis-Hastings sampling. We employ a subsampling technique to reduce the computational cost, and illustrate with three examples the computational effectiveness of the likelihood method under general Metropolis-Hastings sampling.

蒙特卡洛方法接受-拒绝抽样Metropolis-Hastings算法方差缩减计算统计