gp regression.py-20170902

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Fri Sep  1 16:27:03 2017

@author: vicky
"""

import numpy as np  
import matplotlib.pyplot as plt  
from sklearn.gaussian_process import GaussianProcessRegressor as GPR
from sklearn.gaussian_process.kernels import RBF, ConstantKernel as C  
from mpl_toolkits.mplot3d import Axes3D  

pylab
pb.ion()

k = GPy.kern.rbf(d,var,theta)

 X = np.linspace(0.05,0.95,10)[:,None]
 Y = -np.cos(np.pi*X) +np.sin(4*np.pi*X) + np.random.randn(10,1)*0.05
 pb.figure()
 pb.plot(X,Y,'kx',mew=1.5)
 
 k = GPy.kern.rbf(input_dim=1, variance=1., lengthscale=1.)
 m = GPy.models.GPRegression(X,Y,k)