728x90

dist(data, method)

mahalanobis(data, center, cov)

hclust(data, method)

agnes(data, metric, method)

daisy(data, metric)

 

#유클리드 거리 측정
US.dist_euclidean = dist(USArrests, "euclidean")
US.dist_euclidean

#맨하탄 거리 측정
US.dist_manhattan = dist(USArrests, "manhattan")
US.dist_manhattan

#마할라노비스 거리 측정
mahalanobis(USArrest, colMeans(USArrests), cov(USArrests))

#계층적 군집 분석
US.single = hclust(US.dist_euclidean^2, method="single") #최단거리법
plot(US.single)

US.complete = hclust(US.dist_euclidean^2, method="complete") #최장거리법
plot(US.complete)

#군집 분석을 통한 그룹 확인
group = cutree(US.single, k=6) #6개의 그룹으로 분할
group

rect.hclust(US.single, k=6, border="blue")

 

728x90