Для расчета сегментации на основе этого метода пользователю нужно задать вероятность, с которой максимум статистики на каждом сегменте должен превысить порог. Чем больше мы зададим вероятность, тем больше точек нестационарности будет найдено и тем менее заметные изменения в последовательности будут выявлены.
Как и в случае с методом фрактальной размерности, для просмотра результатов каждому сегменту последовательно будем присваивать стадии от 1 до 4, в результате получим гипнограмму, показывающую границы найденных сегментов. Результат будем сравнивать с гипнограммой, построенной экспертом.
Пример 1. Вероятность, с которой максимум статистики превысит порог, равна 0,02. Количество сегментов в этом случае равно 9190. Результаты сегментации приведены на рисунке 21. Ярко выражены сегменты соответствующие стадии бодрствования (выделены красном на рис.21), а также некоторые из сегментов 4 стадии.
Рисунок 21. Верхний рисунок – сегментация при вероятности 0,02. Нижний рисунок – гипнограмма эксперта
|
|
Приведем процентное соотношение длин сегментов. Посчитаем количество сегментов длиной до 5 секунд, до 10, 30 и 60. Результаты приведены в таблице 3.
Табл.3. Длины сегментов
Длина, сек | ||||
Проценты | 88 | 8.2 | 2.3 | 0.3 |
Эти данные подтверждают результаты исследования ученых в том, что длина большинства стационарных сегментов ЭЭГ не превышает 5 секунд.
Увеличим вероятность обнаружения точек нестационарности до 0,05. Число сегментов возрастает до 12600 (рис.22).
Рисунок 22. Верхний рисунок – сегментация при вероятности 0,05. Нижний рисунок – гипнограмма эксперта
По-прежнему хорошо выделены сегменты бодрствования. Заметить соответствие сегментов сна стадиям гинограммы уже сложнее, так как задавая большую вероятность, мы заставляем алгоритм находить более мелкие изменения в характере ЭЭГ, тем самым создавая сегменты маленькой длины. Количество сегментов длиной до 5 секунд увеличилось до 93% (табл.4).
Табл.4. Длины сегментов
Длина, сек | ||||
Проценты | 93 | 5.5 | 1.1 | 0.2 |
Рисунок 23. Верхний рисунок – сегментация при вероятности 0,1. Нижний рисунок – гипнограмма эксперта
Пример 2. Вероятность превышения порога максимумом статистики равна 0,1. Количество сегментов в данном случае равно 13700. Результаты сегментации приведены на рисунке 24. В данном примере по сравнению с предыдущим нам пришлось увеличить вероятность обнаружения нестационарностей, так как алгоритм некоторые участки ЭЭГ определял в один сегмент, в то время как эксперт разделял его на несколько стадий. Видимо это происходило из-за того, что изменения ЭЭГ на этом участке слабо выделялось на фоне остальных нестационарностей. Процентное соотношение длин сегмент, как и в предыдущем примере, подтверждает, что средняя длина стационарного сегмента не превышает 5 секунд (табл.5).
|
|
Табл.5. Длины сегментов
Длина, сек | ||||
Проценты | 89,8 | 7,2 | 2 | 0,4 |
Увеличим величину вероятности обнаружения до 0,2. Количество сегментов возросло до 15800 (рис.24). Также увеличилось количество сегментов длиной до 5 секунд до 92% (табл.6).
Рисунок 24. Верхний рисунок – сегментация при вероятности 0,2. Нижний рисунок – гипнограмма эксперта
Табл.6. Длины сегментов
Длина, сек | ||||
Проценты | 91,9 | 6 | 1,4 | 0,2 |