Abstract
We present a hierarchical decomposition scheme for computing the n-dimensional integral of multivariate normal probabilities that appear frequently in statistics. The scheme exploits the fact that the formally dense covariance matrix can be approximated by a matrix with a hierarchical low rank structure. It allows the reduction of the computational complexity per Monte Carlo sample from O(n2) to O(mn+knlog(n/m)), where k is the numerical rank of off-diagonal matrix blocks and m is the size of small diagonal blocks in the matrix that are not well-approximated by low rank factorizations and treated as dense submatrices. This hierarchical decomposition leads to substantial efficiencies in multivariate normal probability computations and allows integrations in thousands of dimensions to be practical on modern workstations.
Original language | English (US) |
---|---|
Pages (from-to) | 268-277 |
Number of pages | 10 |
Journal | Journal of Computational and Graphical Statistics |
Volume | 27 |
Issue number | 2 |
DOIs | |
State | Published - May 17 2018 |
Bibliographical note
KAUST Repository Item: Exported on 2020-10-01Fingerprint
Dive into the research topics of 'Hierarchical Decompositions for the Computation of High-Dimensional Multivariate Normal Probabilities'. Together they form a unique fingerprint.Datasets
-
Hierarchical Decompositions for the Computation of High-Dimensional Multivariate Normal Probabilities
Genton, M. G. (Creator), Keyes, D. E. (Creator), Turkiyyah, G. (Creator) & Turkiyyah, G. (Creator), figshare, 2017
DOI: 10.6084/m9.figshare.5386996, http://hdl.handle.net/10754/663816
Dataset