dendrogram​関数の結果から、閾値​以下のグループ数をカ​ウントする方法につい​て

以下のコードを実行すると閾値0.4以下のノードが4色に塗り分けられますが、塗り分けられた色の数（グループの数）をカウントする方法はないでしょうか？

rng('default') 
X = rand(20,2);
tree = linkage(X,'average');
dendrogram(tree,0,'ColorThreshold',0.4);

Answer 1

[michio]

かなり不格好ですが・・ラインハンドル H から全ての色を確認し、何種類の色が存在するかをチェックすることでグループ数を確認できます。

rng('default') 
X = rand(20,2);
tree = linkage(X,'average');
dendrogram(tree,0,'ColorThreshold',0.4);
[H,T,outperm] = dendrogram(tree,0,'ColorThreshold',0.4);

lineColors = vertcat(H.Color); % Lineオブジェクトの色を行列に結合
uniqueColors = unique(lineColors, 'rows'); % ユニークな色（黒も合わせて出てきます）
size(uniqueColors,1) - 1 % 黒以外の色の数

[denpika]

いただいたコードで解決しました。続けての質問で恐縮なのですが、分かれたグループ毎に番号（データ点）は取得可能でしょうか？

[michio]

上で得られたグループ数を使って

[H,T,outperm] = dendrogram(tree,4);

と実行すれば、少なくとも各グループの番号は取得できますがいかがでしょうか。

find(T==1)

[denpika]

解決しました、ありがとうございました。