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.

Holonomic Gradient Method for Two Way Contingency Tables

Expectile regression for analyzing heteroscedasticity in high dimension

Semiparametric Expectile Regression for High-dimensional Heavy-tailed and Heterogeneous Data

Pathwise Derivatives for Multivariate Distributions

Analysis Of A Sequence Of Dependent 2 X 2 Tables

An Application of the Holonomic Gradient Method to the Neural Tangent Kernel

Two-Way Contingency Tables under Conditional Hot Deck Imputation

Gradient and Hessian of Joint Probability Function with Applications on Chance-Constrained Programs

BOOSTED NONPARAMETRIC HAZARDS WITH TIME-DEPENDENT COVARIATES

A Nonparametric Measure of Local Association for two-way Contingency Tables

A Method Expanding 2 by 2 Contingency Table by Obtaining Tendencies of Boolean Operators: Boolean Monte Carlo Method

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

Random Transpositions on Contingency Tables

Multinomial Logistic Regression: Asymptotic Normality on Null Covariates in High-Dimensions

Two-Sample Test for Sparse High Dimensional Multinomial Distributions

Using a monotonic density ratio model to increase the power of the goodness-of-fit test for logistic regression models with case-control data

Beyond Log-Concavity: Theory and Algorithm for Sum-Log-Concave Optimization

The Geometry of Statistical Models for Two-Way Contingency Tables with Fixed Odds Ratios

A Cramér Moderate Deviation Theorem for Hotelling’s $T^{2}$-Statistic with Applications to Global Tests

Technical Note-Knowledge Gradient for Selection with Covariates: Consistency and Computation

Adaptive generalized conditional gradient method for multiobjective optimization