pith. sign in

arxiv: 2411.08338 · v1 · pith:PGN33AOHnew · submitted 2024-11-13 · 🧮 math.NA · cs.NA· physics.data-an· stat.ME

Quantifying uncertainty in the numerical integration of evolution equations based on Bayesian isotonic regression

classification 🧮 math.NA cs.NAphysics.data-anstat.ME
keywords variancesbayesiandiscretizationdistributionequationserrorsisotonicnumerical
0
0 comments X
read the original abstract

This paper presents a new Bayesian framework for quantifying discretization errors in numerical solutions of ordinary differential equations. By modelling the errors as random variables, we impose a monotonicity constraint on the variances, referred to as discretization error variances. The key to our approach is the use of a shrinkage prior for the variances coupled with variable transformations. This methodology extends existing Bayesian isotonic regression techniques to tackle the challenge of estimating the variances of a normal distribution. An additional key feature is the use of a Gaussian mixture model for the $\log$-$\chi^2_1$ distribution, enabling the development of an efficient Gibbs sampling algorithm for the corresponding posterior.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.