$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Демографические данные
Всего 665 пациентов с переломами бедра соответствовали критериям включения и были зачислены в это исследование (рисунок 1). Среди них 44 пациента (6,6%) были диагностированы с ранним МКВТ, выявленным на УЗИ в течение 48 часов после травмы. Подробное сравнение демографических и клинических характеристик между группами с MCVT и не-группами изложено в таблице 2.
Регрессия LASSO и результаты RCS
Кандидатские переменные, выявленные как статистически значимые (p-значение < 0,01) в одномерном анализе, вместе с клинически значимыми параметрами, дополнительно открывались с помощью регрессии LASSO (см. рисунок 2). Было оценено одиннадцать кандидатных признаков, включая девять непрерывных переменных и две переменные типа «фиктивные переломы». Окончательная модель прогноза сохранила TFIA, общий холестерин (TC), время активированного частичного тромбопластина (APTT), NLR, PLR и TG. PLR показал пограничную статистическую значимость в многомерной модели (p-значение = 0,068), но был сохранён на основе скрининга LASSO и клинической правдоподобности. Для ослабления предположения о линейности был проведён анализ RCS для шести непрерывных предикторов с учётом остальных выявленных переменных (Таблица 3). Наблюдалась значимая нелинейная зависимость дозы и ответа между ТГ и риском MCVT (общее p-значение < 0,001). Как показано на рисунке 3, риск MCVT резко вырос на нижних и средних значениях TG, а затем увеличился более постепенно, с большей неопределённостью на более высоких уровнях TG. Анализ RCS проводился с использованием 3 узлов. Для TG узлы были размещены на 0,62 ммоль/л, 1,03 ммоль/л и 1,64 ммоль/л, что соответствует 10-му, 50-му и 90-му перцентилям распределения TG. Кажущаяся точка отсчёта около 1,03-1,04 ммоль/л следует рассматривать как исследовательскую, а не как валидированную клиническую пороговую границу.
Многомерный логистический регрессионный анализ
С учётом переменных, выявленных с помощью анализа LASSO и RCS, была построена окончательная многопеременная логистическая регрессионная модель. В этой модели TFIA, TC, APTT, NLR и PLR были включены как линейные термины, а TG был включен как ограниченный кубический сплайн для учёта его нелинейного эффекта (Таблица 4). В окончательной модели прогноза были сохранены шесть переменных: TFIA (OR = 1,026, 95% CI: 1,004–1,049, p-значение = 0,021), TC (OR = 0,244, 95% CI: 0,146–0,406, p-значение < 0,001), APTT (OR = 0,856, 95% CI: 0,764–0,959, p-значение = 0,007), NLR (OR = 1,189, 95% CI: 1,103–1,280, p-значение <0,001), PLR (OR = 1,004, 95% CI: 1.000–1.008, p-значение = 0.068), и TG, смоделированный с использованием RCS. Общий эффект TG в модели был очень значительным (p-значение < 0,001). На основе этих переменных была разработана клиническая номограмма для облегчения индивидуальной оценки риска (рисунок 4). Полные коэффициенты регрессии, перехват моделей, коэффициенты TG-сплайна и расположение узлов приведены в Таблице 1. Определения предикторов приведены в таблице 5.
Производительность модели
Дискриминативная производительность модели оценивалась с помощью анализа операционных характеристик приёмника (ROC) (рисунок 5), что дало AUC 0,894 (95% CI: 0,829–0,941). Калибровочная кривая показала общее согласование между предсказанными вероятностями и наблюдаемыми исходами (рисунок 6), что подтверждается баллом Бриера 0,041 (bootstrap 95% ДИ: 0,031–0,053). DCA предположила, что использование этой модели может дать более высокую чистую выгоду, чем стратегии «лечить всё» и «лечить ничего» при диапазоне пороговых вероятностей (рисунок 7). Были рассчитаны дополнительные калибровочные метрики. Видимый перехват калибровки был примерно 0, видимый наклон калибровки — 1.000, а калибровочный наклон — примерно 0. Наклон калибровки, скорректированный bootstrap, составил 0,896, что свидетельствует о умеренном оптимизме после внутренней валидации. Внутренняя валидация с 1 000 перевыборок bootstrap показала ограниченный оптимизм, с оптимизмом-скорректированным C-индексом 0,884 (bootstrap 95% CI: 0,837-0,946) по сравнению с видимым AUC 0,894. Калибровочная кривая с корректировкой смещения аппроксимировала идеальную линию со средней абсолютной ошибкой 0,013. Bootstrap 95% доверительные интервалы для оценки чистой выгоды DCA показаны на рисунке 7. Тест Хосмера-Лемешоу показал χ2 = 21,77, df = 8, p-значение = 0,005, что указывает на некоторое несоответствие между предсказанными и наблюдаемыми вероятностями между децилями риска. Поэтому калибровку следует интерпретировать осторожно и переоценивать во внешних когортах.
Доступность данных:
Де-идентификационный набор данных участников, заполненный чек-лист TRIPOD+AI и полный статистический код R, используемый для обработки данных, разработки модели, валидации и оценки производительности, предоставлены в виде дополнительных файлов 1–3. Эти материалы содержат все необходимые данные и аналитические ресурсы для воспроизведения результатов, представленных в данном исследовании.

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

Рисунок 2: Регрессия LASSO для выбора предиктора. График показывает кросс-валидированное биномиальное отклонение как функцию log(lambda) в 10-кратной перекрёстной валидации. Было оценено одиннадцать кандидатных признаков, включающих девять непрерывных переменных и две фиктивные переменные типа перелома. Верхние числа обозначают количество ненулевых коэффициентов при каждом лямбда-значении. Полоски ошибок указывают на стандартные ошибки. Вертикальные пунктирные линии обозначают lambda.min (0.00526) и lambda.1se (0.02558). Сокращения: CCI = Индекс сопутствующих заболеваний Чарлсона; TFIA = Время от травмы до госпитализации; WBC = лейкоциты; N = нейтрофил; L = лимфоциты; M = моноцит; HB = Гемоглобин; HCT = гематокрит; PLT = тромбоциты; NLR = соотношение нейтрофилов к лимфоцитам; PLR = соотношение тромбоцитов к лимфоцитам; TP = общий белок; Alb = альбумин; Tbil = общий билирубин; Креа = креатинин; TC = Общий холестерин; GFR = скорость фильтрации клубочков; IP = неорганический фосфат; INR = Международное нормированное отношение; Fbg = фибриноген; APTT = Активированное частичное время тромбопластина; TT = Время тромбина; PT = Протромбиновое время. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 3: RCS-анализ TG и раннего MCVT. Кривая показывает скорректированные OR для раннего MCVT в зависимости от уровня TG (ммоль/л). Модель была скорректирована для TC, NLR, TFIA, APTT и PLR и использовала 3 узла при 0,62 ммоль/л, 1,03 ммоль/л и 1,64 ммоль/л. Эталонным значением был медианный уровень TG (1,03 ммоль/л). Затемнённая область составляет 95% ДИ. ТГ был смоделирован как непрерывный нелинейный предиктор, и кривая не должна интерпретироваться как определяющая валидированный клинический порог. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 4: Номограмма для прогнозирования раннего MCVT, выявленная в течение 48 часов после травмы. Номограмма оценивает предсказанную вероятность раннего MCVT с использованием TFIA (h), TG (ммоль/л), TC (ммоль/л), APTT (секунды), NLR и PLR. Определения предикторов и единицы измерения приведены в Таблице 5. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 5: Кривая ROC окончательной модели предсказания. Кривая ROC показывает способность модели предсказания для ранних MCVT. AUC составил 0,894 (95% CI: 0,829-0,941). Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 6: Калибровочная кривая окончательной модели предсказания. График калибровки сравнивает предсказанные и наблюдаемые вероятности ранних MCVT. Идеальная линия представляет собой идеальную калибровку; Видимая кривая отражает производительность модели на наборе данных разработки; а кривая, скорректированная смещением, представляет собой калибровку, скорректированную bootstrap, после 1000 повторных выборок. Балл Брайера составил 0,041, видимый калибровочный перехват — примерно 0, видимый наклон калибровки — 1,000, а наклон калибровки, скорректированный bootstrap, — 0,896. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 7: DCA окончательной модели прогнозирования. DCA показывает чистую пользу номограммы по пороговым вероятностям. Кривая модели отражает чистую выгоду от использования номограммы; линия лечения всех предполагает, что у всех пациентов развивается MCVT; а линия «не лечить» предполагает, что у пациентов не развивается MCVT. Кривая следует интерпретировать как доказательство потенциальной клинической полезности только во время внутренней валидации, поскольку валидированный порог лечения не был установлен. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.
| Предиктор/член | Преобразование / определение | Коэффициент бета |
| Перехват | Перехват моделей | -1.621145541 |
| TFIA | Линейные; Часы от травмы до госпитализации | 0.025977241 |
| NLR | Линейные; Соотношение нейтрофилов к лимфоцитам | 0.172703853 |
| PLR | Линейные; Соотношение тромбоцитов к лимфоцитам | 0.004034446 |
| TG_linear | Линейный член RCS; TG в mmol/l | 6.24150855 |
| TG_rcs1 | Нелинейный член RCS; узлы — 0,62, 1,03, 1,64 ммоль/л | -4.814103739 |
| TC | Линейные; Общий холестерин в ммоль/л | -1.411933879 |
| APTT | Линейные; секунды | -0.155260603 |
| Узлы TG RCS | 0,62, 1,03, 1,64 ммоль/л | |
| Предсказанная вероятность | P = 1 / (1 + exp(-LP)) | |
ТАБЛИЦА 1: Полная спецификация модели предсказания. В этой таблице представлена полная модель прогнозирования, включая перехват, коэффициенты регрессии, коэффициенты ограниченного кубического сплайна для TG и расположение узлов, необходимые для реализации модели и независимой валидации.
| Переменная | Группа MCVT | Группа без MCVT | Значение p |
| В среднем ± SD или n (%) | (n = 44) | (n = 621) |
| Пол (мужский/женский) | 15/29 (34.1/65.9) | 224/397 (36.1/63.9) | 0,872a |
| Возраст | 79.68 ± 12.03 | 75.66 ± 15.65 | 0,041b |
| ИМТ | 22.05 ± 2.74 | 21.94 ± 3.12 | 0,796b |
| Гипертония (да/нет) | 15/29 (34.1/65.9) | 196/425 (31.6/68.4) | 0,739a |
| Диабет (да/нет) | 6/38 (13.6/86.4) | 91/530 (14.7/85.3) | 1.000a |
| Коронарная болезнь (да/нет) | 1/43 (2.3/97.3) | 29/592 (4.7/95.3) | 0,713a |
| Тип перелома (шейка бедренной кости/интертрохантерная/субтрохантерская) | 22/19/3 (50.0/43.2/6.8) | 462/153/6 (74.4/24.6/1.0) | <0,001c |
| CCI | 2.55 ± 2.12 | 2.31 ± 2.09 | 0,472b |
| TFIA (часы) | 25.58 ± 18.89 | 16.22 ± 14.60 | 0,002b |
| Количество лейкоцитов (x 10⁹/L) | 9.19 ± 2.92 | 8.90 ± 3.18 | 0,529b |
| N количество (*10⁹/L) | 7.48 ± 2.99 | 4.20 ± 3.11 | <0.001b |
| L (x 10⁹/L) | 0.99 (0.78-1.21) | 1.15 (0.92-1.52) | 0,003d |
| M (*10⁹/L) | 0.57 (0.44-0.71) | 0.59 (0.43-0.75) | 0,941d |
| Количество ГБ (г/л) | 105.13 ± 22.33 | 110.02 ± 21.16 | 0,165b |
| HCT | 32.07 ± 6.46 | 33.62 ± 6.09 | 0,127b |
| Количество PLT (x 10⁹/L) | 195.06 ± 69.21 | 165.82 ± 65.07 | 0,009b |
| NLR | 6.16 (4.30-12.05) | 3.21 (1.56-5.20) | <0.001d |
| PLR | 181.51 (126.19-260.71) | 131.76 (96.00-176.23) | <0.001d |
| TP (г/л) | 61.79 ± 7.25 | 62.33 ± 7.31 | 0,637b |
| Альб (г/л) | 35.77 ± 5.90 | 36.58 ± 4.91 | 0,381b |
| Тбил (мкмоль/л) | 11.20 (9.23-14.38) | 12.70 (9.00-16.70) | 0,117d |
| Креа (мкмоль/л) | 70.50 (58.80-89.50) | 67.00 (57.47-82.00) | 0,280d |
| Мочевина (mmol/L) | 6.57 (5.28-9.15) | 6.18 (4.92-8.14) | 0,121d |
| TG (mmol/L) | 1.24 (1.03-1.62) | 1.01 (0.77-1.27) | <0.001d |
| TC (mmol/L) | 3.11 ± 0.76 | 3,84 ± 0,95 | <0.001b |
| GFR (мл/мин) | 70.1 ± 25.32 | 78.04 ± 24.45 | 0,050b |
| IP (mmol/L) | 1.04 ± 0.24 | 1.06 ± 0.24 | 0,787b |
| INR | 1.04 (0.99-1.07) | 1.05 (1.01-1.10) | 0,094d |
| FBG (g/L) | 3.7 ± 1.01 | 4.04 ± 1.06 | 0,034b |
| APTT (вторая) | 28.16 ± 3.23 | 29,48 ± 3,88 | 0,012b |
| TT (второй) | 15.95 (15.30-16.95) | 16.10 (15.60-16.80) | 0,533d |
| Физиотерапия (второй) | 11.90 (11.40-12.50) | 12.27 (11.70-12.80) | 0,064d |
| D-димер (мкг/л) | 7186.00 (4,447.25-17,346.75) | 4780.90 (1,980.00-12,140.00) | 0,017d |
| Данные представлены в виде среднего ± SD, медианы (IQR) или n (%), в зависимости от необходимости. aТест Уэлча; bтест на сумму ранга Уилкоксона; cТест хи-квадрата; dТочный тест Файшера (точный тест Фишера-Фримана-Халтона использовался для многокатегорийных переменных, когда ожидаемое количество ячеек было незначительным). |
ТАБЛИЦА 2: Исходные характеристики пациентов с MCVT и без него. Непрерывные переменные представлены как среднее ± стандартное отклонение или медиана с межквартильным диапазоном, в зависимости от необходимости. Категориальные переменные представлены в виде чисел и процентов. В итоговой аналитической когорте не было отсутствующих значений кандидатов-предикторов или переменных исхода; 12 обследованных пациентов были исключены из-за отсутствия ключевых клинических или лабораторных данных.
| Переменная | Узлы (10-й, 50-й, 90-й перцентили) | Нелинейное p-значение | Общее p-значение |
| TC | 2.61, 3.77, 4.92 | 0.369 | <0.001 |
| NLR | 0.708, 3.35, 8.68 | 0.058 | <0.001 |
| TG | 0.62, 1.03, 1.64 | <0.001 | <0.001 |
| TFIA | 2, 12, 48 | 0.327 | 0.072 |
| APTT | 25.5, 28.9, 33.9 | 0.554 | 0.003 |
| PLR | 74, 134, 250 | 0.923 | 0.163 |
ТАБЛИЦА 3: Нелинейные ассоциации, оцениваемые с помощью анализа RCS. В таблице представлены расположения узлов, общие p-значения и нелинейные p-значения для выбранных предикторов. Модели RCS использовали 3 узла, расположенные на 10-м, 50-м и 90-м процентилях каждого распределения предиктора.
| Переменные | ИЛИ | CI-lower | CI-аппер | p-значение |
| TFIA | 1.026 | 1.004 | 1.049 | 0.021 |
| TC | 0.244 | 0.146 | 0.406 | <0.001 |
| APTT | 0.856 | 0.764 | 0.959 | 0.007 |
| NLR | 1.189 | 1.103 | 1.28 | <0.001 |
| PLR | 1.004 | 1 | 1.008 | 0.068 |
| TG (моделируется с использованием RCS) | - | - | - | <0.001 |
ТАБЛИЦА 4: Многомерная логистическая регрессионная модель для ранних MCVT. Для линейных членов показаны OR, 95% CI и p-значения. Поскольку TG моделировалось с использованием RCS, один OR и 95% CI не были зафиксированы. Скорректированное соотношение дозы и ответа показано на рисунке 3, а коэффициенты сплайна и расположение узлов приведены в таблице 1.
| Предиктор | Определение | Подразделение | Время |
| TFIA | Время от травмы до госпитализации | Часы работы | Анализ крови при поступлении |
| TG | Триглицериды | mmol/l | Анализ крови при поступлении |
| TC | Общий холестерин | mmol/l | Анализ крови при поступлении |
| APTT | Время активированного частичного тромбопластина | секунды | Анализ крови при поступлении |
| NLR | Количество нейтрофилов / лимфоцитов | Соотношение | Анализ крови при поступлении |
| PLR | Количество тромбоцитов / лимфоцитов | Соотношение | Анализ крови при поступлении |
ТАБЛИЦА 5: Определения переменных модели предсказания. В этой таблице приведены предикторы, включённые в итоговую модель, а также их определения, единицы измерения и, где применимо, клинические интерпретации.
Дополнительный файл 1: Заполненный чек-лист отчётности TRIPOD+AI. Этот файл содержит завершённый контрольный список «Прозрачная отчетность многопеременной модели предсказания для индивидуального прогноза или диагностики плюс искусственного интеллекта» (TRIPOD+AI), документирующий соответствие рекомендациям по отчётности для исследований моделей прогнозирования. Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 2: Результаты статистического анализа для разработки и валидации моделей прогнозирования. Эта рабочая тетрадь содержит результаты статистического анализа, полученные в ходе разработки и валидации модели, включая одномерные анализы, результаты регрессии LASSO, анализы ограниченных кубических сплайнов, коэффициенты многомерной логистической регрессии, метрики производительности модели, калибровочные анализы, результаты тестов Hosmer-Lemeshow на соответствие, анализ кривой принятия решений и связанные сводные таблицы. Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 3: скрипт R для обработки данных и статистического анализа. Этот файл содержит полный R-скрипт, используемый для предварительной обработки данных, выбора переменных, разработки модели, внутренней валидации, калибровки, анализа дискриминации, анализа кривой принятия решений и генерации статистических результатов, представленных в этом исследовании. Пожалуйста, нажмите здесь, чтобы скачать этот файл.