Skip to Main content Skip to Navigation
Journal articles

Estimating orthant probabilities of high dimensional Gaussian vectors with an application to set estimation

Abstract : The computation of Gaussian orthant probabilities has been extensively studied for low-dimensional vectors. Here, we focus on the high-dimensional case and we present a two-step procedure relying on both deterministic and stochastic techniques. The proposed estimator relies indeed on splitting the probability into a low-dimensional term and a remainder. While the low-dimensional probability can be estimated by fast and accurate quadrature, the remainder requires Monte Carlo sampling. We further refine the estimation by using a novel asymmetric nested Monte Carlo (anMC) algorithm for the remainder and we highlight cases where this approximation brings substantial efficiency gains. The proposed methods are compared against state-of-the-art techniques in a numerical study, which also calls attention to the advantages and drawbacks of the procedure. Finally, the proposed method is applied to derive conservative estimates of excursion sets of expensive to evaluate deterministic functions under a Gaussian random field prior, without requiring a Markov assumption. Supplementary material for this article is available online.
Complete list of metadata

https://hal.archives-ouvertes.fr/hal-01289126
Contributor : Azzimonti Dario Connect in order to contact the contributor
Submitted on : Friday, November 30, 2018 - 11:03:33 AM
Last modification on : Friday, November 30, 2018 - 1:50:58 PM
Long-term archiving on: : Friday, March 1, 2019 - 1:57:43 PM

Files

mainArticle.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : hal-01289126, version 3
  • ARXIV : 1603.05031

Citation

Dario Azzimonti, David Ginsbourger. Estimating orthant probabilities of high dimensional Gaussian vectors with an application to set estimation. Journal of Computational and Graphical Statistics, Taylor & Francis, 2018, 27 (2), pp.255-267. ⟨hal-01289126v3⟩

Share

Metrics

Les métriques sont temporairement indisponibles