Sparse Gaussian Process Regression Beyond Variational Inference