apunkt/GPy
1
0
Fork 0
mirror of https://github.com/SheffieldML/GPy.git synced 2026-05-01 07:46:22 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

FITC example added

Browse source
This commit is contained in:
Ricardo 2013-06-05 16:30:57 +01:00
parent 7fb3466fba
commit d3f84816df
1 changed files with 10 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

13
GPy/examples/classification.py
View file

@ -114,7 +114,8 @@ def sparse_toy_linear_1d_classification(seed=default_seed):
return m
def sparse_crescent_data(inducing=10, seed=default_seed):
"""Run a Gaussian process classification on the crescent data. The demonstration calls the basic GP classification model and uses EP to approximate the likelihood.
"""
Run a Gaussian process classification with DTC approxiamtion on the crescent data. The demonstration calls the basic GP classification model and uses EP to approximate the likelihood.
:param model_type: type of model to fit ['Full', 'FITC', 'DTC'].
:param seed : seed value for data generation.
@ -137,7 +138,8 @@ def sparse_crescent_data(inducing=10, seed=default_seed):
return m
def FITC_crescent_data(inducing=10, seed=default_seed):
"""Run a Gaussian process classification on the crescent data. The demonstration calls the basic GP classification model and uses EP to approximate the likelihood.
"""
Run a Gaussian process classification with FITC approximation on the crescent data. The demonstration uses EP to approximate the likelihood.
:param model_type: type of model to fit ['Full', 'FITC', 'DTC'].
:param seed : seed value for data generation.
@ -146,13 +148,18 @@ def FITC_crescent_data(inducing=10, seed=default_seed):
:type inducing: int
"""
data = GPy.util.datasets.crescent_data(seed=seed)
Y = data['Y']
Y[Y.flatten()==-1]=0
data = GPy.util.datasets.crescent_data(seed=seed)
Y = data['Y']
Y[Y.flatten()==-1]=0
m = GPy.models.FITCClassification(data['X'], Y)
m.ensure_default_constraints()
m['.*len'] = 10.
m['.*len'] = 3.
m.update_likelihood_approximation()
m.optimize()
print(m)