Efficient Marginal Likelihood Computation for Gaussian Process Regression

Details Efficient Marginal Likelihood Computation for Gaussian Process Regression Efficient Marginal Likelihood Computation for Gaussian Process Regression

2011-10-29

arxiv.org/abs/1110.6546v1

Statistics

Machine Learning

20 pages, 3 figures

Scientific paper

In a Bayesian learning setting, the posterior distribution of a predictive model arises from a trade-off between its prior distribution and the conditional likelihood of observed data. Such distribution functions usually rely on additional hyperparameters which need to be tuned in order to achieve optimum predictive performance; this operation can be efficiently performed in an Empirical Bayes fashion by maximizing the posterior marginal likelihood of the observed data. Since the score function of this optimization problem is in general characterized by the presence of local optima, it is necessary to resort to global optimization strategies, which require a large number of function evaluations. Given that the evaluation is usually computationally intensive and badly scaled with respect to the dataset size, the maximum number of observations that can be treated simultaneously is quite limited. In this paper, we consider the case of hyperparameter tuning in Gaussian process regression. A straightforward implementation of the posterior log-likelihood for this model requires O(N^3) operations for every iteration of the optimization procedure, where N is the number of examples in the input dataset. We derive a novel set of identities that allow, after an initial overhead of O(N^3), the evaluation of the score function, as well as the Jacobian and Hessian matrices, in O(N) operations. We prove how the proposed identities, that follow from the eigendecomposition of the kernel matrix, yield a reduction of several orders of magnitude in the computation time for the hyperparameter optimization problem. Notably, the proposed solution provides computational advantages even with respect to state of the art approximations that rely on sparse kernel matrices.

Affiliated with

Also associated with

No associations

Rating

Efficient Marginal Likelihood Computation for Gaussian Process Regression does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Efficient Marginal Likelihood Computation for Gaussian Process Regression, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Efficient Marginal Likelihood Computation for Gaussian Process Regression will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-200544