Cette étude est une analyse rétrospective approuvée par le Comité d’éthique institutionnelle de l’Hôpital First People’s de Nantong. L’essai a été enregistré sur ClinicalTrials.gov (numéro d’enregistrement : NCT06833099). Comme tous les renseignements sur la santé des participants étaient anonymisés, le consentement éclairé n’était pas requis. Les consommables et équipements utilisés sont répertoriés dans la table des matériaux.
1. Population étudiée
L’étude comprenait des données cliniques et des dossiers d’imagerie préopératoire de 436 patients ayant subi une discectomie interlaminaire endoscopique percutanée (PEID) pour des douleurs lombaires et des douleurs à la jambe dues à une hernie discale L5-S1 à l’Hôpital First People’s de Nantong entre le 1er janvier 2020 et le 30 juin 2024. Sur la base des critères d’exclusion, 309 patients ont finalement été inclus, parmi lesquels la hernie discale lombaire récurrente (rLDH) a été confirmée par les scores postopératoires de l’échelle visuelle analogique (EVA) et l’imagerie de suivi. Sur les 309 patients, 33 ont présenté une rLDH après la chirurgie, tandis que les 276 autres ont eu un soulagement significatif des douleurs lombaires et des jambes, avec des scores EVA réduits de plus de 60 %. Les variables d’imagerie de tous les participants ont été dérivées des radiographies, des tomodensitogrammes et des examens IRM préopératoires, combinées à des informations cliniques détaillées, notamment le sexe, l’âge, la taille, le poids, l’IMC, les scores EVA et d’autres variables pertinentes (figure 1).
2. Critères d’inclusion et d’exclusion
Critères d’inclusion pour la rLDH : (1) Patients atteints d’une hernie discale lombaire L5-S1 ayant subi une PEID à segment unique. (2) Examens d’imagerie complets effectués dans un délai d’un mois avant l’opération. (3) Réduction du score VAS postopératoire ≥60 %, suivie d’une augmentation du score et d’une confirmation par imagerie. (4) Aucune autre anomalie détectée à l’imagerie. (5) Période de suivi minimale de 6 mois.
Critères d’inclusion pour la non-rLDH : (1) Patients atteints d’une hernie discale lombaire L5-S1 ayant subi une PEID à segment unique. (2) Examens d’imagerie complets effectués dans un délai d’un mois avant l’opération. (3) Réduction du score EVA postopératoire ≥60 %, sans récidive. (4) Aucune autre anomalie détectée à l’imagerie. (5) Période de suivi minimale de 6 mois.
Critère d’exclusion : (1) Présence d’autres conditions pathologiques causant des douleurs lombaires, telles que l’infection discale, les tumeurs de la colonne vertébrale, les maladies osseuses métaboliques ou l’ostéoporose. (2) Antécédents de chirurgies discales lombaires ou d’autres chirurgies de la colonne vertébrale. (3) Mauvaise qualité d’imagerie ou données d’examen incomplètes. (4) Perte de suivi.
3. Variables catégorielles et continues
Les caractéristiques cliniques préopératoires et les données d’imagerie des patients ont été analysées et mesurées statistiquement (tableau 1 et tableau 2). Des variables catégorielles ont été utilisées pour distinguer les caractéristiques de base de la maladie, les facteurs de mode de vie et d’autres variables, tandis que les variables continues représentaient des mesures spécifiques décrivant l’état physiologique des patients et les changements d’imagerie préopératoires. Pour réduire les biais, des mesures strictes de contrôle de la qualité ont été mises en place. Deux radiologues et chirurgiens de la colonne vertébrale avec plus de 10 ans d’expérience clinique étaient responsables de l’analyse statistique et de la mesure des données d’imagerie. Pour les cas complexes, les deux médecins ont résolu les problèmes par le biais d’une consultation conjointe afin d’assurer l’exactitude et la cohérence du traitement des données.
Les variables catégorielles comprenaient : le sexe, le diabète, l’hypertension, les maladies cardiovasculaires et cérébrovasculaires (DCC), la scoliose, la sténose spinale, les facteurs déclenchants (p. ex., activité intense, entorse, exposition au froid, impact ou aucun déclencheur clair), la durée de la maladie (plus de 6 mois ou moins de 6 mois), l’engourdissement ou la faiblesse, le type de protrusion (contenue ou non contenue), la dégénérescence discale (grades I, II, III et IV, V), la dégénérescence discale adjacente (grades I, II, III et IV, V), modifications modiques et calcification discale.
Les variables continues comprenaient : l’âge, la durée chirurgicale, la taille, le poids, l’indice de masse corporelle (IMC), le diamètre maximal de la hernie discale (MDHD), la hauteur du disque postérieur (PDH), la distance entre deux centres de vertèbres (DBTVC), l’angle d’orientation de l’articulation fatictale (FJOA), l’angle du disque (DA), l’angle de pente sacrée (SSA), l’angle de lordose lombaire (LLA) et l’indice de hauteur du disque postérieur (PDHI, PDHI = PDH / DBTVC).
4. Nettoyage des données et sélection des variables
Tout d’abord, le fichier de données « zjkj3 » a été lu dans R (version 4.3.1, https://www.r-project.org/, Plate-forme : x86_64-w64-mingw32/x64 (64 bits), Copyright (C) 2023 The R Foundation for Statistical Computing) à l’aide de la fonction read_excel et stocké dans une variable de données. L’encodage de caractères sélectionné était UTF-8 (système par défaut), qui est une méthode d’encodage largement utilisée capable de prendre en charge les langues écrites de la plupart des pays du monde. La fonction select a ensuite été utilisée pour séparer la variable cible des variables de fonctionnalité. Pour garantir la reproductibilité de la répartition des données, une valeur de départ aléatoire fixe de 3 a été définie et l’ensemble de données a été divisé en un ensemble d’apprentissage (80 %) et un ensemble de test (20 %). Enfin, une fonction de régression LASSO a été définie, et une partie des données d’apprentissage a été échantillonnée de manière aléatoire à l’aide de la fonction glmnet pour s’adapter à la régression Lasso. Le paramètre de régularisation optimal λ a été déterminé à l’aide de la régularisation L1 et de la validation croisée (par défaut 10 fois), tout en veillant à ce qu’aucune donnée de patient en double ne soit mélangée entre les ensembles d’entraînement et de test. L’erreur de validation croisée pour chaque valeur λ a été calculée pour déterminer la valeur λ optimale.
5. Élaboration et évaluation du modèle
SVM (Support Vector Machine) : Ce modèle mappe les données d’entraînement à un hyperplan, maximisant la marge entre deux classes pour prédire la cible. Les paramètres d’hyperparamètres spécifiques au modèle SVM dans cette étude sont les suivants : noyau = « linéaire », coût = 1. La figure supplémentaire 1 illustre les matrices de confusion pour la SVM.
DT (Decision Tree) : Ce modèle divise les données de manière récursive pour créer une structure arborescente permettant de prédire la variable cible. Le paramètre d’hyperparamètre spécifique pour le modèle DT dans cette étude est : max_depth = 3. La figure supplémentaire 2 illustre les matrices de confusion pour DT.
ADA (AdaBoost) : Ce modèle combine plusieurs apprenants faibles (généralement des arbres de décision) pour créer un apprenant fort, améliorant ainsi les performances de classification. Les paramètres d’hyperparamètres spécifiques pour le modèle ADA dans cette étude sont les suivants : n_estimators = 150, learning_rate = 0,1, graine = 80. La figure supplémentaire 3 illustre les matrices de confusion pour l’ADA.
LGBM (Light Gradient Boosting Machine) : Ce modèle utilise un algorithme d’histogramme pour diviser les caractéristiques continues, accélérant le processus d’entraînement et réduisant la charge de la mémoire, ce qui en fait un cadre de boosting de gradient efficace et rapide, particulièrement adapté aux ensembles de données à grande échelle. Les paramètres d’hyperparamètres spécifiques pour le modèle LGBM dans cette étude sont : num_leaves = 5, learning_rate = 0,05, n_estimators = 50. La figure supplémentaire 4 illustre les matrices de confusion pour LGBM.
RF (Random Forest) : Ce modèle construit plusieurs arbres de décision et combine leurs prédictions pour améliorer la précision et contrôler le surapprentissage. Les paramètres d’hyperparamètres spécifiques pour le modèle RF dans cette étude sont les suivants : ntree = 310, mtry = 1, maxnodes = 10, max_depth = 1, seed = 80. La figure supplémentaire 5 illustre les matrices de confusion pour les RF.
XGB (Extreme Gradient Boosting) : Ce modèle est un algorithme amélioré basé sur des arbres de décision, utilisant un nouvel algorithme d’arbre de décision généralisé pour accélérer la construction du modèle, avec une forte applicabilité pour les tâches de classification et de régression. Les paramètres d’hyperparamètres spécifiques pour le modèle XGB dans cette étude sont les suivants : nrounds = 100, max_depth = 2, eta = 0,39, gamma = 0, colsample_bytree = 0,88, seed = 80. La figure supplémentaire 6 représente les matrices de confusion pour XGB.
Les performances des six modèles ont été évaluées à l’aide des paramètres suivants : ensemble de tests ROC AUC, précision, sensibilité, spécificité, valeur prédictive positive (PPV), valeur prédictive négative (NPV) et score F1. La formule de calcul est la suivante22 :

TP : Vrai positif ; TN : Vrai négatif ; FP : Faux positif ; FN : Faux négatif.
6. Recherche de grille et réglage des hyperparamètres
Pour optimiser les performances du modèle, la recherche par grille a été utilisée pour explorer diverses combinaisons d’hyperparamètres. Grâce à un réglage et à une évaluation itératifs, la combinaison optimale d’hyperparamètres a été identifiée afin de maximiser les valeurs AUC du ROC pour le train et les ensembles de test, améliorant ainsi les performances globales du modèle. Les plages des hyperparamètres ajustés sont présentées dans le tableau 3. La recherche par grille est fournie dans la figure supplémentaire 7, la figure supplémentaire 8 et le tableau supplémentaire 1, le tableau supplémentaire 2, le tableau supplémentaire 3 et le tableau supplémentaire 4.
7. Classement de l’importance variable
Les modèles les plus performants, Random Forest (RF) et Extreme Gradient Boosting (XGB), ont été sélectionnés pour le classement d’importance variable. Ces classements aident à identifier les prédicteurs les plus critiques de la rLDH L5-S1 après PEID.
8. Analyse statistique
Toutes les analyses statistiques ont été effectuées à l’aide du logiciel R. Les variables catégorielles ont été déclarées en pourcentage, tandis que les variables continues ont été exprimées en moyenne ± en écart-type. Les caractéristiques cliniques et les paramètres d’imagerie ont été comparés entre le groupe non récurrent (non-rLDH, n = 276) et le groupe récurrent (rLDH, n = 33).