Данное исследование представляет собой ретроспективный анализ, одобренный Комитетом по институциональной этике Первой народной больницы Наньтун. Исследование зарегистрировано на ClinicalTrials.gov (регистрационный номер: NCT06833099). Поскольку информация о состоянии здоровья всех участников была анонимизирована, информированное согласие не требовалось. Используемые расходные материалы и оборудование указаны в Таблице материалов.
1. Исследуемая популяция
В исследование были включены клинические данные и предоперационные визуализационные записи 436 пациентов, которым была выполнена чрескожная эндоскопическая интерламинарная дискэктомия (PEID) по поводу боли в пояснице и боли в ногах из-за грыжи диска L5-S1 в Первой народной больнице Наньтун в период с 1 января 2020 года по 30 июня 2024 года. На основании критериев исключения в исследование было включено 309 пациентов, среди которых рецидивирующая грыжа межпозвоночного диска поясничного отдела позвоночника (рЛДГ) была подтверждена с помощью послеоперационной шкалы визуальной аналоговой шкалы (ВАШ) и последующей визуализации. Из 309 пациентов у 33 наблюдался рLDH после операции, в то время как у остальных 276 наблюдалось значительное облегчение боли в пояснице и ногах, при этом баллы по шкале VAS снизились более чем на 60%. Визуализирующие переменные для всех участников были получены на основе предоперационных рентгенологических, компьютерных и магнитно-резонансных томографий в сочетании с подробной клинической информацией, включая пол, возраст, рост, вес, ИМТ, баллы по шкале ВАШ и другие соответствующие переменные (рис. 1).
2. Критерии включения и исключения
Критерии включения рЛДГ: (1) Пациенты с грыжей межпозвоночного диска поясничного отдела L5-S1, которым была выполнена односегментная ПЭИД. (2) Комплексные визуализирующие обследования, проводимые в течение одного месяца до операции. (3) Послеоперационное снижение балла по ВАШ ≥60% с последующим повышением балла и подтверждением с помощью визуализации. (4) При визуализации не обнаружено других отклонений. (5) Минимальный период наблюдения 6 месяцев.
Критерии включения для пациентов без рЛДГ: (1) Пациенты с грыжей межпозвоночного диска поясничного отдела L5-S1, которым была выполнена односегментная ПЭД. (2) Комплексные визуализирующие обследования, проводимые в течение одного месяца до операции. (3) Послеоперационное снижение балла по ВАШ ≥60%, без рецидивов. (4) При визуализации не обнаружено других отклонений. (5) Минимальный период наблюдения 6 месяцев.
Критерии исключения: (1) Наличие других патологических состояний, вызывающих боль в пояснице, таких как инфекция межпозвоночного диска, опухоли позвоночника, метаболические заболевания костей или остеопороз. (2) Операции на межпозвоночном диске поясничного отдела позвоночника или другие операции на позвоночнике в анамнезе. (3) Низкое качество визуализации или неполные данные обследования. (4) Потеря для последующего наблюдения.
3. Категориальные и непрерывные переменные
Предоперационные клинические характеристики и данные визуализации пациентов были статистически проанализированы и измерены (Таблица 1 и Таблица 2). Категориальные переменные использовались для различения основных характеристик заболевания, факторов образа жизни и других переменных, в то время как непрерывные переменные представляли собой конкретные измерения, описывающие физиологическое состояние пациентов и предоперационные изменения визуализации. Для снижения предвзятости были внедрены строгие меры контроля качества. Два радиолога и спинальные хирурги с более чем 10-летним клиническим опытом отвечали за статистический анализ и измерение данных визуализации. В сложных случаях два врача решали вопросы путем совместной консультации, чтобы обеспечить точность и последовательность обработки данных.
Категориальные переменные включали: пол, диабет, гипертонию, сердечно-сосудистые и цереброваскулярные заболевания (CCD), сколиоз, стеноз позвоночного канала, провоцирующие факторы (например, напряженная активность, растяжение связок, воздействие холода, удар или отсутствие явного триггера), продолжительность заболевания (более 6 месяцев или менее 6 месяцев), онемение или слабость, тип протрузии (локализованная или неограниченная), дегенерация диска (степени I, II, III и IV, V), дегенерация соседнего диска (степени I, II, III и IV, V), Модические изменения и кальцификация диска.
Непрерывные переменные включали: возраст, продолжительность операции, рост, вес, индекс массы тела (ИМТ), максимальный диаметр грыжи диска (MDHD), высоту заднего диска (PDH), расстояние между двумя центрами позвонков (DBTVC), угол ориентации фасеточных суставов (FJOA), угол диска (DA), угол крестцового наклона (SSA), угол поясничного лордоза (LLA) и индекс высоты заднего диска (PDHI, PDHI = PDH / DBTVC).
4. Очистка данных и выбор переменных
Сначала файл данных "zjkj3" был прочитан в R (версия 4.3.1, https://www.r-project.org/, платформа: x86_64-w64-mingw32/x64 (64-бит), Copyright (C) 2023 The R Foundation for Statistical Computing) с помощью функции read_excel и сохранен в переменной данных. В качестве кодировки символов была выбрана кодировка UTF-8 (по умолчанию в системе), которая является широко используемым методом кодирования, способным поддерживать письменные языки большинства стран мира. Затем функция select использовалась для отделения целевой переменной от переменных функции. Чтобы обеспечить воспроизводимость разбиения данных, было установлено фиксированное случайное начальное значение 3, а набор данных был разделен на обучающий набор (80%) и тестовый набор (20%). Наконец, была определена функция регрессии LASSO, и часть обучающих данных была случайным образом отобрана с помощью функции glmnet для соответствия регрессии LASSO. Оптимальный параметр регуляризации λ был определен с помощью регуляризации L1 и перекрестной проверки (по умолчанию 10-кратная), при этом гарантировалось, что между обучающим и тестовым наборами не было смешивания дубликатов данных о пациентах. Для определения оптимального значения λ была рассчитана ошибка перекрестной проверки для каждого значения λ.
5. Разработка и оценка модели
SVM (метод опорных векторов): эта модель сопоставляет обучающие данные с гиперплоскостью, максимизируя запас между двумя классами для прогнозирования цели. Конкретные настройки гиперпараметров для модели SVM в этом исследовании: ядро = "линейный", стоимость = 1. На дополнительном рисунке 1 изображены матрицы неточностей для SVM.
DT (дерево принятия решений): Эта модель рекурсивно разделяет данные для создания древовидной структуры для прогнозирования целевой переменной. Конкретная настройка гиперпараметров для модели DT в этом исследовании: max_depth = 3. На дополнительном рисунке 2 изображены матрицы несоответствий для DT.
ADA (AdaBoost): Эта модель объединяет несколько слабых учеников (обычно деревьев решений) для создания сильного ученика, улучшая производительность классификации. Конкретные настройки гиперпараметров для модели ADA в этом исследовании: n_estimators = 150, learning_rate = 0,1, начальное значение = 80. На дополнительном рисунке 3 показаны матрицы несоответствий для ADA.
LGBM (Light Gradient Boosting Machine): Эта модель использует алгоритм гистограммы для разделения непрерывных признаков, ускоряя процесс обучения и снижая нагрузку на память, что делает ее эффективной и быстрой средой градиентного бустинга, особенно подходящей для крупномасштабных наборов данных. Конкретные настройки гиперпараметров для модели LGBM в этом исследовании: num_leaves = 5, learning_rate = 0,05, n_estimators = 50. На дополнительном рисунке 4 изображены матрицы неточностей для LGBM.
RF (Random Forest): Эта модель строит несколько деревьев решений и объединяет их прогнозы для повышения точности и контроля переобучения. Конкретные настройки гиперпараметров для радиочастотной модели в этом исследовании: ntree = 310, mtry = 1, maxnodes = 10, max_depth = 1, seed = 80. На дополнительном рисунке 5 изображены матрицы неточностей для RF.
XGB (Extreme Gradient Boosting): Эта модель представляет собой усовершенствованный алгоритм, основанный на деревьях решений, использующий новый алгоритм обобщенного дерева решений градиентного бустинга для ускорения построения модели, с высокой применимостью для задач классификации и регрессии. Конкретные настройки гиперпараметров для модели XGB в этом исследовании: nrounds = 100, max_depth = 2, eta = 0,39, gamma = 0, colsample_bytree = 0,88, seed = 80. На дополнительном рисунке 6 показаны матрицы неточностей для XGB.
Производительность шести моделей оценивалась с использованием следующих показателей: тестовый набор ROC AUC, точность, чувствительность, специфичность, положительная прогностическая ценность (PPV), отрицательная прогностическая ценность (NPV) и оценка F1. Формула расчета следующая:22:

ТП: Истинный позитив; ТН: Истинный негатив; FP: ложное срабатывание; ФН: Ложноотрицательный.
6. Поиск по сетке и настройка гиперпараметров
Для оптимизации производительности модели был использован поиск по сетке для изучения различных комбинаций гиперпараметров. С помощью итерационной настройки и оценки была определена оптимальная комбинация гиперпараметров, позволяющая максимизировать значения ROC AUC как для поезда, так и для тестовых наборов, тем самым повышая общую производительность модели. Диапазоны скорректированных гиперпараметров представлены в таблице 3. Поиск по сетке представлен на Дополнительном рисунке 7, Дополнительном рисунке 8 и Дополнительной таблице 1, Дополнительной таблице 2, Дополнительной таблице 3 и Дополнительной таблице 4.
7. Ранжирование переменной важности
Наиболее эффективные модели — Random Forest (RF) и Extreme Gradient Boosting (XGB) — были выбраны для ранжирования переменной важности. Эти рейтинги помогают определить наиболее важные предикторы L5-S1 рЛДГ после ПЭЙД.
8. Статистический анализ
Все статистические анализы были выполнены с использованием программного обеспечения R. Категориальные переменные были представлены в процентах, в то время как непрерывные переменные были выражены в виде среднего ± стандартного отклонения. Клинические характеристики и параметры визуализации сравнивали в группе с рецидивом (без рЛДГ, n=276) и группой с рецидивом (рЛДГ, n=33).