Inscrição na biblioteca: Guest
Portal Digital Begell Biblioteca digital da Begell eBooks Diários Referências e Anais Coleções de pesquisa
International Journal for Uncertainty Quantification
Fator do impacto: 3.259 FI de cinco anos: 2.547 SJR: 0.417 SNIP: 0.8 CiteScore™: 1.52

ISSN Imprimir: 2152-5080
ISSN On-line: 2152-5099

Open Access

International Journal for Uncertainty Quantification

DOI: 10.1615/Int.J.UncertaintyQuantification.2013006330
pages 133-149

GAUSSIAN PROCESS ADAPTIVE IMPORTANCE SAMPLING

Keith R. Dalbey
Department of Optimization and Uncertainty Quantification, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
Laura P. Swiler
Optimization and Uncertainty Quantification Department, Sandia National Laboratories, P.O. Box 5800, MS 1318, Albuquerque, New Mexico 87185, USA

RESUMO

The objective is to calculate the probability, PF, that a device will fail when its inputs, x, are randomly distributed with probability density, p (x), e.g., the probability that a device will fracture when subject to varying loads. Here failure is defined as some scalar function, y (x), exceeding a threshold, T. If evaluating y (x) via physical or numerical experiments is sufficiently expensive or PF is sufficiently small, then Monte Carlo (MC) methods to estimate PF will be unfeasible due to the large number of function evaluations required for a specified accuracy. Importance sampling (IS), i.e., preferentially sampling from "important" regions in the input space and appropriately down-weighting to obtain an unbiased estimate, is one approach to assess PF more efficiently. The inputs are sampled from an importance density, p' (x). We present an adaptive importance sampling (AIS) approach which endeavors to adaptively improve the estimate of the ideal importance density, p* (x), during the sampling process. Our approach uses a mixture of component probability densities that each approximate p* (x). An iterative process is used to construct the sequence of improving component probability densities. At each iteration, a Gaussian process (GP) surrogate is used to help identify areas in the space where failure is likely to occur. The GPs are not used to directly calculate the failure probability; they are only used to approximate the importance density. Thus, our Gaussian process adaptive importance sampling (GPAIS) algorithm overcomes limitations involving using a potentially inaccurate surrogate model directly in IS calculations. This robust GPAIS algorithm performs surprisingly well on a pathological test function.


Articles with similar content:

ADAPTIVE SAMPLING WITH TOPOLOGICAL SCORES
International Journal for Uncertainty Quantification, Vol.3, 2013, issue 2
Dan Maljovec, Valerio Pascucci, Bei Wang, Ana Kupresanin, Gardar Johannesson, Peer-Timo Bremer
A MULTI-INDEX MARKOV CHAIN MONTE CARLO METHOD
International Journal for Uncertainty Quantification, Vol.8, 2018, issue 1
Ajay Jasra, Yan Zhou, Kengo Kamatani, Kody J. H. Law
IMPROVEMENTS TO GRADIENT-ENHANCED KRIGING USING A BAYESIAN INTERPRETATION
International Journal for Uncertainty Quantification, Vol.4, 2014, issue 3
Hester Bijl, Richard P. Dwight, Jouke H.S. de Baar
WASSERSTEIN METRIC-DRIVEN BAYESIAN INVERSION WITH APPLICATIONS TO SIGNAL PROCESSING
International Journal for Uncertainty Quantification, Vol.9, 2019, issue 4
Mohammad Motamed, Daniel Appelo
FAST AND FLEXIBLE UNCERTAINTY QUANTIFICATION THROUGH A DATA-DRIVEN SURROGATE MODEL
International Journal for Uncertainty Quantification, Vol.8, 2018, issue 2
Gerta Köster, Hans-Joachim Bungartz, Felix Dietrich, Tobias Neckel, Florian Künzner