This project implements a Bayesian Gaussian Process regression model which is fit using an adaptive Metropolis within Gibbs algorithm.