Deze studie is een retrospectieve analyse die is goedgekeurd door de Institutional Ethics Committee van het Nantong First People's Hospital. De studie is geregistreerd op ClinicalTrials.gov (registratienummer: NCT06833099). Omdat de gezondheidsinformatie van alle deelnemers werd geanonimiseerd, was geïnformeerde toestemming niet vereist. De gebruikte verbruiksartikelen en apparatuur staan vermeld in de Materiaaltabel.
1. Studiepopulatie
De studie omvatte klinische gegevens en preoperatieve beeldvormingsdossiers van 436 patiënten die tussen 1 januari 2020 en 30 juni 2024 een percutane endoscopische interlaminaire discectomie (PEID) ondergingen voor lage rugpijn en beenpijn als gevolg van een L5-S1-hernia, in het Nantong First People's Hospital. Op basis van uitsluitingscriteria werden uiteindelijk 309 patiënten geïncludeerd, bij wie recidiverende lumbale discushernia (rLDH) werd bevestigd door middel van postoperatieve Visual Analog Scale (VAS)-scores en follow-up beeldvorming. Van de 309 patiënten ervoeren er 33 rLDH na de operatie, terwijl de overige 276 een significante verlichting van lage rug- en beenpijn hadden, waarbij de VAS-scores met meer dan 60% daalden. Beeldvormingsvariabelen voor alle deelnemers werden afgeleid van preoperatieve röntgen-, CT- en MRI-onderzoeken, gecombineerd met gedetailleerde klinische informatie, waaronder geslacht, leeftijd, lengte, gewicht, BMI, VAS-scores en andere relevante variabelen (Figuur 1).
2. In- en uitsluitingscriteria
Inclusiecriteria voor rLDH: (1) Patiënten met L5-S1 lumbale discushernia die een PEID met één segment hebben ondergaan. (2) Uitgebreide beeldvormende onderzoeken die binnen een maand preoperatief zijn voltooid. (3) Postoperatieve VAS-scoreverlaging ≥60%, gevolgd door een toename van de score en bevestiging door beeldvorming. (4) Geen andere afwijkingen gedetecteerd op beeldvorming. (5) Minimale follow-up periode van 6 maanden.
Inclusiecriteria voor niet-rLDH: (1) Patiënten met L5-S1 lumbale discushernia die een PEID met één segment hebben ondergaan. (2) Uitgebreide beeldvormende onderzoeken die binnen een maand preoperatief zijn voltooid. (3) Postoperatieve VAS-scoreverlaging ≥60%, zonder recidief. (4) Geen andere afwijkingen gedetecteerd op beeldvorming. (5) Minimale follow-up periode van 6 maanden.
Uitsluitingscriteria: (1) Aanwezigheid van andere pathologische aandoeningen die lage rugpijn veroorzaken, zoals schijfinfectie, spinale tumoren, metabole botziekten of osteoporose. (2) Geschiedenis van lumbale schijf of andere spinale operaties. (3) Slechte beeldkwaliteit of onvolledige onderzoeksgegevens. (4) Verlies voor follow-up.
3. Categorische en continue variabelen
Preoperatieve klinische kenmerken en beeldvormingsgegevens van patiënten werden statistisch geanalyseerd en gemeten (Tabel 1 en Tabel 2). Categorische variabelen werden gebruikt om basisziektekenmerken, leefstijlfactoren en andere variabelen te onderscheiden, terwijl continue variabelen specifieke metingen vertegenwoordigden die de fysiologische status van patiënten en preoperatieve beeldvormingsveranderingen beschreven. Om vooroordelen te verminderen, werden strikte kwaliteitscontrolemaatregelen geïmplementeerd. Twee radiologen en wervelkolomchirurgen met meer dan 10 jaar klinische ervaring waren verantwoordelijk voor de statistische analyse en meting van beeldvormingsgegevens. Voor complexe gevallen losten de twee artsen problemen op door middel van gezamenlijk overleg om de nauwkeurigheid en consistentie van de gegevensverwerking te waarborgen.
Categorische variabelen waren onder meer: geslacht, diabetes, hypertensie, cardiovasculaire en cerebrovasculaire aandoeningen (CCD), scoliose, spinale stenose, uitlokkende factoren (bijv. inspannende activiteit, verstuiking, blootstelling aan kou, impact of geen duidelijke trigger), ziekteduur (meer dan 6 maanden of minder dan 6 maanden), gevoelloosheid of zwakte, type uitsteeksel (ingesloten of niet-ingeperkt), schijfdegeneratie (graad I, II, III en IV, V), aangrenzende schijfdegeneratie (graad I, II, III en IV, V), Modische veranderingen en schijfverkalking.
Continue variabelen waren onder meer: leeftijd, operatieduur, lengte, gewicht, body mass index (BMI), maximale diameter van de hernia (MDHD), hoogte van de achterste schijf (PDH), afstand tussen twee wervelcentra (DBTVC), oriëntatiehoek van het facetgewricht (FJOA), schijfhoek (DA), sacrale hellingshoek (SSA), lumbale lordosehoek (LLA) en hoogte-index van de achterste schijf (PDHI, PDHI = PDH / DBTVC).
4. Opschonen van gegevens en selectie van variabelen
Eerst werd het gegevensbestand "zjkj3" ingelezen in R (versie 4.3.1, https://www.r-project.org/, Platform: x86_64-w64-mingw32/x64 (64-bit), Copyright (C) 2023 The R Foundation for Statistical Computing) met behulp van de read_excel functie en opgeslagen in een gegevensvariabele. De geselecteerde tekencodering was UTF-8 (systeemstandaard), een veelgebruikte coderingsmethode die de geschreven talen van de meeste landen over de hele wereld kan ondersteunen. De selectiefunctie werd vervolgens gebruikt om de doelvariabele te scheiden van de functievariabelen. Om de reproduceerbaarheid van de datasplitsing te garanderen, werd een vaste willekeurige seed van 3 ingesteld en werd de dataset opgedeeld in een trainingsset (80%) en een testset (20%). Ten slotte werd een LASSO-regressiefunctie gedefinieerd en werd een deel van de trainingsgegevens willekeurig bemonsterd met behulp van de glmnet-functie om in de LASSO-regressie te passen. De optimale regularisatieparameter λ werd bepaald met behulp van L1-regularisatie en kruisvalidatie (standaard 10-voudig), terwijl ervoor werd gezorgd dat er geen dubbele patiëntgegevens werden gemengd tussen de trainings- en testsets. De kruisvalidatiefout voor elke λ-waarde werd berekend om de optimale λ-waarde te bepalen.
5. Ontwikkeling en evaluatie van modellen
SVM (Support Vector Machine): Dit model wijst trainingsgegevens toe aan een hypervlak, waarbij de marge tussen twee klassen wordt gemaximaliseerd om het doel te voorspellen. De specifieke hyperparameterinstellingen voor het SVM-model in dit onderzoek zijn: kernel = "lineair", kosten = 1. Aanvullende figuur 1 geeft de verwarringsmatrices voor SVM weer.
DT (beslissingsboom): Dit model splitst gegevens recursief om een boomstructuur te creëren voor het voorspellen van de doelvariabele. De specifieke hyperparameterinstelling voor het DT-model in dit onderzoek is: max_depth = 3. Aanvullende figuur 2 toont de verwarringsmatrices voor DT.
ADA (AdaBoost): Dit model combineert meerdere zwakke studenten (meestal beslissingsbomen) om een sterke student te creëren, waardoor de classificatieprestaties worden verbeterd. De specifieke hyperparameterinstellingen voor het ADA-model in dit onderzoek zijn: n_estimators = 150, learning_rate = 0,1, seed = 80. Aanvullende figuur 3 toont de verwarringsmatrices voor ADA.
LGBM (Light Gradient Boosting Machine): Dit model maakt gebruik van een histogramalgoritme om continue kenmerken te splitsen, het trainingsproces te versnellen en de geheugenbelasting te verminderen, waardoor het een efficiënt en snel gradiëntversterkend raamwerk is, bijzonder geschikt voor grootschalige datasets. De specifieke hyperparameterinstellingen voor het LGBM-model in dit onderzoek zijn: num_leaves = 5, learning_rate = 0,05, n_estimators = 50. Aanvullende figuur 4 toont de verwarringsmatrices voor LGBM.
RF (Random Forest): Dit model construeert meerdere beslissingsbomen en combineert hun voorspellingen om de nauwkeurigheid en overfitting te verbeteren. De specifieke hyperparameterinstellingen voor het RF-model in dit onderzoek zijn: ntree = 310, mtry = 1, maxnodes = 10, max_depth = 1, seed = 80. Aanvullende figuur 5 toont de verwarringsmatrices voor RF.
XGB (Extreme Gradient Boosting): Dit model is een verbeterd algoritme op basis van beslissingsbomen, waarbij gebruik wordt gemaakt van een nieuw gegeneraliseerd algoritme voor het stimuleren van gradiënten om de modelconstructie te versnellen, met een sterke toepasbaarheid voor classificatie- en regressietaken. De specifieke hyperparameterinstellingen voor het XGB-model in dit onderzoek zijn: nrounds = 100, max_depth = 2, eta = 0,39, gamma = 0, colsample_bytree = 0,88, seed = 80. Aanvullende figuur 6 toont de verwarringsmatrices voor XGB.
De prestaties van de zes modellen werden geëvalueerd aan de hand van de volgende statistieken: ROC AUC van de testset, nauwkeurigheid, gevoeligheid, specificiteit, positief voorspellende waarde (PPV), negatief voorspellende waarde (NPV) en F1-score. De berekeningsformule is als volgt22:

TP: Waar positief; TN: Waar Negatief; FP: Vals positief; FN: Vals negatief.
6. Raster zoeken en hyperparameterafstemming
Om de modelprestaties te optimaliseren, werd rasterzoeken gebruikt om verschillende combinaties van hyperparameters te verkennen. Door middel van iteratieve afstemming en evaluatie werd de optimale hyperparametercombinatie geïdentificeerd om de ROC AUC-waarden voor zowel de trein- als de testsets te maximaliseren, waardoor de algehele modelprestaties werden verbeterd. De bereiken van aangepaste hyperparameters worden weergegeven in tabel 3. De rasterzoekopdracht wordt gegeven in Aanvullende Figuur 7, Aanvullende Figuur 8 en Aanvullende Tabel 1, Aanvullende Tabel 2, Aanvullende Tabel 3 en Aanvullende Tabel 4.
7. Rangschikking van variabel belang
De best presterende modellen - Random Forest (RF) en Extreme Gradient Boosting (XGB) - werden geselecteerd voor een rangschikking van variabel belang. Deze ranglijsten helpen bij het identificeren van de meest kritische voorspellers van L5-S1 rLDH na PEID.
8. Statistische analyse
Alle statistische analyses zijn uitgevoerd met behulp van R-software. Categorische variabelen werden gerapporteerd als percentages, terwijl continue variabelen werden uitgedrukt als gemiddelde ± standaarddeviatie. Klinische kenmerken en beeldvormingsparameters werden vergeleken tussen de niet-recidiverende groep (niet-rLDH, n=276) en de recidiverende groep (rLDH, n=33).