Toggle Main Menu Toggle Search

Open Access padlockePrints

The Computational Cost of Blocking for Sampling Discretely Observed Diffusions

Lookup NU author(s): Dr Murray Pollock

Downloads


Licence

This work is licensed under a Creative Commons Attribution 4.0 International License (CC BY 4.0).


Abstract

© 2022, The Author(s). Many approaches for conducting Bayesian inference on discretely observed diffusions involve imputing diffusion bridges between observations. This can be computationally challenging in settings in which the temporal horizon between subsequent observations is large, due to the poor scaling of algorithms for simulating bridges as observation distance increases. It is common in practical settings to use a blocking scheme, in which the path is split into a (user-specified) number of overlapping segments and a Gibbs sampler is employed to update segments in turn. Substituting the independent simulation of diffusion bridges for one obtained using blocking introduces an inherent trade-off: we are now imputing shorter bridges at the cost of introducing a dependency between subsequent iterations of the bridge sampler. This is further complicated by the fact that there are a number of possible ways to implement the blocking scheme, each of which introduces a different dependency structure between iterations. Although blocking schemes have had considerable empirical success in practice, there has been no analysis of this trade-off nor guidance to practitioners on the particular specifications that should be used to obtain a computationally efficient implementation. In this article we conduct this analysis and demonstrate that the expected computational cost of a blocked path-space rejection sampler applied to Brownian bridges scales asymptotically at a cubic rate with respect to the observation distance and that this rate is linear in the case of the Ornstein–Uhlenbeck process. Numerical experiments suggest applicability both of the results of our paper and of the guidance we provide beyond the class of linear diffusions considered.


Publication metadata

Author(s): Mider M, Jenkins PA, Pollock M, Roberts GO

Publication type: Article

Publication status: Published

Journal: Methodology and Computing in Applied Probability

Year: 2022

Volume: 24

Pages: 3007-3027

Print publication date: 01/12/2022

Online publication date: 14/07/2022

Acceptance date: 10/03/2022

Date deposited: 12/08/2022

ISSN (print): 1387-5841

ISSN (electronic): 1573-7713

Publisher: Springer Nature

URL: https://doi.org/10.1007/s11009-022-09949-y

DOI: 10.1007/s11009-022-09949-y


Altmetrics

Altmetrics provided by Altmetric


Funding

Funder referenceFunder name
EP/K034154/1
EP/R018561/1
EP/K014463/1
EP/L016710/1
EP/N510129/1

Share