פרוטוקול המחקר נבדק ואושר על ידי ועדת האתיקה של בית החולים הכללי של Northern Theater Command (מספר אישור Ethics Y (2026) 75). למרות שהנתונים הקליניים נגזרו מחולים שטופלו בין ינואר 2024 לדצמבר 2025, המחקר הנוכחי נערך כניתוח רטרוספקטיבי של רשומות קליניות קיימות. האישור האתי שהתקבל בשנת 2026 כיסה חילוץ נתונים רטרוספקטיבי, הסרת זיהוי, ניתוח ופרסום של נתונים קליניים אלו שנאספו בעבר. לא בוצעו התערבויות פרוספקטיביות או גיוס מטופלים לפני קבלת האישור האתי. הדרישה להסכמה מדעת בכתב בוטלה בשל עיצוב המחקר כניתוח תצפיתי רטרוספקטיבי. כל ההליכים בוצעו בהתאם להצהרת הלסינקי, וכל הנתונים הקליניים עברו הסרת זיהוי לפני הניתוח. הניתוח הנוכחי נערך כחלק מפרויקט מחקר מאושר של נתונים קליניים נוירו-אינטרוונציוניים רטרוספקטיביים. החומרים והציוד הנדרשים להליכים המתוארים להלן מסוכמים בטבלת החומרים.
גיוס מטופלים וקביעת CIE
מטופלים שעשויים להיות מתאימים זוהו באמצעות שאילתה במערכת הרשומות הרפואיות האלקטרוניות עבור מטופלים שטופלו במחלקה לנוירוכירורגיה בבית החולים הכללי של Northern Theater Command בין ינואר 2024 לדצמבר 2025. החיפוש הראשוני הוגבל לפי תאריך אשפוז, מחלקה, אבחנה של מחלה צרברו-וסקולרית ורישומים של הליכים נוירו-אינטרוונציוניים. זרימת העבודה של סינון המטופלים והקצאת когоורט מוצגת ב-איור 2. מטופלים היו זכאים להכללה אם עמדו בכל הקריטריונים הבאים: (1) גיל בין 30 ל-80 שנים; (2) אבחנה של מחלה צרברו-וסקולרית הדורשת הערכה או טיפול נוירו-אינטרוונציוני; (3) זמינות של רשומות קליניות מלאות ותוצאות בדיקות מעבדה; ו-(4) בדיקת CT ראש לאחר ניתוח שבוצעה תוך 3 ימים לאחר הניתוח.
מטופלים הוצאו מהמחקר אם היו להם נתונים קליניים או נתוני דימות חסרים, אם לא עברו התערבות אנדווסקולרית, אם סבלו מגסרים נוירולוגיים חמורים לפני הניתוח (ציון של ≥4 במד Rankin Scale המעודכן), אם הודגם אוטם מוחי חריף בבדיקת דימוי שקלול דיפוזיה (diffusion-weighted imaging) לפני הניתוח, אם סבלו ממחלת כליות כרונית בשלב G4 ומעלה לפי הקריטריונים של Kidney Disease: Improving Global Outcomes (KDIGO), או אם סבלו ממחלה קרדיווסקולרית חמורה. קריטריוני ההכללה וההדרה הוחלו באופן רצופי. לאחר החיפוש האלקטרוני הראשוני, כל רשומה של מועמד נבדקה ידנית כדי לאשר את התאמתו ולתעד את סיבת ההדרה במידת הצורך. לאחר הסינון, מטופלים מתאימים הוקצו באופן אקראי לקבוצות אימון ובדיקה ביחס של 8:2, כולל קבוצת אימון (n = 128) וקבוצת בדיקה (n = 3). נעשה שימוש בדגימה אקראית מרובדת כדי לשמור על התפלגות עקבית של מקרי CIE ומקרי non-CIE בין שתי הקבוצות.
האבחנה של CIE התבססה בעיקר על הקשר הכרונולוגי בין הופעת התסמינים לבין פרוצדורות התערבות צרברו-וסקולרית, בשילוב עם הערכה דימויית לשלילת גורמים אחרים. שני מומחים עצמאיים להתערבות צרברו-וסקולרית סקרו את הרשומות הקליניות לאחר הניתוח, את זמן הופעת התסמינים, את הביטויים הנוירולוגיים ואת ממצאי ההדמיה לאחר הניתוח עבור כל המקרים החשודים. מחלוקות נפתרו באמצעות דיון עד להשגת הסכמה. מטופלים נחשבו תואמים ל-CIE אם התפתחו עיוורון קורטיקלי, שינוי במצב ההכרה (כולל נמנום, דליריום או קומה), פרכוס, או גרעונות נוירולוגיים מוקדיים כגון חמיפרזיס או אפזיה, בתוך 24 h מסיום הפרוצדורה הנוירו-אינטרוונציונית.
לאחר הופעת התסמינים, סיבות חלופיות העשויות לגרום למניפסטציות נוירולוגיות דומות נשללו באמצעות סקירה של תמונות CT או MRI לאחר הניתוח, יחד עם המעקב הקליני. CT ראש לאחר הניתוח בוצע לפי הפרוטוקול המוסדי הסטנדרטי ל-CT קרניאלי ללא חומר ניגוד. פרמטרי הרכישה המהותיים כללו מתח צינור של 80 kVp, זרם צינור של 260 mA, עובי חתך של 0.5 mm ושחזור אקסיאלי סטנדרטי. התמונות נסקרו בחלונות מוח ועצם על ידי קלינאים מנוסים. במהלך האבחנה המבדלת נשקלו דימום תו-גולגולתי חריף, דימום תת-עכבישי, דימום תוך-מוחי, אוטם מוחי נרחב שהתפתח לאחרונה, שינויים הקשורים להתקפים, זיהום והפרעות מטבוליות. בדיקות CT לאחר הופעת התסמינים הראו בדרך כלל בצקת מוחית מוקדית או דיפוזית וצלליים בצפיפות גבוהה המערבים את הקורטקס, אזורים תת-קורטיקליים או את החלל התת-עכבישי, המדמים נגעים דימומיים. ממצאי MRI מראים בדרך כלל נפיחות קורטיקלית, במיוחד באזור הקורטקס הטמפורופריאטלי אוקסיפיטלי. עם זאת, בחלק מהמטופלים, בדיקות CT או MRI לא הראו חריגות ברורות למרות מניפסטציות קליניות תואמות.
איסוף משתנים פרי-אופרטיביים ובניית CGR
נתונים דמוגרפיים כלליים, כולל גיל, מין ומשקל גוף, הופקו משדות מובנים במערכת הרשומות הרפואיות האלקטרוניות. משתנים של היסטוריה קלינית, כולל היסטוריית עישון, היסטוריית צריכת אלכוהול, יתר לחץ דם, סוכרת ומחלה כלילית של הלב, הופקו מהערות קבלה, מרשומות היסטוריה רפואית קודמת ומאבחנות שחרור, ונבדקו להצלבה לצורך עקביות. נתונים מעבדתיים הופקו ממערכת המידע המעבדתית. עבור כל מטופל, נעשה שימוש בדגימת דם ורידית ראשונה שנלקחה לאחר הקבלה ולפני הפרוצדורה הנוירו-אינטרוונציונית. המשתנים המעבדתיים שהופקו כללו קריאטינין בנסיב, קצב סינון גלומרולרי מוערך (eGFR), כולסטרול כולל ורמות טריגליצרידים. eGFR חושב באמצעות המשוואה המבוססת על קריאטינין של Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) והתקבל ממערכת המידע המעבדתית של בית החולים.
מכיוון שמורכבות הפרוצדורה, החשיפה לחומר ניגוד ומיקום הנגע עשויים להשפיע על פגיעה במחסום הדם-מוח ועל שימור חומר הניגוד במהלך פרוצדורות נוירו-אינטרוונציוניות, משתנים פרוצדורליים פרי-אופרטיביים שולבו גם הם במסגרת החיזויית. המשתנים הקשורים לפרוצדורה כללו את מיקום הנגע, סוג הפרוצדורה, משך הפרוצדורה, סוג חומר הניגוד ונפח הניגוד הכולל. מיקום הנגע וסוג הפרוצדורה נקבעו מתוך דוחות ניתוח ורישומים אנגיוגרפיים. משך הפרוצדורה הוגדר כזמן מרגע הניקור העורקי ועד לסיום הפרוצדורה הנוירו-אינטרוונציונית. נפח הניגוד הכולל הוגדר כנפח המצטבר של חומר ניגוד יודנני שניתן מתחילת הפרוצדורה ועד סופה. חומרי הניגוד ששימשו במחקר זה היו ICA לא-יוניים, כולל iodixanol (320 mg של יוד/mL), שניתנו בהתאם לפרקטיקה הקלינית המוסדית. כל הפרוצדורות הנוירו-אינטרונבנציוניות שנכללו במחקר זה בוצעו תחת הנחיית אנגיוגרפיית חיסור דיגיטלית על ידי צוות נוירו-אינטרוונציוני מנוסה במרכז המחקר. כל הפרוצדורות בוצעו באמצעות גישה עורקית טרנס-פמורלית. במהלך ההתערבות, החולים קיבלו נוגדי קרישה תוך-ניתוחיים עם heparin, וזמן קרישה מופעל (activated clotting time) נשמר בטווח היעד של 250–30 s. סימנים חיוניים נוטרו באופן רציף לאורך כל הפרוצדורה.
התערבויות טיפוליות, כולל אמבוליזציה של מפרצת, אנגיופלסטיה והשמה של תומך (סטנט), בוצעו בהתאם למאפייני הנגע ולשיקול דעת המנתח. מיקרוקתטרים, סלילים, תומכים או בלונים מתאימים נבחרו בהתאם לדרישות הפרוצדורה. במהלך ההדמיה האנגיוגרפית וההתערבות, חומרי ניגוד הומסקו באמצעות הזרקה בלחץ גבוה או עירוי ידני בקצב קבוע. סוג חומר הניגוד ונפח הניגוד הכולל הופקו מרישומי הפרוצדורה ונבדקו מול רישומי ההרדמה או הסיעות, כאשר אלו היו זמינים. רשומות עם מידע לא עקבי או חסר לגבי נפח חומר הניגוד נסקרו ידנית לפני הכללתן במערך הנתונים האנליטי הסופי.
כל המשתנים הקשורים להליך, כולל משך הפרוצדורה, סוג ההתערבות הטיפולית, שימוש במכשירים ורישומים של מתן חומר ניגוד, נסקרו באופן שיטתי ונבדקו לשלימות ועקביות לפני הניתוח. כדי להעריך באופן כמותי את חוסר האיזון בין עומס חומר הניגוד לבין קיבולת הפינוי הכליתי, נקבע ה-CGR כמשתנה ניבוי מרכזי וחושב לפי המשוואה הבאה:
(1)
נפח חומר הניגוד הכולל נרשם ב-mL, וה-eGFR נרשם ב-mL/min/1.73 m2. ערך ה-eGFR ששימש לחישוב ה-CGR היה ערך ה-eGFR הפרה-אופרטיבי שהתקבל מדגימת הדם הוורידית הראשונה שנלקחה לאחר הקבלה ולפני ההליך הנוירו-אינטרוונציוני. עבור כל מטופל, ה-CGR חושב לאחר אימות של נפח חומר הניגוד הכולל ושל ה-eGFR. המונה היה נפח חומר הניגוד הכולל ב-mL, והמכנה היה ה-eGFR ב-mL/min/1.73 m2. כלל חישוב זה הוחל על כל המטופלים לפני פיתוח המודל.
בחירת מאפיינים וזרימת עבודה של למידה חישובית (ML)
על מנת להבטיח רובסטות אנליטית, כל המשתנים הקליניים והמעבדתיים עברו הליכי בקרת איכות שיטתיים לפני האנליזה. שיעור הנתונים החסרים חושב עבור כל משתנה מועמד לפני פיתוח המודל. במערך הנתונים האנליטי הסופי, לא נצפו ערכים חסרים בקרב המנבאים הנכללים או תגי התוצאה; לפיכך, אף משתנה לא הוצא מהניתוח בשל חסר בנתונים, ולא נדרשה השלמה מרובה (multiple imputation). סיכום מפורט של הנתונים החסרים מופיע ב- טבלה משלימה 1.
כל המשתנים המנבאים הפוטנציאליים הוערכו בתחילה באמצעות ניתוחים חד-משתניים כדי לבחון את קשריהם להתרחשות של CIE. שיטות סטטיסטיות נבחרה בהתאם למאפייני ההתפלגות של הנתונים. משתנים בעלי התפלגות נורמלית הוצגו כממוצע ± סטיית תקן והושוו באמצעות מבחן t למדגמים בלתי תלויים. משתנים בעלי התפלגות לא נורמלית הוצגו כחציונית (טווח בין-רבעוני, IQR) ונותחו באמצעות מבחן Mann–Whitney U. משתנים קטגוריאליים הוצגו כשכיחויות ואחוזים ונותחו באמצעות מבחן חי-ריבוע או מבחן Fisher המדויק. מובהקות סטטיסטית הוגדרה כ- p < 0.05.
כל מערך הנתונים חולק באופן אקראי לקבוצות אימון ובדיקה ביחס של 8:2 לפני פיתוח המודל. כדי למזער התאמת יתר (overfitting) ולהפחית מולטיקוליניאריות, הוחלה רגרסיית LASSO (least absolute shrinkage and selection operator) בשילוב עם אימות צולב (cross-validation) של 10 קיפולים לצורך בחירת מאפיינים. בחירת המאפיינים, כוונון ההיפר-פרמטרים ופיתוח המודל בוצעו באמצעות קבוצת האימון בלבד כדי למנוע זליגת נתונים. קבוצת הבדיקה הפנימית לא שימשה במהלך הערכת פרמטרי השלמה (imputation), בחירת מאפיינים, כוונון היפר-פרמטרים או אימון המודל, והוצמדה לשימוש חד-פעמי לצורך הערכת הביצועים הסופית. אופטימיזציה של ההיפר-פרמטרים בוצעה באמצעות אסטרטגיית חיפוש רשת (grid search) עם אימות צולב של 10 קיפולים בתוך קבוצת האימון. באופן ספציפי, קבוצת האימון חולקה באופן אקראי ל-10 תת-קבוצות מוexcluded-ly. במהלך כל חזרה, תשע תת-קבוצות שימשו לאימון המודל, והתת-קבוצה הנותרת שימשה לאימות. תהליך זה חזר 10 פעמים כדי להבטיח שכל תת-קבוצה שימשה כקבוצת אימות בדיוק פעם אחת.
כל האנליזות החישוביות בוצעו בסביבת Python. פיתוח והערכה של מודלי ML בוצעו באמצעות scikit-learn. סקריפטים של אנליזה המשמשים לעיבוד מקדים של נתונים, בחירת תכונות, אימון מודלים, הערכת ביצועים ו-DCA סופקו כקבצי קוד משלימים. לצורך יישום המודל, התכונות שנבחרו מרגרסיית LASSO שימשו כמשתני קלט עבור כל המסווגים המועמדים. משתנים רציפים עברו סטנדרטיזציה באמצעות נורמליזציית z-score בתוך קבוצת האימון, ואותם פרמטרי סקיילינג הוחלו על קבוצת הבדיקה הפנימית. משתנים קטגוריאליים קודדו לפי סכמת הקידוד שהוגדרה מראש, והגדרות המשתנים וסכמת הקידוד מופיעות ב-Supplementary Table 2. חמישה מסווגים יושמו: Gaussian Naive Bayes, SVM, KNN, LightGBM ו-MLP. עבור SVM, כווננו סוג הגרעין (kernel type), פרמטר הרגולריזציה C ומקדם הגרעין gamma. עבור KNN, כווננו מספר השכנים ואסטרטגיית שקלול המרחק. עבור LightGBM, כווננו מספר האומדנים (estimators), קצב הלמידה, עומק העץ המקסימלי ומספר העלים. עבור MLP, כווננו מבנה השכבה החבויה, פונקציית ההפעלה, פרמטר הרגולריזציה וקצב הלמידה. Gaussian Naive Bayes יושם תוך שימוש בהגדרת variance smoothing שצוינה. שילוב ההיפר-פרמטרים האופטימלי עבור כל מודל נבחר בהתאם לביצועי ולידציה צולבת (cross-validation) בתוך קבוצת האימון. שילובי ההיפר-פרמטרים האופטימליים הסופיים שנבחרו עבור כל מודל מסוכמים ב-Supplementary Table 3. לאחר בחירת ההיפר-פרמטרים, כל מודל סופי הותאם מחדש (refitted) על קבוצת האימון המלאה והוערך פעם אחת בקבוצת הבדיקה הפנימית.
השוואות של ביצועי המודל
על בסיס המנבאים שנבחרו, פותחו והושוו חמישה מודלים של למידה חישובית (ML). ביצועי המודלים הוערכו הן בקבוצת האימון והן בקבוצת הבדיקה באמצעות השטח תחת עקומת מאפייני תגובת המקלט (AUC), רווח ברחב信頼 95% (95% CI), רגישות, סגוליות ודיוק. עקומות מאפייני תגובת המקלט הופקו כדי להשוות את ביצועי ההבחנה של המודלים השונים. סף הסיווג האופטימלי עבור כל מודל נקבע באופן בלעדי בתוך קבוצת האימון באמצעות מדד Youden. הסף שהתקבל נקבע לאחר מכן כקבוע והוחל ללא שינוי על קבוצת הבדיקה הפנימית לצורך חישוב מדדי ביצועים התלויים בסף.
ניתוח עקומת החלטה
ניתוח עקומת תועלת (DCA) בוצע על ידי חישוב התועלת נטו של כל מודל על פני טווח של הסתברויות סף, והשוואת האסטרטגיה המבוססת על המודל לשתי אסטרטגיות ברירת מחדל: טיפול בכל המטופלים ואי-טיפול באף מטופל. DCA נבחר מכיוון שהוא מאפשר הערכה סימולטנית של ביצועי ההבחנה והשמישות הקלינית הפוטנציאלית במסגרת קבלת החלטות פרי-אופרטיביות. כדי להקל על השחזור, זרימת העבודה המלאה בוצעה ברצף הבא: זיהוי מטופלים ממערכת הרשומות הרפואיות האלקטרוניות, יישום קריטריוני הכללה והחרגה, קביעת CIE על ידי שני מומחים עצמאיים, חילוץ ובדיקה הצלבית של משתנים דמוגרפיים, היסטוריה קלינית, מעבדה ומשתנים הקשורים להליך, חישוב CGR, הערכת נתונים חסרים ופיצול לקבוצת אימון ובדיקה, בחירת מאפיינים מבוססת LASSO בתוך קבוצת האימון, כוונון היפר-פרמטרים באמצעות ולידציה צולבת של 10 קיפולים, התאמה מחדש של המודל הסופי בקבוצת האימון המלאה, הערכת קבוצת בדיקה פנימית, הערכת ביצועים מבוססת ROC, קביעת סף לפי מדד Youden, ו-DCA.