Title
Probabilistic Solution of Differential Equations for Bayesian Uncertainty Quanitification and Inference
Speaker
Oksana Chkrebtii, Simon Fraser University
Abstract
In many scientific disciplines, the time and space evolution of system states is naturally described by differential equation models, which define states implicitly as functions of their own rates of change. Inference for differential equation models requires an explicit representation of the states. This solution is typically not known in closed form, but can be approximated by a variety of discretization-based numerical methods. However, the associated numerical error analysis cannot currently be propagated through the statistical inverse problem, and is thus ignored in practice.
We resolve this problem by developing a Bayesian framework to characterize discretization uncertainty in models defined by high-dimensional systems of differential equations with unknown solutions. Viewing solution estimation as an inference problem allows us to quantify and propagate discretization error into uncertainty in the model parameters and subsequent predictions. In this talk, I will introduce our new methodology, discussing efficient computational algorithms and application to a wide range of challenging forward and inverse problems.
