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

Data-dependent density estimation for the Fokker-Planck equation in higher dimensions

Interacting Particle Solutions of Fokker–Planck Equations Through Gradient–Log–Density Estimation

Solving high-dimensional Fokker-Planck equation with functional hierarchical tensor

Fisher information and shape-morphing modes for solving the Fokker-Planck equation in higher dimensions

Beating the curse of dimension with accurate statistics for the Fokker–Planck equation in complex turbulent systems

Deterministic Fokker-Planck Transport -- With Applications to Sampling, Variational Inference, Kernel Mean Embeddings & Sequential Monte Carlo

Data-driven probability density forecast for stochastic dynamical systems

Discretizing the Fokker-Planck equation with second-order accuracy: a dissipation driven approach

Efficient numerical approximation of a non-regular Fokker-Planck equation associated with first-passage time distributions

A BDF2-Approach for the Non-linear Fokker-Planck Equation

A statistical approach for simulating the density solution of a McKean-Vlasov equation

Fractional diffusion for Fokker-Planck equation with heavy tail equilibrium: an à la Koch spectral method in any dimension

Convex Relaxation for Fokker-Planck

Solving Fokker-Planck equations using the zeros of Fokker-Planck operators and the Feynman-Kac formula

Strong Convergence of the JKO Scheme for the Fokker–Planck Equation

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

Dimensional reduction of gradient-like stochastic systems with multiplicative noise via Fokker-Planck diffusion maps

Epsilon-Expansion of the Stationary Probability-Distribution of 2-Dimensional Fokker-Planck Equations

Long-time behavior for discretization schemes of Fokker-Planck equations via couplings

Exponential Stability and Hypoelliptic Regularization for the Kinetic Fokker–Planck Equation with Confining Potential