Доступ предоставлен для: Guest
Портал Begell Электронная Бибилиотека e-Книги Журналы Справочники и Сборники статей Коллекции
International Journal for Uncertainty Quantification
Импакт фактор: 3.259 5-летний Импакт фактор: 2.547 SJR: 0.531 SNIP: 0.8 CiteScore™: 1.52

ISSN Печать: 2152-5080
ISSN Онлайн: 2152-5099

Свободный доступ

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

Краткое описание

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:

LOW-COST MULTI-DIMENSIONAL GAUSSIAN PROCESS WITH APPLICATION TO UNCERTAINTY QUANTIFICATION
International Journal for Uncertainty Quantification, Vol.5, 2015, issue 4
Guang Lin, Bledar A. Konomi
BIVARIATE QUANTILE INTERPOLATION FOR ENSEMBLE DERIVED PROBABILITY DENSITY ESTIMATES
International Journal for Uncertainty Quantification, Vol.5, 2015, issue 2
Brad Eric Hollister, Alex Pang
AN ENSEMBLE KALMAN FILTER USING THE CONJUGATE GRADIENT SAMPLER
International Journal for Uncertainty Quantification, Vol.3, 2013, issue 4
Heikki Haario, Antti Solonen, Albert Parker, Marylesa Howard, Johnathan M. Bardsley
GRADIENT-BASED STOCHASTIC OPTIMIZATION METHODS IN BAYESIAN EXPERIMENTAL DESIGN
International Journal for Uncertainty Quantification, Vol.4, 2014, issue 6
Youssef Marzouk, Xun Huan
CHARACTERIZING UNCERTAINTIES IN A THERMAL STORAGE PACKED BED USING BAYESIAN INFERENCEWITH MCMC AND VARIATIONAL APPROXIMATIONS
Second Thermal and Fluids Engineering Conference, Vol.50, 2017, issue
Ashley F. Emery, Sang Ien