2025-11-24T16:10:25.080119

Using Information Geometry to Characterize Higher-Order Interactions in EEG

Albers, Marriott, Tatsuno
In neuroscience, methods from information geometry (IG) have been successfully applied in the modelling of binary vectors from spike train data, using the orthogonal decomposition of the Kullback-Leibler divergence and mutual information to isolate different orders of interaction between neurons. While spike train data is well-approximated with a binary model, here we apply these IG methods to data from electroencephalography (EEG), a continuous signal requiring appropriate discretization strategies. We developed and compared three different binarization methods and used them to identify third-order interactions in an experiment involving imagined motor movements. The statistical significance of these interactions was assessed using phase-randomized surrogate data that eliminated higher-order dependencies while preserving the spectral characteristics of the original signals. We validated our approach by implementing known second- and third-order dependencies in a forward model and quantified information attenuation at different steps of the analysis. This revealed that the greatest loss in information occurred when going from the idealized binary case to enforcing these dependencies using oscillatory signals. When applied to the real EEG dataset, our analysis detected statistically significant third-order interactions during the task condition despite the relatively sparse data (45 trials per condition). This work demonstrates that IG methods can successfully extract genuine higher-order dependencies from continuous neural recordings when paired with appropriate binarization schemes.
academic

استخدام الهندسة المعلوماتية لتوصيف التفاعلات من الرتبة الأعلى في تخطيط كهربائية الدماغ

المعلومات الأساسية

  • معرّف البحث: 2510.14188
  • العنوان: استخدام الهندسة المعلوماتية لتوصيف التفاعلات من الرتبة الأعلى في تخطيط كهربائية الدماغ
  • المؤلفون: إريك ألبرز، بول ماريوت، ماسامي تاتسونو
  • التصنيف: q-bio.NC (الخلايا العصبية والإدراك)، q-bio.QM (الطرق الكمية)
  • تاريخ النشر: 16 أكتوبر 2025 (نسخة أولية على arXiv)
  • رابط البحث: https://arxiv.org/abs/2510.14188

الملخص

يوسع هذا البحث طرق الهندسة المعلوماتية (Information Geometry, IG) من بيانات تسلسلات الشوكة الثنائية التقليدية إلى تحليل إشارات تخطيط كهربائية الدماغ (EEG) المستمرة. من خلال التحليل المتعامد لتباعد كولباك-لايبلر والمعلومات المتبادلة، يتم تحديد التفاعلات المختلفة من الرتب بين الخلايا العصبية. طور البحث ثلاث طرق ثنائية التقسيم لتحديد التفاعلات من الرتبة الثالثة في تجارب الخيال الحركي، واستخدم بيانات بديلة مع عشوائية الطور لتقييم الأهمية الإحصائية. من خلال التحقق من النموذج الأمامي، تم تحديد كمية تناقص المعلومات في كل خطوة من خطوات التحليل. تظهر النتائج أنه على الرغم من ندرة البيانات نسبياً (45 محاولة لكل حالة)، يمكن للطريقة اكتشاف تفاعلات من الرتبة الثالثة ذات دلالة إحصائية في ظروف المهمة.

الخلفية البحثية والدافع

تعريف المشكلة

يركز البحث العصبي التقليدي بشكل أساسي على العلاقات الثنائية بين مناطق الدماغ (التفاعلات من الرتبة الثانية)، لكن الدماغ كنظام معقد قد يحتوي على تفاعلات من رتبة أعلى تتجاوز العلاقات الثنائية. قد تفشل شبكات الاتصال الوظيفي الحالية المبنية على الارتباطات الثنائية في التقاط التعقيد الكامل لمعالجة المعلومات في الدماغ.

الأهمية

  1. الأهمية النظرية: فهم ما إذا كان الدماغ يحتاج إلى تفاعلات من الرتبة الثالثة أو أعلى لإكمال الوظائف الإدراكية
  2. الأهمية المنهجية: توسيع طرق الهندسة المعلوماتية من بيانات الشوكة المنفصلة إلى إشارات تخطيط كهربائية الدماغ المستمرة
  3. القيمة التطبيقية: توفير أدوات تحليل جديدة لواجهات الدماغ والحاسوب وتشخيص الأمراض العصبية

حدود الطرق الموجودة

  1. طرق الهندسة المعلوماتية: تطبق بشكل أساسي على بيانات الشوكة الثنائية، وتفتقر إلى استراتيجيات تقسيم فعالة للإشارات المستمرة
  2. تحليل تخطيط كهربائية الدماغ التقليدي: يعتمد بشكل أساسي على الارتباطات الثنائية، متجاهلاً العلاقات الاعتماد من الرتبة الأعلى
  3. الاستدلال الإحصائي: في ظروف البيانات النادرة، قد لا تكون الأدوات المقاربة القياسية (مثل توزيع χ²) مناسبة

الدافع البحثي

توسيع طرق الهندسة المعلوماتية الناجحة المطبقة على تحليل الشوكة إلى بيانات تخطيط كهربائية الدماغ، وتطوير استراتيجيات تقسيم مناسبة لالتقاط التفاعلات الحقيقية من الرتبة الأعلى في التسجيلات العصبية المستمرة.

المساهمات الأساسية

  1. الابتكار المنهجي: تطوير ثلاث طرق تقسيم (Sign و Diff و Power) لتحويل إشارات تخطيط كهربائية الدماغ المستمرة إلى تمثيلات ثنائية مناسبة لتحليل الهندسة المعلوماتية
  2. إطار التحقق: إنشاء طريقة اختبار الأهمية الإحصائية بناءً على بيانات بديلة مع عشوائية الطور
  3. النمذجة الأمامية: تنفيذ نماذج أمامية بعلاقات اعتماد معروفة من الرتبة الثانية والثالثة، مع تحديد كمية تناقص المعلومات خلال عملية التحليل
  4. النتائج التجريبية: اكتشاف تفاعلات من الرتبة الثالثة ذات دلالة إحصائية في بيانات تخطيط كهربائية الدماغ للخيال الحركي
  5. الرؤى النظرية: الكشف عن أقصى خسارة معلومات تحدث عند الانتقال من الحالات الثنائية المثالية إلى تنفيذ العلاقات في الإشارات المتذبذبة

شرح الطريقة

تعريف المهمة

الإدخال: إشارات تخطيط كهربائية الدماغ متعددة القنوات المستمرة الإخراج: مكونات المعلومات المتبادلة من الرتبة الأولى والثانية والثالثة بين ثلاثيات القنوات القيود: التعامل مع البيانات النادرة (45 محاولة/حالة) وتحديات تقسيم الإشارات المستمرة

الأساس النظري للهندسة المعلوماتية

بالنسبة لثلاثة متغيرات ثنائية X₁ و X₂ و X₃، يمكن تمثيل التوزيع الاحتمالي المشترك كمتجه من ثماني احتمالات:

p = (p₀₀₀, p₀₀₁, p₀₁₀, p₀₁₁, p₁₀₀, p₁₀₁, p₁₁₀, p₁₁₁)

نظام الإحداثيات المعاملة المتوقعة η:

  • η₁ و η₂ و η₃: معدلات التنشيط الهامشية
  • η₁₂ و η₁₃ و η₂₃: معدلات التنشيط الثنائية
  • η₁₂₃: معدل التنشيط الثلاثي

نظام الإحداثيات المعاملة الطبيعية θ يُعرّف من خلال النسب اللوغاريتمية، مثل:

θ₁₂₃ = log(p₀₀₁p₀₁₀p₁₀₀p₁₁₁)/(p₁₁₀p₁₀₁p₀₁₁p₀₀₀)

التحليل المتعامد لتباعد كولباك-لايبلر

باستخدام نظام إحداثيات مختلط، يمكن تحليل تباعد كولباك-لايبلر بشكل متعامد إلى:

D[p : q] = D[p : p̄] + D[p̄ : p̃] + D[p̃ : q]

حيث:

  • Dp : p̄: معلومات التفاعل الثلاثي
  • Dp̄ : p̃: معلومات التفاعل الثنائي
  • Dp̃ : q: معلومات تعديل معدل التنشيط

طرق التقسيم

1. طريقة Sign

binary_signal = 1 if EEG_signal > 0 else 0

تلتقط معلومات الطور الخشنة، متجاهلة السعة.

2. طريقة Diff

diff_signal = diff(EEG_signal)
binary_signal = 1 if diff_signal > 0 else 0

تلتقط أنماط تحولات الطور.

3. طريقة Power

power = EEG_signal²
envelope = moving_average(power, 30_samples)
z_scores = (envelope - mean) / std
binary_signal = 1 if z_scores > 1 else 0

تلتقط فترات السعة العالية، بغض النظر عن الطور.

اختبار الأهمية الإحصائية

استخدام إحصائية الاختبار:

λ = 2N·D[p : p̄] ~ χ²(1)

نظراً لندرة البيانات، فإن التقريب χ² ليس جيداً، لذا يتم استخدام اختبار غير معاملي بناءً على بيانات بديلة IAAFT (تحويل فورييه المعدل التكراري للسعة).

إعداد التجربة

مجموعة البيانات

مجموعة بيانات الخيال الحركي من OpenNeuro (Triana-Guzman et al., 2022):

  • المشاركون: 32 متطوعاً صحياً (16 امرأة)
  • الأقطاب الكهربائية: 17 قطباً كهربائياً، موضوعة وفقاً لنظام 10-20 الدولي
  • معدل العينات: 250 هرتز
  • تصميم التجربة:
    • 6 كتل (3 جالسة، 3 واقفة)
    • 30 محاولة لكل كتلة (15 خيال حركي، 15 حالة خمول)
    • إجمالي 45 محاولة لكل حالة

هيكل المحاولة:

  1. التثبيت (4 ثوان): النظر إلى الصليب على الشاشة
  2. الملاحظة (3 ثوان): عرض المهمة المراد تنفيذها
  3. الخيال (4 ثوان): تنفيذ المهمة العقلية (الخيال الحركي أو حالة الخمول)
  4. الراحة (4 ثوان): النشاط الحر

معالجة البيانات المسبقة

  1. التصفية: تصفية عالية التمرير 0.5 هرتز، تصفية شق 58-62 هرتز
  2. إزالة الأثر: استخدام طريقة ASR (إعادة بناء الفضاء الفرعي للأثر)
  3. تصفية النطاق الترددي: تقسيم إلى Delta (0.5-4 هرتز) و Theta (4-8 هرتز) و Alpha (8-12 هرتز) و Beta (12-30 هرتز) و Gamma (30-60 هرتز)
  4. استخراج الفترة: فترة 11 ثانية من 7 ثوان قبل بدء مهمة الخيال إلى 4 ثوان بعد البدء

مقاييس التقييم

  • المعلومات المتبادلة من الرتبة الأولى (I₁): معلومات تعديل معدل التنشيط
  • المعلومات المتبادلة من الرتبة الثانية (I₂): معلومات التفاعل الثنائي
  • المعلومات المتبادلة من الرتبة الثالثة (I₃): معلومات التفاعل الثلاثي
  • الأهمية الإحصائية: p < 0.01 (بناءً على بيانات بديلة IAAFT)

الطرق المقارنة

  1. بيانات بديلة بضوضاء بيضاء: تأثير البنية العشوائية البحتة
  2. بيانات بديلة IAAFT: الحفاظ على الطيف الترددي وتوزيع السعة، مع عشوائية الطور
  3. طرق تقسيم مختلفة: مقارنة Sign مقابل Diff مقابل Power

نتائج التجربة

النتائج الرئيسية

التحقق من البيانات البديلة

  1. بيانات IAAFT تنتج قيماً معلوماتية أعلى من الضوضاء البيضاء، وهذا متوقع لأن IAAFT تحافظ على الاختلافات الطيفية بين مراحل المحاولة
  2. طريقة Power تظهر اتجاهاً تناقصياً I₁ > I₂ > I₃ في جميع النطاقات الترددية
  3. طرق Sign و Diff تظهر انحيازاً نحو I₂، خاصة في النطاقات عالية التردد، مما يحد من قدرتها على التقاط التفاعلات الثلاثية

نتائج بيانات الخيال الحركي

  1. فشل التقريب χ²: نظراً لندرة البيانات (45 محاولة)، التوزيع المقارب القياسي غير مناسب
  2. تفاعلات ثلاثية ذات دلالة إحصائية: اكتشاف I₃ ذات دلالة إحصائية في مراحل الملاحظة والخيال
  3. التحكم في الإيجابيات الكاذبة: حوالي 1% من الثلاثيات ذات دلالة إحصائية في مرحلة التثبيت، مما يتحقق من ملاءمة الفرضية الصفرية
  4. الديناميكيات الزمنية: ديناميكيات زمنية مختلفة لمعلومات الرتبة الثالثة عبر النطاقات الترددية والثلاثيات المختلفة

نتائج التحقق من النموذج الأمامي

تحديد كمية خسارة المعلومات

  1. أقصى خسارة معلومات: تحدث عند الانتقال من الإشارات الثنائية المثالية إلى الإشارات المتذبذبة (حوالي 50% خسارة I₃)
  2. تأثير النقل الحجمي أقل: خسارة المعلومات من إشارات المصدر إلى الأقطاب الكهربائية على فروة الرأس نسبياً صغيرة
  3. حساسية الضوضاء: تحت نسبة الإشارة إلى الضوضاء المتوسطة، ينخفض كل من I₂ و I₃ بشكل كبير

تنفيذ العلاقات الاعتماد

تم تنفيذ بنجاح علاقات اعتماد معروفة من الرتبة الثانية والثالثة في إشارات المصدر المتذبذبة:

  • الحالة الثنائية: الإشارة المستهدفة مرتبطة، بغض النظر عن حالة الإشارة الضابطة
  • الحالة الثلاثية: الإشارة المستهدفة مرتبطة عندما تكون الإشارة الضابطة عالية (1)، وفي ارتباط معاكس عندما تكون منخفضة (0)

النتائج الخاصة بالنطاق الترددي

  • Delta و Theta: الانحياز نحو I₂ في طرق Sign و Diff أقل وضوحاً
  • Alpha والنطاقات الأعلى: طرق Sign و Diff تظهر انحيازاً كبيراً نحو I₂، مما يحد من اكتشاف I₃
  • جميع النطاقات: طريقة Power تحافظ على هيكل هرمي معقول من I₁ > I₂ > I₃

الأعمال ذات الصلة

تطبيقات الهندسة المعلوماتية في علم الأعصاب

  • Amari & Nagaoka (2000): النظرية الأساسية للهندسة المعلوماتية
  • Nakahara & Amari (2002): قياسات الهندسة المعلوماتية لشوكة الأعصاب
  • Tatsuno et al. (2009): التقدير القوي لقوة الاتصال والمدخلات الخارجية

طرق تحليل تخطيط كهربائية الدماغ

  • الطرق التقليدية: تركز بشكل أساسي على تحليل الطيف الترددي والارتباطات الثنائية
  • الاتصال الوظيفي: العلاقات بين المناطق بناءً على الاعتماد الإحصائي كما عرّفه Friston (1995)
  • تحليل الشبكات: تحليل شبكات الدماغ المعقدة بواسطة Bullmore & Sporns (2009)

دراسات التفاعلات من الرتبة الأعلى

  • Battiston et al. (2020, 2021): بنى الشبكات والديناميكيات التي تتجاوز التفاعلات الثنائية
  • مساهمة هذا البحث: أول تطبيق منهجي لطرق الهندسة المعلوماتية على تحليل التفاعلات من الرتبة الأعلى في بيانات تخطيط كهربائية الدماغ

الخلاصة والمناقشة

الاستنتاجات الرئيسية

  1. جدوى الطريقة: يمكن توسيع طرق الهندسة المعلوماتية بنجاح لتحليل إشارات تخطيط كهربائية الدماغ المستمرة
  2. أهمية استراتيجية التقسيم: طريقة Power هي الأنسب لاكتشاف التفاعلات من الرتبة الأعلى
  3. التفاعلات الحقيقية من الرتبة الأعلى: اكتشاف تفاعلات من الرتبة الثالثة ذات دلالة إحصائية في مهام الخيال الحركي
  4. آليات خسارة المعلومات: تحدث خسارة المعلومات الرئيسية أثناء عملية التحويل من الثنائي إلى الإشارات المتذبذبة

القيود

  1. التعقيد الحسابي: 17 قناة قريبة بالفعل من حد الجدوى، وقد تكون المصفوفات عالية الكثافة (128-256 قناة) صعبة حسابياً
  2. الدقة الزمنية: نافذة غير متداخلة بمدة ثانية واحدة توفر ديناميكيات زمنية أقل دقة
  3. التحليل داخل النطاق: يأخذ في الاعتبار فقط التفاعلات داخل نفس النطاق الترددي، لم يتناول التحليل عبر النطاقات
  4. قيود التقسيم: قد يفوت أنماط تفاعل غير خطية أكثر تعقيداً

الاتجاهات المستقبلية

  1. التقسيم الهجين: دمج طرق تقسيم مختلفة لاكتشاف ظواهر مثل الاقتران بين الطور والسعة
  2. النوافذ التكيفية: استخدام نوافذ متداخلة متكيفة مع التردد لتحسين الدقة الزمنية
  3. التحليل عبر النطاقات: توسيع إلى التفاعلات من الرتبة الأعلى بين النطاقات الترددية المختلفة
  4. التفاعلات من رتبة أعلى: استكشاف أنماط التفاعل من الرتبة الرابعة وما فوق

التقييم المتعمق

المزايا

  1. الابتكار المنهجي: توسيع ناجح لطرق الهندسة المعلوماتية من المجال المنفصل إلى المجال المستمر
  2. التحقق الصارم: توفير إطار تحقق شامل من خلال النماذج الأمامية والبيانات البديلة
  3. القيمة العملية: توفير أدوات قابلة للتطبيق لتحليل التفاعلات من الرتبة الأعلى في بيانات تخطيط كهربائية الدماغ
  4. المساهمة النظرية: تحديد كمية خسارة المعلومات في كل خطوة من خطوات التحليل

أوجه القصور

  1. قيود حجم العينة: 45 محاولة نسبياً قليلة، قد تؤثر على القوة الإحصائية
  2. تبسيط التقسيم: تحويل الإشارات المستمرة المعقدة إلى ثنائية قد يفقد معلومات مهمة
  3. قابلية التوسع الحسابية: لم يتم حل التحديات الحسابية لمصفوفات تخطيط كهربائية الدماغ عالية الكثافة بشكل كافٍ
  4. التفسير البيولوجي: النقاش حول الأهمية العصبية البيولوجية للتفاعلات الثلاثية المكتشفة غير كافٍ

التأثير

  1. التأثير المنهجي: توفير أدوات رياضية جديدة لتحليل الرتبة الأعلى للإشارات العصبية
  2. الآفاق التطبيقية: يمكن تطبيقها على واجهات الدماغ والحاسوب وتشخيص الأمراض العصبية وغيرها
  3. القيمة النظرية: تعزيز فهمنا للمبادئ التنظيمية للشبكات الدماغية المعقدة
  4. قابلية الاستنساخ: توفير رموز مفتوحة المصدر ومجموعات بيانات عامة تدعم استنساخ النتائج

السيناريوهات المطبقة

  1. البحث الأساسي في علم الأعصاب: استكشاف مبادئ التنظيم من الرتبة الأعلى لشبكات الدماغ
  2. التطبيقات السريرية: تحليل أنماط الاتصال من الرتبة الأعلى في الأمراض العصبية
  3. واجهات الدماغ والحاسوب: استخراج ميزات إشارة عصبية أكثر ثراءً للتحكم
  4. علم الأعصاب الإدراكي: دراسة التفاعلات العصبية المعقدة في المهام الإدراكية

المراجع

يستشهد هذا البحث بـ 28 مرجعاً مهماً، تشمل بشكل أساسي:

  1. أساسيات الهندسة المعلوماتية: Amari & Nagaoka (2000)، Amari (2001)
  2. التطبيقات العصبية: Nakahara & Amari (2002)، Tatsuno et al. (2009)
  3. منهجية تخطيط كهربائية الدماغ: Delorme & Makeig (2004)، Oostenveld et al. (2011)
  4. الشبكات من الرتبة الأعلى: Battiston et al. (2020, 2021)
  5. مصدر البيانات: Triana-Guzman et al. (2022)

التقييم العام: هذا بحث عالي الجودة في المنهجية، نجح في توسيع نظرية الهندسة المعلوماتية إلى مجال تحليل إشارات تخطيط كهربائية الدماغ. على الرغم من وجود بعض القيود في قابلية التوسع الحسابية والتفسير البيولوجي، فإن إطار التحقق الصارم واستراتيجية التقسيم المبتكرة توفر مساهمات نظرية وعملية مهمة لتحليل التفاعلات من الرتبة الأعلى في الإشارات العصبية.