يظهر مخطط عام لسير عمل المعالجة المسبقة لتخطيط كهربية الدماغ (EEG)، واستخراج الميزات، والتصنيف في الشكل 1. والمواد والبرمجيات المستخدمة في الدراسة مدرجة في جدول المواد.

الشكل 1: سير عمل خط معالجة تخطيط كهربائية الدماغ (EEG) من حيث المعالجة الأولية، واستخراج الميزات، والتصنيف. نظرة عامة تخطيطية لسير عمل الدراسة، بما في ذلك التحقق من بيانات EEG، والترشيح بنطاق تمرير (band-pass filtering) وإعادة المرجعية، واستخراج ميزات جذر متوسط المربعات (RMS)، وكثافة القدرة الطيفية (PSD)، وميزات الإنتروبيا، وإضافة العمر، وتقسيم البيانات إلى مجموعات تدريب واختبار، وتطوير نموذج الغابة العشوائية (Random Forest)، وتحليل استقرار الميزات، وتحليل SHAP، وتحليل حجم التأثير باستخدام معامل Cohen’s d، والاختيار المتكامل للميزات، واستئصال ميزات الإنتروبيا، وإعادة تدريب النموذج، وتقييم الأداء. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
1. الحصول على البيانات
تم الحصول على مجموعة بيانات تخطيط كهربائية الدماغ (EEG) في حالة الراحة المتاحة للعموم من مستودع OpenNeuro (مجموعة البيانات ds004504، الإصدار 1.0.9)33. وتم التحقق من تسميات المشاركين والمعلومات الديموغرافية. تكونت مجموعة البيانات من 88 مشاركاً، من بينهم 36 مشاركاً مصاباً بمرض الزهايمر (AD)، و23 مصاباً بالخرف الجبهي الصدغي (FTD)، و29 من الضوابط الصحية (HC).
2. استيراد تسجيلات تخطيط كهربائية الدماغ (EEG)
بالنسبة لكل مشارك، تم تحميل تسجيل تخطيط كهربائية الدماغ باستخدام مكتبة MNE-Python. تم فحص كل ملف EEG للتأكد من إمكانية الوصول إليه، وإمكانية استيراده بنجاح، واحتوائه على معرف صالح للمشارك. وقد تم استبعاد التسجيلات التي كانت ملفاتها مفقودة أو تالفة. تم الإبقاء على تردد أخذ العينات الأصلي لـ EEG وهو 500 Hz في هذه المرحلة. ثم خُفض معدل أخذ العينات للإشارات لتقليل الحمل الحسابي مع الحفاظ على جميع معلومات EEG ذات الصلة سريريًا بهذه الدراسة.
3. الترشيح الممرر للنطاق
تم ترشيح تسجيلات EEG المستمرة باستخدام مرشح ممرر للنطاق من الدرجة الرابعة، مع تردد قطع منخفض قدره 0.5 Hz وتردد قطع مرتفع قدره 40 Hz.
4. المرجعية المتوسطة الشائعة وتقسيم الفترات الزمنية
تم تطبيق المرجعية المتوسطة الشائعة (CAR) على تسجيلات EEG المُرشحة. وعند كل عينة زمنية، تم حساب متوسط الإشارة عبر جميع قنوات EEG المتاحة وطرحه من كل قناة. وقد تم الاحتفاظ بإشارات EEG المرجعية للتحليلات اللاحقة. ثم قُسّم كل تسجيل EEG مرجعي إلى فترات زمنية (epochs) ثابتة الطول باستخدام دالة make_fixed_length_epochs() في مكتبة MNE-Python، حيث استُخدمت فترة زمنية مدتها 10 s مع تداخل قدره 0 s.
5. التحقق من صحة الحقبات الزمنية
بالنسبة لكل مشارك، تم الاحتفاظ بالحقبات الزمنية الكاملة التي تبلغ مدتها 10s فقط. وتم استبعاد أي جزء متبقٍ من تخطيط كهربائية الدماغ (EEG) تقل مدته عن 10s في نهاية التسجيل. واستُخدمت كل حقبة زمنية تم الاحتفاظ بها لاحقاً كعينة فردية لاستخراج الميزات.
6. تحضير الميزات
على الرغم من توفر بيانات فحص الحالة العقلية المصغر (MMSE) في مجموعة البيانات المصدر، إلا أنه تم استبعادها من التحليل الحالي. دُمجت ميزات تخطيط كهربائية الدماغ (EEG) المستخلصة مع المعلومات الديموغرافية للمشاركين، بما في ذلك العمر والتصنيفات التشخيصية، والتي تم الحصول عليها من ملف البيانات الوصفية participants.csv. استُخدم معرف المشارك كفتاح مشترك، وتم الاحتفاظ فقط بالمشاركين الذين لديهم سجلات ميزات تخطيط كهربائية الدماغ والمعلومات الديموغرافية المقابلة لها.
تضمنت مصفوفة سمات النطاقات المتعددة خمس سمات لجذر متوسط المربعات (RMS) وهي (delta_rms وtheta_rms وalpha_rms وbeta_rms وgamma_rms)، وخمس سمات لكثافة الطيف القدرية (PSD) وهي (delta_psd وtheta_psd وalpha_psd وbeta_psd وgamma_psd)، وخمس سمات للإنتروبيا وهي (delta_entropy وtheta_entropy وalpha_entropy وbeta_entropy وgamma_entropy). كما أُدرج العمر كمتغير ديموغرافي تكميلي. ونظراً لأن توزيعات الأعمار قد تختلف عبر المجموعات التشخيصية، فلم يكن من الممكن استبعاد التداخل الديموغرافي بشكل كامل. وقد تم تعيين المجموعة التشخيصية، التي عُرِّفت بأنها AD أو HC أو FTD، كملصق مستهدف.
7. المعالجة المسبقة للبيانات وتقسيم بيانات التدريب والاختبار
خلال مرحلة المعالجة المسبقة، تم خفض معدل أخذ العينات للتسجيلات من 500 Hz إلى 250 Hz لتقليل المتطلبات الحسابية مع الاحتفاظ بمعلومات تردد EEG ذات الأهمية.
تم تقسيم مجموعة البيانات إلى مجموعات فرعية للتدريب (80%) والاختبار (20%) على مستوى الشخص باستخدام استراتيجية تقسيم مجمعة. بعد ذلك، تم الحصول على التنبؤات على مستوى الشخص عن طريق التصويت بالأغلبية عبر العصور الزمنية المتوقعة التي تنتمي إلى كل مشارك. اعتمد التقييم الأساسي على تقسيم التدريب والاختبار المجمع على مستوى الشخص لمنع ظهور العصور الزمنية لنفس المشارك في كلتا المجموعتين الفرعيتين. تم الحصول على التنبؤات على مستوى الشخص عن طريق التصويت بالأغلبية عبر العصور الزمنية المتوقعة التي تنتمي إلى كل مشارك. كما تم فحص مجموعة البيانات بحثًا عن القيم المفقودة، وتمت إزالة الملاحظات المفقودة أو تعويضها حسب الاقتضاء. وتم ترميز التسميات التشخيصية. ثم تم ضبط StandardScaler على بيانات التدريب وتطبيقه بعد ذلك على كل من مجموعتي بيانات التدريب والاختبار.
8. تطوير نموذج الغابة العشوائية
تم تهيئة مصنف الغابة العشوائية (Random Forest) باستخدام 80 شجرة، وبحد أقصى لعمق الشجرة يبلغ 10، وبحد أقصى أربع سمات، وبحد أدنى خمس عينات لكل ورقة، وحالة عشوائية قدرها 42. وقد تم تدريب المصنف باستخدام مجموعة بيانات التدريب الموحدة.
9. تقييم النموذج
تَمّ التنبؤ بتصنيفات الفئات لكل من مجموعتي بيانات التدريب والاختبار. كما أُنشئت مصفوفة ارتباك (confusion matrix)، وحُسبت الدقة (accuracy)، والضبط (precision)، والاسترجاع (recall)، ومقياس F1-score، وذلك جنبًا إلى جنب مع تقرير التصنيف. وسُجلت دقة التدريب ودقة الاختبار.
بالنسبة للتقييم الخاص بكل مشارك، تأكد من تخصيص جميع الحقبات (epochs) الخاصة بمشارك معين حصرياً إما لمجموعة التدريب أو لمجموعة الاختبار. وقد تم تدريب مصنف الغابة العشوائية (Random Forest) باستخدام نفس المعلمات الفائقة (hyperparameters) المستخدمة في التحليل الأساسي.
10. التحقق المتبادل
كإجراء إضافي لتقييم النموذج، تم إجراء تحقق متبادل طبقي خماسي الطيات (five-fold stratified cross-validation) باستخدام shuffle = True و random_state = 30. وقد تم حساب متوسط الدقة والانحراف المعياري عبر الطيات الخمس. واعتُبر هذا التحليل منفصلاً عن تقييم الاحتجاز الأساسي على مستوى الأشخاص.
11. تحليل استقرار الميزات
تم تقييم استقرار الميزات عن طريق تكرار تدريب الغابة العشوائية (Random Forest) 10 مرات باستخدام بذور عشوائية من 0 إلى 9. وفي كل عملية تشغيل، تم تسجيل دقة الاختبار ودرجات أهمية الميزة. وحُسب المتوسط والانحراف المعياري لدرجة أهمية الميزة لكل متنبئ عبر عمليات التشغيل العشر، ثم رتبت المتنبئات وفقاً لاستقرارها. استُخدم هذا التحليل لتقييم اتساق تصنيفات الميزات بدلاً من استبدال التقييم الأساسي على مستوى العينة.
12. تحليل الذكاء الاصطناعي القابل للتفسير
طُبق SHAP TreeExplainer على نموذج Random Forest المدرب. حُسبت قيم SHAP لتقدير مساهمة كل ميزة في تنبؤات النموذج. كما حُسب متوسط القيمة المطلقة لـ SHAP لكل ميزة، ورُتبت الميزات وفقاً لمساهماتها في قيم SHAP. تم تحديد الميزات التي أظهرت قيم SHAP منخفضة باستمرار ومقارنتها بنتائج تحليلات استقرار الميزات وCohen’s d. واختيرت الميزات ذات المساهمات المنخفضة باستمرار لإجراء عملية الاستئصال (ablation) وإعادة تدريب النموذج لاحقاً.
13. تحليل حجم التأثير الإحصائي
تم حساب قيمة Cohen’s d لكل مؤشر حيوي من مؤشرات EEG للمقارنات التالية: AD مقابل HC، وAD مقابل FTD، وFTD مقابل HC. وقد تم تفسير مقادير حجم التأثير باستخدام عتبات قدرها 0.20 للتأثير الصغير، و0.50 للتأثير المتوسط، و0.80 للتأثير الكبير.
14. اختيار الميزات المتكامل
تمت مقارنة نتائج تحليل أهمية الميزات باستخدام الغابة العشوائية (Random Forest)، وتحليل SHAP، وتحليل معامل Cohen’s d. وقد تم تحديد المتنبئات التي أظهرت باستمرار أهمية منخفضة للميزات، ومساهمة منخفضة في قيم SHAP، وأحجام تأثير صغيرة كمرشحين للحذف.
15. استئصال السمات
تم تقييم سمات الإنتروبيا في تحليل استئصالي استكشافي، وتم إنشاء مصفوفة سمات مختزلة تحتوي على RMS و PSD والعمر. أُعيد تدريب مصنف الغابة العشوائية (Random Forest) باستخدام نفس المعلمات الفائقة. كما كُررت عمليات تدريب النموذج، والاختبار، والتحقق المتقاطع، وتحليل خصائص تشغيل المستقبل (ROC) باستخدام مجموعة السمات المختزلة.
16. تحليل منحنى خصائص تشغيل المستقبِل (ROC)
تم الحصول على احتمالات الفئة من مصنف الغابة العشوائية (Random Forest) المُحسَّن. كما تم إنشاء منحنيات خصائص تشغيل المستقبل (ROC) متعددة الفئات باستخدام نهج الـ one-مقابل- استراتيجية الراحة. تم حساب قيم المساحة تحت المنحنى (AUC) الخاصة بكل فئة، بالإضافة إلى متوسط المساحة تحت المنحنى.
17. مقارنة الأداء
تمت مقارنة أداء نموذج الميزات الكامل بأداء نموذج الميزات المختزل الذي تم الحصول عليه بعد استئصال ميزات الإنتروبيا.