Inschrijving van deelnemers, voltooiing en basiskenmerken op T0
Figuur 1A illustreert de stroom van deelnemers vanaf de inschrijving en screening op geschiktheid tot de nulmeting en randomisatie. Figuur 1B illustreert het schema van de beoordelingen tijdens de nulmeting (T0), de interventieperiode van 6 weken, de post-interventiebeoordeling uitgevoerd binnen 48 u na week 6 (T1) en de follow-upbeoordeling na 4 weken (T2). Negentig studenten werden in een verhouding van 1:1:1 gerandomiseerd naar PME (n = 30), een actieve controlegroep voor muziekappreciatie (n = 30) of een wachtlijstcontrolegroep (n = 30). Demografische kenmerken bij aanvang, eerdere muziekopleiding, status van belangrijke stressoren in de afgelopen 7 dagen, blootstelling aan muziek bij aanvang, aanwezigheid en naleving van het protocol zijn samengevat in Tabel 1. De uitkomstmaten bij aanvang zijn samengevat in Tabel 2. De dataset op deelnemerniveau die is gebruikt om de tabellen en figuren te genereren, is opgenomen in Aanvullend Databestand 1.
De drie studiegroepen waren bij aanvang grotendeels vergelijkbaar. De gemiddelde leeftijd was vergelijkbaar tussen de groepen (PME, 20,34 ± 1,50 jaar; actieve controlegroep, 20,45 ± 1,70 jaar; wachtlijstcontrolegroep, 20,48 ± 1,54 jaar). Vrouwelijke deelnemers maakten 18/30 (60,0%) van de PME-groep uit, 19/30 (63,3%) van de actieve controlegroep en 24/30 (80,0%) van de wachtlijstcontrolegroep. Geen of minimale eerdere muziektraining werd respectievelijk gerapporteerd door 18/30 (60,0%), 18/30 (60,0%) en 20/30 (66,7%) van de deelnemers (Tabel 1). Uitkomstmaten bij aanvang vertoonden eveneens geen opvallende verschillen tussen de groepen vóór de interventie (Tabel 2). De PANAS-scores voor positieve affect waren 29,23 ± 6,13 in de PME-groep, 28,17 ± 5,75 in de actieve controlegroep en 29,13 ± 6,58 in de wachtlijstcontrolegroep. De PANAS-scores voor negatieve affect waren respectievelijk 22,33 ± 4,33, 20,00 ± 6,58 en 21,50 ± 5,12. De gemiddelde PSS-10-scores waren 18,87 ± 5,02, 19,23 ± 5,13 en 17,87 ± 4,07, terwijl de SCS-R-scores respectievelijk 77,67 ± 11,40, 74,27 ± 11,47 en 77,53 ± 13,76 waren.
Interventieblootstelling, getrouwheid en veiligheid tijdens week 1–6
Figuur 2 illustreert de gestandaardiseerde structuur van een PME-sessie van 90 min, inclusief de aardingsopwarming, ritmecirkel, vastgesteld pauzemoment, improvisatie, co-creatie en de gestructureerde reflectiecomponenten. De blootstelling aan de interventie was hoog in beide actieve studiearmen. De gemiddelde aanwezigheid was 0,87 ± 0,09 in de PME-groep en 0,84 ± 0,09 in de actieve controlegroep. Protocolnaleving, gedefinieerd als aanwezigheid bij vijf of meer van de zes sessies, werd bereikt door 25/30 (83,3%) deelnemers in de PME-groep en 23/30 (76,7%) deelnemers in de actieve controlegroep.
Tabel 3 vat de volledige indicatoren voor implementatie, getrouwheid, afwijking, correctieve acties en veiligheid over de periode van 6 weken samen. De getrouwheid bleef gedurende de gehele interventie consistent hoog. De getrouwheidsscores voor PME varieerden van 8,92 ± 0,58 tot 9,33 ± 0,44, met een gemiddelde over 6 weken van 9,20 ± 0,49. De getrouwheidsscores voor de actieve controlegroep varieerden van 8,86 ± 0,61 tot 9,09 ± 0,66, met een overeenkomstig gemiddelde van 8,99 ± 0,56. Het aangewezen pauzemoment werd gebruikt in 2 van de 6 PME-sessies (33,3%), wat wijst op een selectieve implementatie van de resetprocedure in plaats van een routinematige onderbreking van de sessieuitvoering. Implementatieafwijkingen kwamen zelden voor. Geen enkele PME-sessie bereikte de vooraf gedefinieerde drempelwaarde voor lage getrouwheid. Eén sessie in de actieve controlegroep (1/6; 16,7%) werd geclassificeerd als lage getrouwheid vanwege een schending van de grenzen van de actieve controle. Facilitator-gerelateerde afwijkingen kwamen één keer voor in de PME-groep en twee keer in de actieve controlegroep, terwijl deelnemer-gerelateerde afwijkingen respectievelijk vier en vijf keer voorkwamen. Correctieve acties werden één keer geïmplementeerd in de PME-groep en twee keer in de actieve controlegroep. Hertraining van de facilitator was één keer nodig in de actieve controlegroep en was niet nodig in de PME-groep. Er werden geen bijwerkingen gemeld die gerelateerd waren aan deelname aan het onderzoek, en geen enkele deelnemer stopte vanwege een stressvolle gebeurtenis. Kortstondig ongemak dat geen stopzetting van de deelname vereiste, trad op bij 1/30 (3,3%) van de deelnemers in zowel de PME- als de actieve controlegroep.
Resultaten na interventie en follow-up op T1 en T2
Tabel 4 vat de uitkomstscores samen op T1 en T2, de waarden van de verandering ten opzichte van de nullijn (ΔT1 = T1 − T0 en ΔT2 = T2 − T0) en de schattingen van de attenuatie van T1 naar T2 (ΔT2 − ΔT1). Figuur 3 illustreert de resultaten van de verandering ten opzichte van de nullijn voor PANAS positieve affect (Figuur 3A), PANAS negatieve affect (Figuur 3B), PSS-10 (Figuur 3C) en SCS-R (Figuur 3D). Op T1 vertoonde de PME-groep de grootste verbeteringen over alle vier de uitkomstmaten. De PANAS positieve affect nam toe, terwijl de PANAS negatieve affect en PSS-10 scores sterker afnamen dan in beide vergelijkingsgroepen. De SCS-R scores namen eveneens meer toe in de PME-groep dan in de actieve controlegroep of de wachtlijstcontrolegroep (Tabel 4; Figuur 3). Dit gunstige patroon bleef behouden bij de follow-up na 4 weken (T2). De PME-groep bleef grotere verbeteringen ten opzichte van de nullijn vertonen dan beide vergelijkingsgroepen voor alle uitkomstmaten, terwijl de actieve controlegroep kleinere verbeteringen vertoonde en de wachtlijstcontrolegroep minimale of inconsistente veranderingen liet zien. De attenuatieanalyse van T1 naar T2 wees op behoud of voortdurende verbetering in de PME-groep. Positieve affect, ervaren stress en sociale verbondenheid bleven in de gunstige richting veranderen, terwijl reducties in negatieve affect grotendeels werden behouden. Ter vergelijking waren de attenuatieschattingen kleiner in de actieve controlegroep en was er weinig consistente verbetering zichtbaar in de wachtlijstcontrolegroep (Tabel 4).
Modelgebaseerde schattingen
Tabel 5 vat de vooraf gedefinieerde contrasten van het lineaire mixed-effectsmodel samen, waarbij de PME-groep werd vergeleken met de actieve controlegroep en de wachtlijstcontrolegroep op T1 en T2. De modellen bevatten vaste effecten voor studiegroep, beoordelingstijdstip, de interactie tussen groep × tijd, de status van belangrijke stressoren in de afgelopen 7 dagen bij baseline, en eventuele significant ongebalanceerde covariabelen bij baseline, samen met een random intercept op deelnemerniveau. De geschatte marginale gemiddelde verschillen, 95% BI's, P-waarden en Hedges' g worden gepresenteerd in Tabel 5. De op het model gebaseerde resultaten waren consistent met de descriptieve analyses. Op zowel T1 als T2 vertoonde de PME-groep hogere scores voor PANAS positieve affect en SCS-R en lagere PSS-10 scores dan beide vergelijkingsgroepen. Voor PANAS negatieve affect verschilde de PME-groep significant van de wachtlijstcontrolegroep op beide tijdstippen, terwijl de vergelijking met de actieve controlegroep niet statistisch significant was. Gestandaardiseerde effectgroottes waren over het algemeen gunstig voor de PME-groep en bleven in dezelfde richting bij beide beoordelingen na de interventie, wat de consistentie van de waargenomen groepsverschillen over de uitkomstmaten ondersteunt. Sensitiviteitsanalyses waarin multiple imputatie voor missende gegevens en clustering op sessieniveau waren opgenomen, veranderden de richting van de vooraf gedefinieerde groep × tijd-contrasten niet. Exploratieve subgroepanalyses werden geïnterpreteerd na correctie voor de false discovery rate.
De aanvullende materialen bij dit artikel omvatten de PME-facilitatorscripts (Supplementary File 1), de muziekbibliotheek en discussiegids voor de actieve controlelaag (Supplementary File 2), de online screeningvragenlijst (Supplementary File 3), het R-script voor randomisatie (Supplementary File 4), het sjabloon voor bijwerkingen en stressreacties (Supplementary File 5), Supplementary Figures 1–2, en Supplementary Data File 1 met de dataset op deelnemerniveau die is gebruikt om de tabellen en figuren in deze studie te genereren.
Aanvullend bestand 1. PME-facilitatorscripts. Facilitatorscripts voor alle zes de sessies participatieve muziekeducatie (PME), inclusief timing van de sessies, aanwijzingen voor de facilitator, overgangen, begeleidende reflectievragen, bewoordingen van pausemomenten en controlepunten voor de getrouwheid van de interventie.Klik hier om dit bestand te downloaden.
Aanvullend bestand 2. Muziekbibliotheek voor actieve controlegroep en discussiegids. Wekelijkse luisterthema's, gecureerde muziekfragmenten, instructies voor het maken van aantekeningen, gestructureerde discussievragen en verboden actieve muziekactiviteiten voor de actieve controleconditie.Klik hier om dit bestand te downloaden.
Aanvullend bestand 3. Online screeningsvragenlijst. Online screeningsvragenlijst voor deelnemers, inclusief geschiktheidscriteria, beschikbaarheid voor deelname aan de studie, eerdere muziekopleiding, uitsluitingscriteria en screening van muziekactiviteiten buiten de sessies.Klik hier om dit bestand te downloaden.
Aanvullend bestand 4. R-script voor randomisatie.R-script voor sekse-gestratificeerde 1:1:1 blokrandomisatie met een blokgrootte van 6 en een random seed van 20260302.Klik hier om dit bestand te downloaden.
Aanvullend bestand 5. Sjabloon voor bijwerkingen en distressreacties. Veiligheidsmonitoringsformulier waarin de ernst van het event, de ondernomen acties, aangeboden verwijzingen, de status van de terugkeer van de deelnemer en de beoordeling van de relatie tot het onderzoek worden gedocumenteerd.Klik hier om dit bestand te downloaden.
Aanvullende tabel 1. Getrouwheidsdrempels en codering van afwijkingen. Drempels voor de getrouwheid van de interventie, timingcriteria, afwijkingscodes, definities van lage getrouwheid en de bijbehorende correctieve acties.Klik hier om dit bestand te downloaden.
Aanvullende Tabel 2. Variabelenwoordenlijst en scoringsregels. Variabelenamen, coderingsconventies, scoringsprocedures van de vragenlijst, regels voor omgekeerde scoren, verwerking van ontbrekende gegevens en afgeleide veranderingsscore-variabelen.Klik hier om dit bestand te downloaden.
Aanvullend gegevensbestand 1. Geanonimiseerde dataset op deelnemersniveau. Geanonimiseerde dataset op deelnemersniveau in lang formaat, gebruikt voor alle statistische analyses en voor het genereren van de tabellen en figuren die in deze studie worden gerapporteerd.Klik hier om dit bestand te downloaden.