Compute and Visualize Profile Extrema Functions

Computes profile extrema functions for arbitrary functions. If the function is expensive-to-evaluate it computes profile extrema by emulating the function with a Gaussian process (using package 'DiceKriging'). In this case uncertainty quantification on the profile extrema can also be computed. The different plotting functions for profile extrema give the user a tool to better locate excursion sets.


Computation and plots of profile extrema functions. The package main functions are:

  • coordinateProfiles: Given a km objects computes the coordinate profile extrema function for the posterior mean and its quantiles.

  • coordProf_UQ: UQ part of coordinateProfiles.

  • obliqueProfiles: Given a km objects computes the profile extrema functions for a generic list of matrices Psi for the posterior mean and its quantiles.

  • obliqueProf_UQ: The UQ part of obliqueProfiles.

  • getAllMaxMin: computes coordinate profile extrema with full optimization for a deterministic function.

  • approxMaxMin: approximates coordinate profile extrema for a deterministic function.

  • getProfileExtrema: computes profile extrema given a list of matrices Psi for a deterministic function.

  • approxProfileExtrema: approximates profile extrema given a list of matrices Psi for a deterministic function.

Plotting:

  • plot_univariate_profiles_UQ: plots for the results of coordProf_UQ or obliqueProf_UQ. Note that this function only works for univariate profiles.

  • plotBivariateProfiles: plots the bivariate maps results of a call to obliqueProfiles with a two dimensional projection matrix Psi.

  • plotMaxMin: simple plotting function for univariate profile extrema.

  • plotOneBivProfile: simple plotting function for bivariate profile extrema.

Note

This work was supported in part the Hasler Foundation, grant number 16065 and by the Swiss National Science Foundation, grant number 167199. The author warmly thanks David Ginsbourger, Jérémy Rohmer and Déborah Idier for fruitful discussions and accurate, thought provoking suggestions.

References

Azzimonti, D., Bect, J., Chevalier, C., and Ginsbourger, D. (2016). Quantifying uncertainties on excursion sets under a Gaussian random field prior. SIAM/ASA Journal on Uncertainty Quantification, 4(1):850–874.

Azzimonti, D., Ginsbourger, D., Rohmer, J. and Idier, D. (2017+). Profile extrema for visualizing and quantifying uncertainties on excursion regions. Application to coastal flooding. arXiv:1710.00688.

Chevalier, C. (2013). Fast uncertainty reduction strategies relying on Gaussian process models. PhD thesis, University of Bern.

Chevalier, C., Picheny, V., Ginsbourger, D. (2014). An efficient and user-friendly implementation of batch-sequential inversion strategies based on kriging. Computational Statistics & Data Analysis, 71: 1021-1034.

Johnson, S. G. The NLopt nonlinear-optimization package, http://ab-initio.mit.edu/nlopt

Koenker, R. (2017). quantreg: Quantile Regression. R package version 5.33.

Nocedal, J. and Wright, S. J. (2006). Numerical Optimization, second edition. Springer- Verlag, New York.

Neuwirth, E. (2014). RColorBrewer: ColorBrewer Palettes. R package version 1.1-2.

Roustant, O., Ginsbourger, D., Deville, Y. (2012). DiceKriging, DiceOptim: Two R Packages for the Analysis of Computer Experiments by Kriging-Based Metamodeling and Optimization. Journal of Statistical Software, 51(1): 1-55.

News

profExtrema 0.2.0

Changes with respect to profExtrema 0.1.0

  • New plotting functions such as plotBivariateProfiles and plotOneBivProfile give a more principled way to visualize bivariate profile extrema functions.

  • updated univariate profile extrema plotting functions plot_univariate_profiles_UQ for more informative plots.

  • non-coordinate profile extrema functions are now computed with a faster null-space method.

  • Fixed many minor errors and bugs.

Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.

install.packages("profExtrema")

0.2.0 by Dario Azzimonti, a year ago


Browse source code at https://github.com/cran/profExtrema


Authors: Dario Azzimonti [aut, cre, cph]


Documentation:   PDF Manual  


GPL-3 license


Imports microbenchmark, quantreg, lhs, splines, methods, RColorBrewer, utils, MASS, rcdd

Depends on DiceKriging, KrigInv, pGPx


See at CRAN