Variational inference for robust Gaussian process regression