Consistency of Bayesian inference with Gaussian process priors for a parabolic inverse problem

We consider the statistical nonlinear inverse problem of recovering the absorption term in the heat equation \partial_tu-\frac{1}{2}\Delta u+fu=0 \quad \text{on }\quad u = g \quad \text{on }\quad u(\cdot,0)=u_0 \quad \text{on }, where is a bounded domain, is a fixed time, and are given sufficiently smooth functions describing boundary and initial values respectively. The data consists of discrete noisy point evaluations of the solution on . We study the statistical performance of Bayesian nonparametric procedures based on a large class of Gaussian process priors. We show that, as the number of measurements increases, the resulting posterior distributions concentrate around the true parameter generating the data, and derive a convergence rate for the reconstruction error of the associated posterior means. We also consider the optimality of the contraction rates and prove a lower bound for the minimax convergence rate for inferring from the data, and show that optimal rates can be achieved with truncated Gaussian priors.
View on arXiv