$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
המחקר אושר על ידי ועדת האתיקה הביו-רפואית של אוניברסיטת פקין (IRB00001052-11015). כל הנבדקים קיבלו הסכמה מדעת.
אוכלוסיית המחקר
הנתונים למחקר זה התקבלו מ-CHARLS 2020, הכולל 150 יחידות ברמת מחוז ו-450 יחידות ברמת כפר ברחבי המדינה, הכולל כ-10,000 משקי בית ו-19,395 אנשים בגיל הביניים וקשישים בגיל 45ומעלה. הסקר השתמש בגישת דגימה הסתברותית רב-שלבית. הפרטים נכללו אם עמדו בקריטריונים הבאים: גיל 45 ומעלה; בית נבחר בתגובה לפריט השאלון BA007: מהו סוג המגורים בכתובת המגורים?; וסיפק תשובה ברורה לסולם הדיכאון של מרכז ללימודים אפידמיולוגיים (CES-D-10) בן 10 פריטים. אנשים הוצאו מכלל אם עמדו באחד או יותר מהקריטריונים הבאים: תגובות המצביעות על מגורים מחוץ לבית (כולל קהילות, בתי חולים ואחרים או נעדרים); חוסר מידע על דיכאון; חוסר מידע על משתנים רלוונטיים. בסופו של דבר, 14,466 מבוגרים שעמדו בקריטריונים הראשוניים נכללו במחקר זה. תהליך בחירת המשתתפים המפורט מוצג באיור 1.

איור 1: תרשים זרימה לבחירת אוכלוסיית המחקר. אנא לחצו כאן כדי לצפות בגרסה גדולה יותר של איור זה.
הגדרת המשתנים
תסמיני דיכאון הוערכו באמצעות גרסה מפושטת של CES-D-10, שיש לה אמינות ותוקף גבוהים במחקרים קודמים22,23. ל-CES-D-10 יש 10 פריטים עם ארבעה סולמות: 0 = אף פעם, 1 = לפעמים או לעיתים רחוקות, 2 = לעיתים קרובות, ו-3 = תמיד, והפריטים החמישי והשמיני מדורגים הפוך24. הציון הכולל בסולם הדיכאון מתקבל על ידי סכימת דירוגי כל 10 הפריטים, הנעים בין 0 ל-30; ככל שהציון גבוה יותר, כך תסמיני הדיכאון של האדם חמורים יותר. בהתבסס על מחקרים קודמים באוכלוסיות בגיל הביניים וקשישים25,26, משתתפים עם ציוני CES-D-10 ≥10 סווגו כבעלי תסמינים דיכאוניים, בעוד שציונים <10 סווגו כלא מדוכאים.
המידע שנאסף מכל משתתפי המחקר כולל מידע דמוגרפי (גיל, מין, מצב משפחתי, מגורים ורמת השכלה), הרגלי חיים ומצב בריאותי פיזי (פעילות גופנית, עישון, שתייה, אינטראקציות חברתיות, בריאות מדורגת עצמית, מצב מחלות כרוניות וזמן תנומה), ומידע הקשור לתנאי מגורים וילדים (שביעות רצון מתנאי המחיה, מצב משפחתי של הצאצאים, ומצב הבריאות של הצאצא). מחלות כרוניות מאובחנות על ידי רופאים וכוללות לחץ דם גבוה, דיסליפידמיה, גידולים ממאירים, מחלות ריאה כרוניות, מחלות לב, שבץ, דלקת מפרקים ושיגרון.
ניתוח סטטיסטי
ניתוח רגרסיה לוגיסטית בינארית
משתנים קטגוריים מוצגים כתדירויות ואחוזים. נערכו מבחני חי-ריבוע כדי לבחון הבדלים בשכיחות דיכאון בין מאפיינים דמוגרפיים שונים באוכלוסיית המחקר. בוצע ניתוח רגרסיה לוגיסטית בינארית לזיהוי גורמים משמעותיים הקשורים לתסמיני דיכאון. ניתוחים סטטיסטיים אלו בוצעו באמצעות תוכנת SPSS (גרסה 21.0).
רגרסיית LASSO
בגרסה 4.3.2 של R, מערך הנתונים חולק באקראי לקבוצת אימון וקבוצת בדיקה ביחס 7:3 באמצעות זרע אקראי קבוע להבטחת שחזוריות. באופן ספציפי, 70% מהנתונים הוקצו לסט האימון לפיתוח המודלים, ו-30% לסט הבדיקה לאימות עצמאי. ראשית, משתנים שהיו מובהקים סטטיסטית (p < 0.05) ברגרסיה לוגיסטית בינארית נכללו בניתוח רגרסיית LASSO כדי למנוע התאמה יתר על ידי הקטנת המקדמים. ערך הלמדה האופטימלי נקבע באמצעות אימות צולב של 10 פעמים. ה-lambda.1se נבחר על בסיס כלל השגיאה הסטנדרטי כדי להשיג את ביצועי בחירת המשתנה האופטימליים. כיווץ LASSO שימש למניעת התאמת יתר על ידי הפחתת מקדמי משתנים פחות חשובים לאפס. לפרטים נוספים, אנא עיינו בתיק המשלים 1 (סעיף 1).
בניית נומוגרמה
החשיבות היחסית של המנבאים שנבחרו דורגה לפי הגודל המוחלט של המקדמים הסטנדרטיים שלהם והוצגה באמצעות תרשים יער. נבנתה נומוגרמה חיזוי באמצעות חבילת rms, בהתבסס על משוואת הרגרסיה הלוגיסטית הסופית. לפרטים נוספים, אנא עיינו בתיק המשלים 1 (סעיף 2).
אימות מודל
ההבחנה הוערכה על ידי חישוב השטח מתחת לעקומת מאפיין ההפעלה של המקלט (AUC) באמצעות חבילת pROC. ה-AUC חושב בנפרד עבור סט האימון וסט הניסוי. בדרך כלל, AUC גבוה מ-0.75 מציין אפליה טובה. לפרטים נוספים, אנא עיינו בתיק המשלים 1 (סעיף 3).
כיול: כיול המודל הוערך באמצעות מבחן ההתאמה הטובה של Hosmer-Lemeshow והוצג באמצעות עקומות כיול שנוצרו על ידי חבילת ה-rms, תוך שרטט הסתברויות חזויות מול תדרים נצפים. מבחן Hosmer-Lemeshow לא משמעותי (p > 0.05) מצביע על כיול טוב. לפרטים נוספים, אנא עיינו בתיק המשלים 1 (סעיף 4).
תועלת קלינית: התועלת הקלינית הוערכה באמצעות ניתוח עקומת החלטה (DCA) עם חבילת ה-rmda, תוך חישוב התועלת נטו בין הסתברויות סף מ-0% עד 100%. לפרטים נוספים, אנא עיינו בתיק המשלים 1 (סעיף 5).
אימות יער אקראי: כאימות משלימה, יושם מודל יער אקראי באמצעות חבילת randomForest להערכת יציבות והכללה של המודל. עקומת ההתכנסות של שיעור השגיאה הוצגה כדי להעריך את הקשר בין מספר העצים לדיוק החיזוי הן בקבוצות אימון והן במערכות מבחן. ערך p דו-צדדי < 0.05 נחשב כמצביע על הבדל מובהק סטטיסטית. לפרטים נוספים, אנא עיינו בתיק המשלים 1 (סעיף 6).