Abstract:Gaussian processes provide probabilistic surrogates for various applications including classification, uncertainty quantification, and optimization. Using a gradient-enhanced covariance matrix can be beneficial since it provides a more accurate surrogate relative to its gradient-free counterpart. An acute problem for Gaussian processes, particularly those that use gradients, is the ill-conditioning of their covariance matrices. Several methods have been developed to address this problem for gradient-enhanced Gaussian processes but they have various drawbacks such as limiting the data that can be used, imposing a minimum distance between evaluation points in the parameter space, or constraining the hyperparameters. In this paper a new method is presented that applies a diagonal preconditioner to the covariance matrix along with a modest nugget to ensure that the condition number of the covariance matrix is bounded, while avoiding the drawbacks listed above. Optimization results for a gradient-enhanced Bayesian optimizer with the Gaussian kernel are compared with the use of the new method, a baseline method that constrains the hyperparameters, and a rescaling method that increases the distance between evaluation points. The Bayesian optimizer with the new method converges the optimality, ie the $\ell_2$ norm of the gradient, an additional 5 to 9 orders of magnitude relative to when the baseline method is used and it does so in fewer iterations than with the rescaling method. The new method is available in the open source python library GpGradPy, which can be found at <a class="link-external link-https" href="https://github.com/marchildon/gpgradpy/tree/paper_precon" rel="external noopener nofollow">this https URL</a>. All of the figures in this paper can be reproduced with this library.

Optimization of the generalized covariance estimator in noncausal processes

Optimization-based Causal Estimation from Heterogenous Environments

Evaluation of combinatorial optimisation algorithms for c-optimal experimental designs with correlated observations

Efficient Algorithms for Covariate Analysis with Dynamic Data Using Nonlinear Mixed-Effects Model.

Generalized covariance‐based inference for models set‐identified from independence restrictions

Inferring Parsimonious Coupling Statistics in Nonlinear Dynamics with Variational Gaussian Processes

Causal Discovery of Linear Non-Gaussian Causal Models with Unobserved Confounding

High Precision Causal Model Evaluation with Conditional Randomization

Portfolio optimization with CVaR under VG process

MODE: Minimax Optimal Deterministic Experiments for Causal Inference in the Presence of Covariates

Testing Generalizability in Causal Inference

Relaxed covariate overlap and margin-based causal effect estimation

Causal Fine-Tuning and Effect Calibration of Non-Causal Predictive Models

A Solution to the Ill-Conditioning of Gradient-Enhanced Covariance Matrices for Gaussian Processes

Cross-Entropy Method Variants for Optimization

Matching as Nonparametric Preprocessing for Reducing Model Dependence in Parametric Causal Inference

High-Dimensional Causal Bayesian Optimization

Dynamic Large Spatial Covariance Matrix Estimation in Application to Semiparametric Model Construction via Variable Clustering: the SCE approach

Covariance estimation for multivariate conditionally Gaussian dynamic linear models

Optimization and Testing in Linear Non-Gaussian Component Analysis

Monte Carlo Integration with adaptive variance selection for improved stochastic Efficient Global Optimization