Нейросетевая модель прогнозирования компенсации сахарного диабета на основе индексов вариабельности гликемии
Самойлова Ю.Г.1,
Кошмелева М.В.1,
Кобякова О.С.1,
Бразовский К.С.1,
Назаренко В.В.1,
Толмачев И.В.1,
Олейник О.А.1,
Филиппова Т.А.1,
Подчиненова Д.В.1,
Сиволобова Т.В.1
'Федеральное государственное бюджетное образовательное учреждение высшего образования "Сибирский государственный медицинский университет" Министерства здравоохранения Российской Федерации (ФГБОУ ВО "СибГМУ" Минздрава России), Томск
Аннотация
Основным показателем компенсации сахарного диабета (СД) является определение уровня гликированного гемоглобина (HbA1c), который отображает лишь средний уровень гликемии. Непрерывный мониторинг гликемии (НМГ) выявляет суточные экскурсии гликемии - вариабельность гликемии (ВГ) и ее специфические индексы, что позволяет практическому врачу принять верное и правильное решение относительно тактики лечения пациента для достижения компенсации заболевания. Цель: провести сравнительный анализ гликированного гемоглобина и индексов вариабельности гликемии для прогнозирования степени компенсации сахарного диабета 1 типа. В исследование было включено 120 пациентов с СД 1 типа. Благодаря технологиям НМГ все пациенты отправляли гликемические отчеты врачу для получения рекомендаций. Обследуемым был проведен анализ HbA1c. С помощью калькулятора EasyGV рассчитывались индексы вариабельности гликемии, на основании которых была построена регрессионная нейросетевая модель в среде статистических вычислений типа R с помощью программного пакета Neuralnet. При использовании НМГ у пациентов с СД 1 типа отмечалось достоверное улучшение индексов вариабельности гликемии к концу исследования. HbA1c в группе дистанционного наблюдения снизился на 1,95% (р<0,000), в группе очного наблюдения уменьшился на 0,8% (р<0,000), в группе стандартной терапии увеличился на 0,05% при р=0,546. Дистанционный мониторинг подтверждает свою клинико-метаболическую эффективность в сравнении с другими моделями наблюдения пациентов с сахарным диабетом 1 типа. Оптимальной оказалась модель на основе многослойного перцептрона с высоким значением коэффициента детерминации R2=0,987, что свидетельствует о высокой достоверности прогнозирования компенсации заболевания, что позволяет спрогнозировать осложнения и обеспечить персонализированный подход в лечении пациентов с сахарным диабетом 1 типа. гемоглобин гликемия диабет
Ключевые слова: сахарный диабет 1 типа; дистанционный мониторинг; вариабельность гликемии; нейросетевая модель.
NEURAL NETWORK PREDICTION MODEL COMPENSATION OF DIABETES BASED ON THE INDICES GLYCEMIC VARIABILITY
Samoilova I.G.1, Koshmeleva M.V.1, Kobyakova O.S.1, Brazovskii K.S.1, Nazarenko V.V.1, Tolmachev I.V.1, Oleynik O.A.1, Filippova T.A.1, Podchinenova D.V.1, Sivolobova T.V.1
'Siberian State Medical University, Tomsk
The main indicator of compensation for diabetes mellitus is the determination of the level of glycated hemoglobin (HbA1c), which displays only the average level of glycemia. Continuous glucose monitoring (CGM) reveals daily glycemic excursions - glycemic variability (VG) and its specific indices, which allows the physician to make the correct decision of treatment of the patient to achieve compensation of the disease. Aim: To conduct a comparative analysis of glycated hemoglobin (HbA1c) and glycemic variability indexes to predict the degree of compensation for the diabetes mellitus type 1. The study included 120 patients with diabetes type 1. Through the technology of CGM all patients transmitted data to the doctor for recommendations and done analysis of HbA1c. Using the EasyGV calculator, we calculated the glycemic variability indices, on the basis of which the regression neural network model was built in the statistical computing environment R using the neuralnet. Using CGM in patients with diabetes type 1, there was a significant improvement in the glycemic variability indices by the end of the study. HbA1c in the remote monitoring group decreased by 1.95% (p<0,000), in the outpatient observation group decreased by 0.8% (p<0,000), in the standard therapy group increased by 0.05% with p=0.546. Remote monitoring confirms its clinical and metabolic efficacy in comparison with other models of observation of patients with diabetes type 1. The optimal model was based on a multilayer perceptron with a high value of the coefficient of determination R2=0.987, which indicates a high accuracy of the prediction of disease compensation, which allows to predict complications and provide a personalized approach in treating patients with diabetes type 1.
Keywords: diabetes mellitus type 1, remote monitoring, glycemic variability, neural network model.
Сахарный диабет (СД) стал эпидемией мирового масштаба, приводящей к серьезным микро- и макрососудистым осложнениям, что в дальнейшем снижает качество и продолжительность жизни, способствует развитию ранней инвалидизации пациентов. В связи с этим прогнозирование риска осложнений при сахарном диабете является одной из приоритетных задач диабетологии [1].
Многие исследования продемонстрировали, что плохой контроль уровня глюкозы является одним из факторов риска, непосредственно связанных с микрососудистыми осложнениями при СД 1 типа [2-4] и смертностью при СД 1 типа [5].
Основным показателем компенсации СД является определение уровня гликированного гемоглобина (HbA1c), который оценивается один раз в три месяца. Такая частота выполнения теста обусловлена международными стандартами достижения высокой степени контроля заболеванияи предотвращения микро- и макрососудистых осложнений [6; 7]. На каждый 1% улучшения показателя HbA1c - пропорционально снижается скорость развития поздних осложнений СД: инсультов, инфарктов, сердечной недостаточности, катаракты, микрососудистых осложнений и ампутаций, связанных с микроангиопатией [8]. За последние два десятилетия HbA1c стал целевым стандартом для измерения уровня общего гликемического контроля и оценки риска осложнений диабета, определяющим выбор терапевтической тактики врачом [9].
Несмотря на то, что HbA1c является основным маркером риска осложнений и смертности от диабета, он отображает лишь средний уровень гликемии и не дает возможности оценки колебаний в течение различных периодов времени. Эта реальность побудила клиницистов и исследователей изучить, есть ли другие показатели помимо HbA1c, которые могли бы спрогнозировать риск развития осложнений сахарного диабета и содействовать более эффективному принятию клинических решений [10].
Сегодня появились качественно новые технические возможности для анализа и контроля гликемических кривых у пациентов с сахарным диабетом 1 типа: инсулиновые помпы с функцией обратной связи, непрерывный мониторинг гликемии, флеш-мониторинг и др. Непрерывный мониторинг гликемии (НМГ) - наиболее объективный метод диагностики нарушений углеводного обмена и оценки степени его компенсации у больных сахарным диабетом. Система постоянного мониторинга глюкозы и данные, полученные с ее помощью, позволяют выявить определенные тенденции суточных колебаний гликемии вариабельность гликемии (ВГ), с регистрацией всех эпизодов гипо- и гипергликемий, выяснить причины их возникновения, которые невозможно оценить при определении гликемии с помощью глюкометра [11; 12].
В настоящее время пристальное внимание в диабетологии уделяется индексам вариабельности гликемии (ИВГ), и ряд исследований доказывает взаимосвязь данных индексов с уровнем контроля и компенсации сахарного диабета [13; 14].
Так как стандартные методы оценки уровня контроля и компенсации сахарного диабета не всегда точны, то именно анализ индексов вариабельности гликемии позволит практическому врачу принимать верное и правильное решение относительно тактики лечения пациента с СД для достижения компенсации заболевания.
Цель исследования: провести сравнительный анализ гликированного гемоглобина и индексов вариабельности гликемии для прогнозирования компенсации сахарного диабета 1 типа.
Материалы и методы исследования. Исследование являлось проспективным, открытым, контролируемым, клиническим, в параллельных группах и выполнялось на базе Федерального государственного бюджетного образовательного учреждения высшего образования "Сибирский государственный медицинский университет" Министерства здравоохранения Российской Федерации (ФГБОУ ВО "СибГМУ" Минздрава России), г. Томск, при поддержке Фонда развития филантропии (Договор № 57/ТФ-14/315 от 14.05.2014) [15].
В исследование было включено 120 пациентов с сахарным диабетом 1 типа, которым проводилось непрерывное мониторирование гликемии в начале и в конце исследования. Все участники исследования были разделены на три группы:
1. Группа А (n=40 чел.) - в группу включались пациенты, использующие инсулинотерапию в помповом режиме с возможностью дистанционного мониторинга.
2. Группа В (n=40 чел.) - данную группу составили пациенты, получающие инсулинотерапию в помповом режиме без возможности дистанционного мониторинга.
3. Группа С (n=40 чел.) - сравнительная группа - больные с базисно-болюсной инсулинотерапией.
Средний возраст всех обследуемых составил 13,6 (10-17) года. Из них девочек 63 чел. (52,5%), мальчиков 57 чел. (47,5%). Длительность заболевания варьировалась от 3 до 17 лет, и в среднем составила 7,3 года. Все пациенты были сопоставимы по полу и возрасту. Обследованные находились под наблюдением в течение 24 недель, каждые 12 недель проводилось 3 очных визита и 24 дистанционные консультации.
Все участники исследования подписали добровольное информированное согласие.
Протокол исследования одобрен Этическим комитетом ФГБОУ ВО "СибГМУ" Минздрава России (№ 5004 от 28.11.2016).
Больные в ходе исследования отправляли данные своих отчетов об уровне глюкозы крови, полученных благодаря различным устройствам НМГ и специализированного программного обеспечения. В начале и в конце исследования с помощью калькулятора EasyGV [16] рассчитывались индексы вариабельности гликемии, которые были выбраны в качестве независимых параметров прогнозирования уровня HbA1c: стандартное отклонение (SD), инде^ длительного повышения гликемии (CONGA), инде^ лaбильности гликемии (LI), индены ртета гипо- и гипергликемии (LBGI, HBGI), среднее значение общего риска (ADRR), средняя амплитуда колебаний гликемии (MAGE), качество контроля гликемии (M- value), скорость изменения гликемии (MAG). Характеристика ИВГ в начале исследования в каждой группе приведена в таблице 1.
Таблица 1
Характеристика индексов вариабельности гликемии у пациентов в начале исследования
|
Индексы вариабельности гликемии |
Группа дистанционного мониторинга (n=40) |
Группа очного наблюдения (n=40) |
Сравнительная группа (n=40) |
р 1-2 |
р 2-3 |
р 1-3 |
|
|
HbA1c, % |
9,1 [8,2; 9,7] |
9,5 [8,7; 10,9] |
9,05 [8,2; 10] |
0,061 |
0,279 |
0,134 |
|
|
Средний уровень гликемии (Mean), ммоль/л |
8.48±1.89 |
9.8±1.71 |
8.95±1.4 |
0,051 |
0,067 |
0,205 |
|
|
Стандартное отклонение (SD), ммоль/л |
2,57 [1,57; 3,38] |
2,75 [2,1; 3,58] |
2,65 [2,12; 3,29] |
0,367 |
0,378 |
0,820 |
|
|
Индекс длительного повышения гликемии (CONGA), ммоль/л |
9,19 [8,21; 10,75] |
10,21 [9,32; 11,63] |
9,72 [8,94; 11,9] |
0,062 |
0,577 |
0,204 |
|
|
Индекс лабильности гликемии (LI), (ммоль/л)2/час |
6,37 [4,92; 10,72] |
9,37 [7,24; 14,28] |
10,34 [9,21; 12,03] |
0,001* |
0,294 |
0,000* |
|
|
Индекс ртета гипогликемии (LBGI) |
5,88 [3,56; 8,05] |
7,33 [6,12; 9,33] |
3,41 [2,45; 5,08] |
0,004* |
0,000* |
0,000* |
|
|
Индекс риска гипергликемии (HBGI) |
12,17 [8,88; 14,7] |
11,24 [8,94; 14,28] |
17,75 [12,45; 20,15] |
0,476 |
0,000* |
0,000* |
|
|
Среднее значение общего риска (ADRR) |
28,38 [15,96; 36,95] |
24,76 [19,76; 48,84] |
26,26 [18,37; 31,23] |
0,115 |
0,168 |
0,199 |
|
|
Оценка межсуточной вариабельности гликемии (MODD), ммоль/л |
3,97 [3,25; 5,52] |
5,08 [3,59; 6,59] |
6,47 [4,5; 8,55] |
0,108 |
0,061 |
0,001* |
|
|
Средняя амплитуда колебаний гликемии (MAGE), ммоль/л |
5,12 [4,66; 7,15] |
7,32 [6,36; 8,26] |
7,81 [6,24; 9,14] |
0,000* |
0,366 |
0,000* |
|
|
Качество контроля (M- value), ммоль/л |
14,2 [4,54; 12,79] |
15,02 [9,88; 22,28] |
13,52 [10,79; 17,97] |
0,071 |
0,501 |
0,072 |
|
|
Скорость изменения гликемии (MAG), ммоль/л/час |
3,36 [1,5; 2,75] |
3,55 [2,12; 4,85] |
4,3 [2,53; 5,8] |
0,002* |
0,075 |
0,000* |
Примечание: * - значимость различий между группами (Z-критерий Манна-Уитни, различия значимы при p<0,05); Me - медиана, Q1; Q3 - нижний; верхний квартили.
Для оценки метаболической эффективности всем пациентам в начале и в конце исследования был проведен анализ HbA1c методом жидкостной хроматографии на анализаторе DS5 Glycomat (фирма DrewScientific, Нидерланды).
Для исследования взаимосвязи между ИВГ и уровнем компенсации сахарного диабета на основе HbA1c было построено две нейросетевые модели: регрессионная и классифицирующая. Регрессионная модель позволяет прогнозировать риск развития сосудистых осложнений на основе ИВГ. Классифицирующая модель дает возможность оценить уровень компенсации сахарного диабета. Обучение модели проводилось на всех включенных пациентах, 70 из которых были случайным образом выбраны для включения в обучающую выборку, 10 пациентов составили тестовую, данные которой не использовались для обучения нейронной сети.
Регрессионная нейросетевая модель была построена в среде статистических вычислений типа R с помощью программного пакета Neuralnet [17]. Структура модели была выбрана путем сравнения более чем 20 тыс. тестовых моделей, имеющих различное количество скрытых слоев и нейронов в этих слоях. Критерием выбора модели служила минимальная сложность, отсутствие переобучения и высокое значение коэффициента детерминации. При построении нейросетевого классификатора компенсации сахарного диабета на основании данных ВГ и концентрации HbA1c были проведены численные исследования на моделях без внутреннего слоя (линейный классификатор) с количеством внутренних слоев от 1 до 3.
Статистический анализ проводился с помощью программы SPSS 23.0 (IBM SPSS Statistics, США). Для оценки вида распределения количественных показателей использовался критерий Шапиро-Уилка. Описательная статистика для ненормально раотределенных количественных параметров представлена медианой и первым и третьим квартилями Me [Q1; Q3], для нормально распределенных количественных параметров приводились значения Хср±о. Статистическую значимость различий для независимых данных оценивали по U- критерию Манна-Уитни, для зависимых данных - по критерию Уилкоксона. Различия считали значимыми при p<0,05.
Результаты исследования и их обсуждение. В отчетах НМГ у пациентов группы А в конце исследования было отмечено достоверное снижение всех показателей ВГ по сравнению с началом исследования: индекса стандартного отклонения (SD), который в начале исследования составил 2,57 [1,57; 3,38] ммоль/л, а в конце - 2,08 [1,07; 2,59] ммоль/л (р<0,001); индекса длительного повышения гликемии (CONGA), отражающего внутрисуточные колебания глюкозы, имел тенденцию к снижению к концу исследования - 9,19 [8,21; 10,75] ммоль/л и 6,16 [5,62; 6,96] ммоль/л, р=0,001. Значительное снижение показал индекс среднего общего риска (ADRR) с 28,38 [15,96; 36,95] до 19,97 [13,34; 27,02] (р<0,000), в 1,5 раза, что свидетельствует об уменьшении риска возникновения микро- и макрососудистых осложнений у пациентов с сахарным диабетом 1 типа на дистанционном мониторинге. Показатель M-value, характеризующий качество контроля заболевания, на визите 1 составил 14,2 [4,54; 12,79] ммоль/л, а на визите 3 - 6,7 [3,23; 10,02] ммоль/л (р=0,005).
| 05_Холера |
| 10.4. Исследование регистров |
| 1112 |
| 12 |
| 1568 |
| 1673 |
| 17 |
| 1703 |
| 1885 |
| 2. Генераторы постоянного тока |