Unbiased Estimation using a Class of Diffusion Processes

Abstract

We study the problem of unbiased estimation of expectations with respect to (w.r.t.) π a given, general probability measure on (Rd,B(Rd)) that is absolutely continuous with respect to a standard Gaussian measure. We focus on simulation associated to a particular class of diffusion processes, sometimes termed the Schr\"odinger-F\"ollmer Sampler, which is a simulation technique that approximates the law of a particular diffusion bridge process \Xt\t∈ [0,1] on Rd, d∈ N0. This latter process is constructed such that, starting at X0=0, one has X1 π. Typically, the drift of the diffusion is intractable and, even if it were not, exact sampling of the associated diffusion is not possible. As a result, sforig,jiao consider a stochastic Euler-Maruyama scheme that allows the development of biased estimators for expectations w.r.t.~π. We show that for this methodology to achieve a mean square error of O(ε2), for arbitrary ε>0, the associated cost is O(ε-5). We then introduce an alternative approach that provides unbiased estimates of expectations w.r.t.~π, that is, it does not suffer from the time discretization bias or the bias related with the approximation of the drift function. We prove that to achieve a mean square error of O(ε2), the associated cost is, with high probability, O(ε-2|(ε)|2+δ), for any δ>0. We implement our method on several examples including Bayesian inverse problems.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…