X = np.linspace(0,10)[:, None]
F = np.round(X*3-4)
F = np.where(F > 0, F, 0)
eps = np.random.randint(0,4, F.shape[0])[:, None]
Y = F + eps
noise_model = GPy.likelihoods.poisson()
After Change
Run a simple demonstration of a standard Gaussian process fitting it to data sampled from an RBF covariance.
x_len = 400
X = np.linspace(0, 10, x_len)[:, None]
f_true = np.random.multivariate_normal(np.zeros(x_len), GPy.kern.rbf(1).K(X))
Y = np.array([np.random.poisson(np.exp(f)) for f in f_true])[:,None]
noise_model = GPy.likelihoods.poisson()
