מחקר זה, הכולל משתתפים אנושיים, רשומות קליניות ודגימות סרום, נערך בהתאם להצהרת הלסינקי ודרישות המוסד למחקר בבני אדם13. פרוטוקול המחקר נבחן ואושר על ידי בית החולים לרפואה סינית של מחוז ג'יאנגסו תחת מספר אישור CJ28J12-45. עבור נתונים קליניים אנונימיים שנאספו רטרוספקטיבית ודגימות סרום שאריות, הדרישה להסכמה מדעת בכתב נוספת בוטלה בהתאם לתקנות האתיקה המוסדיות. עבור דגימות שנאספו פרוספקטיבית, התקבלה הסכמה מדעת בכתב לפני הגיוס. כל מידע על המשתתפים ודגימות הסרום הגיעו למצב של אנונימיות לפני הניתוח. מזהי מחקר שימשו לקישור נתונים, ולא נכלל מידע מזהה אישי במערך הנתונים האנליטי.
מערך המחקר והקצאת קבוצות
מחקר תצפיתי אמביספקטיבי זה, שנערך בשני מרכזים, העריך את הערך האבחנתי והפרוגנוסטי של פאנל משולב של סממנים ביולוגיים בנסיוב, הכולל את carcinoembryonic antigen (CEA), squamous cell carcinoma antigen (SCC-Ag) ו-carbohydrate antigen 125 (CA125), בקרצינומה של תאי קשקש של הוושט (ESCC). נתונים רטרוספקטיביים נאספו מרשומות קליניות רלוונטיות ומדגמי נסיוב שנותרו לפני הטיפול, וגיוס פרוספקטיבי שימש לאיסוף מדגמי נסיוב נוספים לפני הטיפול ולהשלמת המעקב במסגרת אותו מחקר. המשתתפים גויסו מבית חולים שלישוני אחד וממרכז בריאות קהילתי אחד בין ינואר 2021 לדצמבר 2024, כאשר המעקב נמשך עד 31 בדצמבר 2025. מרכז A היה בית החולים השלישוני ושימש כקוהורט האימון לפיתוח המודל, להערכת המקדמים, לבחירת ערכי הסף (cutoff) ולתיקוף פנימי. מרכז B היה מרכז הבריאות הקהילתי ושימש כקוהורט התיקוף החיצוני. המחקר דווח בהתאם להנחיות ה-Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD)14.
בסך הכל נכללו 277 משתתפים, כאשר 194 היו ממרכז A ו-83 ממרכז B. קבוצת האימון כללה 114 חולי ESCC, 40 בקרי מחלהי וריאליים של הוושט ו-40 בקרי בריאות. קבוצת התיקוף החיצונית כללה 53 חולי ESCC, 15 בקרי מחלה וריאליים של הוושט ו-15 בקרי בריאות. נעשה שימוש בהקצאה מבוססת מרכזים במקום פיצול אקראי כדי להפחית זליגת מידע ולבחון את יציבות המודל בסביבה קלינית עצמאית. הניתוח האבחנאי השווה רמות של CEA, SCC-Ag ו-CA125 בסרום בין חולי ESCC, חולים עם מחלה וריאלית של הוושט ובקרי בריאות. הניתוח הפרוגנוסטי הוגבל לחולי ESCC ובחן האם עומס משולב של ביומרקרים שיפר את סיווג הסיכונים להישרדות כוללת (OS) ולהישרדות ללא התקדמות המחלה (PFS). עיצוב המחקר, סינון המשתתפים, הקצאת הקבוצות, זרימת העבודה של ביומרקרים בסרום, פיתוח המודל האבחנאי, התיקוף החיצוני והשוואת המודל הפרוגנוסטי מסוכמים ב-איור 1.
זרימת המשתתפים כללה סינון, החרגות, חלוקה לקבוצות אבחנתיות והקצאה מבוססת מרכזים. סך הכל 326 אנשים עברו סינון. ארבעים ותשעה אנשים הוחרגו, כולל 12 שקיבלו טיפול נגד גידולים לפני דגימת הדם, 8 עם גידול ממאיר ראשוני אחר, 9 עם זיהום חמור, מחלה אוטואימונית או תחלואה נלווית בלתי מאוזנת שעלולה להשפיע על סמני גידולים בסרום, 11 עם המוליזה, ליפמיה, הפרעה עקב ג'aundice, נפח סרום לא מספיק או פרק זמן עיבוד ממושך, ו-9 עם חסר בנתונים קליניים, פתולוגיים או נתוני מעקב מרכזיים. האוכלוסייה האנליטית הסופית כללה 277 משתתפים: 167 חולי ESCC, 55 ביקורות של מחלות וריד מאכל שפירות, ו-55 ביקורות בריאות.
הסבר לבחירת גודל המדגם
גודל המדגם נקבע על פי הזמינות הרצופית של משתתפים מתאימים, דגימות סרום טרום-טיפול ומידע במעקב במהלך תקופת המחקר שהוגדרה מראש. המודל האבחנאי השתמש בשלושה סממנים ביולוגיים (biomarkers) רציפים בסרום כנבואים מועמדים. בקבוצת האימון היו זמינים 114 מקרי ESCC ו-80 ביקורות ללא ESCC, המקבילים ל-38.0 מקרי ESCC ו-26.7 ביקורות ללא ESCC עבור כל נבוא סמן ביולוגי מועמד. דבר זה תמך ברגרסיה לוגיסטית מממד נמוך והפחית את הסיכון להתאמת יתר (overfitting). קבוצת התיקוף החיצונית כללה 53 מקרי ESCC ו-30 ביקורות ללא ESCC, ושימשה לתיקוף בלבד ללא הערכה מחדש של המקדמים. הרציונל לגודל המדגם תאם את העקרונות הנוכחיים למודלים של ניבוי, המדגישים את שכיחות התוצאה, מספר פרמטרי הנבואים המועמדים, ביצועי המודל הצפויים ותיקוף עצמאי, במקום פיצול אקראי של נתונים15.
גיוס משתתפים, קביעת זכאות ואיסוף נתונים
המשתתפים גויסו באופן רציף במהלך תקופת המחקר. במרכז A, חולי ESCC גויסו ממחלקות לניתוחים בחזה, גסטרואנטרולוגיה ואונקולוגיה. במרכז B, חולי ESCC זוהו באמצעות סריקה מבוססת קהילה, רישומי מעקב או נתיבי הפניה, ונכללו רק כאשר היו זמינים אישור פתולוגי ודגימות סרום שלפני הטיפול. ביקורות של מחלות ושט שפירות גויסו מאנשים שעברו גסטרוסקופיה בשל דיספגיה, תסמיני רפלקס או חשד למומיות בושט, והפגינו ממצאים פתולוגיים שאינם ממאירים. ביקורות בריאות נבחרו מבדיקות בריאות שגרתיות במהלך אותה תקופה, ללא היסטוריה של ממאירות, זיהום פעיל, תפקוד כבדי או כלייתי לק, או עדות בדימות/אנדוסקופיה לנגעים תופסי מקום בושט.
חולי ESCC היו זכאים להשתתפות אם הוכח כי הם סובלים מ-ESCC ראשוני באמצעות ביופסיה גסטרוסקופית או פתולוגיה לאחר ניתוח, לא קיבלו טיפול אנטי-גידולי לפני דגימת הדם, והיו בעלי נתונים קליניים מלאים בקו הבסיס, רישומים פתולוגיים ודגימות סרום לפני טיפול. משתתפים הוצאו מהמחקר אם סבלו מממאירות ראשונית אחרת, טיפול אנטי-גידולי קודם לפני דגימת הדם, זיהום חריף קשה, מחלה אוטואימונית פעילה, תחלואה נלווית לא מאוזנת שעלולה להשפיע על רמות סמני גידול בסרום, המוליזה נראית לעין, ליפמיה, הפרעה עקב צהובון, נפח סרום לא מספיק, היעדר משתנים קליניים מרכזיים או נתוני מעקב לא זמינים. הניתוח הפרוגנוסטי הוגבל לחולי ESCC עם שלב TNM (גידול-קשרית-גרורה) זמין ותוצאות מדידות של הישרדות או התקדמות המחלה. חולים עם מעקב קצר מ-6 חודשים נשמרו במחקר אם התרחשה מוות או התקדמות המחלה בתוך תקופה זו.
משתני הבסיס כללו גיל, מין, היסטוריית עישון, היסטוריית צריכת אלכוהול, מדד מסת גוף (BMI), מרכז השתתפות, קבוצה אבחנתית, סוג הנגע ורמות ביומיות של CEA, SCC-Ag ו-CA125. עבור חולי ESCC, משתנים נוספים כללו מיקום הגידול, אורך הגידול, דיפרנציאציה היסטולוגית, עומק הפלישה, מצב קשרי לימפה, שלב TNM, פלישה לימפו-וסקולרית, שיטת טיפול ראשונית, רצידיב, התקדמות המחלה, מצב הישרדות ותאריך המעקב האחרון. היסטוריית עישון הוגדרה כעישון נוכחי או עישון סדיר בעבר של לפחות 100 סיגריות במהלך החיים. היסטוריית צריכת אלכוהול הוגדרה כצריכת אלכוהול לפחות פעם בשבוע למשך יותר מ-6 חודשים. שלב ה-TNM נקבע בהתאם למהדורה השמינית של מערכת הדירוג של ה-American Joint Committee on Cancer/Union for International Cancer Control לסרטן הוושט16. שקופיות פתולוגיות ממקרי ESCC ונגעים שפירים נבחנו באופן עצמאי על ידי שני פתולוגים בכירים, וכל מחלוקת נפתרה על ידי פתולוג שלישי.
שני חוקרים מיומנים חילצו באופן עצמאי מידע קליני מתוך רשומות רפואיות אלקטרוניות, דוחות פתולוגיים, מערכות מידע של מעבדות, רשומות טיפול, רשומות מעקב ורישומי תמותה, תוך שימוש במילון נתונים שהוגדר מראש. סתירות נפתרו באמצעות בדיקה של הרשומות המקוריות, כאשר חוקר שלישי קיבל את ההחלטה הסופית בעת הצורך. מיקום הגידול סווג כוורסיגוס תורקי עליון, אמצעי או תחתון. הדיפרנציאציה ההיסטולוגית סווגה כדיפרנציאציה טובה, בינונית או דלה. יחידות של משתנים רציפים סטנדרטו בשני המרכזים המשתתפים לפני הניתוח.
דגימת דם, עיבוד סרום ומדידת סממנים ביולוגיים
כל דגימות הדם נאספו לפני הטיפול הראשוני, בתוך 7 ימים לאחר האישור הפתולוגי, ולפחות 3 ימים לפני תחילת הטיפול. דגימות דם ורידי פריפרי של 5–8 mL נאספו בין השעות 7:00 ל-9:00 בבוקר לאחר צום לילה, באמצעות מבחנות להפרדת סרום. הדגימה נדחתה אם המשתתף קיבל עירוי דם, עירוי אלבומין בנפח גדול או עבר הליך פולשני בתוך 72 h לפני האיסוף. הדגימות הושארו לקרישה בטמפרטורת החדר למשך 30 min, ולאחר מכן סובבו בצנטריפוגה ב-1,500 × g למשך 10 min. הסרום הופרד וחולק ל-3–5 מבחנות מסומנות בברקוד, בנפח של 300–500 µL לאליקוט. האליקוטים נשמרו ב-−80°C עד לבדיקה, ולאף אליקוט לא בוצע יותר ממחזור הקפאה-הפשרה אחד. דגימות עם המוליזה נראית לעין, ליפמיה, הפרעה בשל צהבת, נפח לא מספיק או מרווח של יותר מ-2 h בין האיסוף לההקפאה הוצאו מהמחקר.
רמות של CEA, SCC-Ag ו-CA125 בסרום נמדדו באמצעות מערכות אוטומטיות של immunoassay כימילומינצנטי, בהתאם להוראות היצרן. מנות סרום הופשרו בטמפרטורה של 4°C לפני הבדיקה, ערבובו של התמיסה נעשה בעדינות, והניתוח בוצע בריצה אחת במידת האפשר. במקרים בהם נדרשה בדיקה של סדרות דגימות, סדר הדגימות נקבע באופן שמונע קיבוץ סיסטמטי לפי מצב המחלה, מרכז, שלב פתולוגי או תוצאת מעקב. דגימות שחרגו מטווח המדידה האנליטי דוללו ונבדקו שוב בהתאם להוראות הרכיבים. צוות המעבדה לא היה מודע (blinded) לקבוצת המשתתפים, לשלב הפתולוגי, למצב הטיפול או לתוצאות המעקב.
הספים הקליניים היו 5.0 ng/mL עבור CEA, 1.5 ng/mL עבור SCC-Ag, ו-35 U/mL עבור CA125, בהתבסס על סטנדרטים של מעבדה מקומית ועל טווחי ייחוס של היצרנים. ספים אלו שימשו להגדרת עלייה במדדים ועומס גבוה של סממנים ביולוגיים משולבים לצורך ניתוחים של קשר קלינו-פתולוגי וניתוחי סטרטיפיקציה של הישרדות. לצורך סיווג אבחנתי, נבחרו ערכי חתך המבוססים על נתונים רק בקוהורט האימון באמצעות מדד Youden, ולאחר מכן הוחלו ללא שינוי על קוהורט התיקוף החיצוני. כיול המכשיר ובקרת איכות פנימית בוצעו לפני הבדיקה. בקרות ברמה גבוהה ונמוכה נכללו במהלך כל יום בדיקה, והבדיקות המשיכו רק כאשר תוצאות בקרת האיכות היו בטווח המקובל.
השחזוריות הבין-בדיקתית והבין-מרכזית הוערכו לפני ניתוח המודל הסופי באמצעות מנות סרום גישור עיוורות וחומרי בקרת איכות שגרתיים. סט גישור של 30 מנות סרום עיוורות נבדק בשני זרימות העבודה של המעבדה. מקדמי השונות הבין-מרכזיים היו 4.8% עבור CEA, 5.6% עבור SCC-Ag, ו-5.1% עבור CA125. מקדמי המתאם התוך-מעמדיים התואמים היו 0.94 עבור CEA, 0.92 עבור SCC-Ag, ו-0.93 עבור CA125. לא בוצע כיול מחדש של הבדיקה עבור מרכז ספציפי, כיוון ששום סטייה שיטתית לא חרגה מטווח בקרת האיכות המעבדתי שהוגדר מראש. ערכי הביומארקרים הסופיים הותאמו באמצעות סטנדרטיזציה של יחידות וסטנדרטיזציה סטטיסטית המבוססת על קבוצת האימון לפני בניית המודל.
הגדרה משולבת של סממנים ביולוגיים ובניית מודל
פאנל הסמנים הביולוגיים המשולב בנסיוב (serum) כלל את CEA, SCC-Ag ו-CA125. לצורך ניתוח קשר קליני-פתולוגי והסתרטות הישרדות לפי Kaplan-Meier, כל סמן סווג כמעלה או לא מעלה בהתאם לסף הקליני שהוגדר מראש. עומס גבוה של סמנים ביולוגיים משולבים הוגדר כעלייה של לפחות שניים מתוך שלושת הסמנים. לצורך מידול אבחנתי ופרוגנוסטי, שלושת הסמנים הביולוגיים נשמרו כמשתנים רציפים, הומרו באמצעות לוגריתם טבעי, ונסטנדרו באמצעות הממוצע וסטיית התקן מקבוצת האימון. אותם פרמטרי המרה הוחלו על קבוצת התיקוף החיצונית.
הביומארקרים הסטנדרטיים שעברו טרנספורמציה לוגריתמית חושבו באמצעות הפרמטרים הבאים של קהלת האימון:
zlnCEA = (lnCEA - 1.42)/0.71
zlnSCC - Ag = (lnSCC - Ag - 0.28)/0.57
zlnCA125 = (lnCA125 - 3.05)/0.62
הפאנל האבחנתי המשולב נבנה באמצעות רגרסיה לוגיסטית רב-משתנית בקבוצת האימון. המנבא הליניארי האבחנתי חושב באופן הבא:
LPdiagnostic = -0.46 + 0.59 × zlnCEA + 1.05 × zlnSCC - Ag+0.42 × zlnCA125
ההסתברות החזויה האישית ל-ESCC חושבה באופן הבא:
PESCC = 1/1+exp(-LPdiagnostic)
נקודת החיתוך האבחנתית האופטימלית עבור הפאנל המשולב הייתה 0.57, שנבחרה בקבוצת האימון באמצעות מדד Youden והוחלה ללא שינוי בתיקוף חיצוני.
מדד הביומארקר הפרוגנוסטי נבנה באמצעות מודל סיכונים פרופורציונליים של Cox. מודל שלב ה-TNM לבדו שימש כמשווה עבור השלב האנטומי. מודל הביומארקר השתמש בעומס ביומארקרים משולב. המודל הקליני-פתולוגי השתמש בגיל, אורך הגידול, דיפרנציאציה היסטולוגית, שלב pT, מצב קשרי הלימפה ואופן הטיפול הראשוני. המודל המשולב הוסיף עומס ביומארקרים משולב גבוה למודל הקליני-פתולוגי. עבור OS, החזאי הליניארי הפרוגנוסטי המשולב חושב באופן הבא:
LPOS = 0.33 × גיל ≥65 years +0.36 × אורך הגידול ≥5cm + 0.25 × בידול דל + 0.31 × pT3
- שלב 4 + 0.66 × גרורות בקשריות לימפה
-0.58 × טיפול מבוסס ניתוח רדיקלי
+0.81 × עומס ביומרקרים משולב גבוה
עבור PFS, המנבא הליניארי הפרוגנוסטי המשולב חושב באופן הבא:
LPPFS = 0.27 × age ≥ 65years + 0.31 × tumorlength ≥ 5cm + 0.24 × диффереנציאציה לקויה +0.25 × pT3
-4stage + 0.55 × גרורות בקשרי לימפה
- 0.46 × טיפול המבוסס על ניתוח רדיקלי
+0.71 × עומס גבוה של סמנים ביולוגיים משולבים
כל המשתנים הבינאריים במשוואות הפרוגנוסטיות קודדו כ-1 כאשר היו נוכחים וכ-0 כאשר היו חסרים. כל המקדמים, פרמטרי הסטנדרטיזציה ונקודות הסף נאמדו בקוהורט האימון והוחלו ישירות על קוהורט התיקוף החיצוני ללא אופטימיזציה מחדש.
נקודות קצה אבחנאיות ופרוגנוסטיות
האבחנה הפתולוגית שימשה כסטנדרט הייחוס האבחנתי. נקודת הקצה האבחנתית העיקרית הייתה ההבחנה בין ESCC לבקרי non-ESCC, כולל מחלה ורידית של הוושט ובקרים בריאים. נקודות קצה אבחנתיות משניות כללו הבחנה בין ESCC למחלה ורידית של הוושט, בין ESCC לבקרים בריאים, ובין ESCC בשלב מוקדם לבקרים שאינם סרטניים. ESCC בשלב מוקדם הוגדר כמחלה בשלבי TNM I–II. הביצועים האבחנתיים הוערכו באמצעות עקומות מאפיין של מקלט (ROC), שטח תחת העקומה (AUC), רגישות, סגוליות, ערך ניבוי חיובי, ערך ניבוי שלילי, דיוק, מדד F1, עקומות דיוק- recall, עקומות כיול, מדדי Brier, חיתוכי כיול, שיפורי כיול וניתוח עקומות החלטה.
ניתוח פרוגנוסטי בוצע רק בחולים עם ESCC. נקודת הסיום הפרוגנוסטית הראשונית הייתה OS, שהוגדרה כפרק הזמן מהתאריך של הטיפול האנטי-גידולי המכריע הראשון ועד המוות מכל סיבה שהיא. נקודת הסיום הפרוגנוסטית המשנית הייתה PFS, שהוגדרה כפרק הזמן מהטיפול הראשון ועד להתקדמות תיועדה של המחלה, חזרת המחלה או מוות. עבור חולים שעברו ניתוח רדיקלי, חזרה מקומית, גרורות מרוחקות או מוות הקשור לגידול נספרו כאירוע PFS כאשר לא נותר נגע לאחר הניתוח. עבור חולים שקיבלו כימורדיותרפיה מכריעה או טיפול פליאטיבי, התקדמות המחלה נקבעה על סמך ממצאי הדמיה, הערכה אנדוסקופית והערכה קלינית.
נתוני מעקב נאספו באמצעות ביקורי מרפאה, רשומות אשפוז, רשומות מעקב קהילתיות, רשומות הפניות וקשר טלפוני. המטופלים היו במעקב מדי 3 חודשים במהלך השנתיים הראשונות לאחר הטיפול, מדי 6 חודשים מהשנה השלישית עד החמישית, ולאחר מכן פעם בשנה. מטופלים שלא חוו מוות, הישנות או התקדמות המחלה עד לתאריך המעקב הסופי סווגו כנצנרים (censored) במועד הקשר האחרון המאושר שלהם. אירועי תוצאה נסקרו באופן עצמאי על ידי שני חוקרים. כאשר תאריכי האירועים היו שונים בין המקורות, נבדקו דוחות דימות, רשומות שחרור, דוחות אנדוסקופיה, מסמכי הפניה ורישומי תעודות פטירה כדי לקבוע את תאריך האירוע הסופי. תאריך הסיום של המעקב היה ה-31 בדצמבר 2025, עבור שני המרכזים המשתתפים.
ניתוח סטטיסטי
הניתוח הסטטיסטי בוצע באמצעות R 4.3.2 ו-SPSS 27.0. משתנים רציפים סוכמו כממוצע ± סטיית תקן או חציון (טווח בין-רבעוני), בהתאם להתפלגות. משתנים קטגוריאליים סוכמו כמספרים ואחוזים. השוואות בין-קבוצתיות בוצעו באמצעות t-test למדגמים בלתי תלויים, מבחן Mann–Whitney U, ניתוח שונות חד-כיווני (one-way ANOVA), מבחן Kruskal–Wallis, מבחן chi-square או מבחן Fisher's exact, לפי העניין. נתונים חסרים הוערכו לפני בניית המודל. משתתפים שחסרו סטטוס אבחנתי מרכזי, ערכי ביומרקרים לפני טיפול, שלב TNM או נתוני הישרדות/התקדמות הדרכו הוצאו מהניתוח בהתאם לקריטריוני ההתאמה שהוגדרו מראש. משתנים נלווים שאינם מרכזיים עם שיעור חסר נמוך מ-5.0% טופלו באמצעות imputing מרובה לפי משוואות שרשראתיות עם 20 מערכי נתונים מושלמים ו-20 איטרציות17. מודלי ה-imputation כללו קבוצה אבחנתית, מרכז, גיל, מין, היסטוריית עישון, שימוש באלכוהול, ערכי ביומרקרים, שלב TNM, שיטת טיפול, סטטוס OS, סטטוס PFS וזמן מעקב. ניתוחים של מקרים מלאים (complete-case analyses) בוצעו כניתוחי רגישות.
המודל האבחנתי פותח בקוהורט האימון. עקומות מאפייני תגובה של מקלט (ROC) שימשו לחישוב השטח תחת העקומה, רגישות, סגוליות, ערך ניבוי חיובי, ערך ניבוי שלילי, דיוק ומדד F1. שטחים תחת העקומה הושוו באמצעות מבחן DeLong18. עקומות דיוק-החזרה (Precision-recall) שימשו כהערכה אבחנתית משלימה. כיול נבחן באמצעות עקומות כיול, מדדי Brier, נקודות חיתוך של כיול ושיפורי כיול. תועלת קלינית הוערכה באמצעות ניתוח עקומות החלטה19. תיקוף פנימי בוצע באמצעות 1,000 דגימות חוזרות בשיטת bootstrap להערכת אפ்டימיזם מתוקן של ההבחנה והכיול. תיקוף חיצוני בוצע על ידי החלת המשוואה של קוהורט האימון, פרמטרי הסטנדרטיזציה וערך הסף ישירות על מרכז B ללא התאמה מחדש.
לניתוח פרוגנוסטי, נוצרו עקומות Kaplan–Meier כדי להשוות בין OS ל-PFS בין קבוצות בעלות עומס ביומרקרים שונה, והשתמשו במבחן log-rank כדי להעריך הבדלים בהישרדות. רגרסיית סיכונים פרופורציונלית של Cox שימשה להערכת יחסי סיכון (hazard ratios) ומרווחי בר סמך של 95%. משתנים בעלי רלוונטיות קלינית ומשתנים הקשורים לתוצאות בניתוח חד-משתני נשקלו עבור מודלים רב-משתניים. ביצועי המודל הוערכו באמצעות מדד Harrell’s C-index, עקומות מאפייני תפעול של מקלט (ROC) תלויות-זמן, עקומות כיול, מדדי Brier וניתוח עקומות החלטה20. כיול פרוגנוסטי הוערך בנקודות זמן של 1, 3 ו-5 שנים. הנחת הסיכונים הפרופורציונליים הוערכה באמצעות שאריות Schoenfeld ועלילות הישרדות של log-minus-log21. ערכי P במבחן Schoenfeld הגלובלי היו 0.42 עבור מודל ה-OS ו-0.37 עבור מודל ה-PFS. מולטי-קולנאריות (Multicollinearity) הוערכה באמצעות גורמי הגדלת שונות (variance inflation factors) לפני פרשנות המודל הסופי, וכל המשתנים המועמדים במודלים האבחנתיים והפרוגנוסטיים הסופיים היו בעלי גורמי הגדלת שונות נמוכים מ-2.1. ערך פרוגנוסטי מוסף מעבר לשלב TNM ומודלים קלינו-פתולוגיים קונבנציונליים הוערך באמצעות שיפור בסיווג נקי (net reclassification improvement) ושיפור באפליקציה אינטגרטיבית (integrated discrimination improvement). פרמטרי המודל שהוערכו בקוהורט האימון הוחלו על קוהורט התיקוף החיצוני ללא התאמה מחדש22. כל המבחנים היו דו-זנביים, וערך P < 0.05 נחשב למשמעותי סטטיסטית.