【发布时间】:2018-08-10 15:15:35
【问题描述】:
我对 python 很陌生,对做高斯回归很感兴趣。 我在 py3.6 和 SKlearn 0.19 下。
我有简单的代码,但我得到一个关于 cdist 中由 predict 调用的向量维度的错误。我知道我的输入有问题。但我不明白为什么......
我查看了高斯过程回归器的示例,但它似乎不是最常用的工具。
提前感谢您的帮助。
干杯。
这是我的代码示例:
import pandas as pd
import numpy as np
import numpy as np
from sklearn.gaussian_process import GaussianProcessRegressor as gpr
from sklearn.gaussian_process.kernels import RBF, ConstantKernel as C
....
#X_train are the training samples
X_train= np.column_stack((xc,yc,zc))
print('X_train')
print(X_train.shape)
print(X_train)
这是 X_train 的打印:
X_train (4576, 3)
[[ 0.71958336 -1.12719598 0.47889958]
[ 0.71958336 -1.12719598 0.47889958]
[ 0.71958336 -1.12719598 0.34285071]
...
[ 0.55255508 -1.18817547 -1.63666023]
[ 0.55255508 -1.18817547 -1.70468466]
[ 0.55255508 -1.18817547 -1.77270909]]
这是训练中的目标特征:
print('v1')
print(v1.shape)
print(v1)
它的印刷品
v1
(4576,)
0 10.0
1 14.0
2 13.0
3 19.0
....
4573 39.0
4574 16.0
4575 12.0
这是要预测的样本:
x = np.column_stack((xp,
yp,
zp))
print('x')
print(x.shape)
print(x)
这是印刷品:
x (75, 3) [[-1.41421356 -1.41421356 -1.22474487] [-0.70710678 -1.41421356 -1.22474487] [ 0. -1.41421356 -1.22474487] [ 0.70710678 -1.41421356 -1.22474487] ..... [ 0.70710678 -0.70710678 -1.22474487] [ 1.41421356 -0.70710678 -1.22474487] [-1.41421356 0. -1.22474487] [-0.70710678 0. -1.22474487] [ 0. 0. -1.22474487]
这里是拟合和预测
v1 = v1.ravel()
#default kernel
kernel = C(1.0, (1e-3, 1e3)) * RBF(10, (1e-2, 1e2))
X_train, v1 = make_regression()
model = gpr(kernel=kernel, n_restarts_optimizer=9)
model.fit(X_train,v1)
#Predict v1
v1_pred = model.predict(x)
运行时出现以下错误:
文件“test.py”,第 189 行,在测试中 v1_pred = model.predict(x) 文件“/usr/local/lib/python3.6/site-packages/sklearn/gaussian_process/gpr.py”, 第 315 行,在预测中 K_trans = self.kernel_(X, self.X_train_) 文件“/usr/local/lib/python3.6/site-packages/sklearn/gaussian_process/kernels.py”, 第 758 行,在 调用 返回 self.k1(X, Y) * self.k2(X, Y) 文件“/usr/local/lib/python3.6/site-packages/sklearn/gaussian_process/kernels.py”, 第 1215 行,在 call 中 metric='sqeuclidean') 文件 "/usr/local/lib/python3.6/site-packages/scipy/spatial/distance.py", 第 2373 行,在 cdist 中 raise ValueError('XA and XB must have the same number of columns ' ValueError: XA and XB must have the same number of columns (i.e. 特征维度。)
【问题讨论】:
-
xp、yp和zp是什么? -
抱歉,可能不清楚。为了解释我的情况,我对 3D 空间中的一个组件进行了化学分析。 xc, yc, zc 是我的训练样本上的变量,我知道 3D 空间中 xc,yc,zc 位置的化学分析值。 xp, yp, zp 是相同的,但在数据中我想预测未知的化学分析。简而言之,这些是笛卡尔坐标,应该与化学分析相关联
-
xp, yp, zp 和 xc,yc,zx 是:
type <class 'numpy.ndarray'>
标签: python pandas numpy machine-learning scikit-learn