Быстрое развитие медицины, огромный поток поступающей информации об этиологии, патогенезе заболеваний диктуют правила для поиска новых методов обработки полученных результатов. В таких условиях важно уметь находить нужную информацию и принимать правильные решения, от которых будут зависеть течение и исход заболевания. Широкое внедрение компьютерных технологий, создание программных обеспечений позволяют анализировать задачи диагностики и прогнозирования различных заболеваний. Одним из перспективных направлений медицинской нейроинформатики являются нейронные сети (НС). Актуальность применения искусственных нейронных сетей объясняется их востребованностью и практичностью в использовании.

Впервые о нейронных сетях заговорили в 1943 году, когда У. Маккалог и У. Питтс создали первую математическую модель нейронной сети. Практическая модель при помощи компьютера была создана в 1957 году Ф. Розенблаттом. С тех пор нейронная сеть стала активно применяться в различных сферах жизни общества и в науке [1][2][3][4].

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

В зависимости от количества нейронов НС подразделяются на однослойные, настроенные на выполнение простейших команд, и многослойные, обладающие большей вычислительной способностью. На их основе создаются различные автоматизированные системы диагностики, системы распознавания текста, системы анализа и прогнозирования, системы автоматической классификации и сверки информации.

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

На наш взгляд, наглядным примером возможности изучения вопросов прогнозирования течения заболевания является рожистое воспаление. Согласно общепринятой классификации рожистого воспаления (Черкасов В.Л., 1986) выделяют эритематозную, эритемо-буллезную, эритемо-геморрагическую и буллезно-геморрагическую формы воспаления [2][5][6][7][8][9]. Ряд авторов выделяют еще флегмонозную и некротическую формы рожистого воспаления. При этом у части больных отмечается четкая стадийность течения процесса т.е. переход поверхностных форм воспаления в более глубокие со значительным повреждением мягких тканей конечностей и увеличением числа летальных исходов. У некоторых больных инфекционный процесс останавливается на уровне эритематозной или эритематозно-буллезной формы воспаления. С указанных позиций вопросы прогнозирования течения рожистого воспаления являются актуальными [9][10]. Многолетние исследования в области нейросетевого прогнозирования привели к созданию большого количества публикаций [8][9][10][11][12][13][14][15][16]. Однако эти исследования не содержат детального описания математических моделей прогнозирования различных заболеваний. Данные работы содержат лишь единые рекомендации использования нейронных сетей, а также алгоритмы их обучения, наглядные примеры.

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

Y (t–N+1), Y(t–N+2), ..., Y(t),(1)

где Y — вектор контролируемых параметров двигателя; t — дискретное время.

Задача прогнозирования в нейросетевом базисе сводится к построению нейросетевой модели (предиктора), позволяющей найти значения вектора Y в момент времени t+1 по предшествующим N значениям временного ряда, т.е.

Y(T+1) = f [Y(t), Y(t–1), ..., Y(t–N+1),(2)

где f(•) — некоторая нелинейная вектор — функция, которую требуется оценить с помощью НС.

Достоверность прогноза можно оценить по следующему выражению:

||εt+1|| = ||Ŷt+1 – Yt+1||,(3)

где Ŷt+1 — прогнозируемое значение, вычисленное НС для момента времени t+1; Yt+1 — реальное значение вектора Y в этот же момент времени; εt+1 — ошибка прогноза.

Достоинствами НС являются:

- способность фильтрации посторонних шумов. После обучения НС способны воспринимать только нужную информацию;

- во время работы НС непрерывно находятся в состоянии самообучения, что лежит в основе их адаптации;

- быстрое выполнение команд за короткое время;

- возможность анализировать несколько параметров одновременно;

- возможность работы с большим объемом данных;

- удобны в использовании.

В данной работе произведена обработка данных ретроспективного исследования историй болезней пациентов, получавших стационарное лечение на базе ГБУЗ РБ ГКБ № 8 г. Уфы в период 2006–2015 гг., для решения задачи прогнозирования рожистого воспаления с использованием современных статистических пакетов программ и среды моделирования MATLAB. Использовалась рекуррентная (динамическая) НС, реализующая зависимость вида (2).

На рисунке 1 приведена архитектура рекуррентной НС, где Z–1 — элемент временной задержки на один такт Δt; N — число элементов задержки (размер «временного окна»); yi(t+1) — прогнозируемое (на один шаг вперед) значение параметра yi(t).

Анамнез пациента в обучающей выборке включает в себя (табл. 1): Y1 — пол пациента; Y2 — место проживания; Y3 — возраст; Y4 — время от начала з/б (ч); Y5 — тяжесть состояния; Y6 — форма выявления; Y7 — cезон; Y8 — локализация; Y9 — D-димер, мкг/мл; Y10 — лейкоциты крови; Y11 — местная температура (холодная зона); Y12 — лихорадочный период (сутки); Y13 — сопут. ССЗ; Y14 — органы дыхания и др.

Таблица 1. Фрагмент обучающей выборки по результатам анамнеза пациентаTable 1. Fragment of the training sample based on a patient’s medical history

TY1Y2Y3Y4Y5Y6Y7Y8Y9Y10Y11Y12Y13Y14Y1520061150744441233335112006216872034113323411200611494822142123200200621741922331311231020061128120231131112012006115416834133213500200621831442312222231020061157288351322234102006217896341131113102006211924111422111002006215238243131213102006124996111232112002006213928825113221200200611547224133211300200611511682123321120020061179821243211210200611272821323111111200622834821133211310200611584821123111210200611352411123111101200621662421113111100200622639623113221210

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

Результаты экспериментальных исследований НС приведены на рисунке 2, откуда следует, что минимальная ошибка обучения НС на интервале наблюдения, содержащем 12 отсчетов, обеспечивается при числе нейронов скрытого слоя равном шести.

Лучше всего сходимость процесса обучения сети обеспечивается при использовании квазиньютоновского алгоритма и алгоритма сопряженных градиентов.

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

При этом прогноз по методу скользящего среднего осуществлялся по формуле:

,(5)

где N — число предшествующих периодов, входящих в скользящее среднее; yt — фактическое значение в момент времени t; yt+1 — прогнозируемое значение в момент времени (t+1).

Для прогноза по методу экспоненциального сглаживания применялось выражение:

yt+1 = yt+ α(At –yt) = α • At + (1–α) • yt ,(6)

где yt+1 — предсказанное значение параметра на основе предыдущего значения yt , скорректированного с учетом погрешности прогноза At – yt и весового коэффициента α (0<α<1). Экспериментальные исследования решения задачи прогноза с помощью метода наименьших квадратов показали, что экстраполирующая функция f(t) может быть выбрана в виде полинома девятого порядка. На рисунке 3 показана поверхность Y = f(t), построенная с помощью среды моделирования MATLAB, где Y = (Y10, Y11)T, а время t  [2008; 2011] соответствует годам исследования.

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

Информация о конфликте интересов. Конфликт интересов отсутствует.

Информация о спонсорстве. Данная работа не финансировалась.