In this article, we consider computing expectations w.r.t. probability measures which are subject to discretization error. Examples include partially observed diffusion processes or inverse problems, where one may have to discretize time and/or space, in order to practically work with the probability of interest. Given access only to these discretizations, we consider the construction of unbiased Monte Carlo estimators of expectations w.r.t. such target probability distributions. It is shown how to obtain such estimators using a novel adaptation of randomization schemes and Markov simulation methods. Under appropriate assumptions, these estimators possess finite variance and finite expected cost. There are two important consequences of this approach: (i) unbiased inference is achieved at the canonical complexity rate, and (ii) the resulting estimators can be generated independently, thereby allowing strong scaling to arbitrarily many parallel processors. Several algorithms are presented, and applied to some examples of Bayesian inference problems, with both simulated and real observed data.
翻译:在本篇文章中,我们考虑计算有离散错误的预期概率尺度,例子包括部分观测到的传播过程或反向问题,为了实际工作,人们可能不得不将时间和/或空间分解,以便实际工作。鉴于只能使用这些离散的可能性,我们考虑建造不带偏见的蒙特卡洛估计期望值,这样的目标概率分布。展示了如何利用随机化计划和Markov模拟方法的新调整来获得这种估计值。在适当的假设下,这些估计值具有有限的差异和预期的有限费用。这种方法有两个重要后果:(一) 以罐头复杂率实现公正的推论,以及(二) 由此产生的估计值可以独立生成,从而可以向许多任意平行的处理器大幅缩放。提出了几种算法,并应用到一些有模拟和真实观察数据的Bayesian推论问题的例子。