На правах рекламы:
ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Авторитетность издания

ВАК - К1
RSCI, ядро РИНЦ

Добавить в закладки

Следующий номер на сайте

2
Ожидается:
16 Июня 2024

Прототип диагностической системы поддержки принятия решений на основе интеграции байесовских сетей доверия и метода Демпстера–Шефера

Prototype of the diagnostic decision support system on the basis of integration bayesian belief networks and the Dempster–Shafer method
Статья опубликована в выпуске журнала № 1 за 2013 год. [ на стр. 11-16 ]
Аннотация:Рассматриваются применение и программная реализация формального аппарата, интегрирующего возможности расширенных байесовских сетей доверия (Bayesian belief networks) и метода Демпстера–Шефера (Dempster–Shafer). Данный аппарат ориентирован на использование в интеллектуальных системах поддержки принятия решений для помощи специалистам при диагностике сложных проблемных ситуаций. В отличие от традиционных байесовских сетей доверия предлагается применение перспективных многоуровневых (сложно-структурированных) сетей, более удобных для комплексного исследования проблемы и представления экспертных данных. Интеграция позволяет использовать при диагностике как экспертные данные, так и числовые (вероятности) данные, полученные в результате измерений. Совместное использование байесовских сетей доверия и метода Демпстера–Шефера позволяет оценить соотношения между возможными исходами (ситуациями) и уточнить их вероятности. Предложенный подход реализован в прототипе интеллектуальной диагностической системы поддержки принятия решений для диагностики сложных заболеваний зрения.
Abstract:Possibilities of the application and program implementation of formal means integrating capabilities of expanded Bayesian belief networks and Dempster–Shafer’s method are considered. This means are oriented towards using in intelligent decision support systems for the help to experts at diagnostics of difficult problem situations. Unlike traditional Bayesian belief networks, the application of perspective multilevel (difficult-structured) networks, more convenient for complex problem research and representing expert data, is offered. Integration allows to use at diagnostics both expert data, and numerical (probabilistic) data obtained in the result of measurements. Joining Bayesian belief networks and Dempstera–Shafer's method allows to estimate ratios between possible outcomes (situations) and to specify their probabilities. The offered approach is implemented in the prototype of the intelligent decision support system for diagnostics of difficult diseases of vision.
Авторы: Еремеев А.П. (eremeev@appmat.ru) - Национальный исследовательский университет «Московский энергетический институт» (профессор), г. Москва, Россия, доктор технических наук, Хазиев Р.Р. (ruslan.haziev@gmail.com) - Национальный исследовательский университет «Московский энергетический институт» (аспирант ), Москва, Россия, Зуева М.В. (visionlab@yandex.ru) - МНИИ глазных болезней им. Гельмгольца (проф., зав. лабораторией), Москва, Россия, доктор биологических наук, Цапенко И.В. (sunvision@mail.ru) - МНИИ глазных болезней им. Гельмгольца (с.н.с.), Москва, Россия, кандидат биологических наук
Ключевые слова: проблемная ситуация., диагностика, поддержка принятия решений, интеллектуальная система
Keywords: problem situation, diagnostics, decision support, intellectual system
Количество просмотров: 13037
Версия для печати
Выпуск в формате PDF (5.29Мб)
Скачать обложку в формате PDF (1.21Мб)

Размер шрифта:       Шрифт:

В статье рассматривается совместное использование байесовских сетей доверия (БСД) [1] и метода Демпстера-Шефера (МДШ) [2], предназначенное для помощи ЛПР при диагностике сложных проблемных ситуаций. В отличие от классической БСД в системе используется многоуровневая архитектура, объединяющая несколько методов и позволяющая исследовать проблему в целом. Разбиение на уровни позволяет также сократить область перебора различных ситуаций для поиска наиболее вероятного результата. В системе используются не только рассуждения экспертов, но и анализ числовых данных, полученных при проведении различных диагностик с применением вероятностных методов.

Установлено [3], что при применении аппарата БСД в интеллектуальных системах поддержки принятия решений (ИСППР), особенно в ИСППР реального времени (РВ), недостаточно определить основные составляющие события в сети и их связи. Нужно также учитывать все возможные ситуации, которые может предвидеть эксперт. Для этого предлагается использовать сложноструктурированную БСД. Такой тип БСД был адаптирован и применен в прототипе ИСППР для диагностирования сложных глазных заболеваний на примере патологий сетчатки, когда невозможно с абсолютной уверенностью сказать, какая именно болезнь у пациента и на какой стадии. Даже после выявления симптомов с помощью специальных компьютерных средств диагностики и заключений опытных специалистов пациентам приходилось повторно проходить обследование из-за неточности и недостаточности имеющейся информации. При построении модели и опытной апробации было выявлено, что использование БСД позволяет анализировать только один из возможных исходов (наличие одного из заболеваний), игнорируя другие варианты. Совместное использование БСД и МДШ позволит оценить соотношения между возможными исходами (ситуациями) и уточнить их вероятности.

Диагностирование патологии зрения типа глаукомы

В качестве объекта исследований была выбрана глаукома. Основными признаками данного заболевания являются повышение внутриглазного давления и снижение зрения. В основе заболевания лежит нарушение циркуляции внутриглазной жидкости (чаще всего затруднение ее оттока) [4]. Статистические данные и данные по пациентам были получены в Московском НИИ глазных болезней им. Гельмгольца. Создаваемый прототип ИСППР должен помочь ЛПР (врачу-офтальмоло­гу) в выявлении на начальной стадии признаков глаукомы и в уточнении типа заболевания с целью как можно более раннего назначения соответствующего лечения.

Рассмотрим основные этапы исследований при диагностике глаукомы. Прежде всего пациент дол­жен пройти комплекс клинических обследований у врача-офтальмолога. На этом этапе применяется определенный набор психофизических тестов для выявления патологий. Если имеются основательные (достоверные) данные в пользу определенного заболевания пациента, то дальнейшие исследования нецелесообразны; в противном случае обследование может продолжиться. Дообследование часто необходимо на ранних стадиях заболевания или в сложнодифференцируемых случаях, например при глаукоме низкого давления, когда данные клинического обследования пациента могут почти не отличаться от нормы.

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

На заключительной стадии специалист-физио­лог (или группа специалистов) на основе полученных в результате обследований пациента данных должен подтвердить или опровергнуть предварительный диагноз офтальмолога (в последнем случае надо уточнить диагноз) и рекомендовать соответствующее лечение.

На рисунке 1 представлена БСД, характеризующая процесс диагностирования глаукомы в виде множества определяемых событий: T – результаты клинических исследований у офтальмолога (подозрение на наличие у пациента заболевания); Z – заключение физиолога о наличии болезни в ходе проведения исследований ЭРГ и данных офтальмолога; I – результат исследования ЭРГ; D – результирующее событие, характеризующее решение поставленной задачи (уточненный диагноз, текущая стадия, рекомендации по лечению или дальнейшему обследованию и т.д.). Цифрами обозначаются события, относящиеся к этапам исследования пациента.

Рассмотрим событие, определяющее результат исследования ЭРГ, на примере первичной открытоугольной глаукомы (ПОУГ). Для диагностирования данного заболевания (состояния) и его стадии необходимо проанализировать ряд ЭРГ, сделать по ним заключение о выявленных признаках (симптоматике) заболевания и определить степень (вероятность) уверенности в диагнозе.

На рисунке 2 представлена модель второго уровня для определения значения события «Результат исследования ЭРГ», где I – целевое событие (диагноз); IК, …, IРТ – результаты проведенных исследований: IK – результат, полученный при исследовании колбочковой ЭРГ, IM – максимальной ЭРГ, IP – палочковой ЭРГ, IOP – осцилляторных потенциалов (ОП), IF – диаграммы фотопического негативного ответа (ФНО), IPT – паттерн–ЭРГ (ПЭРГ).

Подпись:  
Рис. 3. Пример общей ЭРГ
В модели второго уровня представляются результаты основных видов ЭРГ, каждая из которых имеет свой собственный вид и изменяется в зависимости от заболевания, то есть определяются ключевые значения, характеризующие ту или иную болезнь. Палочковая, колбочковая и максимальная ЭРГ используются для проверки способности глаза к темновой и световой адаптации. Каждый из компонентов ЭРГ генерируется различными структурами сетчатки. На рисунке 3 представлен пример ЭРГ и ее компонентов [5].

Каждый из компонентов ЭРГ генерируется различными структурами сетчатки. Основными результирующими показателями являются a- и b-волны [5]. В ходе диагностики анализируются значения амплитуды и времени кульминации волн различных видов ЭРГ. В таблице 1 представлены полученные статистические данные амплитуд колебаний по обследуемым пациентам с ПОУГ 1-й и 2-й стадий.

Таблица 1

Амплитуда a- и b-волн на 1-й и 2-й стадиях ПОУГ (M±m)

Виды стандартных ЭРГ

Норма, мкВ

M±2SD

1-я стадия

2-я стадия

мкВ (% от нормы)

мкВ (% от нормы)

Палочковая

b-волна

58,1±21,7

44,3±5,9 (76,2)

35,1±27,2 (60,4)

Максимальная

a-волна

103,7±37,3

87,8±7,9 (84,6)

70,3±30,4 (67,8)

b-волна

203,2±55,0

176,5±13,3 (86,8)

196,6±46,1 (96,7)

Колбочковая

a-волна

18,2±7,5

16,7±2,7 (91,7)

15,1±6,1 (82,5)

b-волна

91,2±21,0

49,0±6,0 (53,7)

49,6±21,8 (54,4)

По данным амплитуд можно проанализировать изменение значения проходящего через сетчатку сигнала. При наличии заболевания происходит угнетение a- и b-волн. Кроме того, может изменяться и время их кульминации: в худшем случае время увеличивается, а при начальной стадии заболевания время не изменяется относительно нормы. Прогрессирование болезни способствует значительному снижению амплитуды и увеличению времени кульминации волн ЭРГ, снижению зрительных функций.

Рассмотрим модель (БСД) для диагностики ПОУГ на основе проведенного анализа максимальной ЭРГ (рис. 4), где Aa, Ab – значения максимальной амплитуды; Ta, Tb – значения времени кульминации соответствующих волн; Va, Vb – заключения на основе введенных данных, определяющие заболевание; IM – целевое событие, соответствующее общему результату анализа ЭРГ (диагнозу).

Предложенная модель учитывает все ключевые значения стандартных видов ЭРГ для ПОУГ. На данном этапе диагностики ПОУГ возможны только 1-я и 2-я стадии, поэтому можно отметить ряд особенностей:

–      Aa и Ab могут принимать различные значения в зависимости от типа и стадии заболевания в некотором диапазоне относительно средних значений амплитуд, указанных в таблице 1 (например, для заболевания ПОУГ 1-й стадии каждое событие имеет три состояния: «более типично для заболевания ПОУГ 1-й стадии», «типично для заболевания ПОУГ 1-й стадии», «менее типично для заболевания ПОУГ 1-й стадии»; зоны состояний могут пересекаться, то есть нужно проанализировать все состояния, если значения максимальных амплитуд a- и b-волн находятся в зоне пересечений);

–      Ta и Tb изменяются только при ПОУГ 2-й стадии, следовательно, можно выделить два состояния: «больше нормы для здорового человека», «норма для здорового человека»;

–      Va и Vb в зависимости от ключевых значений позволяют констатировать наиболее соответствующее состояние (диагноз): «ПОУГ 1-й стадии», «ПОУГ 2-й стадии» или «здоров»;

–      IM обозначает событие, результат которого при анализе всех заключений по данным максимальной ЭРГ имеет состояния, аналогичные Va и Vb.

Рассмотрим пациента с подозрением на ПОУГ 1-й стадии. В таблицах 2 и 3 заданы ключевые значения амплитуд и времени кульминации a- и b-волн, соответствующие состояния событий и коэффициенты уверенности при анализе максимальной ЭРГ глаза.

Значение амплитуды a-волны соответствует двум состояниям событий – «типично для заболевания ПОУГ 1-й стадии» и «менее типично для заболевания ПОУГ 1-й стадии». Нужно рассчитать вероятность выполнения обоих событий. Наиболее вероятным значением состояний событий Va и Vb является «ПОУГ 1-й стадии». Зададим вероятности этих событий: P(VaçAa, Ta)=0,68 при Aa = «типично для заболевания ПОУГ 1-й стадии»; P(VaçAa, Ta)=0,62 при Aa = «менее типично для заболевания ПОУГ 1-й стадии»; P(VbçAb, Tb)=0,68.

Таблица 2

Амплитуда a- и b-волн при анализе максимальной ЭРГ

Амплитуда

Значение амплитуды

Значение коэффициента уверенности, P(x), x=Aa, Ab

Aa – амплитуда a-волны

85,1 мкВ (82 % от нормы)

0,8 – типично для заболевания ПОУГ 1-й стадии; 0,75 – менее типично для заболевания ПОУГ 1-й стадии

Ab – амплитуда b-волны

177,1 мкВ (87,1 % от нормы)

0,8 –типично для заболевания ПОУГ 1-й стадии

Таблица 3

Время кульминации a- и b-волн при анализе максимальной ЭРГ

Время кульминации

Значение амплитуды

Значение коэффициента уверенности, P(x), x=Ta, Tb

Ta – время кульминации a-волны

24 мкВ (135 % от нормы)

0,65 (Aa = «типично для заболевания ПОУГ 1-й стадии»); 0,8 (Aa = «менее типично для заболевания ПОУГ 1-й стадии»)

Tb – время кульминации b-волны

46 мкВ (124,1 % от нормы)

0,65 – больше нормы для здорового человека

Конечным событием для анализа данных максимальной ЭРГ является событие IM – следствие наиболее вероятного заболевания пациента или опровержение подозрения при клиническом исследовании. В данном случае необходимо подтвердить или опровергнуть, что у пациента ПОУГ 1-й стадии. Значение коэффициента уверенности при такой ситуации будет P(IMçVa, Vb)=0,95. Рассчитаем совместную вероятность событий при условии, что Aa = «типично для заболевания ПОУГ 1-й стадии»: P(Aa, Ta, Ab, Tb, Va, Vb, IM)= =P(Aa)×P(Ta)×P(Ab)×P(Tb)×P(VaçAa, Ta)×P(VbçAb, Tb)´ ´P(IMçVa, Vb)≈0,119.

Чтобы полностью удостовериться, что результатом анализа максимальной ЭРГ является наличие заболевания, нужно определить наибольшее значение совместной вероятности событий. Рассчитав совместную вероятность событий при условии, что Aa = «менее типично для заболевания ПОУГ 1-й стадии», получим P(Aa, Ta, Ab, Tb, Va, Vb, IM)≈0,125.

Таким образом, результат (диагноз) Aa = «менее типично для заболевания ПОУГ 1-й стадии» имеет наибольшее значение совместной вероятности событий. Чтобы удостовериться, что результатом анализа максимальной ЭРГ является наличие у пациента ПОУГ 1-й стадии, нужно рассчитать соответствующую условную вероятность при заданных событиях, определяющих значения ключевых точек на ЭРГ.

Расчет производится по формуле

P(IM | Aa, Ta, Ab, Tb)=P(IM, Aa, Ta, Ab, Tb) / P(Aa, Ta, Ab, Tb) = (P(IM, Aa, Ta, Ab, Tb, Va, Vb) + P(IM, Aa, Ta, Ab, Tb, Va, ØVb) + P(IM, Aa, Ta, Ab, Tb, ØVa, Vb) + P(IM, Aa, Ta, Ab, Tb, ØVa, ØVb)) / ((P(Aa, Ta, Ab, Tb, Va, Vb, IM) + P(Aa, Ta, Ab, Tb, Va, ØVb, IM) + P(Aa, Ta, Ab, Tb, ØVa, Vb, IM) + P(Aa, Ta, Ab, Tb, ØVa, ØVb, IM) + P(Aa, Ta, Ab, Tb, Va, Vb, ØIM) + P(Aa, Ta, Ab, Tb, Va, ØVb, ØIM) + P(Aa, Ta, Ab, Tb, ØVa, Vb, ØIM) + P(Aa, Ta, Ab, Tb, ØVa, ØVb, ØIM)).

После соответствующих вычислений получим результат 0,67 (67 %), означающий, что при анализе максимальной ЭРГ пациента наиболее вероятным диагнозом является наличие ПОУГ 1-й стадии, но данной информации недостаточно, чтобы подтвердить (или опровергнуть) соответствие данного диагноза действительности, поскольку, как уже отмечалось, использование БСД позволяет анализировать лишь один из возможных диагнозов, в частности наличие или отсутствие ПОУГ 1-й стадии. Для устранения этого недостатка (для проведения комплексного анализа с возможностью учета и других возможных диагнозов) предлагается, помимо БСД, использовать также МДП, который позволяет сопоставить все имеющиеся значения как множества признаков, характеризующих заболевания.

Например, значение Va может иметь три возможных значения (состояния). В зависимости от введенных значений амплитуды и времени кульминации волн можно задать вероятность принадлежности для всех возможных состояний события Va. Каждое из состояний Va будем рассматривать как заключение одного из экспертов, который утверждает, что при введенных данных его заболевание является наиболее вероятным. Исходя из этого можно определить интервалы доверия и правдоподобия для каждого из заболеваний. В таблице 4 приведены соответствующие значения масс, коэффициентов доверия и правдоподобия, используемые для расчетов.

Таблица 4

Значения коэффициентов доверия и правдоподобия для события Va

Событие

Масса

Доверие

Правдоподобие

Значение NULL

0

0

0

Здоров

0,05

0,1

0,15

ПОУГ 1-й стадии

0,55

0,62

0,65

ПОУГ 2-й стадии

0,3

0,35

0,4

Одно из (здоров, ПОУГ 1-й ст., ПОУГ 2-й ст.)

0,1

1,0

1,0

Выделенные значения в таблице можно интерпретировать как нижние и верхние пределы интервалов, на которых находятся значения соответствующих вероятностей. Согласно полученным данным можно затем корректировать значения вероятностей для БСД. Совместное использование аппарата БСД и МДШ позволит повысить качество диагностики при наращивании БД по заболеваниям и методам их анализа.

ПО прототипа диагностической ИСППР

Для реализации прототипа диагностирующей ИСППР был выбран перспективный язык искусственного интеллекта Clips [6]. В качестве входных данных используются информация о клинических исследованиях пациента, данные различных видов ЭРГ и заключения врачей. Выходной информацией являются диагноз с соответствующей вероятностью (уверенностью) и предлагаемые методы лечения или рекомендуемые действия для уточнения информации о пациенте. По данным, полученным при анализе различных видов ЭРГ, строится БД (база фактов), хранящая словесное исходное описание задачи (сведения о пациенте). Для построения БСД и оперирования данными из базы фактов используется база знаний (правил), включающая методы получения результата.

Представим фрагмент программы на языке Clips с примером правила добавления нового элемента – значения максимальной амплитуды a-волны при попадании в интервал для ПОУГ 1-й стадии:

(defglobal ?*Vvod-Aa-m* = 0)

(defrule opredelenie-Aa-max-ERG-2 ""

(declare (salience 98))

(not (var Aa ?))

(Aa POUG 1 st.)

=>

(printout t "Выберите интервал, в который попадает значение максимальной амплитуды a-волны:" crlf)

(printout t "Варианты ответа:" crlf)

(printout t "[1] более типично для заболевания ПОУГ 1 ст." crlf)

(printout t "[2] типично для заболевания ПОУГ 1 ст." crlf)

(printout t "[3] менее типично для заболевания ПОУГ 1 ст." crlf)

(bind ?*Vvod-Aa-m* (ask-question "Ответ: " 1 2 3))

(if (eq ?*Vvod-Aa-m* 1)

 then

  (assert (var Aa more for POUG 1 st.))

  (printout t "Aa = более типично для заболевания ПОУГ 1 ст." crlf)

 else

  (if (eq ?*Vvod-Aa-m* 2)

  then

   (assert (var Aa norm for POUG 1 st.))

   (printout t "Aa = типично для заболевания ПОУГ 1 ст." crlf)

 else

  (if (eq ?*Vvod-Aa-m* 3)

  then

   (assert (var Aa less for POUG 1 st.))

   (printout t "Aa = менее типично для заболевания ПОУГ 1 ст." crlf))))

(assert (label Aa)))

После задания в системе множества фактов для описания событий Aa, Ab, Ta, Tb на основании коэффициентов уверенности выполняется поиск правил для определения состояний событий Va и Vb. В процессе поиска может возникнуть ситуация, когда выполняются несколько правил, формирующих некоторое множество фактов при одних и тех же условиях. Это объясняется тем, что экстремумы a- и b-волн могут соответствовать нескольким состояниям. В этом случае используется наибольшее значение условной вероятности события Im.

Литература

1.     Бидюк П.И., Терентьев А.Н., Гасанов А.С. Построение и методы обучения байесовских сетей // Кибернетика и системный анализ. 2005. № 4. С. 133–147.

2.     Люгер Джордж Ф. Искусственный интеллект: стратегии и методы решения сложных проблем; [пер. с англ.]. М.: Вильямс, 2003. 864 с.

3.     Хазиев Р.Р. Расширенные байесовские сети доверия и их применение в системах поддержки принятия решений // IS&IT’2012: тр. конгр. по интеллект. системам и информ. технологиям. М.: Физматлит, 2012. Т. 1. С. 224–232.

4.     Нероев В.В., Зуева М.В., Цапенко И.В., Рябина М.В., Киселева О.А., Каламкаров Г.Р. Ишемические аспекты патогенеза заболеваний сетчатки // Российский офтальмологический журнал. 2010. Т. 3. № 1. С. 42–49.

5.     Зуева М.В., Цапенко И.В. Электрофизиологическая характеристика глиально-нейрональных взаимоотношений при ретинальной патологии // Сенсорные системы. 1992. № 3. С. 58–63.

6.     Джарратано Д., Райли Г. Экспертные системы: принципы разработки и программирование; [пер. с англ.]. М.: Вильямс, 2007. 1152 с.


Постоянный адрес статьи:
http://swsys.ru/index.php?page=article&id=3372%E2%8C%A9=%E2%8C%A9=en&like=1
Версия для печати
Выпуск в формате PDF (5.29Мб)
Скачать обложку в формате PDF (1.21Мб)
Статья опубликована в выпуске журнала № 1 за 2013 год. [ на стр. 11-16 ]

Возможно, Вас заинтересуют следующие статьи схожих тематик: