Skip to main content

Approximate implementation of the logarithm of the matrix determinant in Gaussian process regression

Buy Article:

$63.00 plus tax (Refund Policy)


Maximum likelihood estimation of hyperparameters in Gaussian processes (GPs) as well as other spatial regression models usually requires the evaluation of the logarithm of the matrix determinant, in short, log det. When using matrix decomposition techniques, the exact implementation of log det is of O(N3) operations, where N is the matrix dimension. In this paper, a power-series expansion-based framework is presented for approximating the log det of general positive-definite matrices. Three novel compensation schemes are proposed to further improve the approximation accuracy and computational efficiency. The proposed log det approximation requires only 50N2 operations. The theoretical analysis is substantiated by a large number of numerical experiments, including tests on randomly generated positive-definite matrices, randomly generated covariance matrices, and sequences of covariance matrices generated online in two GP regression examples. The average approximation error is ∼9%.

Keywords: Compensation; Gaussian process; Logarithm of matrix determinant; Power-series expansion

Document Type: Research Article


Affiliations: 1: Hamilton Institute, National University of Ireland, Ireland 2: Department of Electronic and Electrical Engineering, University of Strathclyde, Glasgow, UK

Publication date: 2007-01-01

  • Access Key
  • Free ContentFree content
  • Partial Free ContentPartial Free content
  • New ContentNew content
  • Open Access ContentOpen access content
  • Partial Open Access ContentPartial Open access content
  • Subscribed ContentSubscribed content
  • Partial Subscribed ContentPartial Subscribed content
  • Free Trial ContentFree trial content
Cookie Policy
Cookie Policy
Ingenta Connect website makes use of cookies so as to keep track of data that you have filled in. I am Happy with this Find out more