Kryteria kwalifikacji
Badania zostały uwzględnione, jeśli spełniały wszystkie poniższe kryteria, sformułowane zgodnie ze schematem PICOS (Populacja, Interwencja, Komparator, Wyniki, Projekt badania)17. Wszystkie platformy i materiały wykorzystane w niniejszej pracy są wymienione w Tabeli materiałów. Wszystkie wyszukiwania w bazach danych przeprowadzono w listopadzie 2025 roku. Poza programami RevMan 5.3 oraz EndNote 20 nie wykorzystano żadnego dodatkowego komercyjnego oprogramowania (np. SAS, SPSS, R). Ocena GRADE została przeprowadzona ręcznie przy użyciu narzędzia online GRADEpro GDT; wersja została zapisana z datą dostępu.
Populacja (P)
Uczestnicy wykazywali klinicznie istotne objawy depresyjne, definiowane jako formalna diagnoza zaburzeń depresyjnych (MDD) lub dystymii zgodnie z kryteriami DSM (dowolna wersja) lub ICD, potwierdzone w ustrukturyzowanym wywiadzie klinicznym (np. SCID, MINI, CIDI); lub wynik powyżej zwalidowanego punktu odcięcia w standaryzowanym narzędziu przesiewowym do oceny depresji (np. Beck Depression Inventory ≥ 14, Hamilton Depression Rating Scale ≥ 14, CES-D ≥ 16, PHQ-9 ≥ 10, Geriatric Depression Scale ≥ 11). Uczestnicy mogli być w każdym wieku (dzieci, młodzież, dorośli, osoby starsze)18. Badania, w których wzięli udział uczestnicy z chorobami współistniejącymi o podłożu somatycznym lub psychiatrycznym (np. choroby układu krążenia, cukrzyca, zaburzenia lękowe, zaburzenia związane z używkami), zostały włączone tylko wtedy, gdy w analizie podstawowej wyniki dla podgrupy z depresją raportowano oddzielnie lub gdy ≥80% uczestników spełniało kryteria depresji. Wykluczono badania, w których depresja była wtórnym punktem końcowym w populacji osób niezdepresjonowanych (np. pacjenci z niewydolnością serca bez przesiewowej oceny depresji). Nie wprowadzono ograniczeń dotyczących wyjściowego stopnia nasilenia depresji (lekka, umiarkowana, ciężka), jednak stopień nasilenia został wyodrębniony do analiz podgrup.
Interwencja (I)
Interwencja polegała na ustrukturyzowanym, zaplanowanym wysiłku fizycznym, zdefiniowanym jako każda forma ćwiczeń aerobowych, oporowych, mieszanych (aerobowe + oporowe) lub alternatywnych (joga, tai chi, taniec), prowadzonych przez co najmniej 2 tygodnie. Interwencje wysiłkowe mogły być nadzorowane (przez trenera, terapeutę lub instruktora) lub nienadzorowane (wykonywane w domu, samodzielnie). Badania, w których ćwiczenia łączono z inną aktywną interwencją (np. ćwiczenia + farmakoterapia, ćwiczenia + uważność), włączono tylko wtedy, gdy efekt ćwiczeń mógł zostać wyizolowany (np. ćwiczenia + farmakoterapia vs sama farmakoterapia). Wykluczono badania, w których interwencja polegała wyłącznie na poradach dotyczących aktywności fizycznej (bez ustrukturyzowanego planu ćwiczeń) lub poradnictwie w zakresie stylu życia.
Komparator (C)
Grupa kontrolna nie otrzymała ustrukturyzowanej interwencji w postaci ćwiczeń. Dopuszczalne warunki kontrolne sklasyfikowano jako kontrole pasywne (brak aktywnej interwencji terapeutycznej), np. lista oczekujących (leczenie odroczone), opieka standardowa (standardowe postępowanie medyczne bez ustrukturyzowanych ćwiczeń) lub brak interwencji; oraz kontrole aktywne (nie-ćwiczeniowe interwencje terapeutyczne), np. kontrola uwagi (np. edukacja zdrowotna niezwiązana z ćwiczeniami, sesje relaksacyjne lub lekki stretching nie spełniający kryteriów intensywności ćwiczeń), farmakoterapia (wyłącznie leki przeciwdepresyjne) lub psychoterapia (terapia psychologiczna bez ćwiczeń). Wykluczono badania, w których grupa kontrolna otrzymała jakąkolwiek formę ustrukturyzowanych ćwiczeń (w tym minimalne ćwiczenia lub stretching spełniający kryteria intensywności ćwiczeń). Uwaga dotycząca rozbieżności: choć w sekcji innowacji stwierdzono, że wykluczono badania z aktywnymi interwencjami psychologicznymi lub kontrolami w postaci stretchingu/relaksacji, następnie rozszerzono kryteria kwalifikowalności o aktywne grupy porównawcze (farmakoterapię, psychoterapię, kontrolę uwagi), aby odzwierciedlić rzeczywiste pytanie kliniczne dotyczące tego, jak ćwiczenia wypadają w porównaniu z innymi uznanymi metodami leczenia. Rozbieżność ta została omówiona w dyskusji. Wszystkie analizy są stratyfikowane według rodzaju kontroli (pasywna vs. aktywna), aby zbadać wpływ rodzaju kontroli na wielkości efektu.
Wyniki (O)
Badanie zgłaszało co najmniej jedną miarę wyniku specyficzną dla depresji: remisję depresji (dychotomiczną, zdefiniowaną przez autorów badania jako brak spełnienia kryteriów diagnostycznych lub wynik poniżej punktu odcięcia w zwalidowanej skali); LUB (b) zmianę w ciągłym wyniku nasilenia depresji mierzonego w zwalidowanej skali (np. HAM-D, BDI, CES-D, PHQ-9, GDS, MADRS). Wykluczone zostały badania zgłaszające wyłącznie niespecyficzne miary nastroju (np. „dobrostan”, „stan emocjonalny” bez zastosowania zwalidowanej skali depresji).
Projekt badania (S)
Uwzględniono wyłącznie randomizowane badania kontrolowane (RCT). Wykluczono prospektywne badania kohortowe, badania kliniczno-kontrolne, badania retrospektywne, serie przypadków oraz raporty przypadków, ponieważ wiążą się one z wysokim ryzykiem błędów związanych z czynnikami zakłócającymi i błędem selekcji podczas szacowania efektów leczenia. Badania mogły mieć projekt z grupami równoległymi lub naprzemiennymi (w przypadku badań naprzemiennych wyodrębniono tylko pierwszą fazę). Wymagano, aby w momencie rozpoczęcia badania w każdym ramieniu uczestniczyło co najmniej 10 osób, w celu zminimalizowania wpływu małej liczebności prób w badaniach. Nie wprowadzono ograniczeń językowych, jednak badania nieanglojęzyczne włączono tylko wtedy, gdy dostępny był pełny przekład.
Kryteria wykluczenia
Wykluczono badania, w których depresja nie była głównym punktem końcowym lub w których populacja nie była dobierana pod kątem depresji (np. ogólne populacje medyczne, w których depresja była miernikiem wtórnym). Wykluczono badania, w których interwencja w postaci ćwiczeń była realizowana jako element wieloskładnikowej interwencji stylu życia, bez wyizolowanego efektu ćwiczeń. Wykluczono badania, w których grupa kontrolna brała udział w jakiejkolwiek formie ustrukturyzowanych ćwiczeń. Wykluczono badania z niewystarczającymi danymi do obliczenia wielkości efektu (brak raportowanych średnich, odchyleń standardowych lub wyników zmian, których nie udało się uzyskać od autorów). Wykluczono powielone publikacje dotyczące tej samej kohorty badawczej (uwzględniono tylko najbardziej kompletny lub najnowszy raport). Wykluczono abstrakty konferencyjne, rozprawy doktorskie oraz nieopublikowane manuskrypty, do których nie było dostępnego pełnego tekstu.
Źródła informacji
Całe badanie przedstawiono na Rysunku 1 w formie schematu blokowego PRISMA. Literaturę włączono do badania, jeśli spełnione zostały następujące kryteria włączenia: badanie było randomizowanym badaniem kontrolowanym (RCT); wybrani do badania pacjenci cierpieli na depresję; a warunkiem kontrolnym był brak ćwiczeń. Akceptowalne grupy kontrolne obejmowały opiekę standardową (standardowe leczenie medyczne bez ustrukturyzowanych ćwiczeń), listę oczekujących (leczenie odroczone), brak interwencji, kontrolę uwagi (np. edukację zdrowotną niezwiązaną z ćwiczeniami, sesje relaksacyjne lub lekkie rozciąganie niespełniające kryteriów recepty ćwiczeń) lub leczenie farmakologiczne (wyłącznie leki przeciwdepresyjne). Wykluczono badania, w których grupa kontrolna przechodziła jakiekolwiek ustrukturyzowane treningi aerobowe, oporowe lub mieszane. Badania, w których grupy kontrolne otrzymywały psychoterapię lub terapię poznawczo-behawioralną bez komponentu ćwiczeń, włączono tylko wtedy, gdy zostały one jasno zdefiniowane jako grupy niećwiczące.
W badaniu porównano interwencje polegające na ćwiczeniach fizycznych z warunkami kontrolnymi. Grupy kontrolne zdefiniowano jako porównywacze niećwiczące, obejmujące opiekę standardową, listę oczekujących, brak leczenia, kontrolę uwagi (np. edukację zdrowotną, relaksację lub ćwiczenia rozciągające, które nie spełniały kryteriów intensywności ćwiczeń) lub samo leczenie farmakologiczne. Badania wykluczono, jeśli grupa kontrolna uczestniczyła w jakiejkolwiek ustrukturyzowanej interwencji z ćwiczeniami fizycznymi lub jeśli ćwiczenia porównywano wyłącznie z innym schematem ćwiczeń, bez ramienia kontrolnego niećwiczącego. Wykluczono badania, które nie oceniały wpływu interwencji z ćwiczeniami na objawy depresji, badania nad depresją u pacjentów bez ćwiczeń lub grupy kontrolnej oraz badania bez grupy porównawczej.
Postępowanie w przypadku badań wieloramiennych
W badaniach obejmujących więcej niż dwa istotne ramiona (np. wiele interwencji ćwiczeniowych lub wiele grup porównawczych) zastosowano następujące strategie w celu uniknięcia podwójnego liczenia oraz błędów jednostki analizy: w przypadku wielu ramion ćwiczeniowych względem jednego ramienia kontrolnego, gdy badanie porównywało dwie lub więcej różnych interwencji ćwiczeniowych (np. trening aerobowy vs. trening oporowy) z jedną grupą kontrolną, ramiona ćwiczeniowe połączono w jedną zbiorczą grupę ćwiczeń, stosując wzory zalecane w Podręczniku Cochrane (Rozdział 23). Podejście to pozwala uniknąć podwójnego liczenia uczestników grupy kontrolnej przy zachowaniu niezależności statystycznej; w przypadku ćwiczeń względem wielu porównawczych grup niećwiczących, gdy badanie porównywało ćwiczenia z dwiema lub więcej odrębnymi warunkami kontrolnymi (np. standardową opieką i farmakoterapią), wybrano jedynie ten warunek kontrolny, który najlepiej odpowiadał definicji „kontroli niećwiczącej” (kolejność priorytetów: lista oczekujących ≥ standardowa opieka ≥ kontrola uwagi > farmakoterapia > psychoterapia).
W przypadku występowania wielu grup kontrolnych bez ćwiczeń, wybrano najbardziej konserwatywny (najmniej aktywny) punkt odniesienia, aby uniknąć przeszacowania efektu ćwiczeń; ćwiczenia vs. tylko ćwiczenia (brak kontroli bez ćwiczeń). Badania, w których porównywano jedynie różne rodzaje aktywności fizycznej (np. wysoką intensywność vs. niską intensywność) bez grupy kontrolnej niećwiczącej, zostały wykluczone z analizy podstawowej, ponieważ nie spełniały kryterium porównania „ćwiczenia vs. kontrola”. W przypadku interwencji łączonych, gdzie ramię interwencyjne obejmowało ćwiczenia oraz inny aktywny komponent (np. ćwiczenia + farmakoterapia, ćwiczenia + uważność), ramię to włączono tylko wtedy, gdy można było wyizolować efekt ćwiczeń. Jeśli ramię łączone porównywano z samym komponentem niebędącym ćwiczeniami (np. ćwiczenia + farmakoterapia vs. sama farmakoterapia), różnicę interpretowano jako odzwierciedlającą efekt ćwiczeń. Jeśli taka izolacja nie była możliwa, badanie wykluczono. Aby ocenić stabilność strategii agregacji dla badań wieloramiennych, przeprowadzono analizę wrażliwości z wykluczeniem wszystkich badań wieloramiennych, a wyniki porównano z wynikami analizy podstawowej. Wszelkie istotne różnice zostały zaraportowane.
Wybór ramienia kontrolnego w przypadku występowania wielu komparatorów niebędących ćwiczeniami przeprowadzono w następujący sposób: jeśli badanie porównywało ćwiczenia z dwiema lub więcej odrębnymi warunkami kontrolnymi (np. standardową opieką i farmakoterapią), zastosowano hierarchiczną regułę wyboru w celu uniknięcia arbitralnych decyzji, w której priorytet nad kontrolami aktywnymi otrzymały kontrole pasywne (lista oczekujących ≥ standardowa opieka ≥ brak interwencji), aby lepiej wyizolować specyficzny efekt ćwiczeń. Jeśli dostępne były wyłącznie kontrole aktywne, wybrano najbardziej konserwatywny komparator (farmakoterapia ≥ psychoterapia ≥ kontrola uwagi), aby uniknąć przeszacowania efektu ćwiczeń, oraz przeprowadzono analizy wrażliwości obejmujące wszystkie dostępne ramiona kontrolne (poprzez podział grupy ćwiczącej na poszczególne komparatory), aby sprawdzić, czy reguła wyboru wpłynęła na zbiorczą estymację. Wyniki tych analiz wrażliwości przedstawiono w Materiałach uzupełniających.
Strategia wyszukiwania
Przeprowadzono systematyczne przeszukanie następujących baz danych online od momentu ich powstania do listopada 2025 roku: PubMed, Embase (platforma OVID), Cochrane Central Register of Controlled Trials (CENTRAL), MEDLINE (platforma OVID) oraz Google Scholar. Pełne strategie wyszukiwania, obejmujące terminy MeSH, terminy Emtree, słowa kluczowe, operatory logiczne, ucinanie wyrazów oraz tagi pól, przedstawiono w Tabeli 1. W przypadku baz PubMed i Embase wyszukiwanie ograniczono do badań na ludziach w języku angielskim. W przypadku CENTRAL i MEDLINE wyszukiwanie ograniczono do randomizowanych badań kontrolowanych. W Google Scholar przeszukano pierwsze 50 rekordów posortowanych według trafności, stosując uproszczony ciąg wyszukiwania ('exercise AND depression AND (remission OR medication OR treatment)') ze względu na ograniczenia platformy w zakresie złożonego wyszukiwania logicznego19.
Oprócz przeszukiwania baz danych przeprowadzono następujące dodatkowe wyszukiwania: ręczne przeszukiwanie list referencyjnych wszystkich włączonych badań oraz odpowiednich przeglądów systematycznych (metoda kuli śnieżnej), śledzenie cytowań kluczowych włączonych badań za pomocą Google Scholar i Web of Science, ręczne przeszukiwanie kluczowych czasopism (Mental Health and Physical Activity, Journal of Affective Disorders, Depression and Anxiety oraz Psychosomatic Medicine) dla lat 2020–2025 oraz przeszukiwanie rejestrów badań klinicznych (ClinicalTrials.gov, WHO ICTRP) w celu znalezienia nieopublikowanych lub trwających badań z wykorzystaniem terminów wyszukiwania „exercise” AND „depression”. Pełna historia wyszukiwania, w tym liczba rekordów pobranych na każdym etapie dla każdej bazy danych, została przedstawiona w Tabeli uzupełniającej 1. Postępowano zgodnie z listą kontrolną raportowania wyszukiwań PRISMA 2020 (Plik uzupełniający 1).
Proces selekcji
Wszystkie zidentyfikowane rekordy zostały wyeksportowane do programu EndNote 20 w celu usunięcia duplikatów. Po deduplikacji dwóch niezależnych recenzentów przeprowadziło ćwiczenie kalibracyjne na losowej próbie 10 rekordów, aby zapewnić spójne stosowanie kryteriów kwalifikowalności. Po potwierdzeniu akceptowalnej zgodności (κ ≥ 0.80), recenzenci niezależnie przeskanowali tytuły i streszczenia wszystkich pozostałych rekordów pod kątem wcześniej określonych kryteriów kwalifikowalności. Rekordy uznane za potencjalnie istotne przez któregokolwiek z recenzentów zostały przekazane do pełnotekstowej analizy. Pełne teksty artykułów zostały pozyskane i niezależnie ocenione przez tych samych dwóch recenzentów. Przyczyny wykluczenia na etapie analizy pełnego tekstu zostały udokumentowane zgodnie z wytycznymi PRISMA (np. brak interwencji w formie ćwiczeń, brak diagnozy depresji, brak grupy kontrolnej, projekt nieporównawczy, duplikacja publikacji). Rozbieżności na każdym z etapów selekcji były rozstrzygane w drodze dyskusji; w przypadku braku konsensusu ostateczną decyzję podejmował autor korespondencyjny. Zgodność między recenzentami w zakresie włączenia pełnych tekstów została obliczona przy użyciu statystyki kappa Cohena. Chociaż do wstępnego usuwania duplikatów wykorzystano program EndNote 20, dodatkowe duplikaty, których oprogramowanie nie wykryło (np. ze względu na niewielkie różnice w tytułach, nazwiskach autorów lub skrótach czasopism), zostały zidentyfikowane i usunięte podczas ręcznego przeglądu tytułów i streszczeń przez dwóch niezależnych recenzentów.
Proces przesiewania i selekcji został podsumowany na schemacie blokowym PRISMA 2020 (Rysunek 1), który obejmuje liczbę zidentyfikowanych rekordów, usunięte duplikaty, przesiewne analizy rekordów, raporty przeznaczone do pobrania, raporty oceniane pod kątem kwalifikowalności oraz włączone badania, wraz ze szczegółowymi powodami wykluczenia na każdym etapie. W wyszukiwarce Google Scholar zastosowano uproszczony ciąg zapytania 'exercise AND depression AND (remission OR medication OR treatment)' z następującymi ustawieniami: sortowanie według trafności, wykluczenie cytowań oraz uwzględnienie patentów. Przesiewowo przeanalizowano pierwsze 50 rekordów. Wyszukiwanie przeprowadzono 15 listopada 2025 r. o godzinie 10:0 GMT, bez włączonej personalizacji, a następnie powtórzono 16 listopada 2025 r. w celu weryfikacji spójności.
Ekstrakcja danych
Dwóch recenzentów niezależnie wyodrębniło dane z każdego uwzględnionego badania, korzystając ze zestandaryzowanego, przetestowanego w wersji pilotażowej formularza ekstrakcji danych. Zebrano następujące informacje: nazwisko pierwszego autora, rok publikacji, kraj, w którym przeprowadzono badanie, wielkość próby (całkowita, grupa ćwicząca, grupa kontrolna), charakterystyka uczestników (wiek, metoda diagnozy lub przesiewowej depresji, wyjściowy stopień nasilenia depresji), szczegóły interwencji (rodzaj ćwiczeń, intensywność, częstotliwość, czas trwania, nadzór, format), opis grupy kontrolnej, miary wyników depresji (nazwa skali i zakres), moment oceny wyników oraz wyniki statystyczne (średnie, odchylenia standardowe, szacunki efektu, przedziały ufności, wartości p)20. Rozbieżności między recenzentami rozstrzygano w drodze dyskusji lub poprzez konsultację z autorem korespondencyjnym.
Elementy danych
W przypadku uzyskania odmiennych wyników w badaniu, dane gromadzono niezależnie na podstawie oceny wpływu interwencji w postaci ćwiczeń fizycznych na objawy depresji.
Ocena ryzyka błędu systematycznego
Dwoje autorów niezależnie oceniło ryzyko błędu systematycznego dla każdego włączonego badania, korzystając z narzędzia Cochrane RoB 2 dla randomizowanych badań kontrolowanych (wersja z sierpnia 2019 r.). Narzędzie RoB 2 ocenia pięć obszarów błędu wynikającego z procesu randomizacji, odchyleń od zamierzonych interwencji, brakujących danych o wynikach, pomiaru wyników oraz wyboru raportowanego wyniku. Dla każdego obszaru recenzenci przypisali ocenę „niskie ryzyko błędu” (low-risk of bias), „pewne wątpliwości” (some concerns) lub „wysokie ryzyko błędu” (high risk of bias). Następnie, zgodnie z algorytmami RoB 2, sformułowano ogólną ocenę ryzyka błędu dla każdego badania jako „niskie ryzyko” (wszystkie obszary z niskim ryzykiem), „pewne wątpliwości” (przynajmniej jeden obszar z pewnymi wątpliwościami, ale brak obszarów z wysokim ryzykiem) lub „wysokie ryzyko” (jakikolwiek obszar oceniony jako wysokie ryzyko lub wiele obszarów z pewnymi wątpliwościami). Rozbieżności między recenzentami rozstrzygano w drodze dyskusji, a w razie potrzeby arbiterem był autor korespondencyjny. Badania nie były wykluczane na podstawie ocen ryzyka błędu; zamiast tego zaplanowano analizy wrażliwości, aby ograniczyć metaanalizę do badań z ogólną oceną „niskie ryzyko” lub „niskie ryzyko + pewne wątpliwości”, w celu zbadania odporności wyników. Zgodność między sędziami dla ogólnych ocen ryzyka błędu obliczono przy użyciu współczynnika kappa Cohena21.
Analizy podgrup (analizy moderatorów)
W celu zbadania potencjalnych źródeł heterogeniczności oraz zidentyfikowania moderatorów skuteczności ćwiczeń przeprowadzono szereg analiz podgrup a priori, opartych na zmiennych istotnych klinicznie i metodologicznie. Analizy podgrup wykonano wyłącznie dla wyników, w których w każdej podgrupie znajdowało się co najmniej 10 badań (aby zapewnić odpowiednią moc statystyczną). Zdefiniowano wcześniej następujące zmienne podgrup.
Charakterystyka ćwiczeń
Tryb ćwiczeń: aerobowy (np. chodzenie, bieganie, jazda na rowerze) vs. oporowy (np. trening z ciężarami) vs. mieszany (aerobowy + oporowy) vs. inny (joga, taniec, ruch oparty na uważności).
Intensywność ćwiczeń określona jako niska (np. spokojny spacer, rozciąganie), umiarkowana (np. szybki spacer, umiarkowana jazda na rowerze) lub wysoka (np. bieganie, trening interwałowy o wysokiej intensywności) – zaklasyfikowana na podstawie standardowych kryteriów tętna lub subiektywnego odczucia wysiłku opisanych w każdym badaniu.
Czas trwania ćwiczeń (tygodnie) określono jako krótki (≤8 tygodni), średni (9–12 tygodni) lub długi (≥13 tygodni).
Nadzór jako nadzorowany (sesje ćwiczeń prowadzone przez instruktora, trenera lub terapeutę) w przeciwieństwie do nienadzorowanego (domowego lub samokierowanego, bez bezpośredniego nadzoru).
Sformatuj jako grupowe w przeciwieństwie do indywidualnego (domowego lub jeden na jeden).
Charakterystyka populacji
Diagnozę depresji sklasyfikowano jako zaburzenie depresyjne większe (MDD), zdiagnozowane za pomocą ustrukturyzowanego wywiadu klinicznego, takiego jak ustrukturyzowany wywiad kliniczny dla zaburzeń DSM (SCID) lub mini międzynarodowy wywiad neuropsychiatryczny (MINI), lub jako nasilone objawy depresyjne zidentyfikowane za pomocą walidowanych narzędzi przesiewowych w oparciu o ustalone wartości odcięcia, w tym skalę depresji Centrum Badań Epidemiologicznych (CES-D ≥ 16), inwentarz depresji Becka (BDI ≥ 14) lub kwestionariusz zdrowia pacjenta-9 (PHQ-9 ≥ 10). Grupę wiekową podzielono na dorosłych (18–59 lat) oraz osoby starsze (≥60 lat). Wartość wyjściową nasilenia depresji sklasyfikowano jako umiarkowaną lub ciężką zgodnie ze standardowymi wartościami odcięcia dla każdej skali oceny, gdzie depresję umiarkowaną zdefiniowano, na przykład, jako wynik w skali oceny depresji Hamiltona (HAM-D) w przedziale 14–18 lub wynik BDI w przedziale 14–19, a depresję ciężką jako wynik HAM-D ≥19 lub wynik BDI ≥20.
Charakterystyka projektu badania
Typ grupy kontrolnej zaklasyfikowano jako pasywny, obejmujący warunki kontrolne takie jak lista oczekujących, brak leczenia lub standardowa opieka, lub aktywny, obejmujący interwencje kontrolne w postaci uwagi, farmakoterapię lub psychoterapię. Rok publikacji podzielono na trzy okresy — przed 2010 r., 2010–2019 i 2020–2025 — aby zbadać potencjalne trendy czasowe w wynikach badań. Ryzyko błędu oceniono za pomocą narzędzia Cochrane Risk of Bias 2 (RoB 2) i zaklasyfikowano jako niskie, budzące pewne obawy lub wysokie. Wielkość próby zaklasyfikowano jako małą (łącznie n < 50 uczestników), średnią (50–9 uczestników) lub dużą (≥10 uczestników). Typ grupy kontrolnej podzielono na pasywny (lista oczekujących, standardowa opieka, brak interwencji) – grupy te nie zapewniają aktywnej alternatywy terapeutycznej, więc efekt ćwiczeń może być większy ze względu na efekty oczekiwań lub uwagi; oraz aktywny (kontrola uwagi, farmakoterapia, psychoterapia) – grupy te zapewniają interwencję niebędącą ćwiczeniami o wiarygodnej wartości terapeutycznej, co pozwala na bardziej konserwatywną ocenę specyficznego efektu ćwiczeń. Zgodnie z hipotezą, w wynikach potwierdzono istotnie większy efekt dla grup kontrolnych pasywnych w porównaniu z aktywnymi (p dla interakcji < 0,10; patrz Tabela 3).
Metody statystyczne dla analiz podgrup
W przypadku każdej analizy podgrup, badania w obrębie danej podgrupy zostały połączone przy użyciu modelu efektów losowych2. Do porównania podgrup zastosowano meta-regresję z efektami mieszanymi lub testy interakcji podgrup. W szczególności obliczono statystykę Q dla różnic między podgrupami i przedstawiono powiązaną z nią wartość p dla interakcji. Wartość p < 0,10 (zamiast 0,05) uznano za statystycznie istotną dla różnic między podgrupami, uwzględniając niską moc testów interakcji w metaanalizach. Nie przeprowadzono korekty na wielokrotne porównania (12 analiz podgrup); w związku z tym wyniki dla podgrup należy interpretować jako wstępne.
Analizy wrażliwości
Aby ocenić odporność głównych wyników, przeprowadzono następujące z góry określone analizy wrażliwości. Wybór modelu: ponownie przeanalizowano oba wyniki przy użyciu modelu efektów stałych (metoda odwrotności wariancji), aby porównać je z głównymi wynikami modelu efektów losowych; wykluczenie badań o wysokim ryzyku błędu systematycznego: wykluczono badania z ogólną oceną RoB 2 jako „wysokie ryzyko błędu systematycznego” (n = 17) i ponownie przeanalizowano pozostałe badania (niskie ryzyko + pewne wątpliwości); wykluczenie badań wieloramiennych: wykluczono badania, w których połączono wiele ramion ćwiczeń lub w których jedna grupa kontrolna była wykorzystywana do wielu porównań (n = 9) i ponownie przeanalizowano wyłącznie badania dwuramienne; wykluczenie badań z danymi imputowanymi lub założonymi: wykluczono badania, w których odchylenia standardowe dla wyników zmian były założone (r = 0,50), a nie raportowane (n = 6); wykluczenie badań sprzed wprowadzenia standardów CONSORT: wykluczono badania opublikowane przed 201 rokiem (n = 5), aby ocenić, czy poprawione standardy raportowania wpłynęły na szacunki efektu; analiza wpływu (leave-one-out): przeprowadzono metaanalizę z pominięciem jednego badania, sekwencyjnie usuwając każde badanie i przeliczając łączony efekt w celu zidentyfikowania badań wpływowych (zdefiniowanych jako takie, których usunięcie zmienia oszacowanie punktowe o >10% lub przesuwa przedział ufności poza pierwotne granice); analiza długoterminowego obserwowania: dla badań raportujących dane z obserwacji ≥6 miesięcy po interwencji (n = 6) dane te zostały połączone oddzielnie, aby zbadać trwałość efektów ćwiczeń, bez mieszania ich z danymi punktu końcowego (po interwencji). Wszystkie analizy wrażliwości przeprowadzono z użyciem modeli efektów losowych, chyba że wskazano inaczej. Wyniki analiz wrażliwości przedstawiono w sekcji Wyniki oraz w materiałach uzupełniających. Wszelkie odstępstwa od tych z góry określonych analiz zostały raportowane w sposób przejrzysty.
Ekstrakcja ciągłych danych wynikowych (wyniki skali depresji)
Dla każdego włączonego badania wyciągnięto średni wynik depresji oraz odchylenie standardowe (SD) dla grupy ćwiczącej i grupy kontrolnej w głównym punkcie końcowym (zdefiniowanym jako moment najbliższy zakończeniu interwencji w postaci ćwiczeń). Aby zmaksymalizować spójność między badaniami, zastosowano następującą hierarchiczną zasadę ekstrakcji: priorytetowo traktowano skorygowane szacunki efektu (np. różnice między grupami wyznaczone za pomocą ANCOVA z 95% CI lub SE), ponieważ uwzględniają one różnice wyjściowe i zazwyczaj dostarczają najmniej obciążony szacunek efektu interwencji. W przypadku braku raportowanych skorygowanych efektów wyciągnięto wyniki zmiany względem wartości początkowej (wynik po minus przed lub zmiana procentowa) oraz ich SD. Jeśli SD dla wyników zmiany nie zostało podane bezpośrednio, obliczono je przy użyciu następującego wzoru:
(1)
gdzie przyjęto konserwatywny współczynnik korelacji r = 0.5, chyba że w badaniu podano raportowaną lub możliwą do przypisania korelację. W przypadku braku dostępnych efektów skorygowanych lub wyników zmian, wyodrębniano wyłącznie wyniki po interwencji (bez korekty lub wyników różnicowych). Była to opcja najmniej preferowana, ponieważ nie uwzględnia ona potencjalnych różnic wyjściowych między grupami. W sytuacjach, gdy średnie i SD nie zostały podane bezpośrednio, dokonano konwersji innych statystyk (median i IQR, statystyk t, statystyk F, wartości p lub 95% CI) na średnie i SD, korzystając ze standardowych wzorów z podręcznika Cochrane (Rozdział 6). Badania raportujące jedynie mediany i zakresy włączono tylko wtedy, gdy liczebność próby była wystarczająco duża (n ≥ 25 na grupę), aby przyjąć przybliżoną normalność rozkładu, lub gdy autorzy badania dostarczyli średnie i SD na prośbę.
Wiele punktów czasowych: w przypadku badań zgłaszających wiele punktów czasowych obserwacji, wyodrębniono punkt czasowy najbliższy zakończeniu okresu aktywnej interwencji (zazwyczaj 8–12 tygodni). Dane z długoterminowej obserwacji (≥6 miesięcy po interwencji) wyodrębniono osobno do planowanej wtórnej analizy trwałości efektów, jednak dane te nie były łączone z danymi dotyczącymi głównego punktu końcowego. Wiele skal depresji: gdy w badaniu zgłoszono więcej niż jedną skalę depresji (np. zarówno HAM-D, jak i BDI), priorytetowo traktowano skale oceniane przez klinicystów (HAM-D, MADRS) nad skalami samoopisowymi (BDI, CES-D, PHQ-9, GDS), ponieważ skale oceniane przez klinicystów są uznawane za złoty standard w badaniach nad depresją. Jeśli zgłoszono obie, wyodrębniono główny miernik wyniku zdefiniowany przez autorów badania. Imputacja brakujących SD: w przypadku braku odchyleń standardowych, których nie można było obliczyć na podstawie dostępnych statystyk, zastosowano imputację, przyjmując średnie SD z badań o podobnej wielkości próby i z zastosowaniem tej samej skali depresji. Przeprowadzono analizy wrażliwości, wykluczając badania z imputowanymi SD, aby ocenić wpływ tej imputacji. Wszystkie decyzje dotyczące ekstrakcji danych zostały podjęte niezależnie przez dwóch autorów przy użyciu ustandaryzowanego formularza ekstrakcji danych. Rozbieżności rozstrzygano poprzez dyskusję lub przez autora korespondencyjnego. Odnotowano regułę ekstrakcji oraz wszelkie odstępstwa od niej.
Mierniki efektu
W przypadku wyniku dychotomicznego (remisja depresji) obliczono ilorazy szans (OR) wraz z 95-procentowymi przedziałami ufności (CI). W odniesieniu do wyniku ciągłego (wyniki nasilenia depresji) przewidywano, że w analizowanych badaniach zastosowano różne skale oceny depresji (np. skala oceny depresji Hamiltona, inwentarz depresji Becka, PHQ-9, CES-D, geriatryczna skala depresji). Mimo że w analizowanych badaniach wykorzystano różne skale oceny depresji, zgłoszono standaryzowaną różnicę średnich (SMD), ponieważ wszystkie skale posiadały opracowane równania konwersji na skalę HAM-D.
Synteza statystyczna
Dla wyniku dychotomicznego (remisja depresji) obliczono zbiorcze ilorazy szans (OR) wraz z 95% przedziałami ufności (CI). Dla wyniku ciągłego (wyniki nasilenia depresji) obliczono zbiorcze standaryzowane różnice średnich (SMD, g Hedgesa) wraz z 95% CI, ponieważ w uwzględnionych badaniach zastosowano różne skale oceny depresji. Ze względu na przewidywaną różnorodność kliniczną i metodologiczną między badaniami (zróżnicowanie protokołów ćwiczeń, populacji, warunków kontrolnych i miar wyników), do wszystkich analiz głównych a priori wybrano model efektów losowych (metoda DerSimonian i Lairda)2. Heterogeniczność określono za pomocą statystyki I2, interpretując ją zgodnie z kryteriami podręcznika Cochrane w następujący sposób: 0–40% (może nie być istotna), 30–60% (umiarkowana heterogeniczność), 50–90% (istotna heterogeniczność) oraz 75–10% (znaczna heterogeniczność). Oszacowano również wariancję między badaniami (tau2) i w stosownych przypadkach obliczono 95% przedziały predykcji dla efektów zbiorczych.
Analiza statystyczna i heterogeniczność
Z uwagi na oczekiwaną różnorodność kliniczną i metodologiczną w ramach włączonych badań (w tym różnice w rodzaju ćwiczeń, intensywności, czasie trwania, nadzorze, charakterystyce populacji, stopniu nasilenia depresji, warunkach kontrolnych oraz skalach pomiaru depresji), do wszystkich analiz głównych a priori wybrano model efektów losowych (metoda DerSimoniana i Lairda)2. Model efektów losowych zakłada, że rzeczywisty efekt różni się między badaniami i zapewnia bardziej konserwatywną ocenę z szerszymi przedziałami ufności, co jest odpowiednie w obliczu przewidywanej heterogeniczności.
Ocena heterogeniczności
Heterogeniczność określono za pomocą statystyki I2, która reprezentuje procent całkowitej zmienności między badaniami wynikający z rzeczywistej heterogeniczności, a nie z przypadku. Wartości I2 interpretowano w następujący sposób: 0–40% (może nie być istotna); 30–60% (umiarkowana heterogeniczność); 50–90% (znacząca heterogeniczność); 75–10% (wysoka heterogeniczność), zgodnie z kryteriami podręcznika Cochrane Handbook. Oprócz I2 oszacowano wariancję międzybadawczą (tau2) i, w stosownych przypadkach, podano 95% przedziały predykcyjne dla efektu zbiorczego.
Uzasadnienie wyboru modelu
Do żadnej z analiz głównych nie zastosowano modelu efektów stałych, ponieważ badanie nie są funkcjonalnie identyczne (różnią się protokołami ćwiczeń, populacjami i miernikami wyników), założenie o jednym prawdziwym efekcie wspólnym dla wszystkich badań jest nieprawdopodobne w badaniach nad wpływem ćwiczeń na depresję, a nawet gdy wartość I2 wydaje się niska (np. I2 = 48% dla remisji), sama różnorodność kliniczna uzasadnia podejście z efektami losowymi. Dla kompletności przeprowadzono również analizy wrażliwości z użyciem modelu efektów stałych, aby sprawdzić, czy wybór modelu wpłynął na wyciągnięte wnioski; wyniki te zostały przedstawione w materiałach uzupełniających.
Ocena błędu publikacji
Błąd publikacji oraz efekty małych badań oceniano zarówno za pomocą wizualnej analizy wykresów lejkowych, jak i testu regresji liniowej Eggera. Dla każdego punktu końcowego wygenerowano wykresy lejkowe, przedstawiające wielkość efektu (logarytm ilorazu szans dla remisji; standaryzowana różnica średnich dla wyników w skali depresji) w odniesieniu do błędu standardowego. W ocenie ilościowej test Eggera przeprowadza regresję standaryzowanej oceny efektu względem precyzji (odwrotności błędu standardowego). Statystycznie istotny punkt przecięcia (p < 0,10 w metaanalizach obejmujących ≥10 badań) sugeruje obecność asymetrii wykresu lejkowego, co może wskazywać na błąd publikacji lub efekty małych badań. Odwrotnie, p ≥ 0,10 wskazuje na brak statystycznych dowodów na występowanie takiego błędu.
Przedziały predykcji
W przypadku głównej metaanalizy z efektami losowymi obliczono 95% przedziały predykcji, które szacują zakres, w jakim mieściłby się rzeczywisty efekt w przyszłym badaniu. Szerokie przedziały predykcji wskazują, że efekt ćwiczeń może znacznie różnić się w zależności od warunków i populacji, co ma istotne znaczenie kliniczne. W przypadku badań wieloramiennych, w których wiele ram ćwiczeń połączono w jedną grupę, do obliczenia średnich łącznych, odchyleń standardowych oraz wielkości próby zastosowano wzory zawarte w podręczniku Cochrane Handbook. W badaniach, w których jedna grupa kontrolna była wykorzystywana w wielu porównaniach z ćwiczeniami, unikano podwójnego liczenia poprzez dzielenie wielkości próby grupy kontrolnej przez liczbę ram ćwiczeń podczas obliczania wielkości efektu. Wszystkie analizy przeprowadzono przy użyciu programu Reviewer Manager w wersji 5.3, a sposób postępowania z ramami wieloramiennymi został zweryfikowany niezależnie przez dwóch autorów.
Obliczanie liczby pacjentów koniecznych do leczenia (NNT)
Liczbę pacjentów wymaganych do leczenia (NNT) obliczono z połączonego ilorazu szans (OR) dla remisji, korzystając ze wzoru

gdzie CER (control event rate) stanowiło skumulowany wskaźnik remisji we wszystkich analizowanych badaniach. NNT obliczono również przy użyciu alternatywnej metody zaproponowanej przez Furukawę i Leuchta (201), która przelicza d Cohena (z SMD) na NNT, wykorzystując pole pod krzywą rozkładu normalnego. Wartości NNT obliczono dla analizy głównej (wszystkie badania), dla podgrup (badania o niskim ryzyku błędu systematycznego, badania dotyczące wyłącznie MDD, badania nad ćwiczeniami pod nadzorem) oraz w celu porównania z opublikowanymi wartościami NNT dla psychoterapii i leków przeciwdepresyjnych. Za klinicznie istotną uznano wartość NNT ≤ 10.
Konwersja na klinicznie znane mierniki (BDI oraz HAM-D)
Aby ułatwić interpretację kliniczną, połączona standaryzowana różnica średnich (SMD) dla wyników depresji została przeliczona na szacowane różnice średnich w dwóch powszechnie stosowanych skalach depresji: Inwentarzu Depresji Becka (BDI, zakres 0–63) oraz Skali Oceńania Depresji Hamiltona (HAM-D, zakres 0–52). Do przeliczenia wykorzystano następujący wzór

gdzie SDpooled, reference oznacza połączone odchylenie standardowe wyjściowe z badań, w których wykorzystano odpowiednio skalę BDI (n = 8 badań) lub HAM-D (n = 16 badań). Podejście to zakłada, że wielkość efektu SMD jest spójna dla poszczególnych skal. Obserwowane różnice średnich z tych badań zostały również przedstawione w ramach analizy wrażliwości. Zmiana o ≥3 punkty w skali BDI lub HAM-D została uznana za istotną klinicznie zgodnie z wytycznymi NICE5.
Porównanie z uznanymi metodami leczenia depresji
Aby osadzić efekt ćwiczeń w kontekście innych metod leczenia pierwszego rzutu w depresji, skumulowany wskaźnik NNT dla ćwiczeń porównano z opublikowanymi wartościami NNT z aktualnych, wysokiej jakości metaanaliz. Psychoterapia, jak zgłosili Cuijpers i wsp. (2020), ma NNT wynoszące 2,5 dla wszystkich grup wiekowych. Leki przeciwdepresyjne, według doniesień Cipriani i wsp. (2018), mają NNT wynoszące 4,3 dla leków w porównaniu z placebo. Porównania te mają charakter opisowy i nie opierają się na badaniach bezpośrednich (head-to-head). Dane dotyczące ćwiczeń nie były łączone z danymi dotyczącymi psychoterapii lub farmakoterapii, a różnice w populacjach badanych, definicjach wyników i punktach czasowych ograniczają możliwość bezpośredniego porównania.
Ocena błędu raportowania (błędu publikacji)
Wygenerowano wykresy lejkowe dla każdego punktu końcowego (remisja oraz wynik depresji) – wykresy te przedstawiają wielkość efektu (logarytm ilorazu szans dla remisji; standaryzowana różnica średnich dla wyników depresji) w odniesieniu do błędu standardowego. Asymetria wykresu lejkowego może wskazywać na błąd publikacji, efekty małych badań lub heterogeniczność. W celu formalnej oceny asymetrii wykresu lejkowego przeprowadzono test regresji Eggera, znany również jako liniowy test regresji Eggera. Test Eggera analizuje regresję standaryzowanej oceny efektu względem jej precyzji (odwrotności błędu standardowego). Statystycznie istotny punkt przecięcia (p < 0,10 dla testu Eggera, zgodnie z zaleceniami podręcznika Cochrane dla punktów końcowych obejmujących mniej niż 10 badań, lub p < 0,05 dla większych metaanaliz) sugeruje występowanie efektów małych badań lub błędu publikacji. Odwrotnie, p ≥ 0,05 (lub p ≥ 0,10) wskazuje na brak statystycznie istotnych dowodów na taki błąd. Należy zachować ostrożność w interpretacji, ponieważ asymetria wykresu lejkowego nie jest tożsama z błędem publikacji; może ona również wynikać z rzeczywistej heterogeniczności, różnic w jakości badań lub przypadku, a wyniki należy interpretować odpowiednio23.
Stopień pewności dowodów (ocena GRADE)
Ogólną pewność dowodów dla każdego punktu końcowego (remisja oraz wynik w skali depresji) oceniono z wykorzystaniem systemu Grading of Recommendations Assessment, Development and Evaluation (GRADE). Dwóch autorów niezależnie oceniło pewność dowodów, a rozbieżności rozstrzygano w drodze dyskusji. Podejście GRADE uwzględnia pięć domen, które mogą obniżyć pewność: ryzyko błędu (na podstawie ocen RoB 2, w tym proporcji badań o niskim, pewnym lub wysokim ryzyku), niespójność (na podstawie statystyki I2, przedziałów predykcyjnych i nakładania się przedziałów ufności), pośredniość (na podstawie różnic między populacjami badawczymi, interwencjami, grupami kontrolnymi i punktami końcowymi w odniesieniu do pytania badawczego), nieprecyzyjność (na podstawie optymalnej wielkości informacji oraz tego, czy 95% przedziały ufności przekraczają klinicznie istotne progi, takie jak iloraz szans 1,0 dla remisji lub standaryzowana różnica średnich 0,2 dla małego efektu w wynikach skali depresji) oraz błąd publikacji (na podstawie asymetrii wykresu lejkowego i testu Eggera). Pewność oceniano w następujący sposób: Wysoka (bardzo mało prawdopodobne jest, aby dalsze badania zmieniły poziom ufności co do szacunku efektu); Umiarkowana (dalsze badania prawdopodobnie wywrą istotny wpływ na poziom ufności i mogą zmienić szacunek); Niska (bardzo prawdopodobne jest, aby dalsze badania wywarły istotny wpływ na poziom ufności i zmieniły szacunek); oraz Bardzo niska (każdy szacunek efektu jest bardzo niepewny). Oceny pewności przeprowadzono oddzielnie dla remisji (punkt końcowy dychotomiczny) oraz wyników w skali depresji (punkt końcowy ciągły).