Lo studio è stato approvato dal Comitato Etico per la Biomedicina dell'Università di Jishou (Numero di approvazione: JSDX-2023-0034). Tutti i partecipanti erano studenti universitari a tempo pieno di età compresa tra 18 e 25 anni e hanno fornito un consenso informato scritto prima della partecipazione. È stata ottenuta l'autorizzazione dalle università partecipanti. I dati del sondaggio sono stati utilizzati esclusivamente per ricerca accademica e sono stati analizzati in forma anonima. Lo studio è stato condotto nel rispetto della Dichiarazione di Helsinki e dei requisiti istituzionali locali.
Reclutamento dei partecipanti e progettazione del campionamento
È stato condotto un sondaggio trasversale tra studenti universitari iscritti a corsi di laurea triennale di tre università complete del sud della Cina tra settembre e dicembre 2025, utilizzando un disegno campionario casuale stratificato a grappoli in più stadi. Nella prima fase, tre università complete sono state selezionate da un elenco di università pubbliche complete del sud della Cina mediante campionamento casuale stratificato, con strati definiti in base alla regione amministrativa provinciale al fine di garantire una rappresentatività geografica. Nella seconda fase, all'interno di ciascuna università selezionata, sono state scelte casualmente due o tre classi intere per ogni anno accademico (dal primo all'ultimo anno) come unità campionarie a grappolo, utilizzando numeri casuali generati al computer, ottenendo un totale di 28 classi distribuite tra le tre università. Tutti gli studenti appartenenti alle classi selezionate sono stati invitati a partecipare, come mostrato nella Figura 1. I criteri di inclusione erano i seguenti: (1) status di studente universitario iscritto a tempo pieno; (2) età compresa tra 18 e 25 anni; e (3) partecipazione volontaria con consenso informato scritto. I criteri di esclusione erano: (1) malattia fisica grave recente o controindicazioni all'esercizio fisico, come fratture o gravi malattie cardiache; (2) assunzione attuale di farmaci psichiatrici, ad esempio per il trattamento di gravi depressioni o schizofrenia; e (3) tempo di compilazione del questionario inferiore a 180 s o evidenti incoerenze logiche, come la selezione della stessa opzione per tutti gli elementi. Sono stati distribuiti complessivamente 1.500 questionari. Degli 80 questionari non validi esclusi durante la pulizia dei dati, 24 sono stati esclusi perché compilati in meno di 180 s, 20 per incoerenze logiche (ad esempio risposte identiche a tutti gli elementi o risposte contraddittorie), 18 per età al di fuori dell'intervallo 18–25 anni, 11 per dati mancanti eccessivi (più del 20% degli elementi), e 7 per malattia fisica grave autosegnalata o uso attuale di farmaci psichiatrici. Sono rimasti 1.420 questionari validi per l'analisi, con un tasso di risposta efficace del 94,67%.

Figura 1: Progettazione dello studio e flusso analitico. Schema riassuntivo del reclutamento e della selezione dei partecipanti (1.500 questionari distribuiti; 80 esclusi; N = 1.420 inclusi), valutazione mediante PSQI e IPAQ-SF, analisi dei profili latenti (confronto di modelli da uno a cinque profili utilizzando AIC, BIC, aBIC, entropia, LMRT e BLRT), assegnazione a quattro profili del sonno, test chi-quadrato bivariati, regressione logistica multinomiale con il genere come covariata e analisi di sensibilità (aggiustamento completo per le covariate ed errori standard robusti per cluster). Cliccare qui per visualizzare una versione ingrandita di questa figura.
Valutazione demografica e dello stile di vita
È stato utilizzato un questionario auto-progettato per raccogliere informazioni sulle caratteristiche demografiche e sociologiche, inclusi genere (maschio/femmina), età (in anni), anno di corso (primo, secondo, terzo, quarto anno), luogo di provenienza (urbano/rurale), indice di massa corporea (BMI), storia del fumo, storia dell'assunzione di alcol e tempo medio giornaliero trascorso davanti ai dispositivi schermo. Il BMI è stato calcolato come il peso corporeo in chilogrammi diviso per il quadrato dell'altezza in metri (BMI = peso [kg] / altezza [m]²)2). L'altezza e il peso auto-dichiarati sono stati verificati per valori non plausibili prima dell'analisi; i casi con altezza <140 cm oppure >200 cm, peso <30 kg o >150 kg, o BMI <12 kg/m2 o >40 kg/m2 sono stati segnalati e verificati rispetto ai record originali, con esclusione degli outlier confermati. Lo stato di fumatore è stato classificato in mai, occasionale (fumo alcuni giorni ma non quotidianamente) o regolare (fumo quotidiano), in base al comportamento autodichiarato nel mese precedente. Il consumo di alcol è stato classificato in mai, occasionale (consumo meno di una volta alla settimana) o regolare (consumo almeno una volta alla settimana), riferito al mese precedente. Il tempo medio giornaliero trascorso davanti a schermi è stato valutato chiedendo ai partecipanti di riportare il numero totale di ore al giorno dedicate a smartphone, computer, tablet e televisione nel mese precedente, registrato come variabile continua. Il questionario è stato sottoposto a un test pilota su 30 studenti universitari prima della raccolta formale dei dati, al fine di valutare la chiarezza degli item e il tempo di compilazione; sono state apportate lievi modifiche al testo in base ai feedback ottenuti nel test pilota. Studi precedenti hanno suggerito che queste variabili potrebbero confondere l'associazione tra sonno e attività fisica30,31. Tutti gli strumenti, software e materiali utilizzati in questo studio sono elencati nel Tabella dei materiali.
Valutazione della qualità del sonno
Il PSQI, sviluppato da Buysse et al.9, è stato utilizzato per valutare la qualità del sonno dei partecipanti nell'ultimo mese. La scala comprende 19 voci autovalutate che compongono sette componenti: qualità soggettiva del sonno, latenza dell'addormentamento, durata del sonno, efficienza abituale del sonno, disturbi del sonno, uso di farmaci ipnotici e disfunzione diurna. Ogni componente viene punteggiato da 0 (nessuna difficoltà) a 3 (difficoltà grave). I punteggi di questi sette componenti continui sono stati utilizzati come indicatori osservati per l'LPA. Il PSQI ha dimostrato buona affidabilità e validità in popolazioni di studenti universitari in Cina e a livello internazionale.30Nello studio attuale, l'alpha di Cronbach per la scala è stato 0,83.
Valutazione dell'attività fisica
L'IPAQ-SF, sviluppato da Craig et al.32, è stato utilizzato per valutare l'attività fisica dei partecipanti nei sette giorni precedenti. Il questionario registra la frequenza (giorni/settimana) e la durata (minuti/giorno) dell'attività fisica intensa, dell'attività fisica di intensità moderata e della camminata. I valori dell'equivalente metabolico del compito (MET) sono stati calcolati secondo il protocollo di punteggio IPAQ: MET-min/settimana per la camminata = 3,3 × durata × frequenza; MET-min/settimana per intensità moderata = 4,0 × durata × frequenza; e MET-min/settimana per intensità vigorosa = 8,0 × durata × frequenza. I valori MET per i tre tipi di attività sono stati sommati per ottenere il livello totale di attività fisica. I partecipanti sono stati classificati in alta PA, moderata PA e bassa PA in base ai criteri di punteggio IPAQ-SF validati in studenti universitari33. L'alta PA è definita come attività vigorosa per almeno tre giorni con un livello totale di attività fisica di almeno 1.500 MET-min/settimana, oppure qualsiasi combinazione di camminata, attività di intensità moderata o intensa per sette o più giorni con un livello totale di attività fisica di almeno 3.000 MET-min/settimana. La PA moderata è definita come attività vigorosa per almeno tre giorni per almeno 20 minuti al giorno, oppure attività di intensità moderata e/o camminata per almeno cinque giorni per almeno 30 minuti al giorno, oppure qualsiasi combinazione di attività che raggiunga almeno 600 MET-min/settimana. La bassa PA è definita come mancato raggiungimento dei criteri per la PA moderata o alta. La pulizia dei dati ha seguito il protocollo ufficiale di punteggio IPAQ: sono stati esclusi i partecipanti che riportavano più di 16 h di attività totale al giorno; le durate giornaliere per categoria sono state troncate a 180 min prima del calcolo dei MET-minuti; e i casi con risposte logicamente incoerenti sono stati rimossi durante la pulizia dei dati.
Analisi statistica
Innanzitutto, eseguire l'analisi di profilo latente (LPA) utilizzando un software di modellazione dei profili latenti (vedere Tabella dei materiali)34. Dei 1.420 partecipanti, 33 (2,3%) presentavano dati parzialmente mancanti nel PSQI, con 69 risposte mancanti su 9.940 risposte relative ai singoli componenti (0,69%; cioè 1.420 partecipanti × sette punteggi dei componenti del PSQI). Per gestire questi valori mancanti, è stato utilizzato il metodo della massima verosimiglianza con informazione completa (FIML), assumendo una mancanza casuale dei dati. Sono stati stimati sequenzialmente modelli con da uno a cinque profili. L'adattamento del modello è stato valutato mediante il criterio di informazione di Akaike (AIC), il criterio di informazione bayesiano (BIC), il BIC corretto per la dimensione del campione (aBIC), l'entropia, il test del rapporto di verosimiglianza di Lo-Mendell-Rubin (LMRT) e il test del rapporto di verosimiglianza con bootstrap (BLRT). Valori più bassi di AIC, BIC e aBIC indicano un adattamento relativo migliore, mentre valori di entropia più vicini a 1 indicano una maggiore accuratezza nella classificazione. Il LMRT e il BLRT sono stati utilizzati per determinare se un modello con k profili si adattasse significativamente meglio di un modello con k−1 profili. La soluzione ottimale di profili è stata selezionata considerando congiuntamente gli indici statistici di adattamento, l'entropia, la dimensione delle classi (nessun profilo doveva contenere meno del 5% del campione) e l'interpretabilità sostanziale, seguendo linee guida consolidate per l'analisi di profilo latente11,35. Dopo aver selezionato la soluzione ottimale di profili, l'appartenenza alla classe più probabile di ciascun individuo è stata esportata in un software statistico (vedere Tabella dei materiali) per le analisi successive. Sono stati utilizzati test chi-quadrato di Pearson per esaminare le differenze nelle caratteristiche demografiche e nei livelli di attività fisica tra i diversi profili di sonno. I test chi-quadrato bivariati hanno mostrato che solo il genere e il livello di attività fisica differivano in modo significativo tra i quattro profili di sonno (entrambi p < 0,001), mentre età, classe scolastica, BMI, fumo, consumo di alcol e tempo trascorso davanti agli schermi non differivano (tutti p > 0,05; vedere Tabella Supplementare 1). Pertanto, solo il genere è stato incluso come covariata nel modello di regressione logistica multinomiale, in linea con il principio di parsimonia. Un'analisi di sensibilità che includeva tutte le covariate raccolte ha prodotto lo stesso andamento e significatività dei risultati (Tabella Supplementare 2). È stata eseguita una regressione logistica multinomiale con il livello di attività fisica come variabile dipendente (PA bassa = 0, PA moderata = 1, PA elevata = 2) e il profilo di sonno come variabile indipendente principale. Sono stati calcolati gli odds ratio (OR) e gli intervalli di confidenza al 95% (IC). Tutti i test erano a due code e la significatività statistica è stata fissata a p < 0,05. Poiché i partecipanti sono stati reclutati attraverso classi intere annidate all'interno di università, è stata condotta un'analisi di sensibilità utilizzando il modulo per campioni complessi del software statistico (vedere Tabella dei materiali), specificando università e classe come variabili di raggruppamento, al fine di ottenere errori standard aggiustati per il disegno sperimentale. I risultati sono risultati coerenti con quelli ottenuti dalla regressione logistica multinomiale standard (Tabella Supplementare 3). Il set di dati grezzo anonimizzato, con etichette delle variabili in inglese e un codebook, è fornito come File Supplementare 1.