Fitting Gaussian Process Models in Python