Protokół badania został przejrzany i zatwierdzony przez Komisję Etyczną Szpitala Medycyny Tradycyjnej Chińskiej w hrabstwie Santai w Mianyang, prowincja Sichuan (nr zatwierdzenia 2025003). Badanie przeprowadzono zgodnie z zasadami Deklaracji Helsińskiej. Przed włączeniem do badania od każdego pacjenta lub krewnego pierwszego stopnia uzyskano pisemną świadomą zgodę. W tym retrospektywnym badaniu kohortowym przeanalizowano dokumentację kliniczną osób starszych hospitalizowanych z powodu ostrego zaostrzenia przewlekłej obturacyjnej choroby płuc (AECOPD) w Oddziale Chorób Płuc i Intensywnej Terapii Szpitala Medycyny Tradycyjnej Chińskiej w hrabstwie Santai w okresie od maja 2023 do maja 2025 roku. Badanie zostało zaprojektowane w celu zidentyfikowania rutynowo dostępnych klinicznych predyktorów nadciśnienia płucnego (PH) oraz opracowania nomogramu do indywidualnej oceny ryzyka PH. Instrumenty laboratoryjne, system echokardiograficzny, materiały do pobierania krwi, narzędzia do ekstrakcji danych oraz oprogramowanie statystyczne wykorzystane w tym procesie wymieniono w Tabeli Materiałów.
1. Przesiew pacjentów i ocena kwalifikowalności
Pacjenci hospitalizowani z rozpoznaniem wypisowym AECOPD w okresie badania zostali wyłonieni poprzez system elektronicznej dokumentacji medycznej pacjentów stacjonarnych. Każda dokumentacja została przejrzana w celu potwierdzenia, że pacjent ma co najmniej 60 lat oraz że rozpoznanie przewlekłej obturacyjnej choroby płuc jest zgodne z kryteriami diagnostycznymi Global Initiative for Chronic Obstructive Lung Disease z 2021 roku15. AECOPD zdefiniowano jako ostre pogorszenie objawów oddechowych, które wymagało hospitalizacji i dodatkowego leczenia. W przypadku pacjentów z więcej niż jedną hospitalizacją z powodu AECOPD w okresie badania, do analizy włączono jedynie pierwszą kwalifikującą się hospitalizację.
Pacjentów włączono do badania, gdy spełnione zostały wszystkie poniższe kryteria: wiek co najmniej 60 lat; hospitalizacja z powodu AECOPD między majem 2023 a majem 2025; wykonane badanie echokardiograficzne przezklatkowe podczas indeksowej hospitalizacji; dostępna pierwsza na czczo próbka krwi żylnej po przyjęciu do szpitala; oraz kompletne informacje kliniczne, laboratoryjne i echokardiograficzne niezbędne do budowy modelu. Wykluczono pacjentów z ciężką niewydolnością wątroby, schyłkową niewydolnością nerek, aktywnym procesem nowotworowym, nadciśnieniem płucnym wynikającym z innych znanych przyczyn, pierwotną chorobą lewej komory serca, chorobą tkanki łącznej, taką jak twardzina układowa, aktywną gruźlicą płuc lub w przypadku braku kluczowych zmiennych niezbędnych do klasyfikacji wyników lub budowy modelu.
Przed usunięciem rekordów ze względu na niepełne informacje, sprawdzono wiek, płeć, historię palenia, historię spożywania alkoholu, historię alergii na leki, nadciśnienie tętnicze, cukrzycę, chorobę wieńcową, astmę, rozstrzenie oskrzelowe, infekcję płuc, rozedmę, niewydolność oddechową, niewydolność serca, status nadciśnienia płucnego w echokardiografii, wskaźniki pełnej krwi, albuminy, kreatyninę, wysokoczułe białko C-reaktywne, D-dimer, fibrynogen, peptyd natriuretyczny typu B oraz stosunek neutrofili do limfocytów. Liczbę przesiewowych i wykluczonych rekordów, przyczyny wykluczenia oraz ostateczną liczbę włączonych przypadków odnotowano w celu sporządzenia schematu wyboru pacjentów (Rycina 1).

Rycina 1: Przesiew pacjentów i tworzenie kohorty. Diagram przepływu przedstawia identyfikację rekordów, ocenę kwalifikowalności, wykluczenia oraz ostateczny podział starszych pacjentów hospitalizowanych z powodu AECOPD. Pacjentów zaklasyfikowano do grupy bez PH lub z PH na podstawie wyników echokardiografii przezklatkowej. Skróty: AECOPD, ostre zaostrzenie przewlekłej obturacyjnej choroby płuc; PH, nadciśnienie płucne. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.
2. Grupowanie kliniczne i definicja diagnostyczna nadciśnienia płucnego
Włączeni pacjenci zostali podzieleni na grupę bez nadciśnienia płucnego (non-PH) oraz grupę z nadciśnieniem płucnym (PH) na podstawie wyników echokardiografii przezklatkowej wykonanej podczas hospitalizacji. Schemat diagnostyczny oparto na wytycznych Europejskiego Towarzystwa Kardiologicznego oraz Europejskiego Towarzystwa Oddechowego z 2015 roku dotyczących nadciśnienia płucnego16. W tym retrospektywnym zbiorze danych za podstawę diagnostyczną przyjęto spoczynkową echokardiografię przezklatkową. Badania przeprowadzili przeszkoleni lekarze echokardiografii, stosując standardowe przekroje przyścianowe, wierzchołkowe i podżebrowe. Przed przydzieleniem do grup zweryfikowano interpretowalność obrazów echokardiograficznych, ponieważ słabe okna akustyczne są częste u pacjentów z przewlekłą obturacyjną chorobą płuc. Rekord uznano za odpowiedni do klasyfikacji PH, gdy sygnał Dopplera z niedomykalności zastawki trójdzielnej był mierzalny lub gdy na podstawie interpretowalnych przekrojów serca można było ocenić rozmiar prawej komory, grubość jej ścian oraz funkcję skurczową.
Skurczowe ciśnienie w tętnicy płucnej oszacowano na podstawie maksymalnej prędkości cofania przepływu przez zastawkę trójdzielną oraz szacowanego ciśnienia w prawym przedsionku. Ciśnienie w prawym przedsionku oceniano na podstawie średnicy dolnej żyły głównej i stopnia jej zapadania się podczas wdechu, gdy oba pomiary były dostępne. Nadciśnienie płucne (PH) zdefiniowano jako szacowane skurczowe ciśnienie w tętnicy płucnej wynoszące co najmniej 40 mmHg. W przypadkach, gdy nie można było wiarygodnie oszacować skurczowego ciśnienia w tętnicy płucnej, klasyfikacja PH wymagała pośrednich dowodów echokardiograficznych na zajęcie prawej komory, w tym powiększenia prawej komory, przerostu prawej komory, dysfunkcji skurczowej prawej komory, spłaszczenia przegrody międzykomorowej lub poszerzenia tętnicy płucnej. Dokumentacja echokardiograficzna z nieczytelnymi obrazami, niemożliwymi do prześledzenia sygnałami Dopplera z zastawki trójdzielnej lub niespójnymi znaleziskami w obrębie prawej strony serca została zweryfikowana przez drugiego lekarza echokardiografistę przed ostatecznym przydziałem do grupy. Dokumentacja, która pozostała nieinterpretowalna po weryfikacji, została wykluczona z końcowej analizy. W końcowej kohorcie kryteria kwalifikacji spełniło 230 pacjentów w wieku podeszłym; 120 pacjentów przydzielono do grupy bez PH, a 110 pacjentów do grupy z PH.
3. Ekstrakcja zmiennych demograficznych i klinicznych
Wszystkie dane kliniczne pobrano z elektronicznego systemu dokumentacji medycznej pacjentów stacjonarnych, korzystając ze standaryzowanego formularza gromadzenia danych. Odnotowano płeć, wiek, historię palenia tytoniu, historię spożywania alkoholu, historię alergii na leki, nadciśnienie tętnicze, cukrzycę, chorobę wieńcową, astmę, rozstrzenie oskrzeli, infekcje płuc, rozedmę, niewydolność oddechową oraz niewydolność serca. Historię palenia tytoniu zdefiniowano jako zużycie co najmniej 100 papierosów w ciągu całego życia, a historię spożywania alkoholu jako regularne spożywanie alkoholu co najmniej dwa razy w tygodniu.
Zakażenie płuc diagnozowano na podstawie objawów ze strony układu oddechowego, oznak fizykalnych, wyników obrazowania klatki piersiowej oraz, gdy były dostępne, dowodów mikrobiologicznych. Niewydolność serca diagnozowano w oparciu o objawy kliniczne, oznaki fizykalne, poziom peptydu natriuretycznego typu B oraz wyniki badania echokardiograficznego. Informacje o chorobach współistniejących weryfikowano z rozpoznaniami przy wypisie, wpisami w kartach przebiegu choroby, raportami z badań obrazowych, wynikami badań laboratoryjnych oraz dokumentacją konsultacyjną, w zależności od dostępności.
4. Gromadzenie i przetwarzanie zmiennych laboratoryjnych
Jako źródło danych laboratoryjnych wykorzystano pierwszą próbkę krwi żylnej pobraną na czczo rano, w dniu następującym po przyjęciu do szpitala. Punkt czasowy ten odpowiadał pierwszemu porannemu pobraniu krwi po przyjęciu i przed rutynową dzienną korektą leczenia, o ile informacje te były dostępne w dokumentacji medycznej. Z szpitalnego systemu informacyjnego laboratorium wyekstrahowano wskaźniki morfologii krwi, wskaźniki biochemiczne, markery krzepnięcia, markery stanu zapalnego oraz markery stresu sercowego.
Z raportu z morfologii krwi zapisano liczbę neutrofili, liczbę limfocytów, liczbę płytek krwi, średnią objętość płytki krwi oraz współczynnik zmienności rozkładu szerokości erytrocytów. Stosunek neutrofili do limfocytów obliczono w następujący sposób:

Z odpowiednich raportów laboratoryjnych odnotowano poziomy albumin, kreatyniny, D-dimerów, fibrynogenu, białka C-reaktywnego wysokiej czułości oraz peptydy natriuretyczne typu B. W całym zbiorze danych zastosowano spójne jednostki pomiarowe: albuminy w g/L, kreatyninę w µmol/L, białko C-reaktywne wysokiej czułości w mg/L, D-dimery w µg/mL, fibrynogen w g/L, liczbę neutrofili, liczbę limfocytów oraz liczbę płytek krwi w ×109/L, średnią objętość płytek krwi w fL oraz peptydy natriuretyczne typu B w pg/mL.
5. Czyszczenie i przygotowanie danych
Zbiór danych został sprawdzony przed przeprowadzeniem analizy statystycznej. Każdy włączony rekord zweryfikowano pod kątem kwalifikowalności, przypisania do grupy, kompletności wyników laboratoryjnych, interpretowalności echokardiograficznej oraz dostępności potencjalnych predyktorów. Usunięto bezpośrednie identyfikatory osobowe, a każdemu pacjentowi przypisano numer identyfikacyjny badania. Rekordy powielone, powtarzające się hospitalizacje, niespójne etykiety grup, błędne daty oraz nieprawdopodobne wartości laboratoryjne sprawdzono w odniesieniu do oryginalnej elektronicznej dokumentacji medycznej i raportu laboratoryjnego. Nie przeprowadzano imputacji dla kluczowych zmiennych wynikowych ani predyktorów, ponieważ ostateczny model oparto na analizie przypadków kompletnych.
Zmienne ciągłe przeanalizowano pod kątem charakterystyki rozkładu. Zmienne o rozkładzie normalnym przedstawiono jako średnią ± odchylenie standardowe, zmienne o rozkładzie nienormalnym jako medianę z rozstępem międzykwartylowym, a zmienne kategoryczne jako liczebność i wartości procentowe. Przed modelowaniem wielowymiarowym zbadano zmienne klinicznie powiązane pod kątem redundancji. Liczbę neutrofili oraz liczbę limfocytów przeanalizowano wraz ze stosunkiem neutrofili do limfocytów; infekcję płuc przeanalizowano wraz z liczbą neutrofili, wysokoczułym białkiem C-reaktywnym oraz stosunkiem neutrofili do limfocytów; niewydolność serca przeanalizowano wraz z peptydem natriuretycznym typu B i echokardiograficznymi obrazami serca; D-dimer i fibrynogen przeanalizowano jako zmienne związane z krzepliwością, a liczbę płytek krwi, średnią objętość płytek krwi oraz współczynnik zmienności rozkładu szerokości erytrocytów jako zmienne hematologiczne.
6. Jednowymiarowe porównanie między grupami
Porównano charakterystykę demograficzną, choroby współistniejące, wskaźniki laboratoryjne oraz zmienne dotyczące grupowania echokardiograficznego pomiędzy grupą bez PH a grupą z PH. W przypadku zmiennych ciągłych o rozkładzie normalnym zastosowano test t dla próby niezależnej. Dla zmiennych ciągłych o rozkładzie innym niż normalny wykorzystano test U Manna-Whitneya. W zależności od oczekiwanej liczby obserwacji w komórkach, dla zmiennych kategorycznych zastosowano test χ2 lub dokładny test Fishera. Za istotną statystycznie uznano dwustronną wartość P mniejszą niż 0,05. Zmienne wykazujące istotne różnice między grupami zostały zakwalifikowane do oceny wielowymiarowej wraz z klinicznie istotnymi predyktorami.
7. Wieloczynnikowa analiza regresji logistycznej
Do identyfikacji niezależnych predyktorów związanych z PH u starszych pacjentów z AECOPD zastosowano binarną regresję logistyczną. Jako zmienna zależna posłużył status PH. Predyktory ciągłe wprowadzono w ich oryginalnych jednostkach pomiaru, a predyktory kategoryczne zakodowano jako zmienne binarne. Zmienne kandydackie wybrano na podstawie wyników analizy jednowymiarowej oraz istotności klinicznej. W pierwszej kolejności zidentyfikowano zmienne z P < 0,05 w porównaniu jednowymiarowym, a następnie poddano je analizie pod kątem redundancji, interpretowalności i nakładania się klinicznego.
Pierwotne zmienne kandydackie obejmowały białko C reaktywne wysokiej czułości, D-dimery, współczynnik zmienności szerokości rozkładu erytrocytów, liczbę płytek krwi, średnią objętość płytek krwi, stosunek neutrofilów do limfocytów, peptyd natriuretyczny typu B, chorobę wieńcową oraz astmę w wywiadzie. Liczba neutrofilów i liczba limfocytów nie były wprowadzane jednocześnie ze stosunkiem neutrofilów do limfocytów, ponieważ stanowiły one składowe tego stosunku. Zakażenia płuc i niewydolność serca zostały ocenione podczas analizy redundancji. Aby sprawdzić, czy ich wykluczenie wpłynęło na interpretację modelu, przeprowadzono analizę wrażliwości poprzez dodanie zakażeń płuc i niewydolności serca do głównego zestawu predyktorów. Porównano oszacowania regresji, dyskryminację, kalibrację oraz kliniczną korzyść netto pomiędzy modelem głównym a modelem wrażliwości.
Multikoliniowość wśród potencjalnych predyktorów oceniono za pomocą czynników inflacji wariancji. Zmienne kategoryczne zakodowano przed obliczeniami, a zmienne ciągłe wprowadzono w tych samych jednostkach, których użyto do modelowania regresji. Przyjęto, że czynnik inflacji wariancji poniżej 5 wskazuje na brak silnej multikoliniowości. Dla zmiennych włączonych do modelu wielowymiarowego podano współczynniki regresji, błędy standardowe, wartości χ2 Walda, ilorazy szans, 95% przedziały ufności, wartości P oraz czynniki inflacji wariancji.
8. Konstrukcja nomogramu
Nomogram został opracowany na podstawie niezależnych predyktorów uwzględnionych w wieloczynnikowym modelu regresji logistycznej. Ostatecznymi predyktorami były: D-dimer, stosunek neutrofili do limfocytów, peptyd natriuretyczny typu B oraz historia astmy. Model regresji logistycznej przekształcono w nomogram punktowy przy użyciu oprogramowania R i pakietu rms. Dla zbioru danych modelowych wygenerowano obiekt rozkładu danych, dopasowano model regresji logistycznej, a współczynniki regresji wykorzystano do przypisania punktów każdemu predyktorowi. W przypadku każdego pacjenta zsumowano punkty za D-dimer, stosunek neutrofili do limfocytów, peptyd natriuretyczny typu B oraz historię astmy, aby uzyskać wynik całkowity, który następnie przyporządkowano do szacowanego prawdopodobieństwa wystąpienia PH.
Nomogram zinterpretowano zgodnie ze standardowymi procedurami opartymi na punktacji opisanymi w poprzednich badaniach nad modelami predykcji klinicznej17,18. Wartość D-dimerów pacjenta, stosunek neutrofili do limfocytów, wartość peptydu natriuretycznego typu B oraz status astmatyczny naniesiono na odpowiadające im osie. Od wartości każdego predyktora poprowadzono linię pionową do skali punktowej. Punkty zsumowano w celu uzyskania wyniku całkowitego, a następnie z osi wyniku całkowitego poprowadzono linię pionową do osi przewidywanego prawdopodobieństwa, aby uzyskać zindywidualizowaną szacowaną szansę wystąpienia PH. Ilustracyjny przypadek pacjenta wykorzystany w manuskrypcie został wybrany na podstawie wartości mieszczących się w granicach lub bliskich obserwowanemu rozkładowi klinicznemu w kohorcie.
9. Ocena wydajności modelu
Dyskryminację modelu oceniono za pomocą analizy krzywej charakterystyki operacyjnej odbiornika (ROC) oraz pola pod krzywą (AUC). Pole pod krzywą przedstawiono wraz z 95% przedziałem ufności. Nomogram porównano z poszczególnymi predyktorami oraz z modelem łączącym ilościowe markery laboratoryjne. Kalibrację oceniono poprzez porównanie prawdopodobieństw przewidywanych z zaobserwowanymi wynikami nadciśnienia płucnego (PH). Krzywą kalibracji wygenerowano z wykorzystaniem metody bootstrapowej. Przedstawiono punkt przecięcia kalibracji, nachylenie kalibracji, średni błąd bezwzględny, wynik Briera oraz test Hosmera-Lemeshow. Obliczono również skorygowany metodą bootstrap punkt przecięcia kalibracji oraz skorygowane metodą bootstrap nachylenie kalibracji. Test Hosmera-Lemeshow zinterpretowano łącznie z krzywą kalibracji i ilościowymi wskaźnikami kalibracji.
Użyteczność kliniczną oceniono za pomocą analizy krzywej decyzyjnej. Korzyść netto wykreślono w zależności od prawdopodobieństw progowych i porównano z referencyjnymi strategiami „lecz wszystkich” oraz „nie lecz nikogo”. Zidentyfikowano zakres prawdopodobieństwa progowego, w którym model zapewniał wyższą korzyść netto niż strategie referencyjne.
10. Walidacja wewnętrzna i powtarzalność
Walidację wewnętrzną przeprowadzono za pomocą metody bootstrapowania z 1 000 powtórzeniami. W każdym powtórzeniu bootstrapowym model był ponownie dopasowywany i testowany w celu oszacowania optymizmu w wydajności modelu. Obliczono i przedstawiono wskaźniki dyskryminacji oraz kalibracji skorygowane o optymizm wraz z pozorną wydajnością modelu.
Analizy statystyczne przeprowadzono przy użyciu standardowego oprogramowania statystycznego oraz programu R w wersji 4.5.2. Pakiety R zostały wykorzystane do modelowania regresji logistycznej, konstrukcji nomogramu, analizy krzywych charakterystyki operacyjnej odbiornika (ROC), oceny kalibracji, walidacji wewnętrznej oraz analizy krzywych decyzji. Nomogram został skonstruowany przy użyciu pakietu rms w programie R.