English

Posterior consistency and convergence rates for Bayesian inversion with hypoelliptic operators

Statistics Theory 2016-07-20 v2 Functional Analysis Probability Statistics Theory

Abstract

Bayesian approach to inverse problems is studied in the case where the forward map is a linear hypoelliptic pseudodifferential operator and measurement error is additive white Gaussian noise. The measurement model for an unknown Gaussian random variable U(x,ω)U(x,\omega) is \begin{eqnarray*} M(y,\omega) = A(U(x,\omega) )+ \delta\hspace{.2mm}\mathcal{E}(y,\omega), \end{eqnarray*} where AA is a finitely many times smoothing linear hypoelliptic operator and δ>0\delta>0 is the noise magnitude. The covariance operator CUC_U of UU is 2r2r times smoothing, self-adjoint, injective and elliptic pseudodifferential operator. If E\mathcal{E} was taking values in L2L^2 then in Gaussian case solving the conditional mean (and maximum a posteriori) estimate is linked to solving the minimisation problem \begin{eqnarray*} T_\delta(M) = \text{argmin}_{u\in H^r} \big\{\|A u-m\|_{L^2}^2+ \delta^2\|C_U^{-1/2}u\|_{L^2}^2 \big\}. \end{eqnarray*} However, Gaussian white noise does not take values in L2L^2 but in HsH^{-s} where s>0s>0 is big enough. A modification of the above approach to solve the inverse problem is presented, covering the case of white Gaussian measurement noise. Furthermore, the convergence of conditional mean estimate to the correct solution as δ0\delta\rightarrow 0 is proven in appropriate function spaces using microlocal analysis. Also the contraction of the confidence regions is studied.

Keywords

Cite

@article{arxiv.1507.01772,
  title  = {Posterior consistency and convergence rates for Bayesian inversion with hypoelliptic operators},
  author = {Hanne Kekkonen and Matti Lassas and Samuli Siltanen},
  journal= {arXiv preprint arXiv:1507.01772},
  year   = {2016}
}
R2 v1 2026-06-22T10:07:12.845Z