Исследуйте содействующее вычисление несоответствия на иерархическое кластерное дерево.

Загрузите набор данных examgrades.

load examgrades

Создайте иерархическое кластерное дерево.

Z = linkage(grades);

Создайте матрицу информации о коэффициенте несоответствия с помощью inconsistent. Исследуйте информацию на 84-ю ссылку.

Y = inconsistent(Z);
Y(84,:)

ans = 1×4

    7.2741    0.3624    3.0000    0.5774

Четвертый столбец Y содержит коэффициент несоответствия, который вычисляется с помощью среднего значения в первом столбце Y и стандартного отклонения во втором столбце Y.

Поскольку строки Y соответствуют строкам Z, исследуют 84-ю ссылку в Z.

Z(84,:)

ans = 1×3

  190.0000  203.0000    7.4833

84-я ссылка соединяет 190-е и 203-и кластеры в дереве и имеет высоту 7.4833. 190-й кластер соответствует ссылке индекса $190-120=70$, где 120 количество наблюдений. 203-й кластер соответствует 83-й ссылке.

По умолчанию inconsistent использует два уровня дерева, чтобы вычислить Y. Поэтому это использует только 70-е, 83-и, и 84-е ссылки, чтобы вычислить коэффициент несоответствия для 84-й ссылки. Сравните значения в Y(84,:) с соответствующими вычислениями при помощи высот ссылки в Z.

mean84 = mean([Z(70,3) Z(83,3) Z(84,3)])

mean84 = 7.2741

std84 = std([Z(70,3) Z(83,3) Z(84,3)])

std84 = 0.3624

inconsistent84 = (Z(84,3)-mean84)/std84

inconsistent84 = 0.5774

Создайте выборочные данные.

X = gallery('uniformdata',[10 2],12);
Y = pdist(X);

Сгенерируйте иерархическое кластерное дерево.

Z = linkage(Y,'single');

Сгенерируйте график древовидной схемы иерархического кластерного дерева.

dendrogram(Z)

Вычислите коэффициент несоответствия для каждой ссылки в кластерном дереве Z на глубину 3.

W = inconsistent(Z,3)

W = 9×4

    0.1313         0    1.0000         0
    0.1386         0    1.0000         0
    0.1463    0.0109    2.0000    0.7071
    0.2391         0    1.0000         0
    0.1951    0.0568    4.0000    0.9425
    0.2308    0.0543    4.0000    0.9320
    0.2395    0.0748    4.0000    0.7636
    0.2654    0.0945    4.0000    0.9203
    0.3769    0.0950    3.0000    1.1040