将 scikit-learn 0.17 版代码翻译成 0.20 版答案

【问题标题】：Translating scikit-learn version 0.17 code to verion 0.20将 scikit-learn 0.17 版代码翻译成 0.20 版
【发布时间】：2019-04-11 09:15:01
【问题描述】：

您好，我想知道是否有人知道如何将以下 scikit-learn 0.17 版代码转换为 0.20 版

from sklearn.gaussian_process import GaussianProcess
gp = GaussianProcess(corr='cubic', theta0=1e-2, thetaL=1e-4, thetaU=1E-1,
                     random_start=100)
xfit = np.linspace(0, 10, 1000)
yfit, MSE = gp.predict(xfit[:, np.newaxis], eval_MSE=True)

提前谢谢你！

【问题讨论】：

标签： scikit-learn gaussian

【解决方案1】：

我偶然发现了书中相同的 sn-p 代码 Jake VanderPlas 的“Python 数据科学手册”。我这样做了，

    from sklearn.gaussian_process import GaussianProcessRegressor
    from sklearn.gaussian_process.kernels import RBF, ConstantKernel as C
    
    kernel = C(1.0, (1e-3, 1e3)) * RBF(10, (1e-2, 1e2))
    gp = GaussianProcessRegressor(kernel=kernel, n_restarts_optimizer=10)
    gp.fit(xdata[:,np.newaxis],ydata)
    xfit = np.linspace(0,10,1000)
    yfit, sigma = gp.predict(xfit[:,np.newaxis], return_std=True)
    dyfit = 2*sigma

它和书中的不一样，但它有效，我在官方网站 scikit-learn.org 上关注了这个example。

【讨论】：

【解决方案2】：

GaussianProcess 现在在 GaussianProcessRegressor 和 GaussianProcessClassifier 之间拆分

【讨论】：

是的，我非常清楚这一点，问题是您需要为协方差矩阵选择一个内核，上面的示例显示内核 = corr=cubic，但没有等效内核在gaussian_process.kernels. 的可用选项中。