Radford M. Neal, Dept. of Statistical Sciences and Dept. of Computer Science, University of Toronto

Gaussian Process (GP) models are a powerful and flexible tool for non-parametric regression and classification. Computation for GP models is intensive, since computing the posterior density, $\pi$, for covariance function parameters requires computation of the covariance matrix, $C$, a $pn^2$ operation, where $p$ is the number of covariates and $n$ is the number of training cases, and then inversion of $C$, an $n^3$ operation. We introduce MCMC methods based on the ``temporary mapping and caching'' framework, using a fast approximation, $\pi^*$, as the distribution needed to construct the temporary space. We propose two implementations under this scheme: ``mapping to a discretizing chain'', and ``mapping with tempered transitions'', both of which are exactly correct MCMC methods for sampling $\pi$, even though their transitions are constructed using an approximation. These methods are equivalent when their tuning parameters are set at the simplest values, but differ in general. We compare how well these methods work when using several approximations, finding on synthetic datasets that a $\pi^*$ based on the ``Subset of Data'' (SOD) method is almost always more efficient than standard MCMC using only $\pi$. On some datasets, a more sophisticated $\pi^*$ based on the ``Nystr\"om-Cholesky'' method works better than SOD.

Technical Report, 9 May 2013, 21 pages: pdf.

Also available from arXiv.org.