מחקר זה אושר על ידי ועדת האתיקה הביו-רפואית של אוניברסיטת Jishou (מספר אישור JSDX-2023-0034). כל המשתתפים היו סטודנטים לתואר ראשון במשרה מלאה בגילאי 18–25 שנים, אשר סיפקו הסכמה מדעת בכתב לפני השתתפותם. התקבל אישור מהאוניברסיטאות המשתתפות. נתוני הסקר שימשו למחקר אקדמי בלבד ונולאו באופן אנונימי. המחקר בוצע בהתאם להצהרת הלסינקי ולדרישות המוסדיות המקומיות.
גיוס משתתפים ותכנון דגימה
סקר חתך נערך בקרב סטודנטים לתואר ראשון משלוש אוניברסיטאות מקיפות בדרום סין בין ספטמבר לדצמבר 2025, תוך שימוש במערך דגימה מקראית של אשכולות מרובדית רב-שלבית. בשלב הראשון, נבחרו שלוש אוניברסיטאות מקיפות מתוך רשימה של אוניברסיטאות מקיפות ציבוריות בדרום סין באמצעות דגימה אקראית מרובדת, כאשר הרבדים הוגדרו לפי אזור מנהלי מחוזה כדי להבטיח ייצוג גאוגרפי. בשלב השני, בתוך כל אוניברסיטה שנבחרה, נבחרו באופן אקראי שתיים עד שלוש מחלקות שלמות מכל שנת לימוד (משנה א' ועד שנה ד') כיחידות דגימה של אשכולות באמצעות מספרים אקראים שהופקו במחשב, מה שהניב סך כולל של 28 מחלקות בשלוש האוניברסיטאות. כל הסטודנטים במחלקות שנבחרו הוזמנו להשתתף, כפי שמוצג ב-איור 1. קריטריוני ההכללה היו כדלקמן: (1) סטטוס סטודנט לתואר ראשון במשרה מלאה; (2) גיל בין 18 ל-25 שנים; ו-(3) השתתפות רצונית עם הסכמה מדעת בכתב. קריטריוני ההוצאה היו: (1) מחלה גופנית חמורה לאחרונה או התוויות נגד לפעילות גופנית, כגון שברים או מחלת לב חמורה; (2) שימוש נוכחי בתרופות פסיכיאטריות, כגון טיפול בדיכאון חמור או בסכיזופרניה; ו-(3) זמן מילוי של השאלון הפחות מ-180 s או סתירות לוגיות בולטות, כגון בחירה באותה אפשרות עבור כל הפריטים. סך כולל של 1,500 שאלונים הופצו. מתוך 80 השאלונים הבלתי תקפים שהוצאו במהלך ניקוי הנתונים, 24 הוצאו בשל זמן מילוי נמוך מ-180 s, 20 בשל סתירות לוגיות (למשל, תגובות זהות בכל הפריטים או תשובות סותרות), 18 בשל גיל שאינו בטווח של 18–25 שנים, 11 בשל חוסר בנתונים באופן מוגזם (יותר מ-20% מהפריטים), ו-7 בשל דיווח עצמי על מחלה גופנית חמורה או שימוש נוכחי בתרופות פסיכיאטריות. נותרו 1,420 שאלונים תקפים לניתוח, עם שיעור תגובה אפקטיבי של 94.67%.

איור 1: מערך המחקר וזרימת העבודה האנליטית. סכימה כללית של גיוס וסינון משתתפים (1,500 שאלונים הופצו; 80 הוצאו; N = 1,420 נכללו), הערכה באמצעות PSQI ו-IPAQ-SF, ניתוח פרופיל חבוי (השוואה בין מודלים של פרופיל אחד עד חמישה באמצעות AIC, BIC, aBIC, אנטרופיה, LMRT ו-BLRT), שיוך לארבעה פרופילי שינה, מבחני chi-square דו-משתניים, רגרסיה לוגיסטית רב-נומית עם מגדר כמשתנה תואם (covariate), ואנליזות רגישות (התאמה מלאה של משתנים תואמים ושגיאות תקן חסונות אשכולות). אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של איור זה.
הערכת דמוגרפיה ואורח חיים
שאלון שתוכנן באופן עצמאי שימש להשגת מאפיינים דמוגרפיים וסוציולוגיים, כולל מגדר (זכר/נקבה), גיל (בשנים), שנה בלימודים (שנה א'/שנה ב'/שנה ג'/שנה ד'), מקום מוצא (עירוני/כפרי), מדד מסת הגוף (BMI), היסטוריית עישון, היסטוריית שתייה וממוצע זמן מסך יומי. ה-BMI חושב כמשקל הגוף בקילוגרמים חלקי ריבוע הגובה במטרים (BMI = weight [kg] / height [m]2). גובה ומשקל שדווחו על ידי המשתתפים נסרקו לאיתור ערכים לא סבירים לפני הניתוח; מקרים שבהם הגובה היה <140 cm או >200 cm, המשקל <30 kg או >150 kg, או ה-BMI <12 kg/m2 או >40 kg/m2 סומנו ונבדקו מול הרשומות המקוריות, כאשר ערכים חריגים מאושרים הוצאו מהניתוח. סטטוס העישון סווג כמי שמעולם לא עישן, מעשן מזדמן (עישון בימים מסוימים אך לא מדי יום), או מעשן קבוע (עישון יומי), בהתבסס על התנהגות עישון שדווחה על ידי המשתתפים במהלך החודש האחרון. צריכת אלכוהול סווגה כמי שמעולם לא שתה, שתייה מזדמנת (שתייה פחות מפעם בשבוע), או שתייה קבועה (שתייה לפחות פעם בשבוע), במהלך החודש האחרון. ממוצע זמן מסך יומי הוערך על ידי בקשה מהמשתתפים לדווח על סך השעות ליום שהושקעו בסמארטפונים, מחשבים, טאבלטים וטלוויזיה במהלך החודש האחרון, ונרשם כמשתנה רציף. השאלון נבדק בניסוי ראשוני (pilot) בקרב 30 סטודנטים לתואר ראשון לפני איסוף הנתונים הרשמי כדי להעריך את בהירות הפריטים ואת זמן ההשלמה; התבצעו התאמות קלות בניסוח על סמך המשוב מהניסוי הראשוני. מחקרים קודמים העלו כי משתנים אלה עשויים להוות גורמים מתערבים (confounders) בקשר שבין שינה לפעילות גופנית30,31. כל הכלים, התוכנות והחומרים ששימשו במחקר זה מופיעים בטבלת החומרים.
הערכת איכות השינה
מדד ה-PSQI, שפותח על ידי Buysse et al.9, שימש להערכת איכות השינה של המשתתפים במהלך החודש האחרון. הסולם כולל 19 פריטים בדיווח עצמי המרכיבים שבעה רכיבים: איכות שינה סובייקטיבית, זמן הירדמות, משך השינה, יעילות שינה רגילה, הפרעות שינה, שימוש בתרופות שינה ותפקוד במהלך היום. כל רכיב מקבל ציון מ-0 (ללא קושי) עד 3 (קושי חמור). הציונים של שבעת הרכיבים הרציפים הללו שימשו כאינדיקטורים נצפים עבור LPA. ה-PSQI הפגין מהימנות ותוקף טובים בקרב אוכלוסיות של סטודנטים במכללות בסין ובעולם30. במחקר הנוכחי, אלפא של קרונבך (Cronbach’s alpha) עבור הסולם היה 0.83.
הערכת פעילות גופנית
השאלון IPAQ-SF, שפותח על ידי Craig ואח' 32, שימש להערכת הפעילות הגופנית של המשתתפים במהלך שבעת הימים האחרונים. השאלון מתעד את התדירות (ימים/שבוע) ואת משך הזמן (דקות/יום) של פעילות גופנית מאומצת, פעילות גופנית בעוצמה מתונה והליכה. ערכי שווי מטבולי של פעילות (MET) חושבו על פי פרוטוקול הניקוד של IPAQ: הליכה MET-min/week = 3.3 × משך זמן × תדירות; עוצמה מתונה MET-min/week = 4.0 × משך זמן × תדירות; ועוצמה מאומצת MET-min/week = 8.0 × משך זמן × תדירות. ערכי ה-MET לשלושת סוגי הפעילות סוכמו כדי לקבל את סך הפעילות הגופנית. המשתתפים סווגו לרמות PA גבוהה, PA מתונה ו-PA נמוכה בהתאם לקריטריוני הניקוד של IPAQ-SF שתוקפו בקרב סטודנטים במכללות 33. PA גבוהה הוגדרה כפעילות מאומצת בשלושה ימים לפחות עם רמת פעילות גופנית כוללת של לפחות 1,500 MET-min/week, או כל שילוב של הליכה, פעילות בעוצמה מתונה או פעילות בעוצמה מאומצת בשבעה ימים או יותר עם רמת פעילות גופנית כוללת של לפחות 3,000 MET-min/week. PA מתונה הוגדרה כפעילות מאומצת בשלושה ימים לפחות למשך 20 דקות לפחות ליום, או פעילות בעוצמה מתונה ו/או הליכה בחמישה ימים לפחות למשך 30 דקות לפחות ליום, או כל שילוב של פעילויות המגיעות ללפחות 600 MET-min/week. PA נמוכה הוגדרה כמי שלא עומדת בקריטריונים ל-PA מתונה או גבוהה. ניקוי הנתונים בוצע על פי פרוטוקול הניקוד הרשמי של IPAQ: משתתפים שדיווחו על יותר מ-16 h של פעילות כוללת ליום הוצאו מהמחקר; משכי הזמן היומיים לכל קטגוריה קוצרו ל-180 min לפני חישוב ה-MET-minute; ומקרים עם תגובות שאינן עקביות לוגית הוסרו במהלך ניקוי הנתונים.
ניתוח סטטיסטי
ראשית, ביצעו LPA באמצעות תוכנה למידול פרופילים חבויים (ראו טבלת חומרים)34. מתוך 1,420 המשתתפים, ל-33 (2.3%) היו נתונים חסרים חלקיים ב-PSQI, כאשר 69 מתוך 9,940 תשובות ברמת הרכיב היו חסרות (0.69%; כלומר, 1,420 משתתפים × שבעה ציוני רכיבי PSQI). שיטת maximum likelihood של מידע מלא (FIML) שימשה לטיפול בערכים חסרים אלו תחת הנחת חסר-באקראי. מודלים עם פרופיל אחד עד חמישה נאמדו באופן סדרתי. התאמת המודל הוערכה באמצעות קריטריון המידע של אקאיקה (AIC), קריטריון המידע בייסיאני (BIC), BIC מותאם לגודל המדגם (aBIC), אנטרופיה, מבחן יחס הנראות של Lo-Mendell-Rubin (LMRT), ומבחן יחס הנראות של bootstrap (BLRT). ערכי AIC, BIC ו-aBIC נמוכים יותר מעידים על התאמה יחסית טובה יותר, בעוד שערכי אנטרופיה הקרובים ל-1 מעידים על דיוק סיווג גבוה יותר. LMRT ו-BLRT שימשו כדי לקבוע האם מודל עם k פרופילים התאים באופן מובהק יותר ממודל עם k − 1 פרופילים. הפתרון האופטימלי של הפרופילים נבחר על ידי שקלול משותף של מדדי התאמה סטטיסטיים, אנטרופיה, גודל המחלקה (שום פרופיל לא הכיל פחות מ-5% מהמדגם) וניתנות לפירוש מהותית, בהתאם להנחיות מקובלות לאנליזת פרופילים חבויים11,35. לאחר בחירת הפתרון האופטימלי של הפרופילים, השיוך הסביר ביותר של הפרטים למחלקה יוצא לתוכנה לניתוח סטטיסטי (ראו טבלת חומרים) לצורך ניתוחים המשכיים. מבחני chi-square של פירסון שימשו לבדיקת הבדלים במאפיינים דמוגרפיים וברמות פעילות גופנית בין פרופילי השינה. מבחני chi-square דו-משתני הראו כי רק מגדר ורמת פעילות גופנית נבדלו באופן מובהק בין ארבעת פרופילי השינה (שניהם p < 0.001), בעוד שגיל, שכבת גיל, BMI, עישון, שתייה וזמן מסך לא נבדלו (כולם p > 0.05; ראו טבלה משלימה 1). לכן, רק מגדר נכלל כמשתנה השהיה (covariate) במודל הרגרסיה הלוגיסטית הרב-נומית, בהתאם לעקרון התמציתיות (parsimony). ניתוח רגישות שתיקן לכל המשתנים השהיה שנאספו הניב את אותה תבנית ומובהקות של התוצאות (טבלה משלימה 2). רגרסיה לוגיסטית רב-נומית בוצעה כאשר רמת הפעילות הגופנית היא המשתנה התלוי (PA נמוכה = 0, PA מתונה = 1, PA גבוהה = 2) ופרופיל השינה הוא המשתנה הבלתי תלוי המרכזי. חושבו יחסי סיכויים (ORs) ומרווחי ברחשה של 95% (CIs). כל המבחנים היו דו-זנביים, ומובהקות סטטיסטית נקבעה עבור p < 0.05. מכיוון שהמשתתפים גויסו דרך כיתות שלמות המשובצות בתוך אוניברסיטאות, ביצענו ניתוח רגישות באמצעות מודול המדגמים המורכבים של התוכנה הסטטיסטית (ראו טבלת חומרים), תוך הגדרת האוניברסיטה והכיתה כמשתני קיבוץ (clustering), כדי לקבל שגיאות תקן מותאמות לתכנון המדגם. התוצאות היו עקביות עם אלו של הרגרסיה הלוגיסטית הרב-נומית הסטנדרטית (טבלה משלימה 3). מערך הנתונים הגולמי לאחר הסרת זיהוי, עם תוויות משתנים באנגלית ומסמך קוד (codebook), מסופק כ-קובץ משלים 1.