聚類 - 繪製最大的n個簇

Question

我使用HOPACH聚類 - 是否有一種方法可視化只有最大的n簇（例如，只繪製3個最大）？當前代碼可視化所有羣集。

library(cluster) 
library(hopach) 

distance =distancematrix(DNA[1:30],"cosangle") 
hobpach.DNA =hopach(DNA[1:30],dmat=distance) 

labels = c(hobpach.DNA$clustering$labels) 

table(labels, DNA$class) 

#Plots all clusters 

clusplot(DNA[1:30], hobpach.DNA$clustering$labels, main='Cluster Vis', 
    color=TRUE, shade=TRUE, 
    labels=2, lines=0) 

Answer 1

直接在clusplot沒有這樣的東西。但是，您可以輕鬆找出哪些集羣最大，然後僅將這些集羣提供給clusplot函數。

# find the indicees of the largest clusters 
biggest_indicees <- labels %in% names(sort(table(labels), decreasing = TRUE)[1:3]) 

# plot three largest clusters 
clusplot(df[biggest_indicees, ], labels[biggest_indicees], main = 'Cluster Vis', 
     color = TRUE, shade = TRUE, labels = 2, lines = 0)