המחקר אושר על ידי ועדת האתיקה המוסדית של בית החולים הראשון ביולין, מחוז שאנשי, סין. מכיוון שמדובר בניתוח רטרוספקטיבי של נתונים שנאספו באופן שגרתי בבית החולים, ומכיוון שכל הנתונים אנונימיזו לפני הניתוח, הנדרשת להסכמה מדעת בכתב בוטלה. המחקר נערך בהתאם לדרישות הגנת הנתונים המוסדיות ולעקרונות הצהרת הלסינקי.
נתונים ושיטות
אוכלוסיית המחקר ומערך המחקר
זה היה מחקר עוקבה רטרוספקטיבי במרכז יחיד שנערך בבית החולים Yulin First, במחוז Shaanxi, סין. נסקרו רשומות בית חולים של מטופלים שאושפזו עם aSAH בין יולי 2022 לדצמבר 2024. המטופלים הזכאים היו מבוגרים בני 18 ומעלה עם אבחנה של aSAH שאושרה באמצעות דימות של הגולגולת וכלי הדם התוך-גולגולתיים. נדרש שהמטופלים יאושפזו תוך 72 h מהופעת הסימפטומים ושיעברו טיפול במפרצת במהלך האשפוז הנוכחי.
מטופלים הוצאו מהמחקר אם סבלו מדימום תת-עכבישי טראומטי, דימום תת-עכבישי שאינו אנוריסמלי, מחלה המטולוגית קשה קודמת, מחלה תו-גולגולית משמעותית אחרת המשפיעה על הערכת התוצאות, רישומים לא מספיקים להערכת DCI, או מוות בבית החולים לפני תקופה מספקת להערכת DCI. החרגת מקרי מוות מוקדמים בבית החולים עלולה להוביל להטיית שרידות, מכיוון שהמקרים החמורים ביותר עשויים להיות מוסרים מהניתוח. לכן, יש לפרש את הממצאים ככאלה החלים בעיקר על מטופלים ששרדו מספיק זמן כדי לאפשר הערכת DCI.
איסוף נתונים
משתנים קליניים, מעבדתיים, של דימות ומשתנים הקשורים לטיפול הופקו מהרשומות הרפואיות האלקטרוניות באמצעות טופס הפקת נתונים שהוגדר מראש. המשתנים שהופקו כללו גורמים דמוגרפיים, היסטוריה רפואית, מאפייני מפרצת (aneurysm), דרגה נוירולוגית מוקדמת, ממצאי דימות מוקדמים, תוצאות מעבדה, גישת טיפול וסיבוכים בבית החולים. נקודת הזמן המיועדת לניבוי הייתה לאחר השלמת האשפוז הראשוני וההערכה הפרי-אופרטיבית, אך לפני חלון הסיכון העיקרי ל-DCI.
גיל, מין, מדד מסת גוף (BMI), מצב עישון, שימוש באלכוהול, יתר לחץ דם, סוכרת, גודל המפרצת, מיקום המפרצת, דרגת Fisher משופרת, דרגת Hunt-Hess, דרגת World Federation of Neurological Surgeons, דימום תו-חדרי ודימום חוזר נלקחו מרשומות הקבלה ומבדיקות דימות ראשוניות. משתנים מעבדתיים, כולל המוגלובין, אלבומין בנסיבים ונתרן בנסיבים, הופקו מבדיקת הדם המוקדמת ביותר שזמינה לפני נקודת הזמן של הניבוי. בצקת מוחית הוערכה באמצעות דימות בבסיס או דימות מידי לאחר טיפול שהיה זמין לפני אבחנת DCI. משתנים שנרשמו רק לאחר הופעת DCI לא שימשו לפיתוח המודל כדי להפחית זליגת מידע.
אנמיה הוגדרה כהמוגלובין <110 g/L במטופלות ו-<120 g/L במטופלים. היפואלבומינמיה הוגדרה כאלבומין בנסיוב <35 g/L. היפונטremia הוגדרה כנתרן בנסיוב <130 mmol/L.
המטופלים טופלו בהתאם לנתיב הטיפול המוסדי ל-aSAH. הטיפול הסטנדרטי כלל סגירה מוקדמת של מפרצת, מעקב נוירולוגי, בקרת לחץ דם, ניהול נוזלים ואלקטרוליטים, ומעקב אחר סיבוכים הכוללים הידרוצפלוס, דימום חוזר, פרכוסים, זיהום, ואסו-ספאזם ו-DCI. נעשה שימוש ב-Nimodipine אלא אם כן הייתה התווית נגדית. ניטור של ואסו-ספאזם בוצע באמצעות בדיקה נוירולוגית, דימות כלי דם ואולטרסאונד דופלר טרנס-קרניאלי כאשר הדבר היה מומלץ מבחינה קלינית. הידרוצפלוס טופל באמצעות ניקוז נוזל מידולי כאשר נדרש. אופטימיזציה המודינמית וטיפול הצלה יושמו בהתאם להערכת הצוות המטפל. מכיוון שעצימות הטיפול ופרטי טיפול ההצלה לא תועדו במלואם במאגר הנתונים הרטרוספקטיבי, לא ניתן היה לבצע התאמה מלאה של גורמים אלה במודל, והם נחשבו כמגבלה.
אבחנה של DCI
DCI הוגדר כליקוי נוירולוגי מוקדי חדש, ירידה ברמת ההכרה, או אוטם מוחי חדש בבדיקת CT או MRI במעקב, המתרחשים במהלך חלון הסיכון הצפוי ל-DCI ואינם מוסברים על ידי גורם אחר. חלון הסיכון ל-DCI הוגדר כימים 3–14 לאחר אירוע הדימום (hemorrhage ictus). הסברים חלופיים, הכוללים דימום חוזר, הידרוצפלוס, פרכוס, זיהום, הפרעה מטבולית, השפעת חומרי הרגעה ואוטם הקשור לפרוצדורה, נבחנו לפני סיווג אירוע כ-DCI.
קביעת התוצאות בוצעה על ידי שני קליניקאים בעלי ניסיון בניהול aSAH, הכוללים נוירוכירורג אחד ונוירו-גאולוג/רופא טיפול נמרץ נוירולוגי אחד. השופטים סקרו רשומות קליניות, בדיקות נוירולוגיות, דוחות דימות ודימויי מעקב. הם היו עיוורים לפלט הסופי של המודל במהלך הערכת התוצאות. מחלוקות נפתרו באמצעות דיון להשגת קונצנזוס תוך סקירת הדימויים הרלוונטיים וציר הזמנים הקליני. הסכמה פורמלית בין השופטים לא נמדדה, ומגבלה זו נוספה כיוון שסיווג DCI עשוי לכלול שיקול דעת קליני.
חישוב גודל מדגם
גודל המדגם נאמד באמצעות כלל האגודל של 10 אירועים למשתנה (EPV)11. עם 20 משתני ניבוי מועמדים ושכיחות DCI צפויה של כ-35%, נדרש גודל מדגם מינימלי של 20 × 10 / 0.35 ≈ 572 חולים. הקוהורט הסופית שלנו, שכללה 680 חולים, עלתה על דרישה זו. למרות שכלל האירועים למשתנה (EPV) יושם כקו מנחה כללי, מידול ניבוי מודרני, ובמיוחד אלגוריתמים של למידת מכונה, עשוי לדרוש שיקולים גמישים יותר לגבי גודל המדגם. גודל הקוהורט הסופי נחשב למספיק כדי להבטיח את יציבות המודל, למזער את הסיכון להתאמת יתר (overfitting) ולאפשר הערכה מהימנה של מדדי ביצועי המודל.
פיתוח ותיקוף המודל
המטופלים הוקצו לקבוצת אימון של 544 מטופלים ולקבוצת תיקוף פנימית נפרדת של 136 מטופלים באמצעות חלוקה של 8:2 מותאמת לפי תוצאות. פירוק לפי תוצאות שימש כדי לשמור על פרופורציה דומה של אירועי DCI בשתי הקבוצות. DCI התרחש ב-175 מתוך 544 מטופלים בקבוצת האימון וב-42 מתוך 136 מטופלים בקבוצת התיקוף הפנימי. מבנה ההפרדה של 8:2 נשמר מכיוון שזו הייתה האסטרטגיה שנקבעה מראש לפיתוח המודל, והוא סיפק קבוצה נפרדת להערכה פנימית ראשונית של הביצועים.
בחירת המנבאים בוצעה רק בקוהורט האימון באמצעות רגרסיית LASSO (least absolute shrinkage and selection operator) עם תיקוף צולב של 10-fold. פרמטר הענישה שנבחר היה λ = 0.031. שישה מנבאים נשמרו: גיל, בצקת מוחית, היפואלבומינמיה, דרגת Fisher משופרת, דרגת Hunt-Hess, ודרגת World Federation of Neurological Surgeons.
מודלים של רגרסיה לוגיסטית, extreme gradient boosting, light gradient boosting machine, support vector machine ו-k-nearest neighbor פותחו תוך שימוש באותה קבוצה של Predictors שנבחרו. כל הליכי עיבוד המקדימים, בחירת ה-Predictors ופיתוח המודלים הוגבלו לקבוצת האימון. קבוצת התיקוף הפנימית לא שימשה במהלך בחירת ה-Predictors או כיוונון המודלים, והוערכה רק לאחר פיתוח המודל.
הזרע האקראי המספרי ששימש להקצאה המקורית של הקוהורט לא נשמר ברשומות הניתוח שנארכבו ולא ניתן היה לאמת אותו בדיעבד. תיקון אופטימיות באמצעות Bootstrap וולידציה פנימית חוזרת של k-fold לא היו זמינים אף הם. לפיכך, הניתוח הנוכחי מתואר כוולידציה פנימית בודדת של קבוצת השהיה (hold-out), ולא כהערכת ביצועים מתוקנת לאופטימיות או ככזו שעברה וולידציה צולבת חוזרת.
שחזוריות בלמידת מכונה
לבחירת המנבאים באמצעות LASSO נעשה שימוש באימות צולב של 10 קיפולים (10-fold cross-validation) בתוך קבוצת האימון. מנבאים רציפים עברו סטנדרטיזציה כאשר נדרש על ידי אלגוריתם המידול, ומנבאים קטגוריים יוצגו באמצעות קטגוריות קליניות בינאריות שהוגדרו מראש. קבוצת האימות לא שימשה לכוונון או לבחירת המודל.
חומרי הניתוח שנשמרו לא כללו את רשתות הכוונון הסופיות, את ההיפר-פרמטרים שנבחרו, את גרסאות החבילה המלאות, את הזרע האקראי (random seed) המקורי, או פרוצדורה מאושרת לטיפול בחוסר איזון בין מחלקות (class-imbalance) עבור XGBoost, LightGBM, SVM ו-KNN. לפיכך, הגדרות אלו לא שוחזרו או נאמדו. כתוצאה מכך, יש לפרש את ההשוואות בין מודלי הלמידה המכונית כהשוואות גשמיות (exploratory) של אלגוריתמים ולא כניסויי פיתוח מודלים בעלי שחזור מלא. מידע על יכולת השחזור, הן זה שזמין והן זה שאינו זמין, מסוכם ב-טבלה משלימה 1.
Bias Control & Sensitivity Analysis
חולים זכאים נכללו באופן רציף כדי להפחית הטיות בבחירה. הגדרות משתנים סטנדרטיות וטופס חילוץ נתונים מובנה שימשו להפחתת הטיות מידע. המנבאים הוגבלו למידע שהיה זמין לפני נקודת הזמן המיועדת לחיזוי DCI כדי להפחית דליפת מידע. בחירת המנבאים, עיבוד מקדים ופיתוח המודל בוצעו אך ורק בקוהורט האימון.
דרגת Fisher משוכללת, דרגת Hunt-Hess ודרגת WFNS מייצגות היבטים קשורים של עומס הדימום והחומרה הנוירולוגית. למרות שכל שלושת המשתנים נשמרו לאחר בחירה בשיטת LASSO, ערכי גורם ניפוח השונות (variance inflation factor) המספריים והתוצאות ממודל רגישות שהוציא סולמות חומרה חופפים לא נכללו בפלט הניתוח שנשמר. לא ניתן היה לשחזר ניתוחים אלה מהטבלאות המצטברות מכיוון שנדרשים מתאמים ברמת המטופל ופלטי מודל מותאמים. לפיכך, המקדמים האינדיבידואליים של מדדי חומרה אלה לא פורשו כהשפעות עצמאיות או סיבתיות. הם נשמרו אך ורק כמרכיבים של מודל הניבוי שנבחר בקוהורט האימון.
פרשנות מודלים
הסברה הייתה כי יישומיות קלינית דורשת פרשנות של המודל כאלמנט מכריע. במודל היעיל ביותר, השפעות המנבאים נבחנו כדי להעריך את הסבירות הקלינית והתאימות לידע פתופיזיולוגי קיים בנוגע לאיסכמיה מוחית מאוחרת. שיטה זו אפשרה הערכה ברורה של התנהגות המודל ושיפרה את הפוטנציאל לתרגום קליני.
נתונים חסרים
חוסר בנתונים ברמת המשתנה הוערך במהלך הכנת הנתונים. עם זאת, ספירות ואחוזי החסר המקוריים לפני השלמה (imputation) עבור כל משתנה מועמד לא נשמרו ברשומות הניתוח שנארכבו ולא ניתן היה לשחזרם במדויק מהטבלאות המרוכזות. כתוצאה מכך, לא הניחה המחקר שהנתונים שלמים, ולא דווחו אחוזי חסר שלא אומתו.
יש להפיק טבלת נתונים חסרים משלימה ומאושרת ישירות מתוך מערך הנתונים המקורי ברמת המטופל לאחר הסרת פרטים מזהים. טבלה משלימה 2 צריכה לדווח, עבור כל משתנה מועמד, את מספר ואחוזי התצפיות החסרות לפני הטיפול בנתונים, את השיטה ששימשה לטיפול בחסר, ואת מספר התצפיות שנכללו בניתוח הסופי. חוסר היכולת לשחזר את דפוס החסר המקורי הספציפי למשתנה נחשב כמגבלה של הדיווח הנוכחי.
ניתוח סטטיסטי
בחירת המנבאים בוצעה בקוהורט האימון באמצעות רגרסיית LASSO עם ולידציה צולבת של 10 מקפלים (10-fold cross-validation). היכולת להבחין (Discrimination) הוערכה באמצעות השטח תחת עקומת המאפיין של הפונקציה המקבלת (ROC curve) עם רווחי סמך של 95%. כיול (Calibration) תואר באמצעות תרשימי כיול והערכות נקודתיות של שיפוע הכיול, נקודת חיתוך הכיול ומדד Brier. רווחי סמך בשיטת Bootstrap עבור מדדי הכיול לא היו זמינים.
מדדי סיווג תלויי-סף חושבו מתוך מטריצות בלבול של תיקוף פנימי שהיו זמינות. סף ההסתברות המספרי המדויק ששימש ליצירת מטריצות הבלבול שנשמרו לא נשמר ולא ניתן היה לאמת אותו. לפיכך, מדדים אלו פורשו באופן תיאורי ולא שימשו כבסיס העיקרי לבחירת המודל. ניתוח עקומת החלטה נחשב אף הוא כניסיוני, מכיוון שטווח הסף המוגדר מראש ומסכי תועלת נטו ברמת המטופל לא נשמרו. ערך P דו-צדדי הנמוך מ-0.05 נחשב למשמעותי סטטיסטית.