Для математического анализа отобрано 78 исследуемых признаков: 12 параметров системной гемодинамики, 10 – регионарной гемодинамики, 56 – регионарной микроциркуляции. Выполненный нами математический анализ выявил прогностическую значимость 30 исследуемых признаков у двух групп пациентов: обучающая (1434 человека) и экзаменационная (78 человек). На основе обучающей группы сформирован прогностический алгоритм. Возраст пациентов данной группы составил от 54 до 84 лет (865 женщин, 569 мужчин). После выполнения статистического анализа полученных данных установлены закономерности, позволяющие создать систему выявления риска возникновения тромбоза ретинальных вен при геомагнитных возмущениях у пациентов с ГБ. Для решения этой задачи проанализирована зависимость параметров гемодинамики и случаев возникновения окклюзий ретинальных вен с учетом коэффициента риска:
(1)
где Y – уровень риска, k – номер признака, β – вес признака в модели, x – значение признака, xкр – критическое значение признака. Значения параметров модели, коэффициентов регрессии β оценивали по методу наименьших квадратов (табл. 1).
Таблица 1
Удельный вес изучаемых признаков и их критические значения
| № | Критическое значение признака | Вес фактора в модели | βкxкр |
| Показатели системной гемодинамики, Кр < 4 | |||
| 1 | Систолическое артериальное давление > 150 мм рт. ст. | 1,64 | 246,0 |
| 2 | Диастолическое артериальное давление > 86,5 мм рт. ст. | 1,06 | 91,69 |
| 3 | Частота сердечных сокращений > 83,5 уд./мин. | 1,02 | 85,17 |
| 4 | Индекс Керде > 16 | 0,98 | 15,68 |
| Показатели системной гемодинамики, Кр ≥ 4 | |||
| 5 | Повышение систолического артериального давления ≥ 15% | 1,42 | 21,3 |
| 6 | Повышение диастолического артериального давления ≥ 18% | 1,67 | 30,06 |
| 7 | Число сердечных сокращений < 63 уд./мин. | 1,37 | 86,31 |
| 8 | Уменьшение индекса Керде ≥ 120%, с положительных до отрицательных значений | 1,62 | 19,44 |
| 9 | Уменьшение продолжительности задержки дыхания выдохе ≥ 26,5% | 0,53 | 14,05 |
| Показатели регионарной гемодинамики, Кр < 4 | |||
| 10 | Линейная скорость кровотока в надблоковой артерии ≤ 11,5 см/с | 1,72 | 19,78 |
| 11 | Пульсационный индекс ≥ 4,9 | 1,26 | 6,17 |
| Показатели регионарной гемодинамики, Кр ≥ 4 | |||
| 12 | Снижение линейной скорости кровотока в надблоковой артерии ≥ 30,1% | 1,58 | 47,56 |
| 13 | Снижение объемной скорости кровотока в надблоковой артерии ≥ 46,4% | 0,46 | 21,34 |
| 14 | Увеличение пульсационного индекса в надблоковой артерии ≥ 18,2% | 1,35 | 24,57 |
| Показатели регионарной микроциркуляции, Кр < 4 | |||
| 15 | Снижение показателя микроциркуляции после пробы ≥ 32% | 1,74 | 55,68 |
| 16 | Уменьшение объемного кровенаполнения ткани после пробы ≥ 15,6% | 0,57 | 8,89 |
| 17 | Повышение нейрогенного тонуса после пробы ≥ 16% | 1,04 | 16,64 |
| 18 | Повышение миогенного тонуса после пробы ≥ 14,7% | 0,92 | 13,52 |
| 19 | Уменьшение показателя шунтирования после пробы ≥ 10,5% | 0,48 | 5,04 |
| Показатели регионарной микроциркуляции, Кр ≥ 4 | |||
| 20 | Увеличение показателя микроциркуляции после пробы ≥ 24,0% | 1,72 | 41,28 |
| 21 | Уменьшение показателя микроциркуляции после пробы ≥ 23% | 1,56 | 35,88 |
| 22 | Уменьшение сатурации кислорода после пробы ≥ 12% | 0,63 | 7,56 |
| 23 | Уменьшение амплитуды нейрогенных колебаний после пробы ≥ 20,3% | 0,68 | 13,80 |
| 24 | Увеличение амплитуды нейрогенных колебаний после пробы ≥ 15,6% | 0,94 | 14,66 |
| 25 | Уменьшение амплитуды миогенных колебаний после пробы ≥ 24,6% | 1,32 | 32,47 |
| 26 | Увеличение амплитуды миогенных колебаний после пробы ≥ 10,5% | 1,63 | 17,12 |
| 27 | Повышение нейрогенного тонуса после пробы ≥ 14% | 0,94 | 13,16 |
| 28 | Уменьшение нейрогенного тонуса после пробы ≥ 13% | 0,37 | 4,81 |
| 29 | Повышение миогенного тонуса после пробы ≥ 33% | 0,68 | 22,44 |
| 30 | Уменьшение миогенного тонуса после пробы ≥ 10,8% | 0,79 | 8,53 |
Наибольший удельный вес фактора в модели оказался у следующих параметров: при Кр < 4 – систолическое АД свыше 150 мм рт.ст.; линейная скорость кровотока в надблоковой артерии – от 11,5 см/сек и менее; снижение показателя микроциркуляции после пробы на 32% и более; при Кр ≥ 4 – повышение диастолического АД от 18% и более; уменьшение индекса Керде на 120% и более; увеличение показателя микроциркуляции после пробы от 24% и свыше; увеличение на 10,5% и более амплитуды миогенных колебаний после пробы [4–6].
В результате проведенных расчетов алгоритм прогнозирования можно представить следующим образом:
Для разработки полученной прогностической модели с выяснением количественного показателя отсечения и выяснения ее прогностических качеств использовалась методика ROC-анализа. Применялись показатели чувствительности, специфичности и площадь под ROC-кривой (рис. 1).
Рис. 1. ROC-кривая для обучающей группы пациентов.
Чувствительность модели отражает долю истинно положительных случаев, адекватно выявленных с помощью данного прогностического алгоритма. Специфичность характеризует долю истинно отрицательных случаев, которые также были правильно идентифицированы данным алгоритмом.
С учетом обязательных требований максимальной суммарной чувствительности и специфичности модели нами установлено конкретное количественное значение отсечения, которое оказалось равным 1,05. Чувствительность алгоритма оказалась высокой – 89%, специфичность – 94%. Поскольку площадь под ROC-кривой (AUC) составила 0,95, это означает, что данное прогнозирование характеризуется отличным качеством (табл. 2).
Таблица 2
Площадь под ROC-кривой для обучающей группы пациентов
Примечание: a – в непараметрическом случае; b – нулевая гипотеза: истинная площадь = 0,5.
| Площадь | Стд. ошибкаa | Асимптотическая знач.b | Асимптотический 95% доверительный интервал | |
| нижняя граница | верхняя граница | |||
| 0,948 | 0,011 | 0,000 | 0,927 | 0,970 |
Затем мы выполнили исследование эффективности разработанного прогностического алгоритма. Для этого использована экзаменационная группа из 78 пациентов в возрасте 49–85 лет (36 мужчин и 42 женщины), страдающих артериальной гипертензией свыше трех лет и систематически применяющих гипотензивные препараты. Из их числа у 21 пациента (27%, 13 женщин и 8 мужчин) выявлено, что значения показателя уровня риска составили более 1,05 (от 1,06 до 2,54). Это означало наличие у них высокого риска окклюзионных поражений ретинальных вен при Кр ≥ 4. При динамическом наблюдении обнаружено, что у 17 пациентов в сроки до трех лет развились окклюзии ретинальных вен.
Сформированный нами прогностический алгоритм выяснения риска окклюзионных поражений ретинальных вен характеризуется высокими показателями чувствительности и специфичности (рис. 2).
Рис. 2. ROC-кривая для экзаменационной группы пациентов.
Площадь под ROC-кривой (AUC) составила 0,93, что указывает на прогнозирование отличного качества (табл. 3).
Таблица 3
Площадь под ROC-кривой для экзаменационной группы пациентов
Примечание: a – в непараметрическом случае; b – нулевая гипотеза: истинная площадь = 0,5.
| Площадь | Стд. ошибкаa | Асимптотическая знач.b | Асимптотический 95% доверительный интервал | |
| нижняя граница | верхняя граница | |||
| 0,934 | 0,034 | 0,000 | 0,869 | 0,998 |
Таким образом, методика расчета риска окклюзионных поражений ретинальных вен при геомагнитных возмущениях у пациентов с артериальной гипертензией представляет высокоэффективный алгоритм и рекомендуется для практического применения в офтальмологической клинике.
Конфликт интересов: авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с публикацией настоящей статьи.
Источник финансирования: авторы заявляют о финансировании проведенного исследования из собственных средств.
Участие авторов:
Концепция и дизайн исследования – СЕЛ, ПНВ, КОВ
Сбор и обработка материала – ПНВ, ПЯЕ
Написание текста – ПНВ, СЕЛ
Редактирование – СЕЛ