Questo studio è un'analisi retrospettiva approvata dal Comitato Etico Istituzionale del Nantong First People's Hospital. La sperimentazione è stata registrata il ClinicalTrials.gov (numero di registrazione: NCT06833099). Poiché le informazioni sanitarie di tutti i partecipanti sono state rese anonime, non è stato richiesto il consenso informato. I materiali di consumo e le attrezzature utilizzate sono elencati nella Tabella dei materiali.
1. Popolazione oggetto di studio
Lo studio ha incluso dati clinici e cartelle cliniche di imaging preoperatorio di 436 pazienti sottoposti a discectomia interlaminare endoscopica percutanea (PEID) per lombalgia e dolore alle gambe a causa di ernia del disco L5-S1 presso il Nantong First People's Hospital tra il 1° gennaio 2020 e il 30 giugno 2024. Sulla base dei criteri di esclusione, sono stati infine inclusi 309 pazienti, tra i quali l'ernia del disco lombare ricorrente (rLDH) è stata confermata attraverso i punteggi della Visual Analog Scale (VAS) postoperatoria e l'imaging di follow-up. Dei 309 pazienti, 33 hanno sperimentato rLDH post-operatorio, mentre i restanti 276 hanno avuto un significativo sollievo dal dolore lombare e alle gambe, con punteggi VAS ridotti di oltre il 60%. Le variabili di imaging per tutti i partecipanti sono state derivate da esami radiografici, TC e risonanza magnetica preoperatori, combinati con informazioni cliniche dettagliate, tra cui sesso, età, altezza, peso, BMI, punteggi VAS e altre variabili rilevanti (Figura 1).
2. Criteri di inclusione ed esclusione
Criteri di inclusione per rLDH: (1) Pazienti con ernia del disco lombare L5-S1 sottoposti a PEID a segmento singolo. (2) Esami di imaging completi completati entro un mese prima dell'intervento. (3) Riduzione del punteggio VAS postoperatorio del ≥60%, seguita da un aumento del punteggio e dalla conferma mediante imaging. (4) Non sono state rilevate altre anomalie all'imaging. (5) Periodo minimo di follow-up di 6 mesi.
Criteri di inclusione per non-rLDH: (1) Pazienti con ernia del disco lombare L5-S1 sottoposti a PEID a segmento singolo. (2) Esami di imaging completi completati entro un mese prima dell'intervento. (3) Riduzione del punteggio VAS postoperatorio del ≥60%, senza recidive. (4) Non sono state rilevate altre anomalie all'imaging. (5) Periodo minimo di follow-up di 6 mesi.
Criteri di esclusione: (1) Presenza di altre condizioni patologiche che causano lombalgia, come infezione del disco, tumori spinali, malattie metaboliche ossee o osteoporosi. (2) Storia di disco lombare o altri interventi chirurgici spinali. (3) Scarsa qualità dell'imaging o dati dell'esame incompleti. (4) Perdita al follow-up.
3. Variabili categoriali e continue
Le caratteristiche cliniche preoperatorie e i dati di imaging dei pazienti sono stati analizzati e misurati statisticamente (Tabella 1 e Tabella 2). Le variabili categoriali sono state utilizzate per distinguere le caratteristiche di base della malattia, i fattori dello stile di vita e altre variabili, mentre le variabili continue rappresentavano misurazioni specifiche che descrivevano lo stato fisiologico dei pazienti e i cambiamenti dell'imaging preoperatorio. Per ridurre le distorsioni, sono state implementate rigorose misure di controllo della qualità. Due radiologi e chirurghi della colonna vertebrale con oltre 10 anni di esperienza clinica sono stati responsabili dell'analisi statistica e della misurazione dei dati di imaging. Per i casi complessi, i due medici hanno risolto i problemi attraverso una consultazione congiunta per garantire l'accuratezza e la coerenza nell'elaborazione dei dati.
Le variabili categoriali includevano: sesso, diabete, ipertensione, malattie cardiovascolari e cerebrovascolari (CCD), scoliosi, stenosi spinale, fattori scatenanti (ad esempio, attività intensa, distorsione, esposizione al freddo, impatto o assenza di un fattore scatenante chiaro), durata della malattia (superiore a 6 mesi o inferiore a 6 mesi), intorpidimento o debolezza, tipo di protrusione (contenuta o non contenuta), degenerazione del disco (Gradi I, II, III e IV, V), degenerazione del disco adiacente (Gradi I, II, III e IV, V), alterazioni modiche e calcificazioni del disco.
Le variabili continue includevano: età, durata dell'intervento chirurgico, altezza, peso, indice di massa corporea (BMI), diametro massimo dell'ernia del disco (MDHD), altezza del disco posteriore (PDH), distanza tra due centri delle vertebre (DBTVC), angolo di orientamento delle faccette articolari (FJOA), angolo del disco (DA), angolo di pendenza sacrale (SSA), angolo di lordosi lombare (LLA) e indice di altezza del disco posteriore (PDHI, PDHI = PDH / DBTVC).
4. Pulizia dei dati e selezione delle variabili
Innanzitutto, il file di dati "zjkj3" è stato letto in R (versione 4.3.1, https://www.r-project.org/, Piattaforma: x86_64-w64-mingw32/x64 (64-bit), Copyright (C) 2023 The R Foundation for Statistical Computing) utilizzando la funzione read_excel e memorizzato in una variabile di dati. La codifica dei caratteri selezionata è stata UTF-8 (impostazione predefinita del sistema), che è un metodo di codifica ampiamente utilizzato in grado di supportare le lingue scritte della maggior parte dei paesi del mondo. La funzione di selezione è stata quindi utilizzata per separare la variabile target dalle variabili di funzionalità. Per garantire la riproducibilità della suddivisione dei dati, è stato impostato un seme casuale fisso di 3 e il set di dati è stato diviso in un set di addestramento (80%) e un set di test (20%). Infine, è stata definita una funzione di regressione LASSO e una parte dei dati di addestramento è stata campionata in modo casuale utilizzando la funzione glmnet per adattarsi alla regressione LAZO. Il parametro di regolarizzazione ottimale λ è stato determinato utilizzando la regolarizzazione L1 e la convalida incrociata (impostazione predefinita 10 volte), garantendo al contempo che non venissero mescolati dati duplicati del paziente tra il set di addestramento e quello di test. L'errore di convalida incrociata per ciascun valore λ è stato calcolato per determinare il valore λ ottimale.
5. Sviluppo e valutazione del modello
SVM (Support Vector Machine): questo modello mappa i dati di addestramento su un iperpiano, massimizzando il margine tra due classi per prevedere l'obiettivo. Le impostazioni specifiche degli iperparametri per il modello SVM in questo studio sono: kernel = "linear", cost = 1. La Figura 1 supplementare illustra le matrici di confusione per SVM.
DT (Decision Tree): questo modello suddivide in modo ricorsivo i dati per creare una struttura ad albero per la previsione della variabile target. L'impostazione specifica dell'iperparametro per il modello DT in questo studio è: max_depth = 3. La Figura 2 supplementare illustra le matrici di confusione per DT.
ADA (AdaBoost): questo modello combina più studenti deboli (in genere alberi decisionali) per creare uno studente forte, migliorando le prestazioni di classificazione. Le impostazioni specifiche degli iperparametri per il modello ADA in questo studio sono: n_estimators = 150, learning_rate = 0,1, seme = 80. La Figura 3 supplementare illustra le matrici di confusione per ADA.
LGBM (Light Gradient Boosting Machine): questo modello utilizza un algoritmo di istogramma per suddividere le funzionalità continue, accelerando il processo di addestramento e riducendo il carico di memoria, rendendolo un framework di potenziamento del gradiente efficiente e veloce, particolarmente adatto per set di dati su larga scala. Le impostazioni specifiche degli iperparametri per il modello LGBM in questo studio sono: num_leaves = 5, learning_rate = 0,05, n_estimators = 50. La Figura 4 supplementare illustra le matrici di confusione per LGBM.
RF (Random Forest): questo modello costruisce più alberi decisionali e combina le loro previsioni per migliorare la precisione e controllare l'overfitting. Le impostazioni specifiche degli iperparametri per il modello RF in questo studio sono: ntree = 310, mtry = 1, maxnodes = 10, max_depth = 1, seed = 80. La Figura 5 supplementare illustra le matrici di confusione per RF.
XGB (Extreme Gradient Boosting): questo modello è un algoritmo avanzato basato su alberi decisionali, che utilizza un nuovo algoritmo di albero decisionale generalizzato per accelerare la costruzione del modello, con una forte applicabilità per le attività di classificazione e regressione. Le impostazioni specifiche degli iperparametri per il modello XGB in questo studio sono: nrounds = 100, max_depth = 2, eta = 0,39, gamma = 0, colsample_bytree = 0,88, seed = 80. La Figura 6 supplementare illustra le matrici di confusione per XGB.
Le prestazioni dei sei modelli sono state valutate utilizzando le seguenti metriche: AUC ROC del set di test, accuratezza, sensibilità, specificità, valore predittivo positivo (PPV), valore predittivo negativo (NPV) e punteggio F1. La formula di calcolo è la seguente22:

TP: Vero positivo; TN: Vero negativo; FP: Falso positivo; FN: Falso negativo.
6. Ricerca della griglia e regolazione degli iperparametri
Per ottimizzare le prestazioni del modello, è stata utilizzata la ricerca a griglia per esplorare varie combinazioni di iperparametri. Attraverso l'ottimizzazione e la valutazione iterativa, è stata identificata la combinazione ottimale di iperparametri per massimizzare i valori di AUC ROC sia per il treno che per i set di test, migliorando così le prestazioni complessive del modello. Gli intervalli di iperparametri regolati sono presentati nella Tabella 3. La ricerca nella griglia viene fornita nella Figura 7 supplementare, nella Figura 8 e nella Tabella 1 supplementare, nella Tabella 2 supplementare, nella Tabella 3 e nella Tabella 4 supplementare.
7. Classificazione per importanza variabile
I modelli più performanti - Random Forest (RF) e Extreme Gradient Boosting (XGB) - sono stati selezionati per la classificazione di importanza variabile. Queste classificazioni aiutano a identificare i predittori più critici di L5-S1 rLDH dopo PEID.
8. Analisi statistica
Tutte le analisi statistiche sono state eseguite utilizzando il software R. Le variabili categoriali sono state riportate in percentuale, mentre le variabili continue sono state espresse come media ± deviazione standard. Le caratteristiche cliniche e i parametri di imaging sono stati confrontati tra il gruppo non ricorrente (non-rLDH, n=276) e il gruppo ricorrente (rLDH, n=33).