scikit-GPUPPY: Gaussian Process Uncertainty Propagation with PYthon

This package provides means for modeling functions and simulations using Gaussian processes (aka Kriging, Gaussian random fields, Gaussian random functions). Additionally, uncertainty can be propagated through the Gaussian processes.

Note

The Gaussian process regression and uncertainty propagation are based on Girard’s thesis [1].

An extension to speed up GP regression is based on Snelson’s thesis [2].

Warning

The extension based on Snelson’s work is already usable but not as fast as it should be. Additionally, the uncertainty propagation does not yet work with this extension.

An additional extension for Inverse Uncertainty Propagation is based on my paper (and upcoming PhD thesis) [3].

A simulation is seen as a function \(f(x)+\epsilon\) (\(x \in \mathbb{R}^n\)) with additional random error \(\epsilon \sim \mathcal{N}(0,v)\). This optional error is due to the stochastic nature of most simulations.

The GaussianProcess module uses regression to model the simulation as a Gaussian process. (See Gaussian Process Regression for an explanation)

The UncertaintyPropagation module allows for propagating uncertainty \(x \sim \mathcal{N}(\mu,\Sigma)\) through the Gaussian process to estimate the output uncertainty of the simulation. (See Uncertainty Propagation for an explanation)

The FFNI and TaylorPropagation modules provide classes for propagating uncertainty through deterministic functions.

The InverseUncertaintyPropagation module allows for propagating the desired output uncertainty of the simulation backwards through the Gaussian Process. This assumes that the components of the input \(x\) are estimated from samples using maximum likelihood estimators. Then, the inverse uncertainty propagation calculates the optimal sample sizes for estimating \(x\) that lead to the desired output uncertainty of the simulation. (See Inverse Uncertainty Propagation for an explanation)

Indices and tables

References

[1]Girard, A. Approximate Methods for Propagation of Uncertainty with Gaussian Process Models, University of Glasgow, 2004
[2]Snelson, E. L. Flexible and efficient Gaussian process models for machine learning, Gatsby Computational Neuroscience Unit, University College London, 2007
[3]Baumgaertel, P.; Endler, G.; Wahl, A. M. & Lenz, R. Inverse Uncertainty Propagation for Demand Driven Data Acquisition, Proceedings of the 2014 Winter Simulation Conference, IEEE Press, 2014, 710-721