Stein variational gradient descent on infinite-dimensional space and applications to statistical inverse problems

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Jia, Junxiong, Li, Peijun, Meng, Deyu
Format: Preprint
Veröffentlicht: 2021
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866914309106302976
author Jia, Junxiong
Li, Peijun
Meng, Deyu
author_facet Jia, Junxiong
Li, Peijun
Meng, Deyu
contents In this paper, we propose an infinite-dimensional version of the Stein variational gradient descent (iSVGD) method for solving Bayesian inverse problems. The method can generate approximate samples from posteriors efficiently. Based on the concepts of operator-valued kernels and vector-valued reproducing kernel Hilbert spaces, a rigorous definition is given for the infinite-dimensional objects, e.g., the Stein operator, which are proved to be the limit of finite-dimensional ones. Moreover, a more efficient iSVGD with preconditioning operators is constructed by generalizing the change of variables formula and introducing a regularity parameter. The proposed algorithms are applied to an inverse problem of the steady state Darcy flow equation. Numerical results confirm our theoretical findings and demonstrate the potential applications of the proposed approach in the posterior sampling of large-scale nonlinear statistical inverse problems.
format Preprint
id arxiv_https___arxiv_org_abs_2102_09741
institution arXiv
publishDate 2021
record_format arxiv
spellingShingle Stein variational gradient descent on infinite-dimensional space and applications to statistical inverse problems
Jia, Junxiong
Li, Peijun
Meng, Deyu
Numerical Analysis
65L09, 49N45, 62F15
In this paper, we propose an infinite-dimensional version of the Stein variational gradient descent (iSVGD) method for solving Bayesian inverse problems. The method can generate approximate samples from posteriors efficiently. Based on the concepts of operator-valued kernels and vector-valued reproducing kernel Hilbert spaces, a rigorous definition is given for the infinite-dimensional objects, e.g., the Stein operator, which are proved to be the limit of finite-dimensional ones. Moreover, a more efficient iSVGD with preconditioning operators is constructed by generalizing the change of variables formula and introducing a regularity parameter. The proposed algorithms are applied to an inverse problem of the steady state Darcy flow equation. Numerical results confirm our theoretical findings and demonstrate the potential applications of the proposed approach in the posterior sampling of large-scale nonlinear statistical inverse problems.
title Stein variational gradient descent on infinite-dimensional space and applications to statistical inverse problems
topic Numerical Analysis
65L09, 49N45, 62F15
url https://arxiv.org/abs/2102.09741