6. Методы многопараметрического анализа. Задача кластерного анализа. Алгоритм clope: критерий разделения на кластеры, градиент, функция стоимости.

7. Методы многопараметрического анализа. Задача кластерного анализа. Алгоритм clope: этапы алгоритма кластеризации, задача кластеризации по алгоритму, вычислительная сложность алгоритма.

Пусть имеется база транзакций D, состоящая из множества транзакций {t₁,t₂,…,t_n}. Каждая транзакция есть набор объектов {i₁,…,i_m}. Множество кластеров {C₁,…,C_k} есть разбиение множества {t₁,…,t_n}, такое, что C1 … C_k={t₁,…,t_n} и C_i <> ^С_i C_j = 1<=i, j<=k. Каждый элемент C_i называется кластером, n, m, k – количество транзакций, количество объектов в базе транзакций и число кластеров соответственно.

Каждый кластер C имеет следующие характеристики:

D(C) – множество уникальных объектов; Occ(i,C) – количество вхождений (частота) объекта i в кластер C; W(C) = |D(C)|; H(C) = S(C)/W(C).

Гистограммой кластера C называется графическое изображение его расчетных характеристик: по оси OX откладываются объекты кластера в порядке убывания величины Occ(i,C), а сама величина Occ(i,C) – по оси OY (рис. 2).

Рисунок 2. Иллюстрация гистограммы кластера

На рис. 2 S(C), равное 8, соответствует площади прямоугольника, ограниченного осями координат и пунктирной линией. Очевидно, что чем больше значение H, тем более "похожи" две транзакции. Поэтому алгоритм должен выбирать такие разбиения, которые максимизируют H.

Однако учитывать одно только значение высоты H недостаточно. Возьмем базу, состоящую из 2х транзакций: {abc, def}. Они не содержат общих объектов, но разбиение {{abc, def}} и разбиение {{abc}, {def}} характеризуются одинаковой высотой H=1. Получается, оба варианта разбиения равноценны. Но если для оценки вместо H(C) использовать градиент G(C)=H(C)/W(C)=S(C)/W(C)², то разбиение {{abc},{def}} будет лучше (градиент каждого кластера равен 1/3 против 1/6 у разбиения {{abc, def}}).

Обобщив вышесказанное, запишем формулу для вычисления глобального критерия – функции стоимости Profit(C):

|C_i| – количество объектов в i-том кластере, k – количество кластеров, r – положительное вещественное число большее 1.

С помощью параметра r, названного авторами CLOPE коэффициентом отталкивания (repulsion), регулируется уровень сходства транзакций внутри кластера, и, как следствие, финальное количество кластеров. Этот коэффициент подбирается пользователем. Чем больше r, тем ниже уровень сходства и тем больше кластеров будет сгенерировано.

Формальная постановка задачи кластеризации алгоритмом CLOPE выглядит следующим образом: для заданных D и r найти разбиение C: Profit(C,r) -> max.

8. Методы многопараметрического анализа. Задача кластерного анализа. Меры расстояний между объектами (метрики).

9. Методы многопараметрического анализа. Задача кластерного анализа. Проверка достоверности кластеризации: критерии, дискриминантный анализ.

10. Методы многопараметрического анализа. Задача дискриминантного анализа (да). Ограничения на исходные данные. Применение функции Фишера и расстояния Махаланобиса в да

Дискриминантный анализ позволяет выяснить, действительно ли группы различаются между собой, и каким образом (какие переменные вносят наибольший вклад в имеющиеся различия).

Ограничения: 1- ни одна переменная не может быть линейной комбинацией других переменных. 2-закон распределения для каждой переменной имеет нормальное распределение при фиксированных значениях остальных переменных.

Расстояние Махаланобиса: d(Xi,Xj)=(Xi-Xj)T * S -1 *(Xi-Xj)

где xik , x jk - значения k-й переменной i-го и j-го объекта; Xi , X j - векторы значений переменных у i-го и j-го объектов; S* - общая ковариационная матрица.

линейный дискриминант Фишера

Предположим, что два наблюдаемых класса имеют средние и ковариационные матрицы .Тогда для линейной комбинации признаков средними будут , а ковариационные матрицы будут иметь вид для .Фишер взял за расстояние между этими распределениями величину, равную отношению межклассовой дисперсии к внутриклассовой:

Эта величина в некотором смысле характерезует соотношение сигнал-шум для разметки классов.Можно показать, что наилучшим образом классы разделимы при

.

Если выполняются предположения нормальности и равенства дисперсий, то полученное выше равенство эквивалентно ЛДА.