Abstract:We consider the problem of making nonparametric inference in multi-dimensional diffusion models from low-frequency data. Statistical analysis in this setting is notoriously challenging due to the intractability of the likelihood and its gradient, and computational methods have thus far largely resorted to expensive simulation-based techniques. In this article, we propose a new computational approach which is motivated by PDE theory and is built around the characterisation of the transition densities as solutions of the associated heat (Fokker-Planck) equation. Employing optimal regularity results from the theory of parabolic PDEs, we prove a novel characterisation for the gradient of the likelihood. Using these developments, for the nonlinear inverse problem of recovering the diffusivity (in divergence form models), we then show that the numerical evaluation of the likelihood and its gradient can be reduced to standard elliptic eigenvalue problems, solvable by powerful finite element methods. This enables the efficient implementation of a large class of statistical algorithms, including (i) preconditioned Crank-Nicolson and Langevin-type methods for posterior sampling, and (ii) gradient-based descent optimisation schemes to compute maximum likelihood and maximum-a-posteriori estimates. We showcase the effectiveness of these methods via extensive simulation studies in a nonparametric Bayesian model with Gaussian process priors. Interestingly, the optimisation schemes provided satisfactory numerical recovery while exhibiting rapid convergence towards stationary points despite the problem nonlinearity; thus our approach may lead to significant computational speed-ups. The reproducible code is available online at

Uniqueness and stabilized algorithm for an inverse problem of identifying space-time dependent diffusion coefficients

Simultaneous inversion for the space-dependent diffusion coefficient and the fractional order in the time-fractional diffusion equation

Identification of the initial value for a space-time fractional diffusion equation

Uniqueness, stability and algorithm for an inverse wave-number-dependent source problem

Uniqueness, stability and algorithm for an inverse wave-number-dependent source problems

On the convergence rate of an improved quasi-reversibility method for an inverse source problem of a nonlinear parabolic equation with nonlocal diffusion coefficient

Simultaneous uniqueness and numerical inversion for an inverse problem in the time-domain diffuse optical tomography with fluorescence

A note on uniqueness in the identification of a spacewise dependent source and diffusion coefficient for the heat equation

Simultaneous uniqueness for a coefficient inverse problem in one-dimensional fractional diffusion equation from an interior point measurement

Simultaneous numerical inversion of space-dependent initial condition and source term in multi-order time-fractional diffusion models

Uniqueness and Numerical Method for Determining a Spatial Source Term in a Time-Fractional Diffusion Wave Equation

Identification of Sparsely Representable Diffusion Parameters in Elliptic Problems

A dimension reduction method with applications for coefficient inversion of diffusion equations

Numerical Recovery of the Diffusion Coefficient in Diffusion Equations from Terminal Measurement

Conditional well-posedness and data-driven method for identifying the dynamic source in a coupled diffusion system from one single boundary measurement

Uniqueness of the inverse source problem for fractional diffusion-wave equations

Uniqueness from pointwise observations in a multi-parameter inverse problem

Inverse problem of reconstructing source term for a class of non‐divergence parabolic equations

Numerical Reconstruction of Diffusion and Potential Coefficients from Two Observations: Decoupled Recovery and Error Estimates

Statistical algorithms for low-frequency diffusion data: A PDE approach

Solving forward and inverse problems in a non-linear 3D PDE via an asymptotic expansion based approach