stick_bgplvm to original version with rbf

This commit is contained in:
Andreas 2013-07-17 18:18:00 +01:00
parent f4794796a5
commit 0a4deeefcd

View file

@ -378,11 +378,10 @@ def stick():
def stick_bgplvm(model=None): def stick_bgplvm(model=None):
data = GPy.util.datasets.stick() data = GPy.util.datasets.stick()
Q = 6 Q = 6
kernel = GPy.kern.rbf_inv(Q, ARD=True) + GPy.kern.bias(Q, np.exp(-2)) + GPy.kern.white(Q, np.exp(-2)) kernel = GPy.kern.rbf(Q, ARD=True) + GPy.kern.bias(Q, np.exp(-2)) + GPy.kern.white(Q, np.exp(-2))
m = BayesianGPLVM(data['Y'], Q, init="PCA", num_inducing=35,kernel=kernel) m = BayesianGPLVM(data['Y'], Q, init="PCA", num_inducing=20,kernel=kernel)
# optimize # optimize
m.ensure_default_constraints() m.ensure_default_constraints()
m.constrain_bounded('.*rbf_inv',1e-5, 100)
m.optimize(messages=1, max_iters=3000,xtol=1e-300,ftol=1e-300) m.optimize(messages=1, max_iters=3000,xtol=1e-300,ftol=1e-300)
m._set_params(m._get_params()) m._set_params(m._get_params())
plt.clf, (latent_axes, sense_axes) = plt.subplots(1, 2) plt.clf, (latent_axes, sense_axes) = plt.subplots(1, 2)