From 217fa0e70eaad1eecc3ef77f541a03435ad7ef50 Mon Sep 17 00:00:00 2001 From: Ricardo Andrade Date: Tue, 29 Jan 2013 16:44:12 +0000 Subject: [PATCH] Now it works. --- GPy/examples/poisson.py | 10 ++++------ 1 file changed, 4 insertions(+), 6 deletions(-) diff --git a/GPy/examples/poisson.py b/GPy/examples/poisson.py index 71d80b30..e15f310d 100644 --- a/GPy/examples/poisson.py +++ b/GPy/examples/poisson.py @@ -25,16 +25,14 @@ seed=default_seed """ X = np.arange(0,100,5)[:,None] -F = np.round(np.sin(X/18.) + .1*X) -E = np.random.randint(-3,3,20)[:,None] +F = np.round(np.sin(X/18.) + .1*X) + np.arange(5,25)[:,None] +E = np.random.randint(-5,5,20)[:,None] Y = F + E -pb.plot(X,F,'k-') -pb.plot(X,Y,'ro') pb.figure() -likelihood = GPy.inference.likelihoods.poisson(Y,scale=6.) +likelihood = GPy.inference.likelihoods.poisson(Y,scale=1.) m = GPy.models.GP(X,likelihood=likelihood) -#m = GPy.models.GP(data['X'],Y=likelihood.Y) +#m = GPy.models.GP(X,Y=likelihood.Y) m.constrain_positive('var') m.constrain_positive('len')