Deze studie is goedgekeurd door de Biomedische Ethiekcommissie van Jishou University (Goedkeuringsnummer JSDX-2023-034). Alle deelnemers waren voltijdse bachelorstudenten in de leeftijd van 18–25 jaar en gaven voorafgaand aan deelname schriftelijke geïnformeerde toestemming. Toestemming is verkregen van de deelnemende universiteiten. De enquêtegegevens werden uitsluitend gebruikt voor academisch onderzoek en werden anoniem geanalyseerd. De studie is uitgevoerd in overeenstemming met de Verklaring van Helsinki en de lokale institutionele vereisten.
Werving van deelnemers en steekproefontwerp
Tussen september en december 2025 werd een cross-sectioneel onderzoek uitgevoerd onder bachelorstudenten van drie comprehensive universities in Zuid-China, waarbij gebruik werd gemaakt van een multistage stratified cluster random sampling design. In de eerste fase werden drie comprehensive universities geselecteerd uit een lijst van openbare comprehensive universities in Zuid-China via gestratificeerde steekproeftrekking, waarbij de strata werden gedefinieerd op basis van de provinciale administratieve regio om geografische representativiteit te waarborgen. In de tweede fase werden binnen elke geselecteerde universiteit willekeurig twee tot drie intacte klassen per jaarlaag (eerste- tot vierdejaars) gekozen als clustersteekproefeenheden met behulp van door de computer gegenereerde willekeurige getallen, wat resulteerde in totaal 28 klassen over de drie universiteiten. Alle studenten in de geselecteerde klassen werden uitgenodigd om deel te nemen, zoals weergegeven in Figuur 1. De inclusiecriteria waren als volgt: (1) status als fulltime bachelorstudent; (2) een leeftijd tussen 18 en 25 jaar; en (3) vrijwillige deelname met schriftelijke geïnformeerde toestemming. De exclusiecriteria waren: (1) recente ernstige fysieke aandoeningen of contra-indicaties voor lichaamsbeweging, zoals fracturen of ernstige hartaandoeningen; (2) huidig gebruik van psychiatrische medicatie, zoals behandeling voor ernstige depressie of schizofrenie; en (3) een invultijd van de vragenlijst van minder dan 180 s of duidelijke logische inconsistenties, zoals het selecteren van dezelfde optie voor alle items. In totaal werden 1.50 vragenlijsten verspreid. Van de 80 ongeldige vragenlijsten die tijdens de datacleaning werden uitgesloten, werden er 24 uitgesloten vanwege een invultijd onder de 180 s, 20 vanwege logische inconsistenties (bijv. identieke antwoorden op alle items of tegenstrijdige antwoorden), 18 vanwege een leeftijd buiten de range van 18–25 jaar, 1 vanwege overmatig veel ontbrekende gegevens (meer dan 20% van de items) en 7 vanwege zelfgerapporteerde ernstige fysieke aandoeningen of huidig gebruik van psychiatrische medicatie. Dit liet 1.420 geldige vragenlijsten over voor analyse, met een effectieve responsgraad van 94,67%.

Figuur 1Studieopzet en analytische workflow. Schematisch overzicht van de werving en screening van deelnemers (1.50 verspreide vragenlijsten; 80 uitgesloten; N = 1.420 geïncludeerd), beoordeling met de PSQI en IPAQ-SF, latente profielanalyse (vergelijking van modellen met één tot vijf profielen met behulp van AIC, BIC, aBIC, entropie, LMRT en BLRT), toewijzing aan vier slaapprofielen, bivariante chi-kwadraattoetsen, multinomiale logistische regressie met geslacht als covariabele en sensitiviteitsanalyses (volledige correctie voor covariabelen en cluster-robuuste standaardfouten). Klik hier om een grotere versie van deze figuur te bekijken.
Beoordeling van demografie en levensstijl
Een zelfontworpen vragenlijst werd gebruikt om demografische en sociologische kenmerken te verkrijgen, waaronder geslacht (man/vrouw), leeftijd (in jaren), studiejaar (eerste/tweede/derde/vierde jaar), plaats van herkomst (stedelijk/landelijk), bodymassindex (BMI), rookgeschiedenis, drinkgeschiedenis en gemiddelde dagelijkse schermtijd. De BMI werd berekend als het lichaamsgewicht in kilogram gedeeld door het kwadraat van de lengte in meters (BMI = gewicht [kg] / lengte [m]2). Zelfgerapporteerde lengte en gewicht werden vóór de analyse gescreend op onwaarschijnlijke waarden; gevallen met een lengte <140 cm of >20 cm, een gewicht <30 kg of >150 kg, of een BMI <12 kg/m2 of >40 kg/m2 werden gemarkeerd en geverifieerd aan de hand van de originele dossiers, waarbij bevestigde uitschieters werden uitgesloten. De rookstatus werd gecategoriseerd als nooit, incidenteel (roken op sommige dagen, maar niet dagelijks) of regelmatig (dagelijks roken), gebaseerd op het zelfgerapporteerde rookgedrag gedurende de afgelopen maand. Alcoholconsumptie werd over de afgelopen maand gecategoriseerd als nooit, incidenteel (minder dan één keer per week drinken) of regelmatig (minstens één keer per week drinken). De gemiddelde dagelijkse schermtijd werd beoordeeld door deelnemers te vragen het totale aantal uren per dag dat aan smartphones, computers, tablets en televisie werd besteed over de afgelopen maand te rapporteren, wat werd geregistreerd als een continue variabele. De vragenlijst is voorafgaand aan de formele gegevensverzameling voorlopig getest onder 30 bachelorstudenten om de duidelijkheid van de items en de voltooiingstijd te beoordelen; op basis van de feedback uit deze pilottest zijn kleine aanpassingen in de bewoordingen doorgevoerd. Eerdere studies hebben gesuggereerd dat deze variabelen de associatie tussen slaap en fysieke activiteit kunnen beïnvloeden30,31. Alle instrumenten, software en materialen die in deze studie zijn gebruikt, staan vermeld in de Tabel met Materialen.
Beoordeling van de slaapkwaliteit
De PSQI, ontwikkeld door Buysse et al.9, werd gebruikt om de slaapkwaliteit van de deelnemers over de afgelopen maand te beoordelen. De schaal bevat 19 zelfgerapporteerde items die zeven componenten vormen: subjectieve slaapkwaliteit, slaaplatentie, slaapduur, gebruikelijke slaapefficiëntie, slaapstoornissen, gebruik van hypnotica en overdag disfunctie. Elke component wordt gescoord van 0 (geen problemen) tot 3 (ernstige problemen). De scores van deze zeven continue componenten werden gebruikt als geobserveerde indicatoren voor LPA. De PSQI heeft een goede betrouwbaarheid en validiteit aangetoond onder populaties van studenten in China en internationaal30. In de huidige studie was de Cronbachs alpha voor de schaal 0,83.
Beoordeling van fysieke activiteit
De IPAQ-SF, ontwikkeld door Craig et al.32, werd gebruikt om de fysieke activiteit van de deelnemers over de voorafgaande zeven dagen te beoordelen. De vragenlijst registreert de frequentie (dagen/week) en duur (minuten/dag) van intensieve fysieke activiteit, matig-intensieve fysieke activiteit en wandelen. De metabole equivalenten van taken (MET)-waarden werden berekend volgens het IPAQ-scoreprotocol: wandelen MET-min/week = 3,3 × duur × frequentie; matig-intensief MET-min/week = 4,0 × duur × frequentie; en intensief MET-min/week = 8,0 × duur × frequentie. De MET-waarden voor de drie activiteitstypen werden opgeteld om de totale fysieke activiteit te bepalen. Deelnemers werden ingedeeld in hoge PA, matige PA en lage PA volgens de IPAQ-SF-scorecriteria die zijn gevalideerd bij universiteitsstudenten3. Hoge PA werd gedefinieerd als intensieve activiteit op ten minste drie dagen met een totaal fysiek activiteitsniveau van ten minste 1.500 MET-min/week, of elke combinatie van wandelen, matig-intensieve of intensieve activiteiten op zeven of meer dagen met een totaal fysiek activiteitsniveau van ten minste 3.0 MET-min/week. Matige PA werd gedefinieerd als intensieve activiteit op ten minste drie dagen gedurende ten minste 20 minuten per dag, of matig-intensieve activiteit en/of wandelen op ten minste vijf dagen gedurende ten minste 30 minuten per dag, of elke combinatie van activiteiten die ten minste 60 MET-min/week bereikte. Lage PA werd gedefinieerd als het niet voldoen aan de criteria voor matige of hoge PA. De datareiniging volgde het officiële IPAQ-scoreprotocol: deelnemers die meer dan 16 h totale activiteit per dag rapporteerden, werden uitgesloten; de dagelijkse duur per categorie werd afgekapt op 180 min vóór de MET-minuutberekening; en gevallen met logisch inconsistente antwoorden werden tijdens de datareiniging verwijderd.
Statistische analyse
Voer eerst een LPA uit met behulp van software voor latent profielmodelleren (zie Tabel met Materialen)34. Van de 1.420 deelnemers hadden er 3 (2,3%) gedeeltelijk ontbrekende gegevens in de PSQI, waarbij 69 van de 9.940 antwoorden op componentniveau ontbraken (0,69%; d.w.z. 1.420 deelnemers × zeven PSQI-componentscores). De full-information maximum likelihood (FIML)-methode werd gebruikt om deze ontbrekende waarden te behandelen onder de aanname van missing-at-random. Modellen met één tot vijf profielen werden sequentieel geschat. De modelpassing werd geëvalueerd met behulp van het Akaike-informatiecriterium (AIC), het Bayesian-informatiecriterium (BIC), het op steekproefomvang aangepaste BIC (aBIC), entropie, de Lo-Mendell-Rubin likelihood ratio test (LMRT) en de bootstrap likelihood ratio test (BLRT). Lagere AIC-, BIC- en aBIC-waarden duiden op een betere relatieve passing, terwijl entropiewaarden dichter bij 1 duiden op een hogere classificatienauwkeurigheid. De LMRT en BLRT werden gebruikt om te bepalen of een model met k profielen significant beter paste dan een model met k − 1 profielen. De optimale profieloplossing werd geselecteerd door gezamenlijk te kijken naar statistische passing-indices, entropie, klassengrootte (geen enkel profiel mocht minder dan 5% van de steekproef bevatten) en inhoudelijke interpreteerbaarheid, volgens vastgestelde richtlijnen voor latente profielanalyse11,35. Nadat de optimale profieloplossing was geselecteerd, werd het meest waarschijnlijke klaslidmaatschap van individuen geëxporteerd naar statistische analyse-software (zie Tabel met Materialen) voor verdere analyses. Pearson chi-kwadraattoetsen werden gebruikt om verschillen in demografische kenmerken en fysieke activiteitsniveaus tussen de slaapprofielen te onderzoeken. Bivariate chi-kwadraattoetsen toonden aan dat alleen geslacht en fysiek activiteitsniveau significant verschilden tussen de vier slaapprofielen (beide p < 0,01), terwijl leeftijd, leerjaar, BMI, roken, drinken en schermtijd dat niet deden (allemaal p > 0,05; zie Aanvullende Tabel 1). Daarom werd alleen geslacht als covariaat opgenomen in het multinomiaal logistisch regressiemodel, in overeenstemming met het principe van parsimonie. Een sensitiviteitsanalyse waarbij werd gecorrigeerd voor alle verzamelde covariaten leverde hetzelfde patroon en dezelfde significantie van resultaten op (Aanvullende Tabel 2). Multinomiale logistische regressie werd uitgevoerd met het fysieke activiteitsniveau als afhankelijke variabele (lage PA = 0, matige PA = 1, hoge PA = 2) en het slaapprofiel als de belangrijkste onafhankelijke variabele. Odds ratio's (ORs) en 95% betrouwbaarheidsintervallen (CIs) werden berekend. Alle toetsen waren tweezijdig en de statistische significantie werd vastgesteld op p < 0,05. Omdat deelnemers werden geworven via intacte klassen genest binnen universiteiten, hebben we een sensitiviteitsanalyse uitgevoerd met de complex-samples module van de statistische software (zie Tabel met Materialen), waarbij universiteit en klas als clusteringsvariabelen zijn gespecificeerd om ontwerp-gecorrigeerde standaardfouten te verkrijgen. De resultaten waren consistent met die van de standaard multinomiale logistische regressie (Aanvullende Tabel 3). De gedeïdentificeerde ruwe dataset met Engelse variabelenamen en een codeboek is beschikbaar als Aanvullend Bestand 1.