Прогнозирование заболеваний, вызываемых экологическими факторами риска, характерными для Курской области

ПРОГНОЗИРОВАНИЕ ЗАБОЛЕВАНИЙ, ВЫЗЫВАЕМЫХ ЭКОЛОГИЧЕСКИМИ ФАКТОРАМИ РИСКА, ХАРАКТЕРНЫМИ ДЛЯ КУРСКОЙ ОБЛАСТИ

Н.А. Кореневский, Р.А. Крупчатников, Г.В. Чурсин, Н.А. Коптева, В.И. Афанасьев

Аннотация. В статье рассматриваются вопросы получения нечетких решающих правил для решения задач прогнозирования и ранней диагностики заболеваний, которые вызываются вредными факторами различной природы, порождаемыми окружающей средой, включая техногенные факторы. Показывается, что использование нечеткой логики принятия решений в сочетании с информацией об энергегической реакции биологически активных точек позволяет получать уверенность в принимаемых решениях на уровне 0,92, что вполне приемлемо для практического использования.

Ключевые слова: прогнозирование, диагностика, экология, нечеткие решающие правила.

В качестве исходных данных для построения решающих правил прогнозирования и диагностики заболеваний, вызываемых экологическими факторами, могут служить региональные документы о санитарно-эпидемиологической обстановке и о состоянии окружающей среды, а также, по возможности, карты загрязнений. 1еофизичсских полей и заболеваемости, полученные с помощью технологий ГИС.

В частности для Курской области характерными являются 1акие внешние факторы риска, как:

- напряженность постоянного магнитного поля;

- напряженность и частотные диапазоны переменных электромагнитных полей искусственного происхождения;

- пятна радиационного заражения;

- концентрация вредных веществ в атмосфере (пыль, сернистый газ, окись углерода, окислы азота, фенол и его производные, формальдегид);

- концентрация вредных веществ в поверхностных водах (соединения меди, органические вещества, азот нитратный, азот аммонийный, нефтепродукты, железо, фосфаты);

- концентрация вредных веществ в питьевой воде (физико-химические и бактериологические показатели);

- характеристика состояния почвы (пестициды, ядохимикаты, соли тяжёлых металлов).

С приведёнными факторами риска для Курского региона связывают такие основные нозологические формы заболеваний, как :

- болезни органов пищеварения;

- болезни органов дыхания;

- заболевания сердечно-сосудистой системы;

- болезни мочеполовой системы;

- болезни костно-мышечной системы и соединительных тканей;

- болезни эндокринной системы;

- болезни нервной системы;

- психические расстройства;

- инфекционные заболевания;

- новообразования.

В общей постановке задача синтеза правил прогнозирования и ранней диагностики заболеваний, вызываемых экологическими факторами, сводится к поиску взаимосвязей между перечнем нозологических форм и (или) конкретными заболеваниями, например из классификатора ВОЗ, и между всеми существенными факторами риска (информативными признаками) с учётом временных параметров (динамики поведения объекта исследования).

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

Информация об авторах

Кореневский Николай Алексеевич, доктор технических наук, профессор Курского государственного технического университета.

Крупчатников Роман Анатольевич, кандидат сельскохозяйственных наук, доцент кафедры электротехники и механизации животноводства Курской государственной сельскохозяйственной академии имени профессора И.И. Иванова. (4712) 53-13-70, E-mail: roman0406@yandex.ru.

Чурсин Геннадий Викторович, директор филиала ОАО «ВБД» «ЭКДП».

Коптева Наталья Алексеевна, кандидат технических наук, начальник планово-экономического отдела Курской государственной сельскохозяйственной академии имени профессора И.И. Иванова.

Афанасьев Владимир Иванович, председатель СХП.

дополнительные факторы риска, характерные для конкретного индивидуума (наследственность, перенесённые заболевания, употребление алкоголя, табакокурение, иммунный статус и Т.Д.).

Как правило, для перечисленных типов заболеваний такие списки факторов риска известны, по крайней мере, на понятийном уровне с проверкой статистических взаимосвязей.

Проведённый разведочный анализ показал, что каждый из измеряемых факторов риска (признаков) по отношению к задачам прогноза и диагностики носит неполный и нечёткий характер, а структура классов, относительно которых принимается решение, имеет нечёткие границы с зонами пересечения, переходящими из класса в класс. В таких условиях для синтеза соответствующих решающих правил целесообразно использовать теорию нечёткой логики принятия решений, в рамках которой информативные признаки X, и (или) комплексные показатели, получаемые на их основе представляются функциям принадлежности // (2к) к

рассматриваемым классам СО,, а синтез промежуточных и финальных решающих правил осуществляется через формулы расчёта соответствующих коэффициентов уверенности

КУщ [3,6,71.

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

Задачу прогнозирования заболеваний, вызываемых экологическими факторами, будем рассматривать в двух её вариантах.

В первом варианте, описанном в работе [3], будем определять, перейдёт ли относительно здоровый человек (класс сУ0) в класс заболевания с именем I (класс СО,) через определённый

промежуток времени Т0 с уверенностью, определяемой агрегирующим нечётким решающим правилом.

Во втором варианте уверенность в прогнозе определяется по шкале времён, на которой экспертами определяется верхняя граница времени прогнозирования Тв, когда ещё можно делать

приемлемые для практики классификационные выводы. Выбранная временная шкала разбивается экспертами на интервалы времени наблюдения Тф которые могут быть «привязаны» к приемлемой частоте повторяемости обследований. По интервалам наблюдений синтезируются нечёткие решающие правила прогнозирования, которые в общем случае могут быть различны вплоть до использования различных информативных признаков (факторов риска).

Характерной особенностью задач прогнозирования и диагностики заболеваний, вызываемых экологическими факторами, является то, что собираемая экологическими службами и службами здравоохранения информация в основном характеризует общие тенденции заболеваемости без учёта индивидуальных особенностей каждого из обследуемых, находящихся в зоне исследования. Поэтому получаемые на основе этих данных выражения для определения КУ по известным статистическим от-

а>,

чётным данным определяют некоторую усредненную уверенность в прогнозе (донозологиче-ском диагнозе) со, по всей группе людей, находящейся под воздействием наблюдаемых экологических факторов.

Относительно конкретного человека риск появления заболевания со, зависит не только от

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

С учётом сказанного увеличение точности в определении прогноза (диагноза) по классу со, может быть достигнуто путём агрегации дополняющих друг друга составляющих Л Уэа)( (уверенность в со, от действия экологических факторов) и КУН (уверенность от факторов, характеризующих индивидуальный риск человека по диагнозу со,) в финальное решающее правило вида:

КУп =КУ1М +КУИО)

О), э,0)/ и,О), ' Эу0), ' V '

Другой особенностью влияния экологических факторов на организм человека является то, что большинство из них начинает своё

вредное воздействие при условии, что человек находится в зоне воздействия достаточно долгое время.

Интенсивность действия экологического фактора YJ совместно с временем воздействия в частном решающем правиле по фактору с номером ] могут быть учтены при использовании правила вида

О, если / <ТП тгу =} п

где функция принадлежности к

классу со, по шкале интенсивности действия У^

1^0), (0 - функция принадлежности к классу

со, от действия фактора у, с носителем по шкале времени воздействия; Тп - некоторая пороговая величина времени, меньше которой влиянием У; любой практически существующей интенсивности на возникновения заболевания СО,

можно пренебречь.

Если выбрать форму и параметры функции принадлежности, возрастающие с ростом величин их носителей, то при I > Тп КУа1 ^ растёт

по мере роста интенсивности и времени воздействия по закону, определяемому формой и параметрами соответствующих функций принадлежностей.

Известно, что вредному воздействию экологических факторов препятствуют защитные механизмы человеческого организма, снижая риск возникновения и развития соответствующих заболеваний. Учесть влияние защитных механизмов в классификационных решающих правилах можно, синтезировав решающие правила определения уверенности в уровне защитных свойств по диагнозу со, - КУ3 .

Как показали результаты исследований, хорошей информативностью с точки зрения защитных свойств организма, обладают адаптационный потенциал (АП), определяемый через индекс функциональных изменений (ИФИ) [1] и энергетическая сбалансированность (ЭС) меридианных структур организма, которая может быть определена по электрическим характеристикам БАТ «связанных» с общесистемной реакцией организма (Е23, Е36, ЯР6, У40, У60 и УВ20).

Учитывая, что уровни АП и ЭС отражают различные механизмы деятельности человека и каждый из них вносит свой вклад в защитные функции, удобно в качестве меры доверия к

защите от фактора YJ для заболевания со, выбрать параметр

КУ,,^ =м*.ДЛЛ)+д^(ЭОР-/^я(ЛЯ)], (3) где ци (АП) - функция принадлежностей к уровню защитных свойств по фактору ] для класса СО, по значению адаптационного и потенциала; ; (ЭС ) - функция принадлежности к уровню защитных свойств по фактору ] для класса со, по значению энергетической сбалансированности.

С учётом того, что КУ 3 а)/ / уменьшает

общую уверенность в прогнозе (диагнозе) й),,

можно записать

[ 0, если КУа ,<КУ.а , ,Ач

ггуО _ I ’ <“/.7 >Л./

*"* \КУ* /„ если КУа , > КУ,„ ,

^ а,.] °><.) ».<*>< .7

С учётом множества экологических факторов, каждый из которых в той или иной степени приводит к появлению и развитию заболеваний со,, уверенность КУ э 0)/ ) в соответствии с

общими рекомендациями по синтезу нечётких решающих правил разработанными на кафедре биомедицинской инженерии КурскГТУ [3] может быть определена по формуле:

щщи+ъ=кут

Рассматривая информативные признаки х; , характеризующие индивидуальные риски человека по заболеваниям со, и соответствующие комплексные показатели Z^^ как носители функций принадлежностей к классам со, можно синтезировать правила определения КУ и а)/ аналогично (5), заменяя КУ °а>,(у+|>

на ц (у +1) и р (к +1). Тогда финальное

решающее правило принимает вид (1).

В описанном механизме синтеза нечётких решающих правил качество классификации определяется такими субъективными факторами, как компетенция экспертов, их понимаем существа используемых формул, эффективностью взаимодействия экспертов с инженером по знаниям и т.д.

Снизить величину субъективизма можно, если на этапе синтеза или практического использования нечётких прогностических правил удаётся получить контрольные выборки, которые формируются следующим образом. Набираются группы обследуемых, среди которых есть люди с различными факторами риска, ха-

растеризующими возможность появления патологии со(.

Объём контрольной выборки определяется известными в статистике требованиями.

Для каждого из обследуемых для времени I = 0 определяется и запоминается его коэффициент уверенности. Далее в таблице экспериментальных данных (ТЭД) их удобно пересортировать, расставив в порядке увеличения уверенности в прогнозе 0)(. Учитывая непрерывность шкалы КУ Ш/ для анализа качества работы синтезированных решающих правил ТЭД, ее удобно разбить на слои с выбранным интервалом изменения КУ т - АКУ т .

О) I СО (

Пример такой таблицы представлен в таблице 1.

В приведенном примере первая строка соответствует номеру обследуемого. Вторая строка определяет величины интервалов выбранных экспертами для анализа качества прогнозирования. В третьей строке записаны средние значения КУ 0)( в выбранном интервале,

а в четвертой строке при наступлении времени То записывается оценка вероятности попадания обследуемых в класс со( определяемая как количество людей попадающих в класс со( из

конкретного интервала АКУ к общему количеству людей из этого интервала.

Таблица 1 - Распределение обследуемых по уверенности в со1

№ обследуемого (0 1,2... ... ... ...

ЛКУ 0-0,2 0,2-0,4 0,4- 0,6 0,6- 0,8 >0,8

КУср 0,1 0,3 0,5 0,7 0,9

^

После заполнения таблицы типа I легко определить качество построенных прогностических моделей.

Можно считать, что если по каждому из интервалов КУср и Р (Т0) близки друг к другу в смысле заданной пороговой разности между ними (|КУср- Р'^ (Т0)| < еп), то получены удачные прогностические модели.

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

| КУ ср - Р‘а>'

После чего набирается новая контрольная группа, для которой проверяется совпадение

КУср С />;<То).

Анализ таблицы типа 1 позволяет экспертам легко выбрать пороговую величину КУ а (Т0 ) , превышение которой говорит о

высоком риске

появления у обследуемых класса со, с тем,

чтобы своевременно организовать серию лечебно-профилактических мероприятий.

Если наблюдение за обследуемыми проводить в нескольких временных интервалах Т0ч, то, включив в таблицу типа 1 несколько строк с расчётом КУ т (Т0 ч ) по критерию

|*ГУ„ -Р*.,

можно определить, на какой интервал времени можно распространить надёжное прогнозирование по полученному решающему правилу.

Такая модифицированная таблица может быть использована и для контроля качества временных прогностических правил (таблица 2), если наблюдаемый временной интервал разбить на разумно выбранные временные интервалы.

Таблица 2 - Распределение обследуемых по уверенности в со( в зависимости от временного фактора_______________________________________

№ обследуемо-ю 1 1,2,3... ... ...

АКУ АКУ, ДКУ2 АКУ}

КУср КУ,ср ку2ср КУ]ср

^

^,

...

^;,

М„Ж

Рисунок 1 - Функция принадлежности к классу со* с носителем по шкапе Н

w

Рисунок 2 - Функция принадлежности к классу сож по признаку Н в зависимости от временного фактора

Анализ этой таблицы по величине КУср и Р (T0q) позволяет оценить качество прогностических правил по выполнению критериев типа (7) и при необходимости вносить коррективы во временные прогностические правила аналогично рассмотренному выше варианту.

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

Если задача прогнозирования решается как задача разделения на классы не заболеет со0 и заболеет каким-либо из заболеваний со(> обучающие выборки строятся по набору информативных признаков Х=(Х|, х2 , х„... х„), то в ходе разведочного анализа можно определить взаиморасположение и разделимость исследуемых классов (линейная, кусочно-линейная и нелинейная разделимость; наличие и характер

зон пересечения классов, характер признаковых и дистальных гистограмм и т.д.).

Исходя из анализа структуры признакового пространства согласно рекомендациям [3] с использованием объектов обучающей выборки, могут быть выбраны носители У = Р(х) для функций принадлежностей /л, /ло> таким образом, чтобы обеспечить их минимальное пересечение при возможно меньшем числе информативных признаков. В практических приложениях пользователей больше интересует значение функции /Л(0 (У). Тогда задача поиска

формы и параметров функции принадлежностей может заключаться в поиске такого носителя, при котором цШ1 (У) достигает своего

максимального значения при меньшем числе информативных признаков с учетом дополнительных требований на временные и техникоэкономические ограничения на получение значений информативных признаков.

Предложенный метод синтеза апробирован на решении задачи оценки риска заболеваний желудочно-кишечного тракта (ЖКТ) и нервной системы при воздействии постоянного магнитного поля Курской магнитной аномалии и горно-обогатительного комбината на жителей Же-

лезногорского района Курской области и промышленных предприятий на жителей г.Курска. В ходе проведённых исследований были получены системы функций принадлежностей и синтезированы соответствующие нечеткие прогностические и диагностические правила. В качестве примера на рисунке 1 приведена функция принадлежностей к классу риск возникновения ЖКТ (о)ж) с носителем по шкале напряжённости магнитного поля в миллиэрстедах (Н) - *W(H).

Для учёта временного фактора (в годах) по параметру Н получена функция принадлежности, приведённая на рисунке 2.

Уровень адаптационного потенциала определяется по параметру ИФИ, а энергетическая сбалансированность согласно рекомендациям [4] по электрическому сопротивлению БАТ Е23, Е36, RP40, V40, V60 и VB20. В качестве факторов индивидуального риска по заболеваниям ЖКТ были выбраны употребление алкоголя, табакокурение, приём лекарственных средств, раздражающих слизистую желудка, длительное эмоциональное напряжение и реакция БАТ, связанная с заболеваниями ЖКТ (Е21, Е36, V21, V43 и VB24). Общая уверенность в классе риск заболевания ЖКТ рассчитывается в соответствии с выражением (1), достигает величины 0,92, а по заболеваниям нервной системы 0,89, что вполне приемлемо для практического использования.

Список использованных источников

1 Баевский, P.M. Оценка адаптационных возможностей организма и риск развития забо-

леваний / P.M. Баевский, А.Г1. Берсенева.- М.: Медицина, 1997.- 235 с.

2 Гаваа, Лувсан Очерки методов восточной терапии / Гаваа Лувсан.- 3-е изд., перераб. и доп. - Новосибирск: Наука. Сиб. От-ние, 1991.-432с.

3 Кореневский, Н.А. Проектирование нечетких решающих сетей настраиваемых по структуре данных для задач медицинской диагностики / Н.А. Кореневский // Системный анализ и управление в биомедицинских системах.

- 2005. - Т.4.- № 1. - С. 12-20.

4 Кореневский, Н.А. Синтез моделей взаимодействия внутренних органов с проекционными зонами и их использование в рефлексодиагностике и рефлексотерапии: монография / Н.А. Кореневский, В.В. Буняев, В.Н. Гадалов, Н.Д. Тутов.- Курск, 2005.- 224с.

5 Кореневский, Н.А. Энергоинформационные основы рефлексологии / Н.А. Кореневский, М.И. Рудник, Е.М. Рудник.- Курск, 2001. 236с.

6 Круглов, В.В. Искусственные нейронные сети / В.В. Круглов, В.В. Борисов Теория и практика.- 2-е изд., стереотип.- М.: Горячая линия -телеком, 2002.-382с.

7 Осовский, С. Нейронные сети для обработки информации / Пер. с польского И.Д. Ру-динского. - М.: Финансы и статистика, 2002. -344с.