Hostname: page-component-89b8bd64d-72crv Total loading time: 0 Render date: 2026-05-08T23:50:24.891Z 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
Get access

Abstract

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

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

Article purchase

Temporarily unavailable