2016-11-22 5 views
0

私はHOPACHクラスタリングを使用しています - 最大でn個のクラスタ(例:3番目のプロットのみ)を視覚化する方法はありますか?現在のコードはすべてのクラスタを可視化します。クラスタリング - 最大のn個のクラスタをプロットする

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) 

答えて

0

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) 
+0

私は[、biggest_indicees] DFに若干の問題が...エラーを抱えている:タイプの オブジェクトの閉鎖は、 " –

+0

subsettableではありませんあなたは上記のコードで言及してあなたがラベルを定義しましたか?あなたは 'labels = c(hobpach.DNA $ clustering $ labels)'を持っています。 – shadow

+0

はい、まさに - 私はこの問題をオンラインで検索しましたが、運はありません。 –

関連する問題