Hostname: page-component-74656f9884-j7cmm Total loading time: 0 Render date: 2026-11-08T00:05:05.937Z Has data issue: false hasContentIssue false

Optimal uncertainty quantification for legacy data observationsof Lipschitz functions

Published online by Cambridge University Press:  30 August 2013

T.J. Sullivan
Affiliation:
Mathematics Institute, University of Warwick, Coventry, CV4 7AL, UK.. Tim.Sullivan@warwick.ac.uk
M. McKerns
Affiliation:
Center for Advanced Computing Research, California Institute of Technology, 1200 East California Boulevard, Mail Code 158-79, Pasadena, CA 91125, USA.; mmckerns@caltech.edu
D. Meyer
Affiliation:
Mathematics Institute, University of Warwick, Coventry, CV4 7AL, UK.; f.theil@warwick.ac.uk
F. Theil
Affiliation:
Mathematics Institute, University of Warwick, Coventry, CV4 7AL, UK.; f.theil@warwick.ac.uk
H. Owhadi
Affiliation:
Applied & Computational Mathematics and Control & Dynamical Systems, California Institute of Technology, Mail Code 9-94, 1200 East California Boulevard, Pasadena, CA 91125, USA.; owhadi@caltech.edu
M. Ortiz
Affiliation:
Division of Engineering and Applied Science, California Institute of Technology, Mail Code 105-50, 1200 East California Boulevard, Pasadena, CA 91125, USA.; ortiz@caltech.edu

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the 'Save PDF' action button.

We consider the problem of providing optimal uncertainty quantification (UQ) – and hencerigorous certification – for partially-observed functions. We present a UQ frameworkwithin which the observations may be small or large in number, and need not carryinformation about the probability distribution of the system in operation. The UQobjectives are posed as optimization problems, the solutions of which are optimal boundson the quantities of interest; we consider two typical settings, namely parametersensitivities (McDiarmid diameters) and output deviation (or failure) probabilities. Thesolutions of these optimization problems depend non-trivially (even non-monotonically anddiscontinuously) upon the specified legacy data. Furthermore, the extreme values are oftendetermined by only a few members of the data set; in our principal physically-motivatedexample, the bounds are determined by just 2 out of 32 data points, and the remaindercarry no information and could be neglected without changing the final answer. We proposean analogue of the simplex algorithm from linear programming that uses these observationsto offer efficient and rigorous UQ for high-dimensional systems with high-cardinalitylegacy data. These findings suggest natural methods for selecting optimal (maximallyinformative) next experiments.

Information

Type
Research Article
Copyright
© EDP Sciences, SMAI, 2013