Improving Hyperparameter Learning under Approximate Inference in Gaussian Process Models