$$\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łuc. Kohorta składała się z pacjentów leczonych zgodnie z tym protokołem oraz porównywalnej grupy kontrolnej otrzymującej standardową opiekę okołooperacyjną.
Implementacja protokołu i charakterystyka wyjściowa
Wyjściowe cechy kliniczne, demograficzne i chirurgiczne badanej populacji podsumowano w Tabeli 1. spośród 277 pacjentów, u których przeprowadzono anatomiczną operację klatki piersiowej w placówce 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 wdrażania 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 pomyślnie u 46,0% (vs. 15,0%), przestrzeganie zasad wentylacji chroniącej płuca osiągnęło poziom 92,7% (vs. 52,0%), a wczesną mobilizację w 1. dobie wdrożono u 88,7% (vs. 37,0%); 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, bez 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), wartościach progowych wskaźnika masy ciała, historii palenia tytoniu, rozpowszechnieniu 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 lub zastosowanej metodzie chirurgicznej (wszystkie P > 0,05, Tabela 1). Ta porównywalność statystyczna skutecznie wyklucza historyczny błąd selekcji i waliduje późniejsze porównania wyników między obiema 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). Korzystne, choć nieistotne statystycznie, trendy zaobserwowano 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
Wskaźniki zakażeń ran (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 częstości występowania przedłużającego się przecieku powietrza (1.3% vs. 4.7%, P = 0.187), wysięku opłucnej wymagającego drenażu oraz powikłań moczopochodnych i żołądkowo-jelitowych (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ą z zastosowania zstandaryzowanego protokołu była wyraźna poprawa efektywności wczesnej rekonwalescencji oraz doświadczeń zgłaszanych przez pacjentów. Pooperacyjny czas pobytu w szpitalu był statystycznie porównywalny pomiędzy grupą kompleksowego postępowania a grupą opieki konwencjonalnej (7,16 ± 3,32 dni vs. 7,02 ± 3,66 dni, P = 0,731, Tabela 2). Analiza rozkładu czasu pobytu w szpitalu, porównana wizualnie poprzez wykres pudełkowy na Rysunku 3, wykazała wyraźny profil prawostronnie skośny w obu kohortach, z szerokim zakresem rozciągającym się od minimum 4 dni do maksimum 24 dni. Wysoka zmienność wynikała z indywidualnej dynamiki powrotu do zdrowia, a statystyczne wartości odstające wykraczające poza próg 14 dni zaobserwowano w obu grupach przy użyciu kryteriów wykresu pudełkowego Tukeya, co bezpośrednio odpowiadało osobom cierpiącym na opóźnione ustąpienie przecieków powietrza niskiego stopnia. Jednak grupa interwencyjna wykazała znacznie lepsze wskaźniki jakości rekonwalescencji: maksymalny wynik bólu w numerycznej skali oceny (NRS) w 1–3 dobie pooperacyjnej był istotnie niższy w kohorcie kompleksowego postępowania (3,47 ± 1,01 vs. 5,17 ± 1,45, P < 0,001), co reprezentuje poprawę efektywności klinicznej kontroli bólu. Ponadto, globalne oceny satysfakcji pacjentów dokonane przy wypisie były wyraźnie 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 przedstawiono na Rysunku 3.
Wyniki przeżywalności i analiza czynników ryzyka
Wielodyscyplinarny protokół wykazał znaczący wpływ ochronny na wyniki przeżywalności krótkoterminowej.
Czynniki ryzyka
Aby zbudować solidny model prognostyczny i zapobiec pominięciu istotnych parametrów zakłócających, zmienne kliniczne wyjściowe zostały systematycznie przesianie za pomocą jednowymiarowej regresji logistycznej. Zmienne wykazujące próg istotności wyjściowej P < 0,10 w przesiewie jednowymiarowym, a konkretnie: zaawansowany wiek, upośledzona funkcja płuc wyjściowo, przedoperacyjne COPD 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 kontroli zmiennych zakłócających. Analiza wielowymiarowa zidentyfikowała kompleksowy perioperacyjny protokół zarządzania jako 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 końcowym modelu 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), przedoperacyjne COPD (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 z dokładnymi 95% CI zostały przedstawione graficznie w wykresie leśnym w 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żyła się, lecz nie osiągnęła pełnej istotności statystycznej (P = 0,071, Tabela 2). Analiza przeżywalności metodą Kaplana-Meiera, z rozbiciem 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 (number-at-risk) umieszczone poniżej osi czasu; konkretnie, w grupie interwencyjnej pula osób narażonych na ryzyko wynosiła odpowiednio 150, 144 i 135 osób przeżywających w dniach 0, 30 i 90, w porównaniu do odpowiednio 127, 114 i 89 osób w kohorcie opieki konwencjonalnej.
Ponowne przyjęcie
Wskaźnik nieplanowanych readmisji 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 opis protokołu wielodyscyplinarnego postępowania okołooperacyjnego. Diagram przedstawia trójfazową ścieżkę (przedoperacyjną, śródoperacyjną, pooperacyjną) wraz z kluczowymi interwencjami i rolami odpowiedzialnych członków zespołu. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 2: Porównanie powikłań pooperacyjnych pomiędzy grupami. Wykres słupkowy przedstawiający częstość występowania głównych powikłań pooperacyjnych w grupie kompleksowego zarządzania w porównaniu z grupą opieki konwencjonalnej. Każda kategoria powikłań jest 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ę tej ryciny.

Rysunek 3: Rozkład czasu pobytu w szpitalu po operacji. Wykres pudełkowy ilustrujący porównywalny czas pobytu w szpitalu po operacji pomiędzy grupą kompleksowego zarządzania (7,16 3,32 dni) a grupą opieki konwencjonalnej (7,02 3,66 dni, 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). Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

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 wielowymiarowej analizy regresji logistycznej dla 90-dniowej śmiertelności. Linia pionowa przy OR = 1.0 oznacza brak efektu, a dokładne niezależne szacunki ryzyka/czynników ochronnych są wyświetlane obok każdego istotnego predyktora. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 5: Analiza przeżycia całkowitego metodą Kaplana-Meiera w krótkich okresach obserwacji. (A) Krzywe przeżycia całkowitego w ciągu 30 dni, porównujące prawdopodobieństwo przeżycia pooperacyjnego pomiędzy grupą kompleksowego leczenia a grupą opieki konwencjonalnej. (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 wyświetlić powiększoną wersję tej ryciny.
| Zmienna | Konwencjonalna (n=127) | Kompleksowe (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) | |
| Samiec | 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 nikotynowy, n (%) | | | 0.237 |
| Nie | 52 (40.9) | 50 (33.3) | |
| Tak | 75 (59.1) | 100 (66.7) | |
| Choroby współistniejące: POChP, n (%) | | | 0.916 |
| Nie | 104 (81.9) | 121 (80.7) | |
| Tak | 23 (18.1) | 29 (19.3) | |
| Choroby współistniejące: 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 dostępu chirurgicznego, 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 (endobronchialna ultrasonografia z biopsją cienkoigłową) | 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) | |
| Wynik w skali 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 (wymuszona pojemność wydechowa)₁% (ś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ędzykwartylny) lub liczba (%). Wartości P wyznaczono za pomocą testu χ2 lub testu t, w zależności od zastosowania. Skróty: VATS = wideochirurgia klatki piersiowej; COPD = przewlekła obturacyjna choroba płuc; FEV₁ = natężona objętość wydechowa w 1 s; ASA = klasyfikacja American Society of Anesthesiologists; BMI = wskaźnik masy ciała; CVD = choroby układu sercowo-naczyniowego.
| Zmienna wynikowa | Konwencjonalna (n=127) | Kompleksowa (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 protekcyjna dla płuc, 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ólna liczba powikłań, 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 powietrzny | 6 (4.7) | 2 (1.3) | 0.187 |
| Powikłania moczowe lub żołądkowo-jelitowe | 1 (0.8) | 3 (2.0) | 0.736 |
| Wskaźniki rekonwalescencji | | | |
| Czas pobytu w szpitalu (dni), średnia ± SD | 7.02 ± 3.66 | 7.16 ± 3.32 | 0.731 |
| Maksymalny wynik bólu w skali NRS (0–10) | 5.17 ± 1.45 | 3.47 ± 1.01 | <0.001 |
| Wynik satysfakcji pacjenta | 83.06 ± 6.13 | 91.82 ± 3.67 | <0.001 |
| Wyniki obserwacji (follow-up) | | | |
| Wskaźnik rehospitalizacji w ciągu 30 dni, n (%) | 20 (15.7) | 20 (13.3) | 0.691 |
| Śmiertelność 30-dniowa, n (%) | 13 (10.2) | 6 (4.0) | 0.071 |
| Śmiertelność 90-dniowa, 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 przedstawiono jako liczbę (%), 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 t-testu 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) zapisane w dniach pooperacyjnych 1–3, a ogólna satysfakcja pacjenta została zmierzona za pomocą ustrukturyzowanej 100-punktowej wizualnej skali analogowej (VAS) podczas wypisu.
| Zmienna | Jednowymiarowe OR (95% CI) | Wartość P | Wielowymiarowe 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 sercowo-naczyniowego | 1.74 (1.08–2.80) | 0.023 | 1.68 (1.02–2.78) | 0.041 |
| Przedsoperacyjne FEV₁ < 60% | 1.62 (1.01–2.60) | 0.045 | 1.62 (1.01–2.60) | 0.048 |
| Kompleksowe prowadzenie | 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 | -- | -- |
| Wywiad palenia | 0.78 (0.43–1.45) | 0.432 | -- | -- |
| Stopień zaawansowania 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łuc. OR = iloraz szans; CI = przedział ufności; VATS = torakoskopia wspomagana wideo; 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. Punkt końcowy = 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.