Badanie zostało przeprowadzone zgodnie z Deklaracją Helsińską i zostało zatwierdzone przez instytucjonalną komisję etyki Jinhua Municipal Central Hospital (nr 2022129). Wymóg uzyskania pisemnej świadomej zgody został zniesiony, ponieważ analizowano wyłącznie retrospektywne, zanonimizowane dane. Narzędzia badawcze wykorzystane w niniejszym protokole są wymienione w Tabeli materiałów.
1. Plan badania i populacja źródłowa
Retrospektywne obserwacyjne badanie kohortowe przeprowadzono w kompleksowym centrum udarowym w Jinhua Municipal Central Hospital w prowincji Zhejiang w Chinach. Przesiewowo przeanalizowano kolejne przyjęcia pacjentów od czerwca 2024 do września 2025 roku. Do badania włączono pacjentów, którzy zostali przyjęci w ciągu 24 h od wystąpienia objawów, mieli wyjściową wartość skali NIHSS równą 5 lub niższą oraz otrzymali diagnozę zespołu leczniczego dotyczącą łagodnego udaru niedokrwiennego, zgodną z procesem w obrębie małych naczyń podkorowych lub naczyń perforujących. Pacjentów wykluczono w przypadku innych udokumentowanych podtypów etiologicznych, wyjściowej wartości NIHSS powyżej 5, czasu od wystąpienia objawów do przyjęcia przekraczającego 24 h, ciężkiej istniejącej wcześniej dysfunkcji neurologicznej lub ciężkiej choroby ogólnoustrojowej, a także niekompletnych podstawowych danych klinicznych lub obrazowych. Każde wykluczenie odnotowano w schemacie blokowym.
2. Klasyfikacja etiologiczna i obrazowanie mózgu
Zarchiwizowano klasyfikację etiologiczną przypisaną podczas rutynowej opieki przez zespół leczący udar. Ponieważ było to retrospektywne badanie oparte na elektronicznej dokumentacji medycznej, przypisanie to nie było uznane za prospektywne, scentralizowane, zaślepione ani niezależnie zweryfikowane. W przypadku przypadków potwierdzonych za pomocą MRI wymagana była obecność ostrej zmiany podkorowej zgodnej z zespołem neurologicznym w obrazowaniu zależnym od dyfuzji, z odpowiadającą nieprawidłowością współczynnika dyfuzji pozornej oraz brakiem zawału korowego. W przypadku przypadków opartych wyłącznie na CT wymagane było wykluczenie krwotoku i terytorialnego zawału korowego, wraz z wystąpieniem zgodnego klinicznego zespołu lakunarnego. Przyznano, że samo CT może być potencjalnie niewystarczające do bezpośredniego potwierdzenia małego ostrego zawału podkorowego.
Zapisano następujące parametry: metodę obrazowania, czas od wystąpienia objawów do wykonania badania, maksymalną średnicę zmiany, grubość warstwy, liczbę objętych warstw osiowych, zajęcie brzusznej powierzchni mostu, stenozę tętnicy doprowadzającej oraz stenozę w naczyniach nieodpowiedzialnych za zmianę. Podejrzenie BAD zdefiniowano jako wystąpienie któregokolwiek z poniższych kryteriów: maksymalną średnicę zmiany wynoszącą co najmniej 15 mm, zajęcie co najmniej trzech warstw osiowych lub zajęcie brzusznej powierzchni mostu. Szczegółowa strona zmiany, lokalizacja anatomiczna, segment tętnicy odpowiedzialnej za zmianę oraz możliwość niezależnej, zaślepionej ponownej analizy obrazów nie były dostępne i nie zostały poddane imputacji.
3. Ocena angiograficzna zwężenia
Pomiary zwężenia w CTA lub MRA zostały pobrane z klinicznych raportów obrazowych. Procentowe zwężenie wewnątrzczaszkowe obliczono przy użyciu wzoru WASID:
[1 - (średnica w najwęższym miejscu światła
/średnica sąsiedniej normalnej tętnicy referencyjnej)] x 10015
Stenozę tętnicy doprowadzającej zdefiniowano jako zwężenie wynoszące co najmniej 50% w tętnicy zaopatrującej obszar zawału. Stenozę nieodpowiedzialną za zawał zdefiniowano jako zwężenie wynoszące co najmniej 50% poza tym obszarem. Główną zmienną modelu, jaką była jakakolwiek stenoza naczyniowa wynosząca co najmniej 50%, uznano za dodatnią, jeśli którykolwiek z powyższych komponentów był dodatni. Ze względu na to, że odpowiedzialny odcinek tętnicy nie został zachowany, pomiarów specyficznych dla danego odcinka nie można było zweryfikować niezależnie.
4. Wyjściowe zmienne kliniczne, laboratoryjne i terapeutyczne
Zebrano dane dotyczące wieku, płci, wywiadu naczyniowego, palenia tytoniu, spożycia alkoholu, stosowania leków przeciwpłytkowych przed udarem, wyniku w zmodyfikowanej skali Rankin (mRS) przed udarem, czasu od wystąpienia objawów do przyjęcia do szpitala, skurczowego ciśnienia tętniczego (SBP) i rozkurczowego ciśnienia tętniczego (DBP) przy przyjęciu, wyjściowej punktacji NIHSS, liczby płytek krwi, poziomu glukozy we krwi, liczby neutrofili, liczby limfocytów oraz czasu pobrania krwi. NLR obliczono jako stosunek bezwzględnej liczby neutrofili do bezwzględnej liczby limfocytów.
Zgłoszono następujące parametry: trombolizę dożylną, strategię pojedynczego versus podwójnego leczenia przeciwpłytkowego, czas rozpoczęcia leczenia po przyjęciu do szpitala, intensywność stosowania statyn, szpitalne leczenie obniżające ciśnienie tętnicze, kontrolę poziomu glukozy, przeciwkrzepliwość oraz odnotowany czas ekspozycji na leki przeciwtrombotyczne w stosunku do momentu pobrania krwi. Żaden z uczestników nie otrzymał przeciwkrzepliwości w warunkach szpitalnych. Objętość płynów podanych dożylnie nie została uwzględniona, ponieważ nie można było jej wiarygodnie odtworzyć. Interwencje te traktowano jako potencjalne czynniki zakłócające lub następstwa stanu klinicznego, a nie jako zrandomizowane ekspozycje.
5. Oceny neurologiczne i END
Wykorzystano wyniki w skali NIHSS udokumentowane przy przyjęciu oraz w przybliżeniu po 24 i 72 h, a także uwzględniono dodatkowe oceny wywołane pogorszeniem stanu klinicznego. Wyniki były rejestrowane przez rutynowy zespół ds. udarów. Nie udokumentowano zaś maskowania predyktorów ani scentralizowanego podwójnego punktowania. Główny wskaźnik END zdefiniowano jako wzrost całkowitej liczby punktów w skali NIHSS o co najmniej 2 punkty w ciągu 72 h od przyjęcia w stosunku do wartości wyjściowej, zgodnie ze współczesną literaturą dotyczącą udarów lakunarnych2,3. Odnotowano czas oraz udokumentowaną przyczynę pogorszenia. W głównym punkcie końcowym uwzględniono kwalifikujące się pogorszenia o każdej przyczynie, a model powtórzono, uwzględniając wyłącznie END o przyczynie niedokrwiennej.
6. Obsługa brakujących danych
Przed analizą przeprowadzono audyt brakujących danych. W przypadku trzech rekordów z brakującą surową wartością NLR, ale dostępnymi liczbami neutrofili i limfocytów, wartość NLR została obliczona deterministycznie. W przypadku czterech rekordów z brakującym surowym SD SBP z 24 h, ale dostępną średnią SBP i współczynnikiem zmienności, SD obliczono jako współczynnik zmienności × średnia SBP. Po tych deterministycznych wyprowadzeniach dane dla czterech głównych predyktorów oraz wynik były kompletne, a wszystkich 240 pacjentów włączono do modelu głównego bez zastosowania imputacji wielokrotnej. Pojedynczy rekord z brakującym statusem przeciwprzekrzepowym przed pobraniem krwi został wykluczony wyłącznie z eksploracyjnego modelu skorygowanego o leczenie.
7. Analiza statystyczna i opracowanie modelu
Zmienne ciągłe o rozkładzie zbliżonym do symetrycznego przedstawiono jako średnia ± SD, natomiast zmienne o rozkładzie skośnym lub zmienne porządkowe jako mediana (rozstęp międzykwartylowy). W zależności od potrzeb stosowano testy t Welcha lub testy U Manna–Whitneya. Zmienne kategoryczne przedstawiono jako n (%) i porównano za pomocą testu chi-kwadrat lub dokładnego testu Fishera.
Podstawowy model regresji logistycznej obejmował ciągłe SBP przy przyjęciu przeskalowane na każde 10 mmHg, wyjściową punktację NIHSS na każdy punkt, NLR na jednostkę oraz jakąkolwiek stenozę naczyniową wynoszącą co najmniej 50%. Zgłoszono wyraz wolny, współczynniki, błędy standardowe, OR, 95% CI oraz dwustronne wartości P. Nieliniowość SBP zbadano przy użyciu ograniczonego splotu sześciennego (restricted cubic spline) z trzema węzłami na 10., 50. i 90. percentylu (130,9, 157,0 i 194,0 mmHg), po skorygowaniu o NIHSS, NLR i stenozę. Do oceny składnika nieliniowego zastosowano test ilorazu wiarygodności. W przypadku braku potwierdzenia nieliniowości zastosowano liniowy składnik SBP.
Próg SBP wynoszący 169 mmHg został zachowany jedynie w ramach wtórnej uproszczonej analizy. Próg ten wybrano przy użyciu maksymalnego indeksu Youdena w danych rozwojowych, a niepewność punktu odcięcia ilościowo określono za pomocą 1 000 prób bootstrapowych. Nie przywrócono niezweryfikowanej skali 0–12 oraz kategorii niskiego, pośredniego i wysokiego ryzyka.
Zdolność dyskryminacji oszacowano przy użyciu AUC oraz 95% przedziału ufności (CI) bootstrapowego. Model łączony porównano z każdym z komponentów za pomocą sparowanych różnic bootstrapowych; nie stwierdzano wyższości, gdy przedział ufności różnicy obejmował zero. Optymizm oszacowano przy użyciu 1000 bootstrapowych prób ponownych. Przewidywania bootstrapowe out-of-bag wygenerowano w celu kalibracji i obliczenia wskaźnika Briera. Analiza krzywej decyzyjnej została przeprowadzona wyłącznie dla prawdopodobieństw progowych od 5% do 30%, dla których modelowanym działaniem była wzmożona obserwacja neurologiczna i ponowna ocena kliniczna, a nie automatyczne leczenie.
Analizy wrażliwości przeprowadzono w kohorcie z samym MRI, po wykluczeniu zwężeń tętnicy doprowadzającej o co najmniej 50%, po wykluczeniu podejrzewanego BAD oraz przy uwzględnieniu wyłącznie END o przyczynie niedokrwiennej. W celu zmniejszenia niestabilności dopasowano eksploracyjne rozszerzone modele z regularyzacją grzbietową (ridge-penalized), zawierające zmienne obrazowe, czas wystąpienia objawów oraz zmienne dotyczące leczenia. W analizie punktowej (landmark analysis) z 24-godzinnym oknem usunięto pacjentów, u których stan pogorszył się przed upływem 24 h, a model punktowy porównano z modelami uzupełnionymi o SD SBP, współczynnik zmienności (CV) lub średnią rzeczywistą zmienność (ARV). Modele te zinterpretowano jako dynamiczne predykcje po przyjęciu do szpitala. Przestrzegano zasad przejrzystego raportowania modeli predykcyjnych16. Analizę rewizyjną zapewniającą odtwarzalność przeprowadzono przy użyciu Python 3.13.5, pandas 2.3.3, NumPy 2.3.3 oraz SciPy 1.16.2.