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