Fast matrix algebra for Bayesian model calibration
- Univ. of New Mexico, Albuquerque, NM (United States); Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
In Bayesian model calibration, evaluation of the likelihood function usually involves finding the inverse and determinant of a covariance matrix. When Markov Chain Monte Carlo (MCMC) methods are used to sample from the posterior, hundreds of thousands of likelihood evaluations may be required. In this paper, we demonstrate that the structure of the covariance matrix can be exploited, leading to substantial time savings in practice. Here, we also derive two simple equations for approximating the inverse of the covariance matrix in this setting, which can be computed in near-quadratic time. The practical implications of these strategies are demonstrated using a simple numerical case study and the "quack" R package. For a covariance matrix with 1000 rows, application of these strategies for a million likelihood evaluations leads to a speedup of roughly 4000 compared to the naive implementation
- Research Organization:
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA); USDOE Laboratory Directed Research and Development (LDRD) Program
- Grant/Contract Number:
- NA0003525
- OSTI ID:
- 1828021
- Report Number(s):
- SAND-2020-12827J; 697276
- Journal Information:
- Journal of Statistical Computation and Simulation, Vol. 91, Issue 7; ISSN 0094-9655
- Publisher:
- Taylor & FrancisCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Analyzing rotational bands in odd-mass nuclei using effective field theory and Bayesian methods
Computationally efficient Bayesian inference for inverse problems.