全面了解R语言中的k-means如何聚类？

下面将在iris数据集上演示k-means聚类的过程。
先从iris数据集中移除Species属性，然后再对数据集iris调用函数kmeans，并将聚类结果存储在变kmeans.result中。
在下面的代码中，簇的数目设置为3。
iris2 <- iris
iris2$Species <- NULL
(kmeans.result <- kmeans(iris2, 3))  


将聚类结果与类标号(Species)进行比较，查看相似的对象是否被划分到同一个簇中。
table(iris$Species, kmeans.result$cluster)
 从上面的聚类结果可以看出，setosa类与其他的两个类很容易就划分开了，而versicolor类与virginica类之间存在小范围的重叠。


k-means聚类的结果图


然后，绘制所有的簇以及簇中心，如图所示。
注意，数据集有4个维度，而绘图只使用了前两个维度的数据。



图中所示的一些靠近绿色中心(用星号标记)的黑点实际上在4维空间中更靠近黑色中心。
要注意，多次运行得到的k-means聚类结果可能不同，这是因为初始的簇中心是随机选择的。
plot(iris2[c("Sepal.Length", "Sepal.Width")], col = kmeans.result$cluster)
# plot cluster centers

points(kmeans.result$centers[,c("Sepal.Length", "Sepal.Width")], col = 1:3, pch = 8, cex=2)