Laplace approximation and the natural gradient for Gaussian process regression with the heteroscedastic Student-t model