تصميم الدراسة ومجتمع المرضى
حللت هذه الدراسة الاسترجاعية فحوصات الموجات فوق الصوتية للغدة الدرقية التي أجريت في الفترة ما بين 13 يونيو 2024 و13 يناير 2025. تمت الموافقة على بروتوكول الدراسة من قبل لجنة الأخلاقيات في مستشفى بكين للصداقة، جامعة العاصمة الطبية (رقم الموافقة BFHHZS20240300)، وأجريت وفقاً للمبادئ الأخلاقية الموضحة في إعلان هلسنكي، كما تم التنازل عن الموافقة المستنيرة نظراً للتصميم الاسترجاعي للدراسة. تضمنت مجموعة البيانات 68 عقدة درقية من 63 مريضاً (30 حميدة و38 خبيثة)، وكانت العقدة هي وحدة التحليل. اعتمدت التشخيصات المرجعية على نتائج سحب الخلايا بالإبرة الرفيعة الموجه بالموجات فوق الصوتية (FNAC). وشملت المجموعة الخبيثة 38 عقدة تم تشخيصها على أنها سرطان الغدة الدرقية الحليمي بناءً على نتائج FNAC.
معايير الاشتمال
تطلب من المرضى المؤهلين استيفاء جميع المعايير التالية: (1) الخضوع لتصوير بالموجات فوق الصوتية التقليدية بتدرج الرمادي ودوبلر، يليه تصوير بالموجات فوق الصوتية المعززة بالتباين (CEUS) بجودة صور مرضية تسمح بإجراء تحليل لاحق لتصوير التدفق الوعائي الدقيق؛ (2) وجود تشخيص خلوي مرضي لسرطان الغدة الدرقية الحليمي (PTC) مؤكد عن طريق خزعة بالإبرة الدقيقة، مع أو بدون اختبار إيجابي مصاحب لطفرة BRAFV600E؛ أو وجود نتيجة خلوية من الفئة الثالثة حسب تصنيف بيثيسدا (خلايا غير نمطية ذات أهمية غير محددة)، ولكن مع اختبار إيجابي متزامن لطفرة BRAFV600E؛ (3) وجود تشخيص خلوي مرضي لعقيدات حميدة تكاثرية، أو عقيدات غدية، أو عقيدات جريبية حميدة ضمن الفئة الثانية من تصنيف بيثيسدا دون وجود اختبار إيجابي مصاحب لطفرة BRAFV600E.
معايير الاستبعاد
تم استبعاد المرضى من الدراسة في حال وجود أي من الشروط التالية: (1) صور CEUS ذات جودة رديئة أو إشارة فقاعات دقيقة غير كافية تمنع التقييم الموثوق للتدفق الوعائي الدقيق؛ (2) غياب النتائج الخلوية أو المرضية؛ (3) تشخيص نسيجي مرضي لنوع فرعي نادر أو خاص من سرطان الغدة الدرقية (على سبيل المثال، المتغيرات غير المتمايزة)؛ أو (4) نتائج خلوية تشير إلى ورم غدي جريبي (فئة Bethesda IV) أو أي آفة ذات أصل جريبي غير محدد، بغض النظر عن حالة الطفرات؛ (5) التهاب الغدة الدرقية هاشيموتو
جمع بيانات الموجات فوق الصوتية، والتصوير بالموجات فوق الصوتية المعزز بالتباين (CEUS)، والتصوير بالموجات فوق الصوتية ذات الدقة المكانية العالية (SRUS)
أُجريت جميع الفحوصات باستخدام نظام الموجات فوق الصوتية المذكور ومحول طاقة ذو مصفوفة خطية. وُضع المرضى في وضعية الاستلقاء مع تمديد الرقبة. تم تحديد موقع العقدة المستهدفة وقياسها باستخدام الموجات فوق الصوتية بنمط B-mode؛ كما تم تسجيل التكلسات الدقيقة، واستُخدم مستوى التصوير نفسه المتمركز حول الآفة لتقييم التوعية الدموية داخل العقدة باستخدام Doppler الملون.
بعد إنشاء الوصول الوريدي، تم تحويل النظام إلى وضع CEUS/URM ذي المؤشر الميكانيكي المنخفض (يرمز URM إلى تصوير المجهر فائق الدقة). تم إجراء CEUS عن طريق حقن جرعة بولس وريدية من SonoVue بحجم 1.2 mL، تلاها مباشرةً غسل بمحلول ملحي بحجم 5ml؛ وقد تم تشغيل المؤقت على الشاشة وبدء الاستحواذ السينمائي المستمر مع إعطاء جرعة البولس. تم الحفاظ على المسبار في مستوى ثابت مع ضغط أدنى، وطُلب من المريض تجنب البلع أثناء مرحلتي التدفق الداخل والتدفق الخارج.
بالنسبة لتصوير الموجات فوق الصوتية المعزز بالتباين (CEUS) الكمي، تم استخدام منطقة اهتمام مقتصرة على الآفة لاستخلاص معاملات الوقت والكثافة. أما بالنسبة لتصوير الموجات فوق الصوتية فائقة الدقة (SRUS)، فقد قام سير عمل مجهر الموجات فوق الصوتية (URM) بتحديد وتتبع إشارات الفقاعات الدقيقة بعد التحكم في الحركة، ونتج عن ذلك قياسات نسبة الأوعية، والتعقيد، وكثافة الأوعية الدقيقة، ومؤشر التروية، وسرعة التدفق. وأظهرت صادرات الصور التمثيلية إعدادات VSP 4 وRES 2 وCTR 3 وSM 2 وVEN 3 وCPT 10 s؛ بينما لم تتوفر الإعدادات المقابلة لبقية الفحوصات.
تُدرج متغيرات مدخلات النموذج الـ 25 في الجدول 1 ومصنفة حسب طريقة الحصول عليها: العمر والجنس؛ التكلسات الدقيقة في وضع B-mode؛ التروية الوعائية داخل العقدة باستخدام Doppler الملون؛ التصوير بالموجات فوق الصوتية المعزز بالتباين (CEUS) النوعي؛ التصوير بالموجات فوق الصوتية المعزز بالتباين (CEUS) الكمي؛ و11 قياساً للأوعية الدقيقة باستخدام SRUS.
اختيار الميزات
تم الاحتفاظ بجميع السمات الكمية الـ 25 في نماذج تعلم الآلة. واستُخدمت السمات الرقمية فقط؛ حيث استُبعدت أسماء المرضى، وأرقام التسجيل، وحقول حجم الآفة. لم يتم إجراء أي اختيار للسمات بناءً على البيانات، وأُدخلت نفس المتنبئات المحددة مسبقاً في كل مصنف.
تقييس الميزات
لم يتم تطبيق أي تحويل أو تعويض أو تغيير في المقياس العام قبل إجراء التحقق المتقاطع. وقد طُبقت المعايرة فقط على آلة متجه الدعم ذات دالة الأساس الشعاعي (radial-basis-function SVM) عبر خط معالجة StandardScaler. تم ضبط المقياس بناءً على العقد التدريبية لكل طية، ثم طُبق على العقد التحققية لتلك الطية. أما المصنفات القائمة على الأشجار فقد استلمت المقاييس الرقمية الأصلية:

حيث x هي قيمة الميزة الأصلية، وµ هو المتوسط الحسابي لطية التدريب، وσ هو الانحراف المعياري لطية التدريب. وقد منعت المعالجة المسبقة على مستوى الطية ملاحظات التحقق من المساهمة في معاملات قياس SVM.
تقسيم البيانات
استخدم تقييم الأداء الأولي التحقق المتقاطع من نوع StratifiedGroupKFold بخمس طيات مع ضبط shuffle = True و random_state = 42. حدد معرف المريض 63 مجموعة، وتم تعيين جميع العقد من المريض نفسه إلى الطية ذاتها. لم يساهم أي مريض بعقد في كل من مجموعتي التدريب والتحقق ضمن الطية الواحدة.
تدريب النموذج والتحقق من صحته
بروتوكول التدريب
تم تقييم خمسة مصنفات: الغابة العشوائية (100 شجرة؛ random_state = 42)، وآلة ناقلات الدعم ذات الدالة الشعاعية الأساسية (C = 1.0؛ gamma = scale؛ probability = True؛ مسار عمل StandardScaler؛ random_state = 42)، وشجرة القرار (معيار Gini؛ عمق غير مقيد؛ random_state = 42)، وXGBoost (100 مقدر؛ learning_rate = 0.3؛ max_depth = 6؛ subsample = 1.0؛ colsample_bytree = 1.0؛ random_state = 42)، وتعزيز التدرج (100 مقدر؛ learning_rate = 0.1؛ max_depth = 3؛ random_state = 42). لم يتم إجراء أي بحث شبكي، أو تحسين بايزي، أو ضبط للعتبة، أو اختيار نموذج متداخل.
التحقق المتقاطع
داخل كل طية من الطيات الخمس المقسمة حسب المرضى، تم تدريب النماذج على مجموعات المرضى المتبقية وتقييمها على المجموعات المستبعدة. وتم حساب الدقة، والحساسية، والنوعية، والضبط، ومقياس F1-score، ومساحة منحنى خصائص تشغيل المستقبِل (ROC-AUC) لكل طية، وتم تلخيصها في صورة المتوسط ± الانحراف المعياري (SD). كما تم تجميع تنبؤات خارج العينة (OOF) عبر جميع العقد الـ 68 لإنشاء منحنى ROC واحد مصادق عليه تقاطعيًا ومصفوفة ارتباك واحدة لكل نموذج.
تقييم الأداء
مقاييس التقييم
تم تقييم أداء النموذج ضمن كل طية تحقق مجمعة حسب المرضى ومن خلال تنبؤات OOF المجمعة باستخدام المقاييس التالية:
الدقة: النسبة الإجمالية للتنبؤات الصحيحة

الحساسية (الاسترجاع): نسبة العقد الخبيثة الفعلية التي تم تحديدها بشكل صحيح

النوعية: نسبة العقد الحميدة الفعلية التي تم تحديدها بشكل صحيح

الدقة (القيمة التنبؤية الإيجابية): نسبة الحالات التي تم التنبؤ بأنها خبيثة وكانت خبيثة بالفعل

مقياس F1-Score: المتوسط التوافقي للدقة والاسترجاع

المساحة تحت منحنى خصائص تشغيل المستقبل (ROC-AUC): مقياس لقدرة النموذج على التمييز بين العقد الحميدة والخبيثة عبر جميع عتبات التصنيف
حيث أن TP = الإيجابيات الحقيقية (العقيدات الخبيثة التي تم تحديدها بشكل صحيح)، و TN = السلبيات الحقيقية (العقيدات الحميدة التي تم تحديدها بشكل صحيح)، و FP = الإيجابيات الكاذبة (العقيدات الحميدة التي تم تصنيفها خطأً على أنها خبيثة)، و FN = السلبيات الكاذبة (العقيدات الخبيثة التي تم تصنيفها خطأً على أنها حميدة).
تحليل مصفوفة الارتباك
تم إنشاء مصفوفات الارتباك (confusion matrices) بنظام OOF من خلال تجميع التنبؤات التي تم إجراؤها لكل عقدة فقط في الطية (fold) التي استُبعد فيها ذلك المريض. وبالتالي، حصلت كل عقدة على تنبؤ من نموذج تم تدريبه على عقد من مرضى آخرين.
تحليل أهمية الميزات
بالنسبة لنموذج الغابة العشوائية (Random Forest)، تم حساب درجات أهمية الميزات بناءً على متوسط الانخفاض في شوائب جيني (Gini impurity) عبر جميع أشجار القرار. وتم تحديد وترتيب أهم 15 ميزة لتحديد معاملات الأوعية الدموية الدقيقة التي ساهمت بشكل أكبر في أداء التصنيف.
تحليل SHAP الاستكشافي
أُجري تحليل OOF SHAP مركّز باستخدام مفسر قائم على التبديل لنموذج SVM. وبالنسبة لكل عقيدة مستبعدة، استُخدمت الملاحظات المقابلة في طية التدريب فقط كتوزيع خلفي (128 تبديلاً متضاداً؛ البذرة العشوائية = 20260716). لخصت قيم SHAP المطلقة المتوسطة حجم المساهمات، بينما أشارت القيم ذات الإشارة إلى الاتجاه. كان هذا التحليل استكشافياً ولم يُستخدم لاستنتاج السببية، أو تحديد المؤشرات الحيوية المستقلة، أو تحديد العتبات السريرية.
التحليل الإحصائي
التحليل المقارن
تم تلخيص أداء النموذج وصفياً عبر طيات التحقق الخمس المجمعة حسب المرضى وفي تنبؤات OOF المجمعة. لم يتم إجراء اختبار فرضية رسمي بين النماذج أو ترتيب قائم على قيمة P لأن الطيات مرتبطة وحجم المجموعة صغير.
لإجراء مقارنات المجموعات المرجعية، تم تقييم التوزيع الطبيعي للمتغيرات المستمرة ضمن كل مجموعة من النتائج باستخدام اختبار Shapiro-Wilk. واستُخدم اختبار Welch's t عندما كانت كلتا المجموعتين متوافقتين مع التوزيع الطبيعي؛ وخلاف ذلك، استُخدم اختبار Mann-Whitney U ثنائي الطرف. أما المتغيرات الفئوية فقد تم تقييمها باستخدام اختبار Pearson's chi-square، مع استخدام اختبار Fisher's exact للجداول 2 × 2 ذات البيانات المحدودة. وكانت قيم P ثنائية الطرف، واستكشافية، وغير معدلة (alpha = 0.05).
قابلية التكرار
تم دعم قابلية التكرار من خلال تعريف محدد لمجموعة المرضى، وبذرة عشوائية ثابتة (42)، وتكوينات مصنف ثابتة، ومعالجة مسبقة لكل طية. استُخدمت المعرفات لأغراض التجميع فقط، ولم يتم إدخالها كمتنبئات أو تصديرها مع مخرجات النموذج.