Niniejsze badanie jest analizą retrospektywną zatwierdzoną przez Komisję Etyki Instytucjonalnej szpitala Nantong First People's Hospital. Badanie zostało zarejestrowane w serwisie ClinicalTrials.gov (numer rejestracyjny: NCT06833099). Ponieważ dane zdrowotne wszystkich uczestników zostały zanonimizowane, zgoda na udział w badaniu nie była wymagana. Wykorzystane materiały eksploatacyjne i sprzęt wymieniono w Tabeli Materiałów.
1. Populacja badawcza
Badanie objęło dane kliniczne i zapisy obrazowania przedoperacyjnego 436 pacjentów, którzy zostali poddani przezskórnej endoskopowej dyskektomii międzyblaszawej (PEID) z powodu bólu w odcinku lędźwiowym kręgosłupa i bólu kończyn dolnych wynikających z przepukliny tarczy międzykręgowej L5-S1 w Nantong First People's Hospital w okresie od 1 stycznia 2020 r. do 30 czerwca 2024 r. Na podstawie kryteriów wykluczenia ostatecznie zakwalifikowano 309 pacjentów, u których nawrotową przepuklinę krążka międzykręgowego odcinka lędźwiowego (rLDH) potwierdzono za pomocą pooperacyjnych wyników w wizualnej skali analogowej (VAS) oraz obrazowania kontrolnego. Spośród 309 pacjentów u 33 wystąpiła rLDH po operacji, natomiast u pozostałych 276 odnotowano znaczną redukcję bólu pleców i nóg, z obniżeniem wyników VAS o ponad 60%. Zmienne obrazowe dla wszystkich uczestników uzyskano z przedoperacyjnych badań RTG, TK i MRI, łącząc je ze szczegółowymi informacjami klinicznymi, w tym płcią, wiekiem, wzrostem, masą ciała, BMI, wynikami VAS oraz innymi istotnymi zmiennymi (Rysunek 1).
2. Kryteria włączenia i wykluczenia
Kryteria włączenia dla rLDH: (1) Pacjenci z wypukliną krążka międzykręgowego na poziomie L5-S1, którzy zostali poddani jednosegmentowej procedurze PEID. (2) Kompleksowe badania obrazowe wykonane w ciągu jednego miesiąca przed operacją. (3) Pooperacyjna redukcja wyniku VAS ≥60%, po której nastąpił wzrost wyniku i potwierdzenie w badaniu obrazowym. (4) Brak innych nieprawidłowości wykrytych w badaniach obrazowych. (5) Minimalny okres obserwacji wynoszący 6 miesięcy.
Kryteria włączenia dla grupy Non-rLDH: (1) Pacjenci z przepukliną krążka międzykręgowego L5-S1, poddani jednosegmentowej procedurze PEID. (2) Kompleksowe badania obrazowe wykonane w ciągu jednego miesiąca przed operacją. (3) Redukcja wyniku w skali VAS po operacji ≥60%, bez nawrotów. (4) Brak innych nieprawidłowości wykrytych w badaniach obrazowych. (5) Minimalny okres obserwacji wynoszący 6 miesięcy.
Kryteria wykluczenia: (1) Obecność innych stanów patologicznych powodujących ból w odcinku lędźwiowym kręgosłupa, takich jak infekcje krążka międzykręgowego, guzy kręgosłupa, metaboliczne choroby kostne lub osteoporoza. (2) Historia operacji krążków międzykręgowych odcinka lędźwiowego lub innych operacji kręgosłupa. (3) Słaba jakość obrazowania lub niepełne dane badania. (4) Utrata kontaktu z pacjentem podczas obserwacji.
3. Zmienne kategoryczne i ciągłe
Przeprowadzono analizę statystyczną i pomiary przedoperacyjnych cech klinicznych oraz danych obrazowych pacjentów (Tabela 1 i Tabela 2). Zmienne kategoryjne posłużyły do rozróżnienia podstawowych cech choroby, czynników związanych ze stylem życia i innych zmiennych, natomiast zmienne ciągłe reprezentowały konkretne pomiary opisujące stan fizjologiczny pacjentów oraz przedoperacyjne zmiany w obrazowaniu. W celu zminimalizowania błędu wdrożono rygorystyczne środki kontroli jakości. Za analizę statystyczną i pomiary danych obrazowych odpowiadali dwaj radiolodzy i chirurdzy kręgosłupa z ponad 10-letnim doświadczeniem klinicznym. W przypadku złożonych przypadków obaj lekarze rozwiązywali problemy poprzez wspólną konsultację, aby zapewnić dokładność i spójność przetwarzania danych.
Zmienne kategoryczne obejmowały: płeć, cukrzycę, nadciśnienie tętnicze, choroby sercowo-naczyniowe i mózgowo-naczyniowe (CCD), skoliozę, stenozę kanału kręgowego, czynniki wyzwalające (np. intensywna aktywność, skręcenie, ekspozycja na zimno, uderzenie lub brak wyraźnego czynnika), czas trwania choroby (powyżej 6 miesięcy lub poniżej 6 miesięcy), drętwienie lub osłabienie, typ wypukliny (ograniczona lub nieograniczona), degenerację krążka międzykręgowego (stopnie I, II, III oraz IV, V), degenerację sąsiedniego krążka międzykręgowego (stopnie I, II, III oraz IV, V), zmiany Modica oraz zwapnienie krążka międzykręgowego.
Zmienne ciągłe obejmowały: wiek, czas trwania operacji, wzrost, masę ciała, wskaźnik masy ciała (BMI), maksymalną średnicę przepukliny tarczy międzykręgowej (MDHD), tylną wysokość krążka (PDH), odległość między środkami dwóch kręgów (DBTVC), kąt orientacji stawów międzywyrostkowych (FJOA), kąt krążka (DA), kąt nachylenia kości krzyżowej (SSA), kąt lordozy lędźwiowej (LLA) oraz wskaźnik tylnej wysokości krążka (PDHI, PDHI = PDH / DBTVC).
4. Czyszczenie danych i wybór zmiennych
W pierwszej kolejności plik danych „zjkj3” został wczytany do programu R (wersja 4.3.1, https://www.r-project.org/, platforma: x86_64-w64-mingw32/x64 (64-bit), Copyright (C) 2023 The R Foundation for Statistical Computing) za pomocą funkcji read_excel i zapisany w zmiennej data. Wybrano kodowanie znaków UTF-8 (domyślne dla systemu), które jest powszechnie stosowaną metodą kodowania obsługującą języki pisane większości krajów na świecie. Następnie użyto funkcji select, aby oddzielić zmienną docelową od zmiennych cech. Aby zapewnić powtarzalność podziału danych, ustawiono stałe ziarno losowości (seed) wynoszące 3, a zbiór danych podzielono na zbiór treningowy (80%) i zbiór testowy (20%). Na koniec zdefiniowano funkcję regresji LASSO, a część danych treningowych została wylosowana za pomocą funkcji glmnet w celu dopasowania modelu regresji LASSO. Optymalny parametr regularyzacji λ określono przy użyciu regularyzacji L1 i walidacji krzyżowej (domyślnie 10-krotnej), dbając jednocześnie o to, aby dane tych samych pacjentów nie pojawiły się jednocześnie w zbiorze treningowym i testowym. Aby wyznaczyć optymalną wartość λ, obliczono błąd walidacji krzyżowej dla każdej wartości λ.
5. Opracowanie i ocena modelu
SVM (Support Vector Machine): Model ten odwzorowuje dane treningowe na nadpłaszczyznę, maksymalizując margines między dwiema klasami w celu przewidzenia celu. Konkretne ustawienia hiperparametrów dla modelu SVM w niniejszym badaniu to: kernel = "linear", cost = 1. Rycina uzupełniająca 1 przedstawia macierze pomyłek dla SVM.
DT (drzewo decyzyjne): Ten model rekurencyjnie dzieli dane, tworząc strukturę drzewa w celu przewidzenia zmiennej docelowej. Konkretne ustawienie hiperparametrów dla modelu DT w niniejszym badaniu to: max_depth = 3. Rycina uzupełniająca 2 przedstawia macierze pomyłek dla DT.
ADA (AdaBoost): Model ten łączy wiele słabych klasyfikatorów (zazwyczaj drzewa decyzyjne) w celu stworzenia silnego klasyfikatora, co poprawia wydajność klasyfikacji. Konkretne ustawienia hiperparametrów dla modelu ADA w niniejszym badaniu to: n_estimators = 150, learning_rate = 0.1, seed = 80. Rysunek uzupełniający 3 przedstawia macierze pomyłek dla ADA.
LGBM (Light Gradient Boosting Machine): Model ten wykorzystuje algorytm histogramowy do podziału cech ciągłych, co przyspiesza proces uczenia i zmniejsza obciążenie pamięci, czyniąc go wydajnym i szybkim frameworkiem gradient boosting, szczególnie odpowiednim dla dużych zbiorów danych. Specyficzne ustawienia hiperparametrów dla modelu LGBM w niniejszym badaniu to: num_leaves = 5, learning_rate = 0.05, n_estimators = 50. Rysunek uzupełniający 4 przedstawia macierze pomyłek dla modelu LGBM.
RF (Random Forest): Model ten konstruuje wiele drzew decyzyjnych i łączy ich przewidywania, aby zwiększyć dokładność i zapobiec przeuczeniu. Konkretne ustawienia hiperparametrów dla modelu RF w niniejszym badaniu to: ntree = 310, mtry = 1, maxnodes = 10, max_depth = 1, seed = 80. Rysunek dodatkowy 5 przedstawia macierze pomyłek dla modelu RF.
XGB (Extreme Gradient Boosting): Model ten jest ulepszonym algorytmem opartym na drzewach decyzyjnych, wykorzystującym nowy, uogólniony algorytm budowania drzew decyzyjnych metodą wzmacniania gradientowego w celu przyspieszenia konstrukcji modelu, wykazujący wysoką przydatność w zadaniach klasyfikacji i regresji. Konkretne ustawienia hiperparametrów dla modelu XGB w niniejszym badaniu to: nrounds = 100, max_depth = 2, eta = 0.39, gamma = 0, colsample_bytree = 0.88, seed = 80. Rycina uzupełniająca 6 przedstawia macierze pomyłek dla modelu XGB.
Wydajność sześciu modeli oceniono za pomocą następujących miar: ROC AUC zbioru testowego, dokładności, czułości, swoistości, dodatniej wartości predykcyjnej (PPV), ujemnej wartości predykcyjnej (NPV) oraz wyniku F1. Wzór obliczeniowy jest następujący22:

TP: prawdziwie dodatni; TN: prawdziwie ujemny; FP: fałszywie dodatni; FN: fałszywie ujemny.
6. Przeszukiwanie siatki i strojenie hiperparametrów
Aby zoptymalizować wydajność modelu, zastosowano metodę grid search w celu zbadania różnych kombinacji hiperparametrów. Poprzez iteracyjne dostrajanie i ewaluację zidentyfikowano optymalną kombinację hiperparametrów, która maksymalizuje wartości ROC AUC zarówno dla zbioru treningowego, jak i testowego, poprawiając tym samym ogólną wydajność modelu. Zakresy dostrojonych hiperparametrów przedstawiono w Tabeli 3. Wyniki grid search znajdują się w Ryc. uzupełniającej 7, Ryc. uzupełniającej 8, Tabeli uzupełniającej 1, Tabeli uzupełniającej 2, Tabeli uzupełniającej 3 oraz Tabeli uzupełniającej 4.
7. Ranking ważności zmiennych
Do rankowania istotności zmiennych wybrano modele o najwyższej wydajności – Random Forest (RF) oraz Extreme Gradient Boosting (XGB). Rankingi te pomagają w identyfikacji krytycznych predyktorów rLDH na poziomie L5-S1 po PEID.
8. Analiza statystyczna
Wszystkie analizy statystyczne przeprowadzono przy użyciu oprogramowania R. Zmienne kategoryczne przedstawiono jako wartości procentowe, natomiast zmienne ciągłe wyrażono jako średnia ± odchylenie standardowe. Charakterystykę kliniczną i parametry obrazowania porównano pomiędzy grupą bez nawrotów (non-rLDH, n=276) a grupą z nawrotami (rLDH, n=33).