在TruncatedSVD Python之后绘制K-means集群

Question

我正在尝试绘制在我的数据集上运行群集的结果，但我收到错误：

  File "cluster.py", line 93, in <module>
    Z = kmeans.predict(np.c_[xx.ravel(), yy.ravel()])
  File "/usr/local/lib/python2.7/dist-packages/sklearn/cluster/k_means_.py", line 957, in predict
    X = self._check_test_data(X)
  File "/usr/local/lib/python2.7/dist-packages/sklearn/cluster/k_means_.py", line 867, in _check_test_data
    n_features, expected_n_features))
ValueError: Incorrect number of features. Got 2 features, expected 73122

我对fit()调用工作正常，但绘图是出错的地方。

这是我的代码：

reduced_data = TruncatedSVD(n_components=2).fit_transform(X)

kmeans = KMeans(n_clusters=4, init='k-means++', max_iter=100, n_init=1, verbose=False)
kmeans.fit(X)

h = .02     # point in the mesh [x_min, x_max]x[y_min, y_max].

# Plot the decision boundary. For that, we will assign a color to each
x_min, x_max = reduced_data[:, 0].min() - 1, reduced_data[:, 0].max() + 1
y_min, y_max = reduced_data[:, 1].min() - 1, reduced_data[:, 1].max() + 1
xx, yy = np.meshgrid(np.arange(x_min, x_max, h), np.arange(y_min, y_max, h))

# Obtain labels for each point in mesh. Use last trained model.
Z = kmeans.predict(np.c_[xx.ravel(), yy.ravel()])

# Put the result into a color plot
Z = Z.reshape(xx.shape)
plt.figure(1)
plt.clf()
plt.imshow(Z, interpolation='nearest',
           extent=(xx.min(), xx.max(), yy.min(), yy.max()),
           cmap=plt.cm.Paired,
           aspect='auto', origin='lower')

plt.plot(reduced_data[:, 0], reduced_data[:, 1], 'k.', markersize=2)
# Plot the centroids as a white X
centroids = kmeans.cluster_centers_
plt.scatter(centroids[:, 0], centroids[:, 1],
            marker='x', s=169, linewidths=3,
            color='w', zorder=10)
plt.title('K-means clustering on the digits dataset (PCA-reduced data)\n'
          'Centroids are marked with white cross')
plt.xlim(x_min, x_max)
plt.ylim(y_min, y_max)
plt.xticks(())
plt.yticks(())
plt.show()

任何人都可以建议我如何更改我的代码以获得群集图表？

Answer 1

回溯告诉你问题是什么：

ValueError：功能数量不正确。 有2个功能，预计73122

该kmeans分类器适合与73122维训练样本，因此不能使用kmeans来作出预测， 2维的测试样品。

要修复代码，只需将kmeans.fit(X)更改为kmeans.fit(reduced_data) 。