Сетевое издание

Современные проблемы науки и образования

ISSN 2070-7428

"Перечень" ВАК

ИФ РИНЦ = 1,006

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Тырсин А.Н. 1, 2 Улезко Е.А. 3 Остроушко Д.В. 3 Свирская О.Я. 3 Санковец Д.Н. 3

---

1 ФГАОУ ВО Уральский федеральный университет имени первого Президента России Б.Н. Ельцина

2 ФГАОУ ВО Южно-Уральский государственный университет (национальный исследовательский университет)

3 Государственное учреждение «Республиканский научно-практический центр "Мать и дитя"»

Целью исследования является разработка статистически достоверного метода оценки степени тяжести состояния недоношенных новорожденных со сроком гестации менее 30 недель при рождении с использованием многомерного статистического анализа данных. За основу была взята модель бинарной логистической регрессии. В качестве оценки степени тяжести новорожденного была использована вероятность его отнесения к соответствующей группе пациентов. Для решения поставленной задачи проведена оценка анамнестических данных и результатов лечебно-диагностических процедур 29 недоношенных новорожденных. Информативными признаками явились масса тела ребенка при рождении в граммах, метод родоразрешения, pH и уровень лактата крови КОС при рождении, минимальный уровень насыщения крови кислородом при проведении респираторной поддержки, зафиксированный за первые 24 часа жизни, гемодинамический статус, минимальное и максимальное значение общего содержания кислорода крови КОС первых 24 часов жизни, минимальное значение парциального давления кислорода крови КОС первых 24 часов жизни. По имеющейся выборке прецедентов разработано статистически достоверное решающее правило, на основе которого может быть создана методика распознавания степени тяжести состояния новых пациентов. Проведен сравнительный энтропийный анализ двух групп пациентов.

Статья в формате PDF

197 KB

недоношенные новорожденные

прогноз

модель

логистическая регрессия

дискриминантный анализ

энтропия

1. Ely D.M., Driscoll A.K. Infant mortality in the United States, 2017: data from the period linked birth/infant deathfile. Natl Vital Stat Rep. 2019. Vol. 68 (10). P. 1-15.

2. Harrison W., Goodman D. Epidemiologic Trends in Neonatal Intensive Care 2007-2012. JAMA Pediatr. 2015. Vol. 169 (9). P. 855-862.

3. Chaudhari S. Neonatal intensive care practices harmful to the developing brain. Indian Pediatr. 2011. Vol. 48 (6). P. 437-440.

4. Pursley D.M., Zupancic J.A.F. Using Neonatal Intensive Care Units More Wisely for At-Risk Newborns and Their Families. JAMA Netw Open. 2020. Vol. 3 (6). P. 205693.

5. Белкин А.А., Кузнецов Ю.В., Тырсин А.Н., Помазкин Ю.В. Программа для ЭВМ «Автоматизированная информационно-аналитическая экспертная система (АИАЭС) для оценки прогноза больных с острой церебральной недостаточностью» // Свидетельство о государственной регистрации программы для ЭВМ РФ № 2016615536. Правообладатель Общество с ограниченной ответственностью «Клиника Института Мозга». Заявл. 07.04.2016; опубл. 20.06.2016, реестр программ для ЭВМ.

6. Тырсин А.Н., Шалькевич Л.В., Остроушко Д.В., Шалькевич О.В., Геворгян Г.Г. Исследование перинатального поражения центральной нервной системы у детей в неонатальном периоде методами многомерного статистического анализа // Системный анализ и управление в биомедицинских системах. 2017. Т. 16. № 3. С. 595-605.

7. Sood B.J., McLaughlin K., Cortez J. Near-infrared spectroscopy. Applications in neonates. Seminars Fetal Neonatal Med. 2015. Vol. 20. P.164-172.

8. Kusaka T., Isobe K., Yasuda S. Evaluation of cerebral circulation and oxygen metabolism in infants using near-infrared light. Brain Dev. 2014. Vol. 36 (4). P. 277-283.

9. Мандель И.Д. Кластерный анализ. М.: Финансы и статистика. 1988. 176 с.

10. Hardle W.K., Simar L. Applied Multivariate Statistical Analysis. – 5th.ed. Springer. 2019. 558 p.

11. Hosmer D.W., Lemeshow S., Sturdivant R.X. Applied Logistic Regression. – 3d ed. Wiley. 2013. 397 p.

12. Васильев Н.П., Егоров А.А. Опыт расчета параметров логистической регрессии методом Ньютона–Рафсона для оценки зимостойкости растений // Математическая биология и биоинформатика. 2011. Том 6. № 2. С. 190-199. [Электронный ресурс]. URL: http://www.matbio.org/article_pdf.php?id=82 (дата обращения: 25.10.2020).

13. Тырсин А.Н. Энтропийное моделирование многомерных стохастических систем. Воронеж: Научная книга. 2016. 156 с.

14. Тырсин А.Н., Геворгян Г.Г. Программа вычисления энтропии случайного вектора // Свидетельство о государственной регистрации программы для ЭВМ РФ № 2017612851. Правообладатель Тырсин Александр Николаевич. Заявл. 31.10.2016; опубл. 03.03.2017, реестр программ для ЭВМ.

В результате внедрения современных технологий интенсивной терапии и реанимации в неонатологии ежегодно увеличивается число выживших недоношенных детей со сроком гестации менее 30 недель [1, 2]. Разработка новых подходов выхаживания и терапевтического ведения таких пациентов является весьма актуальной задачей, поскольку для недоношенных новорожденных характерен рост заболеваемости в неонатальном периоде. Для повышения эффективности выхаживания важно своевременно распознать тяжесть нарушений, происходящих в организме ребенка. Это позволит выбрать правильный режим лечения, снижающий риск возможных осложнений пациентов. Перинатальная патология у недоношенных детей характеризуется большей тяжестью и сочетанностью нарушений вследствие незрелости всех систем организма, развившейся на фоне отсутствия внутриутробного благополучия плода во время беременности [3, 4]. Это приводит к необходимости совершенствования методов оценивания состояния здоровья у таких детей.

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

Многие модели успешно зарекомендовали себя и внедрены в лечебный процесс больниц, они являются дополнительным подспорьем для специалиста при принятии решения о выборе той или иной тактики лечения, а также позволяют оценивать имеющийся потенциал здоровья пациента [5]. Другие модели позволяют достоверно распознать возможность развития заболевания в будущем и привлечь дополнительное внимание лечащего врача к конкретному пациенту [6]. Поэтому для выбора оптимального метода выхаживания и реабилитации недоношенных детей важно дать объективный прогноз состояния их здоровья.

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

Материалы и методы исследования

Для решения поставленной задачи проведена оценка анамнестических данных и результатов лечебно-диагностических процедур 29 недоношенных новорожденных со сроком гестации менее 30 недель при рождении. Наблюдение за пациентами и их обследование начинались непосредственно в родильном зале, а далее продолжались в отделении анестезиологии и реанимации (с палатами для новорожденных детей) ГУ «РНПЦ “Мать и дитя”», г. Минск. Средний гестационный возраст младенцев составил 28,25±1,03 недели. Масса тела при рождении составила 1094,48±205,63 г. Мальчиков было 15,00 (51,72%), девочек – 14 (48,27%). 20 (68,97%) новорожденных были от одноплодной беременности, 9 (31,03%) – от многоплодной. Через естественные родовые пути родился 1 (1/3,45%) ребенок, при абдоминальном родоразрешении – 28 (96,55%) младенцев (20 (68,97%) – при экстренном и 8 (27,58%) – при плановом родоразрешении). Все дети (29 (100,00%)) в родильном зале получили сурфактантзаместительную терапию в дозе 227,86±72,47 мг/кг. Среди основных клинических диагнозов отмечены следующие состояния: респираторный дистресс-синдром новорожденных, врожденная пневмония, инфекция, специфичная для перинатального периода неуточненная, врожденный сепсис, инфекция, специфичная для перинатального периода.

Проанализированы 20 показателей: анамнестические данные, кислотно-основное состояние (КОС) артериальной крови при рождении и первого дня жизни, спектроскопия в ближней инфракрасной области (Near-infrared spectroscopy, NIRS), респираторная поддержка новорожденного, гемодинамический статус и гидробаланс первых 24 часов жизни ребенка (табл. 1). Метод NIRS успешно применяется для оценки мозгового кровотока и объема крови в сосудах головного мозга как раннего предиктора перинатальных поражений центральной нервной системы у недоношенных новорожденных [7, 8]. Обследуемые пациенты нуждались в проведении респираторной поддержки. Время осмотра: для признаков P0 – P5 – в родильном зале, для признаков P6 – P19 – в первые 24 часа жизни.

Таблица 1

Анамнестические и лабораторные данные

Исследование включало в себя несколько этапов. Вначале необходимо у недоношенных детей найти число групп по степени тяжести состояния здоровья. Предварительная гипотеза состояла в том, что у недоношенных детей имеет место градация на две группы – «тяжелые» и «очень тяжелые». Необходимо было проверить справедливость предположения о двух уровнях степени тяжести. Здесь были использованы методы кластерного анализа данных [9]. На втором этапе была сформирована система информативных диагностических признаков. Эта процедура была выполнена с помощью дискриминантного анализа [10].

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

Дана обучающая выборка значений признаков

, , (1)

где – вектор значений i-го объекта, ; , – бинарная переменная, указывающая на принадлежность i-го объекта соответствующему классу. Например, для определенности первому классу при и второму – при ; m – число признаков у каждого объекта; n – количество наблюдений. Классификация выполняется с помощью логистической функции [11]:

, (2)

принимающей значения в интервале (0; 1). Пороговым значением является . Вектор в (2) задает разделяющую линейную границу, описываемую уравнением гиперплоскости .

Введем функцию . Зададим область D 1 возможных значений x для первого класса как , а для второго класса – как . Тогда и . Если x принадлежит гиперплоскости Õ, то . Следовательно, для произвольного наблюдения x* вероятность его отнесения к первому классу равна , а ко второму – .

Для оценивания вектора коэффициентов b был использован алгоритм Ньютона–Рафсона [12]. Представляет интерес также исследование системных свойств каждой из групп пациентов с помощью энтропийного моделирования [6, 13].

Методы многомерного статистического анализа реализованы с помощью статистического пакета программ Statistica. Энтропийное моделирование выполнялось с помощью программы вычисления энтропии случайного вектора [14].

Результаты исследования и их обсуждение

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

Дендрограмма

Слева на рисунке приводится расстояние между объединяемыми на данном этапе кластерами (объектами), внизу указаны признаки. С учетом результатов кластерного анализа все пациенты (n=29) были разделены на две группы: к группе «тяжелые» (D1) отнесены 17 пациентов, к группе «очень тяжелые» (D2) отнесены 12 пациентов.

Решение задачи состоит в попытке различить между собой группы пациентов D1 и D2 и в случае успеха построить решающее правило классификации. Для этого нужно было решить задачу многомерной классификации (распознавания) двух групп (кластеров) по показателям анамнестических данных и результатов лечебно-диагностических процедур. Суть решения состояла в нахождении такой совокупности показателей из исходного множества, которая позволила бы (если это в принципе возможно) статистически достоверно распознать различия в этих группах. Данная задача решалась в два этапа. На первом этапе с помощью дискриминантного анализа формировалась система информативных признаков, а на втором с помощью логистической регрессии строилось решающее правило классификации. Строилась дискриминантная функция с минимальным p-уровнем, который составил 0,0004. Информативными признаками явились следующие качественные и количественные показатели: X1 – масса тела ребенка при рождении в граммах; X2 – метод родоразрешения (0 – через естественные родовые пути, 1 – плановое кесарево сечение, 2 – экстренное кесарево сечение); X3 – pH КОС при рождении, X4 – уровень лактата крови КОС при рождении, X5 – минимальная SpO2 при проведении искусственной вентиляции легких, зафиксированная за первые 24 часа жизни, X6 – гемодинамический статус (0 – нет терапии; 1 – допамин ≤ 5 мкг/кг мин; 2 – допамин > 5 и ≤ 10 мкг/кг мин; 3 – допамин > 10 мкг/кг мин или допамин+добутамин ≤ 10 мкг/кг мин; 4 – допамин+добутамин > 10 мкг/кг мин; 5 – эпинефрин или норэпинефрин и/или глюкокортикостероиды; X7 – минимальное значение общего содержания кислорода (ctO2) крови КОС первых 24 часов жизни; X8 – максимальное значение ctO2 крови КОС первых 24 часов жизни; X9 – минимальное значение парциального давления кислорода (p50) крови КОС первых 24 часов жизни.

Распознавание двух групп (D1) и (D2) выполнено формально недостаточно достоверно, поскольку признаки pH крови КОС при рождении и p50 крови КОС первых 24 часов жизни имеют недостаточную статистическую надежность (менее 95%). Однако расчетное значение F-статистики статистически значимо со статистической надежностью 0,9996 (выше требуемой величины 0,95). Данные спектроскопии в ближней инфракрасной области (NIRS) оказались неинформативными показателями для решения поставленной задачи. Средняя вероятность правильного распознавания составила 96,552%, что выше требуемой величины 95%. Это означает, что сформированная система признаков в достаточной степени различает «тяжелых» и «очень тяжелых» пациентов.

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

Учитывая существенные признаки, полученные ранее с помощью дискриминантного анализа, мы рассчитали интегративный прогностический индекс (ИПИ) по формуле:

, (4)

где b0 = –11,591, b1 = 9,333 b2 = 6,808, b3 = –3,056, b4 = –14,060, b5 = 13,964, b6 = –11,353, b7 = 19,082, b8 = –18,484, b9 = –5,997.

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

Все пациенты правильно классифицированы. В таблице 2 приведены результаты распознавания. Здесь обозначены: № – номер наблюдения (порядковый номер ребенка в общей выборке; IND – принадлежность наблюдений к одной из групп (IND=0 – группа «тяжелые», IND=1 – группа «очень тяжелые»); P0 – вероятность отнесения наблюдения (ребенка) к группе «очень тяжелые» (IND=1); 1 – P0 – вероятность отнесения наблюдения (ребенка) к группе «тяжелые» (IND=0).

Таблица 2

Вероятности отнесения пациентов к группам «тяжелые» и «очень тяжелые»

Результаты энтропийного анализа по статистически значимым признакам приведены в таблице 3. Видно, что в группе «очень тяжелые» по сравнению с группой «тяжелые» энтропии хаотичности и самоорганизации меньше. Это означает, что в целом у «очень тяжелых» пациентов в целом вариация показателей, их изменчивость ниже, а теснота взаимосвязи между показателями выше.

Таблица 3

Результаты энтропийного анализа групп «тяжелые» и «очень тяжелые»

Заключение

Использование методов многомерного статистического анализа позволило разработать алгоритм оценки степени тяжести недоношенных детей. Достоверность классификации выборки пациентов по полученному разделяющему правилу составила 100%. Это позволяет использовать данное решающее правило для создания методики в целях прогнозирования степени тяжести состояния недоношенных новорожденных со сроком гестации менее 30 недель. Использование доступных анамнестических, а также лечебно-диагностических данных дает возможность имплементировать разработанное правило после его апробации в практическую деятельность неонатальных отделений интенсивной терапии различных уровней оказания перинатальной помощи, что явится дополнительным критерием для врача при оценке степени тяжести новорожденного пациента и дальнейшем выборе персонифицированной тактики терапевтических манипуляций и медикаментозного лечения.

Благодарности. Работа выполнена при поддержке совместного российско-белорусского проекта РФФИ (грант № 20-51-00001) и БРФФИ (грант № М20Р-008).

Библиографическая ссылка