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.
[1]
Radford M. Neal.
Probabilistic Inference Using Markov Chain Monte Carlo Methods
,
2011
.
[2]
A. P. Dawid,et al.
Regression and Classification Using Gaussian Process Priors
,
2009
.
[3]
J. Weston,et al.
Approximation Methods for Gaussian Process Regression
,
2007
.
[4]
Radford M. Neal.
Slice Sampling
,
2003,
The Annals of Statistics.
[5]
Radford M. Neal.
Monte Carlo Implementation of Gaussian Process Models for Bayesian Regression and Classification
,
1997,
physics/9701026.
[6]
Radford M. Neal.
Sampling from multimodal distributions using tempered transitions
,
1996,
Stat. Comput..
[7]
Petros Drineas,et al.
On the Nyström Method for Approximating a Gram Matrix for Improved Kernel-Based Learning
,
2005,
J. Mach. Learn. Res..
[8]
N. Metropolis,et al.
Equation of State Calculations by Fast Computing Machines
,
1953,
Resonance.