Początkowo do badania włączono 150 pacjentów z zaawansowanym nowotworem, którzy spełniali kryteria kwalifikacji, w tym 82 osoby do grupy badawczej i 68 do grupy kontrolnej. Po zastosowaniu PSM w stosunku 1:1 pomyślnie wygenerowano 60 dobrze dopasowanych par, co dało końcową kohortę analityczną obejmującą 120 pacjentów. Po dopasowaniu nie zaobserwowano istotnych różnic między grupami w żadnej z ocenianych zmiennych towarzyszących: demograficznych (wiek, płeć, stan cywilny, wykształcenie), czynników onkologicznych (typ/stadium guza, przerzuty do kości), stanu funkcjonalnego (ECOG), profilu bólu (wyjściowa skala NRS, rodzaj, czas trwania) oraz wartości laboratoryjnych (hemoglobina, albuminy, kreatynina). Wszystkie zmienne towarzyszące osiągnęły równowagę (P > 0,05, SMD < 0,1), co potwierdza doskonałą porównywalność grup po dopasowaniu. Tabela 1 przedstawia wyjściową charakterystykę populacji badawczej po dopasowaniu.
W 2. tygodniu po leczeniu grupa badawcza wykazała znacznie niższe wyniki w skali NRS (2,8 ± 1,1 vs 4,2 ± 1,3, P < 0,001) oraz wyższy odsetek odpowiedzi na leczenie przeciwbólowe (redukcja NRS ≥50%: 58,3% vs 35,0%, P = 0,011) w porównaniu z grupą kontrolną. Ocena PMI wykazała znacznie mniejszą liczbę przypadków niewystarczającej analizji w grupie badawczej niż w grupie kontrolnej (16,7% vs 41,7%, P = 0,006), co przedstawiono w Tabeli 2.
Po 2 tygodniach grupa badawcza wykazała istotnie lepszą jakość życia, z wyższymi wynikami w dziedzinie funkcjonalnej i niższymi wynikami w dziedzinie objawów w porównaniu z grupą kontrolną. Grupa badawcza wykazała wyższą jakość życia w porównaniu z grupą kontrolną. Tabela 3 przedstawia wyniki jakości życia przed i po leczeniu dla obu grup.
Po 2 tygodniach interwencji pacjenci poddani kompleksowemu leczeniu bólu osiągnęli znacznie większą redukcję dyskomfortu psychicznego w porównaniu z grupą kontrolną, co potwierdzono niższymi wynikami w skali HADS dla lęku (6,2 ± 2,8 vs 9,5 ± 3,1) oraz depresji (5,8 ± 2,5 vs 8,9 ± 3,0) (w obu przypadkach P < 0,001; Tabela 4).
Po dwóch tygodniach interwencji pacjenci otrzymujący kompleksowe leczenie przeciwbólowe wymagali znacznie mniejszej eskalacji dawki opioidów (wzrost o 47,26%) w porównaniu do osób leczonych standardowo (wzrost o 65,27%, P < 0,001), co wskazuje na korzystną skuteczność przeciwbólową przy mniejszym obciążeniu opioidami (Tabela 5).
W trakcie trwania badania nie zaobserwowano istotnych różnic międzygrupowych w częstości występowania działań niepożądanych (P > 0.05). Najczęściej zgłaszanymi zdarzeniami były zaparcia oraz nudności/wymioty. Choć częstość występowania zaparć była nieznacznie niższa w kohorcie prowadzonej w ramach kompleksowego zarządzania w porównaniu z grupą kontrolną, różnica międzygrupowa nie osiągnęła istotności statystycznej (Tabela 6).
Wyniki wielokrotnej analizy regresji liniowej (Tabela 7) wykazały, że po skorygowaniu o charakterystykę demograficzną i zmienne kliniczne, kompleksowe leczenie bólu (grupa badawcza) było istotnie powiązane z poprawą jakości życia (β = 14,042, 95% CI: 8,03 do 20,06, P < 0,001), lęku (β = −3,295, 95% CI: −4,37 do −2,22, P < 0,001) oraz depresji (β = −3,078, 95% CI: −4,08 do −2,08, P < 0,001). Wyższy wyjściowy wynik w skali NRS wiązał się z większą poprawą jakości życia (β = 2,694, P = 0,037). Wyższe wyjściowe wyniki w skalach HADS-A i HADS-D były istotnie powiązane z odpowiednio większą poprawą w zakresie lęku (β = −0,993, P < 0,001) i depresji (β = −0,904, P < 0,001), co sugeruje, że pacjenci z gorszym wyjściowym stanem psychicznym mogą odnieść większe korzyści z kompleksowego leczenia. Wiek i płeć nie były istotnie powiązane z żadnym z trzech wyników (wszystkie P > 0,05).
Wielokrotne analizy wrażliwości przyniosły wyniki zgodne z analizą główną, przedstawione w Tabeli 8. Szacunki efektu pozostały stabilne przy zastosowaniu alternatywnych algorytmów dopasowania — dopasowania jądrowego (kernel matching; różnica średnich NRS: −1.5, 95% CI: −1.9 do −1.1; domena funkcjonalna EORTC: +13.1, 95% CI: 8.1 do 18.1) oraz ważenia odwrotnej prawdopodobieństwem leczenia (IPTW: różnica średnich NRS: −1.5, 95% CI: −1.9 do −1.1; domena funkcjonalna EORTC: +12.8, 95% CI: 7.9 do 17.7) — przy wszystkich wartościach P < 0.001, co sugeruje, że wyniki nie były wrażliwe na wybór metody analitycznej. Analizy podgrup w rozwarstwieniu według typu nowotworu i wyjściowej intensywności bólu wykazały spójne efekty kierunkowe przemawiające za grupą kompleksowego zarządzania bólem we wszystkich podgrupach. Najbardziej wyraźne różnice zaobserwowano u pacjentów z umiarkowanym bólem wyjściowym (różnica średnich: −1.8, 95% CI: −2.5 do −1.0) oraz u pacjentów z rakiem płuca (różnica średnich: −1.5, 95% CI: −2.0 do −1.0) lub rakiem przewodu pokarmowego (różnica średnich: −1.5, 95% CI: −2.3 do −0.6). Pacjenci z innymi typami nowotworów również wykazali istotną różnicę (różnica średnich: −0.7, 95% CI: −1.3 do −0.1, P = 0.027), co sugeruje, że zaobserwowane powiązania generalizują się na różnorodne populacje pacjentów onkologicznych. Podsumowując, analizy wrażliwości sugerują, że wyniki związane z kompleksowym zarządzaniem bólem były dość stabilne, a wyniki analizy głównej nie zostały w znacznym stopniu zawpływnięte przez wybór metody analitycznej lub heterogeniczność populacji.
Dostępność danych
Ze względu na wymogi dotyczące poufności danych pacjentów oraz ograniczenia nałożone przez Komisję Etyczną Szpitala Ludowego w powiecie Lezhi, dane wygenerowane i/lub analizowane w toku ninzego badania nie mogą zostać przekazane do publicznie dostępnego repozytorium danych. Zanonimizowane dane mogą zostać udostępnione wykwalifikowanym badaczom na uzasadnioną prośbę, po uprzednim rozpatrzeniu i zatwierdzeniu przez Komisję Etyczną oraz podpisaniu wymaganej umowy o udostępnieniu danych. Prośby należy kierować do autora korespondencyjnego, który będzie koordynował proces weryfikacji instytucjonalnej.

Rycina 1: Schemat projektu badania. Rycina ta ilustruje proces wyboru i grupowania pacjentów. Szczegółowo przedstawiono kryteria włączenia i wyłączenia, wstępną rekrutację 150 pacjentów, przydział do grupy badawczej (kompleksowe leczenie bólu) i grupy kontrolnej (standardowa farmakoterapia) oraz końcową kohortę 120 pacjentów po dopasowaniu wskaźnika skłonności w stosunku 1:1. Skrót: PSM = propensity score matching (dopasowanie wskaźnika skłonności). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Zmienna | Grupa badawcza (n=60) | Grupa kontrolna (n=60) | Wartość p | SMD |
| Wiek (lata, średnia±SD) | 58.9±11.8 | 59.1±11.9 | 0.919 | 0.017 |
| Płeć [n(%)] | | | 1 | 0.022 |
| Mężczyzna | 33 (55.0) | 34 (56.7) | | |
| Kobieta | 27 (45.0) | 26 (43.3) | | |
| Stan cywilny [n(%)] | | | 0.845 | 0.048 |
| Zamężny/żonaty | 42 (70.0) | 44 (73.3) | | |
| Inny (niezamężna/nieżonaty/rozwiedziony/wdowa/wdowiec) | 18 (30.0) | 16 (26.7) | | |
| Poziom wykształcenia [n(%)] | | | 0.856 | 0.048 |
| Średnie lub wyższe | 28 (46.7) | 26 (43.3) | | |
| Podstawowe lub gimnazjalne | 32 (53.3) | 34 (56.7) | | |
| Typ nowotworu [n(%)] | | | 0.999 | 0.063 |
| Rak płuca | 25 (41.7) | 26 (43.3) | | |
| Nowotwór przewodu pokarmowego | 16 (26.7) | 15 (25.0) | | |
| Rak piersi | 11 (18.3) | 12 (20.0) | | |
| Inny | 8 (13.3) | 7 (11.7) | | |
| Stopień zaawansowania nowotworu [n(%)] | | | 1 | 0.026 |
| Stopień III | 23 (38.3) | 24 (40.0) | | |
| Stopień IV | 37 (61.7) | 36 (60.0) | | |
| Stan sprawności według skali ECOG [n(%)] | | | 0.553 | 0.035 |
| 0-1 | 25 (41.7) | 26 (43.3) | | |
| 2 | 28 (46.7) | 27 (45.0) | | |
| 3-4 | 7 (11.7) | 7 (11.7) | | |
| Wyjściowy wynik NRS (punkty, średnia±SD) | 6.4±1.2 | 6.5±1.2 | 0.637 | 0.083 |
| Typ bólu [n(%)] | | | 0.653 | 0.048 |
| Nocyceptywny | 27 (45.0) | 26 (43.3) | | |
| Neuropatyczny | 13 (21.7) | 12 (20.0) | | |
| Mieszany | 20 (33.3) | 22 (36.7) | | |
| Czas trwania bólu [n(%)] | | | 0.839 | 0.048 |
| <3 miesiące | 22 (36.7) | 20 (33.3) | | |
| ≥3 miesiące | 38 (63.3) | 40 (66.7) | | |
| Przerzuty do kości [n(%)] | 35 (58.3) | 34 (56.7) | 1 | 0.033 |
| Parametry laboratoryjne (średnia±SD) | | | | |
| Hemoglobina (g/L) | 112.5±18.3 | 114.2±17.6 | 0.596 | 0.095 |
| Albumina (g/L) | 36.20 [33.10, 38.85] | 35.95 [31.83, 39.45] | 0.768 | 0.099 |
| Kreatynina (μmol/L) | 73.75 [61.70, 91.00] | 74.75 [60.25, 91.40] | 0.568 | 0.096 |
Tabela 1: Charakterystyka wyjściowa obu grup po dopasowaniu wskaźnika skłonności. Zmienne ciągłe: wiek, wyjściowa skala NRS oraz hemoglobina miały rozkład normalny i są przedstawione jako średnia ± SD, a porównano je za pomocą sparowanych testów t; albuminy i kreatynina miały rozkład nienormalny i są przedstawione jako mediana [Q1, Q3], a porównano je za pomocą sparowanych testów znakowych Wilcoxona. Zmienne kategoryczne: sparowane proporcje porównano za pomocą testu McNemara (zmienne binarne) lub testu McNemara-Bowkera (zmienne wielomianowe). Wszystkie SMD były <0,1, co potwierdza odpowiednie zrównoważenie współzmiennych między grupami. Skróty: PSM = dopasowanie wskaźnika skłonności; SMD = standaryzowana różnica średnich; ECOG = stan sprawności według Eastern Cooperative Oncology Group; NRS = numeryczna skala oceny.
| Wynik | Grupa badawcza (n=60) | Grupa kontrolna (n=60) | Wartość statystyczna | wartość p |
| Wynik w skali NRS po leczeniu (średnia±SD) | 2.8±1.1 | 4.2±1.3 | t=-6.015 | <0.001 |
| Odsetek reakcji bólowej [n(%)] | 35 (58.3) | 21(35.0) | χ²=6,43 | 0.011 |
| Niewystarczająca analgezja [n(%)] | 10(16.7) | 25(41.7) | χ²=7,259 | 0.006 |
Tabela 2: Porównanie wyników kontroli bólu po leczeniu.Odpowiedź przeciwbólowa: redukcja NRS o ≥50% względem wartości wyjściowej. Niedostateczna analgezja: ujemny PMI. Metody statystyczne: test t dla danych ciągłych; test McNemara dla danych kategorycznych. Skróty: NRS = Numeric Rating Scale (numeryczna skala oceny bólu); PMI = Pain Management Index (wskaźnik kontroli bólu).
| Dziedzina | Grupa badawcza (n=60) | Grupa kontrolna (n=60) | Różnica między grupami (95% CI) | Wartość P |
| Dziedzina funkcjonalna | | | +13.2 (8.3, 18.1) | <0.001 |
| Przed leczeniem | 48.6±10.2 | 49.1±10.5 | | |
| Po leczeniu | 68.5±12.4 | 55.3±13.6 | | |
| Dziedzina objawów | | | -12.8 (-16.5, -9.1) | <0.001 |
| Przed leczeniem | 42.3±9.8 | 41.8±10.3 | | |
| Po leczeniu | 25.6±8.9 | 38.4±11.2 | | |
Tabela 3: Wyniki EORTC QLQ-C30 przed i po leczeniu: porównanie między grupami (x̄±s). Jakość życia mierzono za pomocą chińskiej wersji kwestionariusza EORTC QLQ-C30. Wyniki przekształcono liniowo na skalę 0–100; wyższe wyniki wskazują na lepsze funkcjonowanie w Domenie Funkcjonalnej i gorsze objawy w Domenie Objawów. Wszystkie dane przedstawiono jako średnia ± SD. „Wartość netto różnicy między grupami” reprezentuje różnicę między grupą badawczą a grupą kontrolną w wynikach po leczeniu (Grupa badawcza − Grupa kontrolna). Porównania między grupami przeprowadzono za pomocą parowych testów t. Skrót: EORTC QLQ-C30 = European Organization for Research and Treatment of Cancer Quality of Life Core Scale.
| Objaw | Grupa badawcza (n=60) | Grupa kontrolna (n=60) | Różnica między grupami (95% CI) | Wartość p |
| Lęk (HADS-A) | | | -3.3 (-4.4, -2.2) | <0.001 |
| Obróbka wstępna | 12.8±3.2 | 12.5±3.4 | | |
| Po opracowaniu | 6.2±2.8 | 9.5±3.1 | | |
| Depresja (HADS-D) | | | -3.1 (-4.1, -2.1) | <0.001 |
| Obróbka wstępna | 11.6±3.0 | 11.9±3.3 | | |
| Po opracowaniu | 5.8±2.5 | 8.9±3.0 | | |
Tabela 4: Wyniki HADS przed i po leczeniu: porównanie między grupami (x̄±s). HADS (podskale: HADS-A lęk, HADS-D depresja; zakres 0–21 dla każdej podskali, wyższe wyniki wskazują na większe nasilenie objawów). „Wartość netto różnicy między grupami” reprezentuje różnicę między grupą badawczą a grupą kontrolną w wynikach po leczeniu (Grupa badawcza − Grupa kontrolna). Wartości ujemne wskazują na większą poprawę w grupie badawczej. Porównania międzygrupowe przeprowadzono za pomocą testów t dla par zależnych, zgodnie z doborem par. Skrót: HADS = Hospital Anxiety and Depression Scale (Szpitalna Skala Lęku i Depresji).
| Parametr | Grupa badawcza (n=60) | Grupa kontrolna (n=60) | Wartość statystyczna | wartość p |
| Wartość wyjściowa MEDD (mg/24 h) | 48.5±15.6 | 46.8±16.2 | t=0,612 | 0.543 |
| MEDD po leczeniu (mg/24 h) | 69.6±22.0 | 78.2±23.5 | t=-2,122 | 0.038 |
| Wzrost MEDD (%) | 47.26 [32.07, 57.67] | 65.27 [59.73, 83.77] | Z = −5,786 | <0.001 |
Tabela 5: Zmiany dawki przeciwbólowej przed i po leczeniu: porównanie między grupami (x̄±s). MEDD (morfina doustna mg/24 h wyliczona za pomocą czynników konwersji). Wzór na zmianę MEDD: [(Po leczeniu − Baseline) / Baseline] × 100%. Baseline oraz MEDD po leczeniu analizowano za pomocą testów t dla prób zależnych. Przyrost MEDD nie miał rozkładu normalnego i był porównywany przy użyciu testu znaków Wilcoxona dla prób zależnych, przedstawiono go jako medianę [pierwszy kwartyl, trzeci kwartyl]. Skrót: MEDD = równoważna dobowa dawka morfiny.
| Zdarzenie niepożądane | Grupa badawcza (n=60) | Grupa kontrolna (n=60) | Wartość p |
| Zaparcia | 32 (53.3) | 38 (63.3) | 0.146 |
| Nudności i wymioty | 24 (40.0) | 28 (46.7) | 0.125 |
| Senność | 12 (20.0) | 15 (25.0) | 0.25 |
| Zatrzymanie moczu | 4 (6.7) | 6 (10.0) | 0.508 |
| Depresja oddechowa | 1 (1.7) | 2 (3.3) | 0.56 |
Tabela 6: Częstość zdarzeń niepożądanych: porównanie między grupami [n(%)]. Dane przedstawiono jako n (%). Porównania przeprowadzono przy użyciu testu McNemara dla binarnych danych parowanych. W przypadku rzadkich zdarzeń (zatrzymanie moczu i depresja oddechowa) dokładne wartości P uzyskano z rozkładu dwumianowego (obustronnego). Nie podano statystyk testowych, ponieważ test McNemara nie dostarcza wartości chi-kwadrat w tabelach 2 × 2.
| Zmienna zależna | Zmienna niezależna | β | Zstandaryzowane β | 95% CI | Wartość P |
| Poprawa jakości życia (ΔEORTC domena funkcjonalna) | Grupa badawcza (vs grupa kontrolna) | 14.042 | 0.387 | 8.03, 20.06 | <0.001 |
| Wiek (na każde 10 lat wzrostu) | 0.146 | 0.095 | -0.12, 0.41 | 0.277 |
| Kobieta (vs mężczyzna) | -5.551 | -0.152 | -11.74, 0.64 | 0.078 |
| Wyjściowy wynik NRS (na każdy 1 punkt wzrostu) | 2.694 | 0.178 | 0.17, 5.22 | 0.037 |
| Przerzuty do kości (vs brak przerzutów do kości) | -0.126 | -0.003 | -6.55 6.30 | 0.969 |
| Wynik ECOG 2-4 (vs 0-1) | 1.708 | 0.047 | -4.44, 7.86 | 0.583 |
| Poprawa lęku (ΔHADS-A) | Grupa badawcza (vs grupa kontrolna) | -3.295 | -0.348 | -4.37, -2.22 | <0.001 |
| Wiek (na każde 10 lat wzrostu) | 0.014 | 0.036 | -0.03, 0.06 | 0.535 |
| Kobieta (vs mężczyzna) | -0.42 | -0.044 | -1.50, 0.66 | 0.445 |
| Wyjściowy wynik HADS-A (na każdy 1 punkt wzrostu) | -0.993 | -0.688 | -1.16, -0.83 | <0.001 |
| Poprawa depresji (ΔHADS-D) | Grupa badawcza (vs grupa kontrolna) | -3.078 | -0.37 | -4.08, -2.08 | <0.001 |
| Wiek (na każde 10 lat wzrostu) | -0.019 | -0.052 | -0.06, 0.02 | 0.391 |
| Kobieta (vs mężczyzna) | 0.067 | 0.008 | -0.96, 1.09 | 0.897 |
| Wyjściowy wynik HADS-D (na każdy 1 punkt wzrostu) | -0.904 | -0.681 | -1.07, -0.74 | <0.001 |
Tabela 7: Regresja liniowa wielokrotna: predyktory poprawy jakości życia i zdrowia psychicznego. ΔEORTC w dziedzinie funkcjonalnej = po leczeniu minus przed leczeniem; wartości dodatnie wskazują na poprawę jakości życia. ΔHADS = po leczeniu minus przed leczeniem; wartości ujemne wskazują na zmniejszenie poziomu lęku/depresji. Modele skorygowano o wszystkie wymienione zmienne. Model QoL: R2 = 0.212, skorygowane R2 = 0.170, F = 5.070, P < 0.001. Model lęku: R2 = 0.621, skorygowane R2 = 0.607, F = 47.043, P < 0.001. Model depresji: R2 = 0.575, skorygowane R2 = 0.560, F = 38.940, P < 0.001. Skróty: EORTC = European Organization for Research and Treatment of Cancer; HADS = Hospital Anxiety and Depression Scale.
| Metoda analityczna | Główny punkt końcowy | Średnia różnica (95% CI) | Wartość P |
| Analiza podstawowa | | | |
| Dopasowanie PSM 1:1 do najbliższego sąsiada | Wynik NRS | −1.4 (−1.9 do −0.9) | <0.001 |
| Domena funkcjonalna EORTC | +13.2 (8.3 do 18.1) | <0.001 |
| Lęk wg HADS | −3.3 (−4.4 do −2.2) | <0.001 |
| Alternatywne metody dopasowania/ważenia | | | |
| Dopasowanie jądrowe (Kernel matching) | Wynik NRS | −1.5 (−1.9 do −1.1) | <0.001 |
| Domena funkcjonalna EORTC | +13.1 (8.1 do 18.1) | <0.001 |
| Ważenie odwrotnością prawdopodobieństwa leczenia (IPTW) | Wynik NRS | −1.5 (−1.9 do −1.1) | <0.001 |
| Domena funkcjonalna EORTC | +12.8 (7.9 do 17.7) | <0.001 |
| Analizy podgrup (wynik NRS) | | | |
| Stratyfikacja według typu nowotworu | | | |
| Rak płuca | Wynik NRS | −1.5 (−2.0 do −1.0) | <0.001 |
| Nowotwór przewodu pokarmowego | Wynik NRS | −1.5 (−2.3 do −0.6) | 0.001 |
| Rak piersi | Wynik NRS | −1.0 (−1.7 do −0.3) | 0.008 |
| Inne | Wynik NRS | −0.7 (−1.3 do −0.1) | 0.027 |
| Stratyfikacja według wyjściowej intensywności bólu | | | |
| Ból umiarkowany (NRS 4–6) | Wynik NRS | −1.8 (−2.5 do −1.0) | <0.001 |
| Ból silny (NRS 7–10) | Wynik NRS | −1.7 (−2.4 do −1.0) | <0.001 |
Tabela 8: Wyniki analizy wrażliwości. Różnica średnich reprezentuje różnicę między grupą badawczą a grupą kontrolną (Badawcza − Kontrolna). Wartości ujemne wskazują na większą poprawę w grupie badawczej w odniesieniu do wyników dotyczących bólu; wartości dodatnie wskazują na większą poprawę w grupie badawczej w odniesieniu do jakości życia. Trendy w domenach funkcjonalnych EORTC i skali HADS były zgodne z analizami podstawowymi (dane HADS nie zostały przedstawione w tej tabeli dla zwięzłości). Skróty: PSM = dopasowanie wskaźnika tendencji (propensity score matching); IPTW = ważenie odwrotnej prawdopodobieństwem leczenia (inverse probability of treatment weighting); EORTC = Europejska Organizacja Badań i Leczenia Nowotworów (European Organization for Research and Treatment of Cancer); HADS = Hospital Anxiety and Depression Scale; NRS = Numeryczna Skala Bólu (Numerical Rating Scale).