Rekrutacja uczestników, ukończenie badania i charakterystyka wyjściowa w punkcie T0
Rycina 1A przedstawia przepływ uczestników od etapu rekrutacji i weryfikacji kwalifikowalności, poprzez ocenę wyjściową, aż po randomizację. Rycina 1B przedstawia harmonogram ocen w punkcie wyjściowym (T0), w trakcie 6-tygodniowego okresu interwencji, podczas oceny po interwencji przeprowadzonej w ciągu 48 h po 6. tygodniu (T1) oraz podczas oceny kontrolnej po 4 tygodniach (T2). Dziewięćdziesięciu studentów zrandomizowano w stosunku 1:1:1 do grupy PME (n = 30), aktywnej grupy kontrolnej analizy muzycznej (n = 30) lub grupy kontrolnej na liście oczekujących (n = 30). Charakterystyka demograficzna w punkcie wyjściowym, wcześniejsze szkolenie muzyczne, status wystąpienia głównego stresora w ciągu ostatnich 7 dni, wyjściowy poziom kontaktu z muzyką, frekwencja oraz przestrzeganie protokołu zostały podsumowane w Tabeli 1. Wyniki pomiarów wyjściowych podsumowano w Tabeli 2. Zbiór danych na poziomie uczestników, wykorzystany do sporządzenia tabel i rycin, znajduje się w Pliku z danymi uzupełniającymi 1.
Trzy grupy badawcze były w dużym stopniu porównywalne w punkcie wyjściowym. Średni wiek był podobny we wszystkich grupach (PME, 20,34 ± 1,50 lat; aktywna grupa kontrolna, 20,45 ± 1,70 lat; grupa kontrolna na liście oczekujących, 20,48 ± 1,54 lat). Kobiety stanowiły 18/30 (60,0%) grupy PME, 19/30 (63,3%) aktywnej grupy kontrolnej oraz 24/30 (80,0%) grupy kontrolnej na liście oczekujących. Brak lub minimalne wcześniejsze wykształcenie muzyczne zgłosiło odpowiednio 18/30 (60,0%), 18/30 (60,0%) i 20/30 (66,7%) uczestników (Tabela 1). Wyjściowe miary wyników również nie wykazały wyraźnych różnic przed interwencją między grupami (Tabela 2). Wyniki afektu pozytywnego w skali PANAS wynosiły 29,23 ± 6,13 w grupie PME, 28,17 ± 5,75 w aktywnej grupie kontrolnej oraz 29,13 ± 6,58 w grupie kontrolnej na liście oczekujących. Wyniki afektu negatywnego w skali PANAS wynosiły odpowiednio 22,33 ± 4,33, 20,00 ± 6,58 oraz 21,50 ± 5,12. Średnie wyniki PSS-10 wynosiły 18,87 ± 5,02, 19,23 ± 5,13 i 17,87 ± 4,07, natomiast wyniki SCS-R wynosiły odpowiednio 77,67 ± 11,40, 74,27 ± 11,47 i 77,53 ± 13,76.
Ekspozycja na interwencję, wierność protokołowi i bezpieczeństwo w tygodniach 1–6
Rycina 2 przedstawia ustandaryzowaną strukturę 90-minutowej sesji PME, obejmującą rozgrzewkę uziemiającą, krąg rytmiczny, wyznaczone miejsce przerwy, improwizację, współtworzenie oraz elementy ustrukturyzowanej refleksji. Ekspozycja na interwencję była wysoka w obu aktywnych ramionach badania. Średnia frekwencja wyniosła 0,87 ± 0,09 w grupie PME i 0,84 ± 0,09 w grupie aktywnej kontrolnej. Przestrzeganie protokołu, zdefiniowane jako obecność na pięciu lub więcej z sześciu sesji, osiągnęło 25/30 (83,3%) uczestników w grupie PME i 23/30 (76,7%) uczestników w grupie aktywnej kontrolnej.
Tabela 3 podsumowuje kompletne wskaźniki wdrożenia, wierności, odchyleń, działań korygujących oraz bezpieczeństwa w ciągu 6 tygodni. Wierność pozostała na konsekwentnie wysokim poziomie przez cały czas trwania interwencji. Wyniki wierności dla PME mieściły się w przedziale od 8,92 ± 0,58 do 9,33 ± 0,44, przy średniej 6-tygodniowej wynoszącej 9,20 ± 0,49. Wyniki wierności dla grupy aktywnej kontrolnej mieściły się w przedziale od 8,86 ± 0,61 do 9,09 ± 0,66, z odpowiadającą im średnią 8,99 ± 0,56. Wyznaczony punkt przerwania został wykorzystany w 2 z 6 sesji PME (33,3%), co wskazuje na selektywne stosowanie procedury resetowania, a nie rutynowe przerywanie prowadzenia sesji. Odchylenia we wdrożeniu występowały rzadko. Żadna sesja PME nie osiągnęła zdefiniowanego progu niskiej wierności. Jedna sesja w grupie aktywnej kontrolnej (1/6; 16,7%) została zaklasyfikowana jako niska wierność z powodu naruszenia granic grupy aktywnej kontrolnej. Odchylenia związane z facylitatorem wystąpiły raz w grupie PME i dwukrotnie w grupie aktywnej kontrolnej, natomiast odchylenia związane z uczestnikami wystąpiły odpowiednio cztery i pięć razy. Działania korygujące zostały wdrożone raz w grupie PME i dwukrotnie w grupie aktywnej kontrolnej. Ponowne przeszkolenie facylitatora było wymagane raz w grupie aktywnej kontrolnej i nie było wymagane w grupie PME. Nie odnotowano żadnych zdarzeń niepożądanych związanych z udziałem w badaniu i żaden uczestnik nie zrezygnował z powodu wystąpienia stanu stresowego. Krótkotrwały dyskomfort, który nie wymagał przerwania udziału, wystąpił u 1/30 (3,3%) uczestników zarówno w grupie PME, jak i w grupie aktywnej kontrolnej.
Wyniki po interwencji i obserwacji w punktach czasowych T1 i T2
Tabela 4 podsumowuje wyniki mierników w punktach czasowych T1 i T2, wartości zmiany względem stanu wyjściowego (ΔT1 = T1 − T0 oraz ΔT2 = T2 − T0) oraz szacunki osłabienia zmian między T1 a T2 (ΔT2 − ΔT1). Rysunek 3 ilustruje wyniki zmiany względem stanu wyjściowego dla afektu pozytywnego w skali PANAS (Rysunek 3A), afektu negatywnego w skali PANAS (Rysunek 3B), PSS-10 (Rysunek 3C) oraz SCS-R (Rysunek 3D). W punkcie T1 grupa PME wykazała największą poprawę we wszystkich czterech miernikach wyników. Afekt pozytywny w skali PANAS wzrósł, natomiast wyniki afektu negatywnego w skali PANAS i PSS-10 spadły w większym stopniu niż w którejkolwiek z grup porównawczych. Wyniki SCS-R również wzrosły bardziej w grupie PME niż w grupach aktywnej kontroli lub kontroli listy oczekujących (Tabela 4; Rysunek 3). Ten korzystny wzorzec utrzymał się podczas 4-tygodniowej obserwacji (T2). Grupa PME nadal wykazywała większą poprawę względem stanu wyjściowego we wszystkich miernikach niż którakolwiek z grup porównawczych, podczas gdy grupa aktywnej kontroli wykazała mniejszą poprawę, a grupa kontroli listy oczekujących wykazała minimalne lub niespójne zmiany. Analiza osłabienia zmian między T1 a T2 wskazała na utrzymanie efektów lub dalszą poprawę w grupie PME. Afekt pozytywny, odczuwany stres oraz poczucie powiązań społecznych nadal zmieniały się w korzystnym kierunku, podczas gdy redukcja afektu negatywnego została w dużej mierze utrzymana. Dla porównania, szacunki osłabienia były mniejsze w grupie aktywnej kontroli i wykazały niewielką, niespójną poprawę w grupie kontroli listy oczekujących (Tabela 4).
Szacunki oparte na modelu
Tabela 5 podsumowuje zdefiniowane kontrasty w liniowym modelu mieszanym, porównujące grupę PME z grupą aktywnej kontroli oraz grupą listy oczekujących w punktach czasowych T1 i T2. Modele obejmowały efekty stałe dla grupy badania, czasu oceny, interakcji grupa × czas, stanu wystąpienia głównego stresora w ciągu ostatnich 7 dni w punkcie wyjściowym oraz wszelkich istotnie niezrównoważonych współzmiennych wyjściowych, wraz z losowym punktem przecięcia na poziomie uczestnika. Szacowane różnice średnich marginalnych, 95% CI, wartości P oraz współczynnik Hedges’ g przedstawiono w Tabeli 5. Wyniki oparte na modelu były zgodne z analizami opisowymi. Zarówno w T1, jak i T2, grupa PME wykazała wyższe wyniki afektu pozytywnego w PANAS oraz wyższe wyniki w SCS-R, a także niższe wyniki w PSS-10 niż obie grupy porównawcze. W przypadku afektu negatywnego w PANAS, grupa PME różniła się znacząco od grupy listy oczekujących w obu punktach czasowych, podczas gdy porównanie z grupą aktywnej kontroli nie było statystycznie istotne. Standaryzowane wielkości efektu zazwyczaj wskazywały na korzyść grupy PME i zachowały ten sam kierunek w obu ocenach po interwencji, co potwierdza spójność obserwowanych różnic grupowych we wszystkich mierzonych efektach. Analizy wrażliwości z zastosowaniem imputacji wielokrotnej dla brakujących danych oraz grupowania na poziomie sesji nie zmieniły kierunku zdefiniowanych kontrastów grupa × czas. Eksploracyjne analizy podgrup zostały zinterpretowane po korekcie wskaźnika fałszywych odkryć (FDR).
Materiały uzupełniające towarzyszące temu artykułowi obejmują skrypty dla moderatorów PME (Supplementary File 1), bibliotekę muzyczną dla grupy kontrolnej oraz przewodnik dyskusyjny (Supplementary File 2), kwestionariusz przesiewowy online (Supplementary File 3), skrypt R do randomizacji (Supplementary File 4), szablon zgłaszania zdarzeń niepożądanych i reakcji na stres (Supplementary File 5), Supplementary Figures 1–2 oraz Supplementary Data File 1 zawierający zestaw danych na poziomie uczestników, który posłużył do wygenerowania tabel i rycin przedstawionych w niniejszym badaniu.
Plik uzupełniający 1. Scenariusze dla prowadzących PME. Scenariusze dla prowadzących dla wszystkich sześciu sesji partycypacyjnej edukacji muzycznej (PME), zawierające harmonogram sesji, wskazówki dla prowadzących, przejścia, pytania do kierowanej refleksji, sformułowania dla punktów przerwy oraz punkty kontrolne wierności interwencji.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2. Biblioteka muzyczna aktywnej grupy kontrolnej i przewodnik dyskusyjny. Tygodniowe motywy słuchania, wyselekcjonowane fragmenty muzyczne, instrukcje robienia notatek, ustrukturyzowane pytania do dyskusji oraz zakazane aktywności związane z aktywnym tworzeniem muzyki dla warunku aktywnej grupy kontrolnej.Proszę kliknąć tutaj, aby pobrać ten plik.
Plik uzupełniający 3. Kwestionariusz przesiewowy online. Kwestionariusz przesiewowy dla uczestników online, zawierający kryteria kwalifikowalności, dostępność w celu udziału w badaniu, wcześniejszą edukację muzyczną, kryteria wykluczenia oraz badanie aktywności muzycznej poza sesjami.Aby pobrać ten plik, kliknij tutaj.
Plik uzupełniający 4. Skrypt R do randomizacji. Skrypt R do randomizacji blokowej 1:1:1 z rozwarstwieniem według płci, przy użyciu wielkości bloku równej 6 i ziarna losowania 20260302.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 5. Szablon zdarzeń niepożądanych i reakcji na stres. Formularz monitorowania bezpieczeństwa dokumentujący stopień nasilenia zdarzenia, podjęte działania, oferowane skierowania, status powrotu uczestnika oraz ocenę związku ze studium.Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca 1. Progi wierności i kodowanie odchyleń. Progi wierności interwencji, kryteria czasowe, kody odchyleń, definicje niskiej wierności oraz odpowiadające im działania korygujące.Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca 2. Słownik zmiennych i zasady punktacji.Nazwy zmiennych, konwencje kodowania, procedury punktacji kwestionariusza, zasady punktacji odwróconej, sposób postępowania z brakującymi danymi oraz pochodne zmienne zmiany wyniku.Aby pobrać ten plik, kliknij tutaj.
Plik z danymi uzupełniającymi 1. Zanonimizowany zbiór danych na poziomie uczestników. Zanonimizowany zbiór danych na poziomie uczestników w formacie długim, wykorzystany do wszystkich analiz statystycznych oraz do wygenerowania tabel i rycin przedstawionych w niniejszym badaniu.Kliknij tutaj, aby pobrać ten plik.