Variational Optimisation for Non-conjugate Likelihood Gaussian Process Models