Refine
Document Type
- Article (22)
- ZIB-Report (5)
- In Proceedings (3)
- Book (1)
- In Collection (1)
Language
- English (32) (remove)
Is part of the Bibliography
- no (32)
Keywords
- uncertainty quantification (2)
- Bayesian inverse problem (1)
- Bayesian inverse problems (1)
- Gaussian Processes (1)
- Inverse Problems (1)
- Karhunen–Loève expansion (1)
- Meshless Methods (1)
- Partial Differential Equations (1)
- Probabilistic Numerics (1)
- Pseudo-Marginal MCMC (1)
- heavy-tailed distribution (1)
- posterior consistency (1)
- probabilistic numerics (1)
- random likelihood (1)
- randomised misfit (1)
- shape space, shape trajectories, geodesic regression, longitudinal analysis, osteoarthritis (1)
- stable distribution (1)
- surrogate model (1)
- well-posedness (1)
This article attempts to place the emergence of probabilistic numerics as a mathematical-statistical research field within its historical context and to explore how its gradual development can be related to modern formal treatments and applications. We highlight in particular the parallel contributions of Sul'din and Larkin in the 1960s and how their pioneering early ideas have reached a degree of maturity in the intervening period, mediated by paradigms such as average-case analysis and information-based complexity. We provide a subjective assessment of the state of research in probabilistic numerics and highlight some difficulties to be addressed by future works.
It is well understood that Bayesian decision theory and average case analysis are essentially identical. However, if one is interested in performing uncertainty quantification for a numerical task, it can be argued that the decision-theoretic framework is neither appropriate nor sufficient. To this end, we consider an alternative optimality criterion from Bayesian experimental design and study its implied optimal information in the numerical context. This information is demonstrated to differ, in general, from the information that would be used in an average-case-optimal numerical method. The explicit connection to Bayesian experimental design suggests several distinct regimes in which optimal probabilistic numerical methods can be developed.
For Kendall’s shape space we determine analytically Jacobi fields and parallel transport, and compute geodesic regression. Using the derived expressions, we can fully leverage the geometry via Riemannian optimization and reduce the computational expense by several orders of magnitude. The methodology is demonstrated by performing a longitudinal statistical analysis of epidemiological shape data.
As application example we have chosen 3D shapes of knee bones, reconstructed from image data of the Osteoarthritis Initiative. Comparing subject groups with incident and developing osteoarthritis versus normal controls, we find clear differences in the temporal development of femur shapes. This paves the way for early prediction of incident knee osteoarthritis, using geometry data only.