Tingwei Meng, Zongren Zou, Jérôme Darbon, George Em Karniadakis
{"title":"HJ-sampler: A Bayesian sampler for inverse problems of a stochastic process by leveraging Hamilton-Jacobi PDEs and score-based generative models","authors":"Tingwei Meng, Zongren Zou, Jérôme Darbon, George Em Karniadakis","doi":"arxiv-2409.09614","DOIUrl":null,"url":null,"abstract":"The interplay between stochastic processes and optimal control has been\nextensively explored in the literature. With the recent surge in the use of\ndiffusion models, stochastic processes have increasingly been applied to sample\ngeneration. This paper builds on the log transform, known as the Cole-Hopf\ntransform in Brownian motion contexts, and extends it within a more abstract\nframework that includes a linear operator. Within this framework, we found that\nthe well-known relationship between the Cole-Hopf transform and optimal\ntransport is a particular instance where the linear operator acts as the\ninfinitesimal generator of a stochastic process. We also introduce a novel\nscenario where the linear operator is the adjoint of the generator, linking to\nBayesian inference under specific initial and terminal conditions. Leveraging\nthis theoretical foundation, we develop a new algorithm, named the HJ-sampler,\nfor Bayesian inference for the inverse problem of a stochastic differential\nequation with given terminal observations. The HJ-sampler involves two stages:\n(1) solving the viscous Hamilton-Jacobi partial differential equations, and (2)\nsampling from the associated stochastic optimal control problem. Our proposed\nalgorithm naturally allows for flexibility in selecting the numerical solver\nfor viscous HJ PDEs. We introduce two variants of the solver: the\nRiccati-HJ-sampler, based on the Riccati method, and the SGM-HJ-sampler, which\nutilizes diffusion models. We demonstrate the effectiveness and flexibility of\nthe proposed methods by applying them to solve Bayesian inverse problems\ninvolving various stochastic processes and prior distributions, including\napplications that address model misspecifications and quantifying model\nuncertainty.","PeriodicalId":501215,"journal":{"name":"arXiv - STAT - Computation","volume":null,"pages":null},"PeriodicalIF":0.0000,"publicationDate":"2024-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - STAT - Computation","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.09614","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
The interplay between stochastic processes and optimal control has been
extensively explored in the literature. With the recent surge in the use of
diffusion models, stochastic processes have increasingly been applied to sample
generation. This paper builds on the log transform, known as the Cole-Hopf
transform in Brownian motion contexts, and extends it within a more abstract
framework that includes a linear operator. Within this framework, we found that
the well-known relationship between the Cole-Hopf transform and optimal
transport is a particular instance where the linear operator acts as the
infinitesimal generator of a stochastic process. We also introduce a novel
scenario where the linear operator is the adjoint of the generator, linking to
Bayesian inference under specific initial and terminal conditions. Leveraging
this theoretical foundation, we develop a new algorithm, named the HJ-sampler,
for Bayesian inference for the inverse problem of a stochastic differential
equation with given terminal observations. The HJ-sampler involves two stages:
(1) solving the viscous Hamilton-Jacobi partial differential equations, and (2)
sampling from the associated stochastic optimal control problem. Our proposed
algorithm naturally allows for flexibility in selecting the numerical solver
for viscous HJ PDEs. We introduce two variants of the solver: the
Riccati-HJ-sampler, based on the Riccati method, and the SGM-HJ-sampler, which
utilizes diffusion models. We demonstrate the effectiveness and flexibility of
the proposed methods by applying them to solve Bayesian inverse problems
involving various stochastic processes and prior distributions, including
applications that address model misspecifications and quantifying model
uncertainty.