Проверьте расчет коэффициента несогласованности для иерархического дерева кластера.

Загрузить 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-й кластер соответствует звену индекса $\mathrm{}\mathrm{190-120}\mathrm{=}$70, где 120 - количество наблюдений. 203-й кластер соответствует 83-му звену.

По умолчанию inconsistent использует два уровня дерева для вычисления Y. Поэтому для вычисления коэффициента несогласованности для 84-й линии связи используется только 70-я, 83-я и 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