$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Aby ocenić efekty zastosowania tego kompleksowego protokołu opieki, w badaniu zebrano i przeanalizowano dane kliniczne pacjentów poddawanych resekcji raka płuca. Kohorta składała się z pacjentów leczonych zgodnie z tym protokołem oraz z porównywalnej grupy kontrolnej otrzymującej konwencjonalną opiekę okołooperacyjną.
Wdrożenie protokołu i charakterystyka wyjściowa
Wyjściowa charakterystyka kliniczna, demograficzna i chirurgiczna populacji badanej została podsumowana w Tabeli 1. spośród 277 pacjentów, u których przeprowadzono anatomiczną operację klatki piersiowej w danej instytucji między styczniem 2021 a grudniem 2022 roku, 150 objęto kompleksowym protokołem opieki (grupa interwencyjna), a 127 otrzymało opiekę konwencjonalną (grupa kontrolna). Wskaźniki wierności wdrożenia i przestrzegania zaleceń, poddane audytowi w całym kontinuum leczenia, szczegółowo opisano w Tabeli 2. Kluczowe wskaźniki procesowe wskazały na bardzo skuteczną realizację protokołu: ukończenie prehabilitacji osiągnięto u 90,7% pacjentów w grupie interwencyjnej (vs 29,9% w grupie kontrolnej), ustrukturyzowane wsparcie żywieniowe zapewniono w 46,0% (vs 15,0%), przestrzeganie zasad wentylacji chroniącej płuca osiągnęło 92,7% (vs 52,0%), a wczesną mobilizację w 1. dobie wdrożono u 88,7% pacjentów (vs 37,0%), przy czym wszystkie wskaźniki procesowe wykazały wysoce istotną poprawę (wszystkie P < 0,001) w kohorcie interwencyjnej.
Charakterystyka wyjściowa obu grup była dobrze zrównoważona, nie stwierdzono statystycznie istotnych różnic w średnim wieku (56,04 ± 12,80 vs. 54,06 ± 13,60 lat, P = 0,214), rozkładzie płci żeńskiej (62,0% vs. 66,9%, P = 0,467), progach wskaźnika masy ciała, historii palenia tytoniu, częstotliwości występowania głównych chorób współistniejących, w tym przewlekłej obturacyjnej choroby płuc (POChP, 19,3% vs. 18,1%, P = 0,916) i chorób układu sercowo-naczyniowego, stopniu zaawansowania nowotworu ani zastosowanej metodzie chirurgicznej (wszystkie P > 0,05, Tabela 1). Ta porównywalność statystyczna skutecznie wyklucza historyczny błąd selekcji i potwierdza zasadność późniejszych porównań wyników między dwiema grupami.
Powikłania pooperacyjne
Zastosowanie kliniczne protokołu kompleksowej opieki wiązało się ze znaczną redukcją krótkoterminowej zachorowalności pooperacyjnej, co zostało szczegółowo zestawione w Tabeli 2.
Powikłania płucne
Częstość występowania ogólnych powikłań pooperacyjnych była znacząco niższa w grupie interwencyjnej (10,0% vs. 29,1%, P < 0,001). Efekt ochronny ten był najbardziej wyraźny w przypadku zapalenia płuc (3,3% vs. 11,0%, P = 0,022). Zaobserwowano korzystne, choć nieistotne statystycznie trendy w odniesieniu do niedodmy (2,0% vs. 1,6%, P = 1,000) oraz niewydolności oddechowej (0,7% vs. 2,4%, P = 0,501).
Powikłania pozapłucne
Odsetki zakażeń rany (0,0% vs. 3,1%, P = 0,092) oraz zakrzepicy żył głębokich (1,3% vs. 1,6%, P = 1,000) były również niższe w grupie interwencyjnej, choć różnice te nie osiągnęły istotności statystycznej. Nie stwierdzono istotnych różnic między grupami w odsetkach przedłużającego się odmy (1,3% vs. 4,7%, P = 0,187), płynu w jamie opłucnej wymagającego drenażu oraz powikłań moczowo-płciowych/gastrycznych (2,0% vs. 0,8%, P = 0,736). Kompleksowy profil porównawczy powikłań pooperacyjnych przedstawiono na Rysunku 2.
Czas pobytu w szpitalu i efektywność rekonwalescencji
Główną korzyścią ze stosowania zestandaryzowanego protokołu była wyraźna poprawa wydajności wczesnej rekonwalescencji oraz doświadczeń zgłaszanych przez pacjentów. Pooperacyjna długość pobytu w szpitalu była statystycznie porównywalna pomiędzy grupą kompleksowego zarządzania a grupą opieki konwencjonalnej (7,16 ± 3,32 dnia vs. 7,02 ± 3,66 dnia, P = 0,731, Tabela 2). Analiza rozkładu długości pobytu w szpitalu, porównana wizualnie za pomocą wykresu pudełkowego na Rysunku 3, wykazała wyraźny prawoskewny profil w obu kohortach, z szerokim zakresem od minimum 4 dni do maksimum 24 dni. Wysoka zmienność wynikała z indywidualnej dynamiki powrotu do zdrowia, a statystyczne wartości odstające przekraczające próg 14 dni zaobserwowano w obu grupach przy zastosowaniu kryteriów wykresu pudełkowego Tukeya, co bezpośrednio odpowiadało osobom z opóźnionym ustąpieniem niskostopniowych przecieków powietrza. Jednak grupa interwencyjna wykazała znacznie lepsze wskaźniki jakości rekonwalescencji: maksymalna ocena bólu w numerycznej skali NRS w 1–3 dobie po operacji była istotnie niższa w kohorcie kompleksowego zarządzania (3,47 ± 1,01 vs. 5,17 ± 1,45, P < 0,001), co wskazuje na poprawę wydajności klinicznej kontroli bólu. Co więcej, ogólne oceny satysfakcji pacjentów dokonane przy wypisie były znacznie wyższe w kohorcie protokołu (91,82 ± 3,67 vs. 83,06 ± 6,13, P < 0,001, Tabela 2), a wzorce rozkładu danych zilustrowano na Rysunku 3.
Wyniki przeżywalności i analiza czynników ryzyka
Wielodyscyplinarny protokół wykazał znaczący wpływ ochronny na krótkoterminowe wyniki przeżywalności.
Czynniki ryzyka
Aby zbudować wiarygodny model prognostyczny i zapobiec pominięciu istotnych parametrów zakłócających, zmienne kliniczne wyjściowe zostały systematycznie przeanalizowane za pomocą jednowymiarowej regresji logistycznej. Zmienne wykazujące w badaniu jednowymiarowym próg istotności P < 0,10, a konkretnie: zaawansowany wiek, upośledzona wyjściowa funkcja płuc, przedoperacyjna POChP oraz choroby układu sercowo-naczyniowego, wraz z główną zmienną ekspozycji (Kompleksowy Protokół Zarządzania), zostały wprowadzone do końcowego wielowymiarowego modelu regresji logistycznej w celu skorygowania zmiennych zakłócających. Analiza wielowymiarowa wykazała, że kompleksowy perioperacyjny protokół zarządzania stanowi silny, niezależny czynnik ochronny przeciwko śmiertelności (Iloraz Szans [OR] = 0,36, 95% Przedział Ufności [CI]: 0,18–0,74, P = 0,005). Inne niezależne predyktory ryzyka zidentyfikowane w modelu końcowym obejmowały zaawansowany wiek (≥70 lat, OR = 1,98, 95% CI: 1,13–3,48, P = 0,019), upośledzoną funkcję płuc (przedoperacyjne FEV₁ <60%, OR = 1,62, 95% CI: 1,01–2,60, P = 0,048), przedoperacyjną POChP (OR = 1,62, 95% CI: 1,02–2,60, P = 0,044) oraz choroby układu sercowo-naczyniowego (OR = 1,68, 95% CI: 1,02–2,78, P = 0,041). Szczegółowe parametry regresji jednowymiarowej i wielowymiarowej przedstawiono w Tabeli 3, a końcowe ilorazy szans wraz z dokładnymi 95% CI zostały przedstawione graficznie w wykresie leśnym na Rysunku 4.
Śmiertelność
Wskaźnik śmiertelności w ciągu 90 dni był znacznie niższy w grupie interwencyjnej (10,0% vs. 29,9%, P < 0,001), co stanowi istotną kliniczną korzyść w zakresie przeżywalności. Zaobserwowano również niższy wskaźnik śmiertelności w ciągu 30 dni (4,0% vs. 10,2%), choć różnica ta zbliżała się do pełnej istotności statystycznej, ale jej nie osiągnęła (P = 0,071, Tabela 2). Analiza przeżywalności metodą Kaplana-Meiera, podzielona na przeżywalność 30-dniową (Rysunek 5A) i 90-dniową (Rysunek 5B), wykazała wyraźne rozdzielenie krzywych na korzyść kompleksowego protokołu. Ta długofalowa przewaga w przeżywalności jest w pełni potwierdzona przez zaznaczone cieniowaniem 95% przedziały ufności oraz odpowiadające im tabele liczby osób narażonych na ryzyko umieszczone poniżej osi czasu; konkretnie, grupa interwencyjna utrzymywała pulę ryzyka wynoszącą odpowiednio 150, 144 i 135 osób przeżyjących w dniach 0, 30 i 90, w porównaniu do 127, 114 i 89 osób w kohorcie leczonej metodami konwencjonalnymi.
Ponowna hospitalizacja
Wskaźnik nieplanowanych ponownych przyjęć do szpitala w ciągu 30 dni był liczbowo niższy w grupie interwencyjnej (13,3% vs. 15,7%), jednak różnica ta nie była istotna statystycznie (P = 0,691, Tabela 2).

Rycina 1: Schematyczny przegląd protokołu wielodyscyplinarnego postępowania w okresie okołooperacyjnym. Diagram przedstawia trójfazową ścieżkę (przedoperacyjną, śródoperacyjną i pooperacyjną) wraz z kluczowymi interwencjami i rolami odpowiedzialnych członków zespołu. Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.

Rysunek 2: Porównanie powikłań pooperacyjnych pomiędzy grupami. Wykres słupkowy przedstawiający częstość występowania głównych powikłań pooperacyjnych w grupach opieki kompleksowej oraz konwencjonalnej. Każda kategoria powikłań została opatrzona odpowiednią istotnością statystyczną (wartość P lub „ns” dla porównań nieistotnych), wyznaczoną na podstawie końcowego, zweryfikowanego zbioru danych. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Rozkład czasu pobytu w szpitalu po operacji. Wykres pudełkowy ilustrujący porównywalny czas pobytu w szpitalu po operacji w grupie opieki kompleksowej (7,16 ± 3,32 dnia) oraz w grupie opieki konwencjonalnej (7,02 ± 3,66 dnia, P = 0,731, oznaczone jako „ns”). Czerwone okrągłe znaczniki wskazują statystyczne wartości odstające wykraczające poza rozstęp międzykwartylny (IQR). Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Analiza wielowymiarowa czynników związanych z 90-dniową śmiertelnością. Wykres leśny (forest plot) przedstawiający ilorazy szans (OR) i 95% przedziały ufności (CI) z analizy wielowymiarowej regresji logistycznej dla 90-dniowej śmiertelności. Pionowa linia przy OR = 1.0 oznacza brak efektu, a dokładne niezależne oszacowania ryzyka/czynników ochronnych są wyświetlane obok każdego istotnego predyktora. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 5: Analiza przeżycia całkowitego metodą Kaplana-Meiera w ramach krótkoterminowej obserwacji. (A) Krzywe przeżycia całkowitego w ciągu 30 dni, porównujące pooperacyjne prawdopodobieństwa przeżycia między grupą prowadzoną kompleksowo a grupą z opieką konwencjonalną. (B) Krzywe przeżycia całkowitego w ciągu 90 dni, porównujące przeżywalność podłużną między obiema kohortami. Oba panele zawierają zacieniowane pasma 95% przedziału ufności oraz odpowiadające im tabele liczby osób narażonych na ryzyko, umieszczone poniżej osi czasu. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
| Zmienna | Konwencjonalna (n=127) | Kompleksowa (n=150) | Wartość p |
| Wiek (średnia ± SD) | 54,06 ± 13,60 | 56,04 ± 12,80 | 0,214 |
| Płeć, n (%) | | | 0,467 |
| Kobieta | 85 (66,9) | 93 (62,0) | |
| Mężczyzna | 42 (33,1) | 57 (38,0) | |
| BMI, n (%) | | | 0,381 |
| < 25 kg/m² | 87 (68,5) | 111 (74,0) | |
| ≥ 25 kg/m² | 40 (31,5) | 39 (26,0) | |
| Wywiad palenia tytoniu, n (%) | | | 0,237 |
| Nie | 52 (40,9) | 50 (33,3) | |
| Tak | 75 (59,1) | 100 (66,7) | |
| Choroba współistniejąca: POChP, n (%) | | | 0,916 |
| Nie | 104 (81,9) | 121 (80,7) | |
| Tak | 23 (18,1) | 29 (19,3) | |
| Choroba współistniejąca: CVD, n (%) | | | 0,69 |
| Nie | 100 (78,7) | 114 (76,0) | |
| Tak | 27 (21,3) | 36 (24,0) | |
| Stopień zaawansowania guza, n (%) | | | 0,346 |
| Stopień I–II | 63 (49,6) | 84 (56,0) | |
| Stopień III–IV | 64 (50,4) | 66 (44,0) | |
| Metoda chirurgiczna, n (%) | | | 0,552 |
| Biopsja pod kontrolą TK | 0 (0,0) | 1 (0,7) | |
| Lobektomia VATS | 31 (24,4) | 35 (23,3) | |
| Segmentektomia VATS | 94 (74,0) | 112 (74,7) | |
| EBUS-TBNA | 0 (0,0) | 1 (0,7) | |
| Bronchoskopia | 0 (0,0) | 1 (0,7) | |
| Torakoskopowa resekcja zmiany | 1 (0,8) | 0 (0,0) | |
| Torakoskopowa resekcja klinowa | 1 (0,8) | 0 (0,0) | |
| Wskaźnik ASA, n (%) | | | 0,135 |
| I | 20 (15,7) | 32 (21,3) | |
| II | 82 (64,6) | 79 (52,7) | |
| III | 25 (19,7) | 39 (26,0) | |
| FEV₁% (średnia ± SD) | 70,23 ± 15,11 | 73,11 ± 13,65 | 0,097 |
Tabela 1: Charakterystyka wyjściowa pacjentów z rakiem płuca w momencie przyjęcia (n = 277). Wartości przedstawiono jako średnia ± odchylenie standardowe, mediana (rozstęp międzykwartylowy) lub liczba (%). Wartości P wyznaczono za pomocą testu χ2 lub testu t, zależnie od przypadku. Skróty: VATS = torakoskopia wspomagana wideo; COPD = przewlekła obturacyjna choroba płuc; FEV₁ = natężona objętość wydechowa pierwszosekundowa; ASA = American Society of Anesthesiologists; BMI = wskaźnik masy ciała; CVD = choroby układu sercowo-naczyniowego.
| Zmienna wynikowa | Konwencjonalna (n=127) | Kompleksowe (n=150) | wartość p |
| Zgodność z protokołem | | | |
| Ukończenie prehabilitacji, n (%) | 38 (29.9) | 136 (90.7) | <0.001 |
| Ustrukturyzowane wsparcie żywieniowe, n (%) | 19 (15.0) | 69 (46.0) | <0.001 |
| Wentylacja chroniąca płuca, n (%) | 66 (52.0) | 139 (92.7) | <0.001 |
| Wczesna mobilizacja w 1. dobie, n (%) | 47 (37.0) | 133 (88.7) | <0.001 |
| Powikłania pooperacyjne | | | |
| Ogólne powikłania, n (%) | 37 (29.1) | 15 (10.0) | <0.001 |
| Zapalenie płuc | 14 (11.0) | 5 (3.3) | 0.022 |
| Atelektaza | 2 (1.6) | 3 (2.0) | 1 |
| Niewydolność oddechowa | 3 (2.4) | 1 (0.7) | 0.501 |
| Zakażenie rany | 4 (3.1) | 0 (0.0) | 0.092 |
| Zakrzepica żył głębokich | 2 (1.6) | 2 (1.3) | 1 |
| Przedłużający się przeciek powietrza | 6 (4.7) | 2 (1.3) | 0.187 |
| Układ moczowy lub pokarmowy | 1 (0.8) | 3 (2.0) | 0.736 |
| Wskaźniki odzysku | | | |
| Czas pobytu w szpitalu (dni), średnia ± SD | 7.02 ± 3.66 | 7.16 ± 3.32 | 0.731 |
| Maksymalna ocena bólu w skali NRS (0–10) | 5.17 ± 1.45 | 3.47 ± 1.01 | <0.001 |
| Wskaźnik satysfakcji pacjenta | 83.06 ± 6.13 | 91.82 ± 3.67 | <0.001 |
| Wyniki obserwacji po badaniu | | | |
| wskaźnik readmisji w ciągu 30 dni, n (%) | 20 (15.7) | 20 (13.3) | 0.691 |
| śmiertelność w ciągu 30 dni, n (%) | 13 (10.2) | 6 (4.0) | 0.071 |
| śmiertelność po 90 dniach, n (%) | 38 (29.9) | 15 (10.0) | <0.001 |
Tabela 2: Porównanie powikłań pooperacyjnych, efektywności rekonwalescencji i przestrzegania protokołu pomiędzy grupą kompleksowego zarządzania a grupą opieki konwencjonalnej (n = 277). Wartości wyrażono jako liczba (%), z wyjątkiem ciągłych wskaźników rekonwalescencji i jakości (średnia ± SD). Tabela systematycznie zestawia zarówno wskaźniki procesu wdrażania (wskaźniki przestrzegania zasad prehabilitacji, żywienia, wentylacji i wczesnej mobilizacji), jak i kompleksowe wyniki kliniczne. Wartości P uzyskano za pomocą testu chi-kwadrat lub dokładnego testu Fishera dla zmiennych kategorycznych oraz testu t dla zmiennych ciągłych. Powikłania oceniano w ciągu 14 dni po operacji. Wskaźniki bólu reprezentują maksymalne wyniki w numerycznej skali oceny (NRS) zarejestrowane w dniach pooperacyjnych 1–3, a ogólną satysfakcję pacjentów zmierzono za pomocą ustrukturyzowanej 100-punktowej wizualnej skali analogowej (VAS) w momencie wypisu.
| Zmienna | Jednowymiarowy OR (95% CI) | Wartość p | Wielowymiarowy OR (95% CI) | Wartość p |
| POChP | 1.66 (1.06–2.59) | 0.026 | 1.62 (1.02–2.60) | 0.044 |
| Choroby układu krążenia | 1.74 (1.08–2.80) | 0.023 | 1.68 (1.02–2.78) | 0.041 |
| Przedsoperacyjny FEV₁ < 60% | 1.62 (1.01–2.60) | 0.045 | 1.62 (1.01–2.60) | 0.048 |
| Kompleksowe leczenie | 0.45 (0.28–0.72) | <0.001 | 0.36 (0.18–0.74) | 0.005 |
| Wiek ≥ 70 lat | 1.79 (1.03–3.10) | 0.039 | 1.98 (1.13–3.48) | 0.019 |
| Płeć (mężczyzna) | 1.11 (0.59–2.05) | 0.736 | -- | -- |
| BMI ≥ 25 kg/m² | 1.53 (0.80–2.86) | 0.191 | -- | -- |
| Historia palenia | 0.78 (0.43–1.45) | 0.432 | -- | -- |
| Stadium guza (III–IV) | 1.47 (0.81–2.70) | 0.208 | -- | -- |
Tabela 3: Jednowymiarowa i wielowymiarowa analiza regresji logistycznej czynników ryzyka związanych z 90-dniową śmiertelnością u wszystkich pacjentów z rakiem płuca. OR = iloraz szans; CI = przedział ufności; VATS = torakoskopia wspomagana wideowizją; COPD = przewlekła obturacyjna choroba płuc; FEV₁ = natężona objętość wydechowa w 1 s; CVD = choroby układu sercowo-naczyniowego. Do modelowania wielowymiarowego włączono zmienne z P < 0,10 w analizie jednowymiarowej. Wynik = 90-dniowa śmiertelność z wszystkich przyczyn. „--“ oznacza zmienne wykluczone z końcowego modelu wielowymiarowego ze względu na brak niezależnej istotności statystycznej (P ≥ 0,05) podczas etapowej eliminacji wstecznej.