1、準備資料集

set.seed(2) #設定隨機數種子

x = matrix(rnorm(50*2),ncol=2) #生成50X2的正態分佈矩陣

x[1:25,1] = x[1:25,1] +3
x[1:25,2] = x[1:25,2] -4 #構建不同的聚類

2、建立聚類模型

km.out = kmeans(x,2,nstart=20)

km.out$cluster

輸出如下：

很明顯，前25條記錄分為了一類，後25條記錄分為了一類。

3、作圖並標註聚類中心

plot(x,col=(km.out$cluster+1)) #做出資料集的點圖

km.out$centers #列印聚類中心

points(km.out$centers[1,1],km.out$centers[1,2],pch=10,col="red",cex=2)
points(km.out$centers[2,1],km.out$centers[2,2],pch=10,col="blue",cex=2) #標註聚類中心
 

輸出如下：