1

在 Python 中，我试图深入研究用于估计高斯过程模型的GPy库，但在早期的简单绘图中遇到了一个绊脚石。

对于我的数据，我生成了一个简单的正弦波，在中途添加了平方增长率，GPy 成功估计了初始模型。

数据生成：

## Generating data for regression
# First, regular sine wave + normal noise
x = np.linspace(0,40, num=300)
noise1 = np.random.normal(0,0.3,300)
y = np.sin(x) + noise1

# Second, an upward trending starting midway, with its own noise as well
temp = x[150:]
noise2 = 0.004*temp**2 + np.random.normal(0,0.1,150)
y[150:] = y[150:] + noise2

plt.plot(x, y)

初始型号：

## Pre-processing
X = np.expand_dims(x, axis=1)
Y = np.expand_dims(y, axis=1)

## Model
kernel = GPy.kern.RBF(input_dim=1, variance=1., lengthscale=1.)
model1 = GPy.models.GPRegression(X, Y, kernel)

## Plotting
fig = model1.plot()
GPy.plotting.show(fig, filename='basic_gp_regression_notebook')

在此处输入图像描述

但是，此模型指定错误，因为数据仅使用 sin(X) 和 X^2 创建，而不仅仅是 X，因此我创建了下一个模型：

X_all = np.hstack((np.sin(X), np.square(X)))

model2 = GPy.models.GPRegression(X_all, Y, kernel)

fig = model2.plot()
GPy.plotting.show(fig, filename='basic_correct_gp_regression_notebook')

但是，现在，我得到了绘图错误，

Invalid value of type 'builtins.str' received for the 'size' property of scatter.marker Received value: '5'

我认为这是因为绘图不知道使用“X”作为 x 轴，只提供了 sin(X) 和 X^2。

我该如何解决这个问题？

4

0 回答 0