Questions about GPy. GPy is a Gaussian Process (GP) framework written in Python. It includes support for basic GP regression, multiple output GPs (using coregionalization), various noise models, sparse GPs, non-parametric regression and latent variables. Use with the [python] tag
Questions tagged [gpy]
32 questions
6
votes
1 answer
Using GPy Multiple-output coregionalized prediction
I have been facing a problem recently where I believe that a multiple-output GP might be a good candidate. I am at the moment applying a single-output GP to my data and as dimensionality increases, my results keep getting worse. I have tried…

fnaos
- 151
- 1
- 12
4
votes
1 answer
Cannot find plot function in GPy library (python)
I am using the GPy library in Python 2.7 to perform Gaussian Process regressions. I started by following the tutorial notebooks provided in the GitHub page.
Sample code :
import numpy as np
import matplotlib.pyplot as plt
f = lambda x :…

saleml
- 105
- 7
3
votes
1 answer
Most significant input dimensions for GPy.GPCoregionalizedRegression?
I have trained successfully a multi-output Gaussian Process model using an GPy.models.GPCoregionalizedRegression model of the GPy package. The model has ~25 inputs and 6 outputs.
The underlying kernel is an GPy.util.multioutput.ICM kernel consisting…

gehbiszumeis
- 3,525
- 4
- 24
- 41
3
votes
2 answers
How to Save/Load Optimized GPy Regression Model
I'm trying to save my optimized Gaussian process model for use in a different script. My current line of thinking is to store the model information in a json file, utilizing GPy's built-in to_dict and from_dict functions. Something along the lines…

Bocephus85
- 55
- 3
3
votes
1 answer
R6034 An application has made an attempt to load the C runtime library incorrectly
I am running a python 2.7 code (containing GPy and GPyOpt, python implementation of gaussian process and Bayesian optimization) from Matlab on Anaconda 64bit on Windows 10 and I am facing with the following error:
warning in stationary: failed to…

Rayan
- 31
- 3
3
votes
1 answer
Get confidence limits from Gaussian Process model in Python GPy
I calcualted a Gaussian Process model in Python using GPy:
ker0 = GPy.kern.Bias(input_dim=1,variance=1e-2)
...
m = GPy.models.GPRegression(x, y, ker0+ker2)
I can plot it with
m.plot()
plt.show
and it visualizes the points, the spline and the…

horseshoe
- 1,437
- 14
- 42
2
votes
1 answer
LinAlgError: not positive definite, even with jitter. When using a conda environment instead of pip
I am trying to fit some random data to a GP with the RBF kernel, using the GPy package. When I change the active dimensions, I get the LinAlgError: not positive definite, even with jitter error. This error is generated only with a conda environment.…

Katerina
- 21
- 3
2
votes
1 answer
How to combine plot of two GPy models?
I have calculated two GP regression models and would like to have them plotted in the same figure.
Model 1
kernel = GPy.kern.RBF(input_dim=1, variance=.1, lengthscale=1.)
m1 = GPy.models.GPRegression(xa, ya,kernel)
m1.optimize_restarts(num_restarts…

hH1sG0n3
- 121
- 4
2
votes
3 answers
How to find and uninstall the numpy duplicate versions
I am trying to install the library GPy. Although the installation is successful, I have a question on my numpy version.
GPy library can be found here https://github.com/SheffieldML/GPy
The current version of my numpy is 1.9.3
>>> import numpy
>>>…

Vinod Prime
- 371
- 1
- 3
- 13
1
vote
0 answers
gpytorch, regression on targets and classification of gradients to negative or positive
i would like to set up the following model in GPYtorch:
i have 4 inputs and i want to predict an output (regression)
at the same time, i want to constrain the gradients of 3 inputs to be positive and of 1 input to be negative (with respect to the…

john
- 373
- 1
- 16
1
vote
0 answers
Prediction with GPytorch DeepGP was poor
[Link to DeepGP definition:] (https://drive.google.com/file/d/1UwU9fz6vqxOQ3F0NSKKmRSgiP6pNnbwo/view?usp=sharing)
I have defined a DeepGP model using GPytorch, trained it, and the predictions are made as:
test_dataset = TensorDataset(test_x,…

Bikal
- 11
- 1
1
vote
0 answers
GP regression using Poisson likelihood
I am trying to implement GP regression using Poisson likelihood.
I followed the example in GPy by doing
poisson_likelihood = GPy.likelihoods.Poisson()
laplace_inf = GPy.inference.latent_function_inference.Laplace()
m = GPy.core.GP(X=X, Y=Y,…

Abhijith
- 11
- 1
1
vote
1 answer
Multitask/multioutput GPy Coregionalized Regression with non-Gaussian Likelihood and Laplace inference function
I want to perform coregionalized regression in GPy, however I am using a Bernoulli likelihood and then to estimate that as a Gaussian, I use Laplace inference. The code below shows how I would usually run a single-output GP with this set up (with my…

gregory
- 188
- 2
- 17
1
vote
2 answers
How to reproduce results of GPy GPRegression using scikit-learn GaussianProcessRegressor?
Both GPRegression (GPy) and GaussianProcessRegressor (scikit-learn) uses similar initial values and the same optimizer (lbfgs). Why results vary significantly?
#!pip -qq install pods
#!pip -qq install GPy
from sklearn.gaussian_process import…

Zeel B Patel
- 691
- 5
- 16
1
vote
0 answers
Making sense of the standard deviation in GPy for different scaling schemes of training data?
I am trying to make predictions for two different outputs (independent) using two separate models.
I am normalising the input and output data using relevant factors from the physics of my problem. But the uncertainty (standard deviation) in the…

cacophony94
- 11
- 1