Данное исследование было одобрено Комитетом по биомедицинской этике Университета Цзишоу (номер одобрения JSDX-2023-034). Все участники были студентами бакалавриата очной формы обучения в возрасте 18–25 лет и предоставили письменное информированное согласие перед участием. Было получено разрешение от участвующих университетов. Данные опроса использовались только для академических исследований и анализировались анонимно. Исследование проводилось в соответствии с Хельсинкской декларацией и требованиями местных учреждений.
Набор участников и дизайн выборки
С сентября по декабрь 2025 года среди студентов бакалавриата трех многопрофильных университетов на юге Китая был проведен одномоментный опрос с использованием многоступенчатого стратифицированного кластерного рандомизированного отбора. На первом этапе из списка государственных многопрофильных университетов южного Китая методом стратифицированного случайного отбора были выбраны три университета, причем страты были определены по провинциальным административным регионам для обеспечения географического представительства. На втором этапе в каждом выбранном университете из каждого курса (с первого по четвертый) с помощью компьютерного генератора случайных чисел были случайным образом выбраны две или три полноценные группы в качестве единиц кластерной выборки, что в итоге составило 28 групп в трех университетах. Все студенты выбранных групп были приглашены к участию, как показано на Рисунке 1. Критерии включения были следующими: (1) статус студента очного отделения бакалавриата; (2) возраст от 18 до 25 лет; (3) добровольное участие с письменным информированным согласием. Критерии исключения: (1) недавно перенесенные тяжелые соматические заболевания или противопоказания к физическим нагрузкам, такие как переломы или тяжелые заболевания сердца; (2) текущий прием психиатрических препаратов, например, лечение тяжелой депрессии или шизофрении; (3) время заполнения анкеты менее 180 s или явные логические противоречия, такие как выбор одного и того же варианта ответа для всех пунктов. Всего было распространено 1 50 анкет. Из 80 недействительных анкет, исключенных в процессе очистки данных, 24 были исключены из-за времени заполнения менее 180 s, 20 — из-за логических противоречий (например, идентичные ответы во всех пунктах или взаимоисключающие ответы), 18 — из-за возраста вне диапазона 18–25 лет, 1 — из-за чрезмерного количества пропущенных данных (более 20% пунктов) и 7 — из-за самостоятельно задекларированных тяжелых соматических заболеваний или приема психиатрических препаратов. Таким образом, для анализа осталось 1 420 действительных анкет, что составило эффективный уровень отклика 94.67%.

Рисунок 1Дизайн исследования и аналитический рабочий процесс. Схематический обзор процесса набора и скрининга участников (распространено 1 50 анкет; исключено 80; N = 1 420 включено), оценки с помощью PSQI и IPAQ-SF, анализа латентных профилей (сравнение моделей от одного до пяти профилей с использованием AIC, BIC, aBIC, энтропии, LMRT и BLRT), распределения по четырем профилям сна, двумерных критериев хи-квадрат, мультиномиальной логистической регрессии с полом в качестве ковариаты и анализа чувствительности (полная поправка на ковариаты и кластерно-устойчивые стандартные ошибки). Пожалуйста, нажмите здесь, чтобы просмотреть эту фигуру в увеличенном размере.
Оценка демографических данных и образа жизни
Для получения демографических и социологических характеристик использовался специально разработанный опросник, включающий пол (мужской/женский), возраст (в годах), курс обучения (первый/второй/третий/четвертый), место происхождения (город/сельская местность), индекс массы тела (ИМТ), историю курения, историю употребления алкоголя и среднее ежедневное экранное время. ИМТ рассчитывали как массу тела в килограммах, деленную на квадрат роста в метрах (ИМТ = вес [kg] / рост [m]2). Самостоятельно указанные данные о росте и весе проверялись на предмет недостоверных значений перед началом анализа; случаи с ростом <140 cm или >20 cm, весом <30 kg или >150 kg, либо ИМТ <12 kg/m2 или >40 kg/m2 помечались и сверялись с исходными записями, а подтвержденные выбросы исключались. Статус курения разделяли на категории: никогда не курил, курит эпизодически (в некоторые дни, но не ежедневно) или курит регулярно (ежедневно), на основании данных о курении за предыдущий месяц. Потребление алкоголя за предыдущий месяц разделяли на категории: никогда не употреблял, употреблял эпизодически (менее одного раза в неделю) или употреблял регулярно (минимум один раз в неделю). Среднее ежедневное экранное время оценивали путем запроса у участников общего количества часов в день, затрачиваемых на использование смартфонов, компьютеров, планшетов и телевизора за предыдущий месяц; данные записывались как непрерывная переменная. Перед основным сбором данных опросник был апробирован на группе из 30 студентов бакалавриата для оценки ясности формулировок и времени заполнения; на основе полученных отзывов были внесены незначительные правки в формулировки. Предыдущие исследования показали, что эти переменные могут выступать в качестве конфаундеров при анализе связи между сном и физической активностью30,31. Все приборы, программное обеспечение и материалы, использованные в данном исследовании, перечислены в Таблице материалов.
Оценка качества сна
Для оценки качества сна участников за последний месяц использовался опросник PSQI, разработанный Buysse et al.9. Шкала включает 19 пунктов самооценки, которые формируют семь компонентов: субъективное качество сна, латентный период засыпания, продолжительность сна, привычная эффективность сна, нарушения сна, использование гипнотических препаратов и дневная дисфункция. Каждый компонент оценивается от 0 (отсутствие трудностей) до 3 (выраженные трудности). Баллы по этим семи непрерывным компонентам использовались в качестве наблюдаемых показателей для LPA. PSQI продемонстрировал хорошую надежность и валидность среди студентов колледжей как в Китае, так и на международном уровне30. В настоящем исследовании коэффициент альфа Кронбаха для данной шкалы составил 0,83.
Оценка физической активности
Для оценки физической активности участников за последние семь дней использовался опросник IPAQ-SF, разработанный Craig et al.32. В анкете фиксируются частота (дней/неделя) и продолжительность (минут/день) интенсивной физической активности, физической активности умеренной интенсивности и ходьбы. Значения метаболического эквивалента задачи (MET) рассчитывались в соответствии с протоколом подсчета баллов IPAQ: MET-мин/неделя при ходьбе = 3.3 × продолжительность × частота; MET-мин/неделя при умеренной интенсивности = 4.0 × продолжительность × частота; MET-мин/неделя при высокой интенсивности = 8.0 × продолжительность × частота. Значения MET для трех типов активности суммировались для получения общего уровня физической активности. Участники были разделены на группы с высоким, умеренным и низким уровнем ФА в соответствии с критериями оценки IPAQ-SF, валидированными для студентов колледжей3. Высокий уровень ФА определялся как интенсивная активность не менее трех дней с общим уровнем физической активности не менее 1 50 MET-min/week или любая комбинация ходьбы, деятельности умеренной или высокой интенсивности в течение семи и более дней с общим уровнем физической активности не менее 3 0 MET-min/week. Умеренный уровень ФА определялся как интенсивная активность не менее трех дней по крайней мере по 20 минут в день, или активность умеренной интенсивности и/или ходьба не менее пяти дней по крайней мере по 30 минут в день, или любая комбинация видов деятельности, достигающая уровня не менее 60 MET-min/week. Низкий уровень ФА определялся как несоответствие критериям умеренного или высокого уровня ФА. Очистка данных проводилась согласно официальному протоколу подсчета баллов IPAQ: исключались участники, сообщившие о более чем 16 h общей активности в день; ежедневная продолжительность по каждой категории была ограничена 180 min перед расчетом MET-минут; случаи с логически противоречивыми ответами были удалены в процессе очистки данных.
Статистический анализ
Сначала проведите LPA с использованием программного обеспечения для моделирования латентных профилей (см. Таблицу материалов)34. Из 1 420 участников у 3 (2,3%) частично отсутствовали данные по PSQI, при этом отсутствовали 69 из 9 940 ответов на уровне компонентов (0,69%; т. е. 1 420 участников × семь показателей компонентов PSQI). Для обработки этих пропущенных значений при допущении о случайном характере пропусков использовался метод максимального правдоподобия с полной информацией (FIML). Последовательно оценивались модели от одного до пяти профилей. Соответствие модели оценивалось с помощью информационного критерия Акаике (AIC), байесовского информационного критерия (BIC), BIC с поправкой на размер выборки (aBIC), энтропии, критерия отношения правдоподобия Ло-Менделла-Рубина (LMRT) и бутстреп-теста отношения правдоподобия (BLRT). Более низкие значения AIC, BIC и aBIC указывают на лучшее относительное соответствие, тогда как значения энтропии, близкие к 1, указывают на более высокую точность классификации. Тесты LMRT и BLRT использовались для определения того, соответствует ли модель с k профилями значимо лучше, чем модель с k − 1 профилями. Оптимальное решение по профилям было выбрано путем совместного учета статистических индексов соответствия, энтропии, размера классов (ни один профиль не должен содержать менее 5% выборки) и содержательной интерпретируемости в соответствии с установленными рекомендациями по анализу латентных профилей1,35. После выбора оптимального решения по профилям наиболее вероятная принадлежность индивидов к классу была экспортирована в программное обеспечение для статистического анализа (см. Таблицу материалов) для последующего анализа. Для изучения различий в демографических характеристиках и уровнях физической активности между профилями сна использовались критерии хи-квадрат Пирсона. Двумерные тесты хи-квадрат показали, что только пол и уровень физической активности значимо различались между четырьмя профилями сна (в обоих случаях p < 0,01), тогда как возраст, курс обучения, BMI, курение, употребление алкоголя и экранное время не имели значимых различий (все p > 0,05; см. Дополнительную таблицу 1). Таким образом, в модель мультиномиальной логистической регрессии в качестве ковариаты был включен только пол, что соответствует принципу парсимонии. Анализ чувствительности с поправкой на все собранные ковариаты дал такую же закономерность и значимость результатов (Дополнительная таблица 2). Мультиномиальная логистическая регрессия проводилась с уровнем физической активности в качестве зависимой переменной (низкая ФА = 0, умеренная ФА = 1, высокая ФА = 2) и профилем сна в качестве основной независимой переменной. Были рассчитаны отношения шансов (ORs) и 95% доверительные интервалы (CIs). Все тесты были двусторонними, статистическая значимость была установлена на уровне p < 0,05. Поскольку участники набирались из сформированных групп внутри университетов, мы провели анализ чувствительности с использованием модуля для сложных выборок статистического ПО (см. Таблицу материалов), указав университет и группу в качестве переменных кластеризации для получения стандартных ошибок с поправкой на дизайн исследования. Результаты совпали с результатами стандартной мультиномиальной логистической регрессии (Дополнительная таблица 3). Обезличенный набор исходных данных с метками переменных на английском языке и кодовая книга представлены в Дополнительном файле 1.