首先读入csv数据:
np <- read.table("D:/0R/np1cn.csv", header=TRUE, sep=",", row.names="乡镇")
np #显示刚才读入的数据
distance <- dist(np) #计算距离
np.hc <- hclust(distance) #聚类分析,最长距离法
plot(np.hc, hang = -1) #绘画系谱图
re <- rect.hclust(np.hc, k = 4) #分为4类
显示分组后的文本:
re
for (i in 1:4) {
print(paste("第",i,"类"))
print(np[re[[i]],])
}
另外几种分类方法:
np.hc <- hclust(distance,method="centroid") #聚类分析,重心法
np.hc <- hclust(distance,method="median") #聚类分析,中间距离法
np.hc <- hclust(distance,method="ward") #聚类分析,离差平方和法