Este estudio es un análisis retrospectivo aprobado por el Comité de Ética Institucional del Nantong First People's Hospital. El ensayo ha sido registrado en ClinicalTrials.gov (Número de registro: NCT06833099). Dado que toda la información de salud de los participantes fue anónima, no se requirió el consentimiento informado. Los consumibles y equipos utilizados se enumeran en la Tabla de Materiales.
1. Población de estudio
El estudio incluyó datos clínicos y registros de imágenes preoperatorias de 436 pacientes que se sometieron a discectomía interlaminar (PEID) endoscópica percutánea por dolor lumbar y dolor en las piernas debido a una hernia discal L5-S1 en el Nantong First People's Hospital entre el 1 de enero de 2020 y el 30 de junio de 2024. Con base en los criterios de exclusión, finalmente se incluyeron 309 pacientes, entre los cuales se confirmó la hernia discal lumbar recurrente (rLDH) mediante las puntuaciones postoperatorias de la Escala Visual Analógica (EVA) y las imágenes de seguimiento. De los 309 pacientes, 33 experimentaron rLDH después de la cirugía, mientras que los 276 restantes tuvieron un alivio significativo del dolor lumbar y de piernas, con puntuaciones de EVA reducidas en más del 60%. Las variables de imagen de todos los participantes se derivaron de exámenes preoperatorios de radiografías, tomografías computarizadas y resonancias magnéticas, combinados con información clínica detallada, como sexo, edad, altura, peso, IMC, puntuaciones de la EVA y otras variables relevantes (Figura 1).
2. Criterios de inclusión y exclusión
Criterios de inclusión para rLDH: (1) Pacientes con hernia discal lumbar L5-S1 sometidos a PEID de un solo segmento. (2) Exámenes de imagen completos completados dentro de un mes antes de la operación. (3) Reducción de la puntuación EVA postoperatoria ≥60%, seguida de un aumento de la puntuación y confirmación por imagen. (4) No se detectaron otras anomalías en las imágenes. (5) Periodo mínimo de seguimiento de 6 meses.
Criterios de inclusión para la no rLDH: (1) Pacientes con hernia discal lumbar L5-S1 sometidos a PEID de un solo segmento. (2) Exámenes de imagen completos completados dentro de un mes antes de la operación. (3) Reducción de la puntuación EVA postoperatoria ≥60%, sin recurrencia. (4) No se detectaron otras anomalías en las imágenes. (5) Periodo mínimo de seguimiento de 6 meses.
Criterios de exclusión: (1) Presencia de otras condiciones patológicas que causan dolor lumbar, como infección discal, tumores espinales, enfermedades óseas metabólicas u osteoporosis. (2) Antecedentes de cirugía de disco lumbar u otras cirugías de columna. (3) Mala calidad de las imágenes o datos de examen incompletos. (4) Pérdida de seguimiento.
3. Variables categóricas y continuas
Se analizaron y midieron estadísticamente las características clínicas preoperatorias y los datos imagenológicos de los pacientes (Tabla 1 y Tabla 2). Las variables categóricas se utilizaron para distinguir las características básicas de la enfermedad, los factores del estilo de vida y otras variables, mientras que las variables continuas representaron mediciones específicas que describían el estado fisiológico de los pacientes y los cambios en las imágenes preoperatorias. Para reducir el sesgo, se implementaron estrictas medidas de control de calidad. Dos radiólogos y cirujanos de columna vertebral con más de 10 años de experiencia clínica fueron responsables del análisis estadístico y la medición de los datos de imagen. En el caso de los casos complejos, los dos médicos resolvieron los problemas mediante consultas conjuntas para garantizar la precisión y la coherencia en el procesamiento de los datos.
Las variables categóricas incluyeron: sexo, diabetes, hipertensión, enfermedades cardiovasculares y cerebrovasculares (ECC), escoliosis, estenosis espinal, factores desencadenantes (p. ej., actividad extenuante, esguince, exposición al frío, impacto o ausencia de un desencadenante claro), duración de la enfermedad (más de 6 meses o menos de 6 meses), entumecimiento o debilidad, tipo de protrusión (contenida o no contenida), degeneración discal (grados I, II, III y IV, V), degeneración discal adyacente (grados I, II, III y IV, V), cambios módicos y calcificación discal.
Las variables continuas fueron: edad, duración quirúrgica, talla, peso, índice de masa corporal (IMC), diámetro máximo de la hernia discal (MDHD), altura del disco posterior (PDH), distancia entre dos centros de vértebras (DBTVC), ángulo de orientación de la articulación facetaria (FJOA), ángulo del disco (DA), ángulo de la pendiente sacra (SSA), ángulo de la lordosis lumbar (LLA) e índice de altura del disco posterior (PDHI, PDHI = PDH / DBTVC).
4. Limpieza de datos y selección de variables
En primer lugar, el archivo de datos "zjkj3" se leyó en R (versión 4.3.1, https://www.r-project.org/, Plataforma: x86_64-w64-mingw32/x64 (64 bits), Copyright (C) 2023 The R Foundation for Statistical Computing) utilizando la función read_excel y se almacenó en una variable de datos. La codificación de caracteres seleccionada fue UTF-8 (predeterminada del sistema), que es un método de codificación ampliamente utilizado capaz de soportar los idiomas escritos de la mayoría de los países del mundo. A continuación, se utilizó la función select para separar la variable objetivo de las variables de característica. Para garantizar la reproducibilidad de la división de datos, se estableció una semilla aleatoria fija de 3 y el conjunto de datos se dividió en un conjunto de entrenamiento (80 %) y un conjunto de prueba (20 %). Finalmente, se definió una función de regresión LASSO, y una parte de los datos de entrenamiento se muestreó aleatoriamente utilizando la función glmnet para ajustar la regresión LASSO. El parámetro de regularización óptimo λ se determinó mediante regularización L1 y validación cruzada (por defecto 10 veces), al tiempo que se garantizó que no se mezclaran datos de pacientes duplicados entre los conjuntos de entrenamiento y prueba. Se calculó el error de validación cruzada para cada valor λ para determinar el valor λ óptimo.
5. Desarrollo y evaluación de modelos
SVM (Support Vector Machine): Este modelo mapea los datos de entrenamiento a un hiperplano, maximizando el margen entre dos clases para predecir el objetivo. La configuración específica de hiperparámetros para el modelo SVM en este estudio es: kernel = "lineal", costo = 1. La Figura 1 complementaria muestra las matrices de confusión para SVM.
DT (Árbol de decisión): Este modelo divide los datos de forma recursiva para crear una estructura de árbol para predecir la variable objetivo. La configuración de hiperparámetros específica para el modelo DT en este estudio es: max_depth = 3. La Figura 2 suplementaria muestra las matrices de confusión para DT.
ADA (AdaBoost): Este modelo combina varios aprendices débiles (generalmente árboles de decisión) para crear un aprendiz fuerte, mejorando el rendimiento de la clasificación. Los ajustes de hiperparámetros específicos para el modelo ADA en este estudio son: n_estimators = 150, learning_rate = 0,1, seed = 80. La Figura 3 suplementaria muestra las matrices de confusión para ADA.
LGBM (Light Gradient Boosting Machine): Este modelo utiliza un algoritmo de histograma para dividir características continuas, lo que acelera el proceso de entrenamiento y reduce la carga de memoria, lo que lo convierte en un marco de aumento de gradiente eficiente y rápido, particularmente adecuado para conjuntos de datos a gran escala. Los ajustes de hiperparámetros específicos para el modelo LGBM en este estudio son: num_leaves = 5, learning_rate = 0,05, n_estimators = 50. La Figura 4 suplementaria muestra las matrices de confusión para LGBM.
RF (Random Forest): Este modelo construye múltiples árboles de decisión y combina sus predicciones para mejorar la precisión y controlar el sobreajuste. Los ajustes de hiperparámetros específicos para el modelo de RF en este estudio son: ntree = 310, mtry = 1, maxnodes = 10, max_depth = 1, seed = 80. La Figura 5 suplementaria muestra las matrices de confusión para RF.
XGB (Extreme Gradient Boosting): Este modelo es un algoritmo mejorado basado en árboles de decisión, que utiliza un nuevo algoritmo de árbol de decisión generalizado de aumento de gradiente para acelerar la construcción del modelo, con una fuerte aplicabilidad para tareas de clasificación y regresión. Los ajustes de hiperparámetros específicos para el modelo XGB en este estudio son: nrounds = 100, max_depth = 2, eta = 0,39, gamma = 0, colsample_bytree = 0,88, seed = 80. La Figura 6 complementaria muestra las matrices de confusión para XGB.
El rendimiento de los seis modelos se evaluó utilizando las siguientes métricas: AUC del ROC del conjunto de pruebas, precisión, sensibilidad, especificidad, valor predictivo positivo (PPV), valor predictivo negativo (NPV) y puntuación F1. La fórmula de cálculo es la siguiente22:

TP: Verdadero Positivo; TN: Verdadero negativo; FP: Falso positivo; FN: Falso negativo.
6. Búsqueda de cuadrícula y ajuste de hiperparámetros
Para optimizar el rendimiento del modelo, se utilizó la búsqueda en cuadrícula para explorar varias combinaciones de hiperparámetros. A través del ajuste y la evaluación iterativos, se identificó la combinación óptima de hiperparámetros para maximizar los valores de AUC de ROC tanto para el conjunto de entrenamiento como para el de prueba, mejorando así el rendimiento general del modelo. Los rangos de hiperparámetros ajustados se presentan en la Tabla 3. La búsqueda en la cuadrícula se proporciona en la Figura complementaria 7, la Figura complementaria 8 y la Tabla complementaria 1, la Tabla complementaria 2, la Tabla complementaria 3 y la Tabla complementaria 4.
7. Clasificación de importancia variable
Los modelos de mejor rendimiento, Random Forest (RF) y Extreme Gradient Boosting (XGB), fueron seleccionados para la clasificación de importancia variable. Estas clasificaciones ayudan a identificar los predictores más críticos de L5-S1 rLDH después de la PEID.
8. Análisis estadístico
Todos los análisis estadísticos se realizaron con el software R. Las variables categóricas se reportaron como porcentajes, mientras que las variables continuas se expresaron como media ± desviación estándar. Se compararon las características clínicas y los parámetros imagenológicos entre el grupo no recurrente (no rLDH, n=276) y el grupo recurrente (rLDH, n=33).