Bayesian analysis of ODEs: solver optimal accuracy and Bayes factors
From MaRDI portal
Publication:3179310
Abstract: In most relevant cases in the Bayesian analysis of ODE inverse problems, a numerical solver needs to be used. Therefore, we cannot work with the exact theoretical posterior distribution but only with an approximate posterior deriving from the error in the numerical solver. To compare a numerical and the theoretical posterior distributions we propose to use Bayes Factors (BF), considering both of them as models for the data at hand. We prove that the theoretical vs a numerical posterior BF tends to 1, in the same order (of the step size used) as the numerical forward map solver does. For higher order solvers (eg. Runge-Kutta) the Bayes Factor is already nearly 1 for step sizes that would take far less computational effort. Considerable CPU time may be saved by using coarser solvers that nevertheless produce practically error free posteriors. Two examples are presented where nearly 90% CPU time is saved while all inference results are identical to using a solver with a much finer time step.
Recommendations
- Efficient Bayesian estimation and uncertainty quantification in ordinary differential equation models
- Bayesian ODE solvers: the maximum a posteriori estimate
- scientific article; zbMATH DE number 802686
- Bayesian inference for higher-order ordinary differential equation models
- Adaptive step-size selection for state-space probabilistic differential equation solvers
Cites work
- scientific article; zbMATH DE number 52120 (Why is no real title available?)
- scientific article; zbMATH DE number 597911 (Why is no real title available?)
- scientific article; zbMATH DE number 2062250 (Why is no real title available?)
- scientific article; zbMATH DE number 2117879 (Why is no real title available?)
- scientific article; zbMATH DE number 3189754 (Why is no real title available?)
- A Bayesian approach to multiscale inverse problems using the sequential Monte Carlo method
- A First Course in the Numerical Analysis of Differential Equations
- A Posteriori Error Estimation and Global Error Control for Ordinary Differential Equations by the Adjoint Method
- A general purpose sampling algorithm for continuous distributions (the t-walk)
- An Introduction to Numerical Analysis
- An adaptive Metropolis algorithm
- Approximation of Bayesian Inverse Problems for PDEs
- Bayes Factors
- Bayesian MLP neural networks for image analysis
- Bayesian model averaging: A tutorial. (with comments and a rejoinder).
- Estimation of parameters in incomplete data models defined by dynamical systems
- Linear multistep methods, particle filtering and sequential Monte Carlo
- Local regularization method applied to estimating oxygen consumption during muscle activities
- Non-stationary magnetoencephalography by Bayesian filtering of dipole models
- Numerical mathematics.
- On Global Error Estimation and Control for Initial Value Problems
- On adaptive Markov chain Monte Carlo algorithms
- Parameter estimation for stiff deterministic dynamical systems via ensemble Kalman filter
- RECONSTRUCTION OF DOMAIN BOUNDARY AND CONDUCTIVITY IN ELECTRICAL IMPEDANCE TOMOGRAPHY USING THE APPROXIMATION ERROR APPROACH
- Sieve estimation of constant and time-varying coefficients in nonlinear ordinary differential equation models by considering both numerical error and measurement error
- Sparse deterministic approximation of Bayesian inverse problems
- Statistical and computational inverse problems.
- Towards uncertainty quantification and inference in the stochastic SIR epidemic model
Cited in
(13)- Efficient Bayesian estimation and uncertainty quantification in ordinary differential equation models
- Bayesian ODE solvers: the maximum a posteriori estimate
- Statistical analysis of differential equations: introducing probability measures on numerical solutions
- An importance sampling approach for reliable and efficient inference in Bayesian ordinary differential equation models
- Penalised t-walk MCMC
- Strong convergence rates of probabilistic integrators for ordinary differential equations
- Bayesian analysis for a fractional population growth model
- Bayesian solution uncertainty quantification for differential equations
- Disentangling the role of virus infectiousness and awareness-based human behavior during the early phase of the COVID-19 pandemic in the European Union
- Bayesian Probabilistic Numerical Methods
- Bayes linear analysis for ordinary differential equations
- Bayesian mesh adaptation for estimating distributed parameters
- Error control of the numerical posterior with Bayes factors in Bayesian uncertainty quantification
This page was built for publication: Bayesian analysis of ODEs: solver optimal accuracy and Bayes factors
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3179310)
