第三章 聚类分析

load("D:\\文件\\2022下（大三上）\\classR\\health.RData")
head(health)

# ------------第三章 聚类分析
#-------------系统聚类法---------------
x=scale(health)                     #1
d=dist(x,method="euclidean")        #2
hc=hclust(d,method="complete")      #3
plot(hc,hang=-1)                    #4

plot(hc$height[30:1],type="o")      #5

plot(hc,hang=-1)                    #6
h=rect.hclust(hc,k=3)              
print(h)     #聚类结果1             #7

n.hc=cutree(hc,k=3)    # 聚类结果2
n.hc
sort(n.hc)            # 聚类结果3

colMeans(x[n.hc==1,])               #8
aggregate(x,by=list(n.hc),mean) # 各类各个变量的均值

#-------------K均值聚类法---------------

x=scale(health)                      #1

install.packages(("factoextra"))     #2
library(factoextra)
fviz_nbclust(x,kmeans,method="wss")

km<-kmeans(scale(health),4)          #3
km
km$cluster
km$centers
 
fviz_cluster(km,data=x)              #4