An Efficient Implementation of Riemannian Manifold Hamiltonian Monte Carlo for Gaussian Process Models

10/28/2018

∙

by   Ulrich Paquet, et al.

∙

0

∙

share

This technical report presents pseudo-code for a Riemannian manifold Hamiltonian Monte Carlo (RMHMC) method to efficiently simulate samples from N-dimensional posterior distributions p(x|y), where x ∈ R^N is drawn from a Gaussian Process (GP) prior, and observations y_n are independent given x_n. Sufficient technical and algorithmic details are provided for the implementation of RMHMC for distributions arising from GP priors.