Scalable Approximate Inference And Model Selection In Gaussian Process Regression