2025-11-11T16:04:09.679573

A Mixed Model Approach for Estimating Regional Functional Connectivity from Voxel-level BOLD Signals

Liu, Zhang, Tran et al.
Resting-state brain functional connectivity quantifies the synchrony between activity patterns of different brain regions. In functional magnetic resonance imaging (fMRI), each region comprises a set of spatially contiguous voxels at which blood-oxygen-level-dependent signals are acquired. The ubiquitous Correlation of Averages (CA) estimator, and other similar metrics, are computed from spatially aggregated signals within each region, and remain the quantifications of inter-regional connectivity most used by neuroscientists despite their bias that stems from intra-regional correlation and measurement error. We leverage the framework of linear mixed-effects models to isolate different sources of variability in the voxel-level signals, including both inter-regional and intra-regional correlation and measurement error. A novel computational pipeline, focused on subject-level inter-regional correlation parameters of interest, is developed to address the challenges of applying maximum (or restricted maximum) likelihood estimation to such structured, high-dimensional spatiotemporal data. Simulation results demonstrate the reliability of correlation estimates and their large sample standard error approximations, and their superiority relative to CA. The proposed method is applied to two public fMRI data sets. First, we analyze scans of a dead rat to assess false positive performance when connectivity is absent. Second, individual human brain networks are constructed for subjects from a Human Connectome Project test-retest database. Concordance between inter-regional correlation estimates for test-retest scans of the same subject are shown to be higher for the proposed method relative to CA.
academic

نهج نموذج مختلط لتقدير الاتصال الوظيفي الإقليمي من إشارات BOLD على مستوى الفوكسل

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

  • معرّف الورقة: 2211.02192
  • العنوان: نهج نموذج مختلط لتقدير الاتصال الوظيفي الإقليمي من إشارات BOLD على مستوى الفوكسل
  • المؤلفون: Ruobin Liu, Chao Zhang, Chau Tran, Sophie Achard, Wendy Meiring, Alexander Petersen
  • التصنيف: stat.ME (الإحصاء - المنهجية)، stat.AP (الإحصاء - التطبيقات)
  • وقت النشر: نوفمبر 2022 (نسخة arXiv مسبقة، محدثة نوفمبر 2025)
  • رابط الورقة: https://arxiv.org/abs/2211.02192

الملخص

تقدم هذه الورقة طريقة جديدة قائمة على نموذج التأثيرات المختلطة الخطية لمعالجة مشكلة تقدير الاتصال الوظيفي للدماغ في حالة الراحة من التصوير بالرنين المغناطيسي الوظيفي (fMRI). على الرغم من أن مقدّر "الارتباط بين المتوسطات" (Correlation of Averages, CA) التقليدي يُستخدم على نطاق واسع، إلا أنه يعاني من مشاكل الانحياز الناجمة عن الارتباط داخل المنطقة وخطأ القياس. تفصل هذه الورقة من خلال إطار نموذج التأثيرات المختلطة الخطية بين مصادر التباين المختلفة في إشارات الفوكسل، بما في ذلك الارتباط بين المناطق والارتباط داخل المنطقة وخطأ القياس. طورت الدراسة خط أنابيب حسابي جديد يركز على تقدير معاملات الارتباط بين المناطق على المستوى الفردي، مع استخدام طريقة الاحتمالية العظمى للتعامل مع تحديات البيانات الزمكانية عالية الأبعاد. أثبتت نتائج المحاكاة موثوقية تقدير الارتباط وتفوقها على طريقة CA.

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

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

  1. المشكلة الأساسية: كيفية تقدير الاتصال الوظيفي بين مناطق الدماغ بدقة من إشارات BOLD على مستوى الفوكسل، مع تجنب مشاكل الانحياز في الطرق التقليدية
  2. التحديات التقنية:
    • تتمتع بيانات fMRI بتبعيات زمكانية معقدة
    • عدد الفوكسلات أكبر بكثير من البعد الزمني، مما يؤدي إلى تحديات حسابية
    • الارتباط المكاني داخل المنطقة وخطأ القياس يؤثران على دقة تقدير الاتصالية

أهمية البحث

  • الاتصالية الوظيفية هي أساس دراسة الأمراض التنكسية العصبية واضطرابات الوعي
  • يعتبر تقدير الاتصالية الدقيق حاسماً للبحث العصبي على المستويات الفردية والجماعية
  • قد يؤثر الانحياز في الطرق الموجودة على تشخيص الأمراض وأبحاث تمييز الخصائص الفردية

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

المشاكل الموجودة في مقدّر CA التقليدي:

  1. مشكلة الانحياز: يتم تقييدها بالارتباط داخل المنطقة αⱼ ونسبة الإشارة الضوضائية βⱼ، مما يؤدي إلى انحياز التقديرات نحو 0
  2. تجاهل التبعية: لا تأخذ في الاعتبار التبعيات الزمكانية لإشارات الفوكسل
  3. قيود المعاملات: تتأثر معاملات الاتصالية بخطة العينة والضوضاء الآلية، وتفتقر إلى الجوهرية

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

  1. اقتراح نموذج إحصائي جديد: إطار عمل لنمذجة إشارات BOLD على مستوى الفوكسل قائم على نموذج التأثيرات المختلطة الخطية، يميز بوضوح بين التباين بين المناطق والتباين داخل المنطقة
  2. تطوير طريقة تقدير فعالة:
    • استراتيجية تقدير على مرحلتين، تجمع بين تقدير الاحتمالية المقيدة (ReML)
    • تطبيق لأول مرة لطريقة تقريب احتمالية Vecchia في نمذجة الاتصال الوظيفي
  3. ضمانات نظرية: توفير الخصائص الكبيرة للعينة والنظرية الإحصائية المقاربة للمقدّرات
  4. التحقق التجريبي: التحقق من تفوق الطريقة على البيانات المحاكاة والبيانات الحقيقية (مسح الفئران الميتة، بيانات الاختبار-إعادة الاختبار من HCP)

شرح الطريقة

تعريف المهمة

الإدخال: معاملات الموجة الصغيرة لإشارات BOLD على مستوى الفوكسل Xⱼₗₘ، حيث j=1,...,J لمناطق الدماغ، l=1,...,Lⱼ للفوكسلات، m=1,...,M لمعاملات الموجة الصغيرة الإخراج: معاملات الارتباط بين المناطق ρⱼⱼ'، المستخدمة لبناء شبكة الاتصال الوظيفي القيود: التعامل مع الجدوى الحسابية لبيانات زمكانية عالية الأبعاد

بنية النموذج

نموذج التأثيرات المختلطة BOLD

النموذج الأساسي هو:

X = Zμ + Uη + γ + ε

حيث:

  • μⱼ: التأثير الثابت الإقليمي (متوسط المنطقة)
  • ηⱼₘ: التأثير العشوائي الإقليمي (يحفز التبعية بين المناطق)
  • γⱼₗₘ: التأثير العشوائي على مستوى الفوكسل (يحفز التبعية داخل المنطقة)
  • εⱼₗₘ: خطأ القياس

معاملات هيكل التباين المشترك

  1. الارتباط بين المناطق: Var(η) = (SRS) ⊗ A، حيث R = {ρⱼⱼ'} هي مصفوفة الارتباط المستهدفة
  2. الهيكل داخل المنطقة: Λⱼ = Cⱼ ⊗ Bⱼ (التباين المشترك المكاني-الزمني القابل للفصل)
  3. تعريف دوال النواة:
    • النواة المكانية: نواة Matérn K(d; ν, φ)
    • النواة الزمنية: النواة الغاوسية H(|m-m'|; τ)

نقاط الابتكار التقني

استراتيجية التقدير على مرحلتين

المرحلة 1: تقدير المعاملات الخاصة بالمنطقة

  • استخدام ReML لتقدير معاملات كل منطقة θⱼ = kᵧⱼ, σ²ᵧⱼ, φᵧⱼ, τᵧⱼ
  • القضاء على تأثير التأثيرات الإقليمية من خلال الاحتمالية المقيدة

المرحلة 2: تقدير المعاملات العامة والبين-إقليمية

  • تقدير معاملات الارتباط بين المناطق θ = τη, kη, ρ₁₂, σ²η
  • تثبيت تقديرات المرحلة 1، والتركيز على معاملات الاتصالية

تقريب Vecchia

لحل مشكلة التعقيد الحسابي (الوقت O(N³)، الذاكرة O(N²))، يتم استخدام تقريب احتمالية Vecchia:

p(X) ≈ p(X_π(1)) ∏ᵢ₌₂ᴺ p(X_π(i) | X_π(j), j ∈ Jᵢ)

من خلال مجموعة شرطية صغيرة |Jᵢ|=100 لتحقيق تحسن الكفاءة الحسابية

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

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

  1. البيانات المحاكاة:
    • J=3 مناطق دماغية، M=60 معامل موجة صغيرة
    • استخدام الإحداثيات المكانية من تجارب الفئران الحية (L₁=41, L₂=25, L₃=77 فوكسل)
    • تباين قوة الإشارة δⱼ ∈ {0.1, 0.5, 0.7} والتباين المشترك المكاني ψⱼ ∈ {0.2, 0.5, 0.8}
  2. البيانات الحقيقية:
    • بيانات مسح الفئران الميتة (التحقق من معدل الإيجابيات الكاذبة)
    • قاعدة بيانات الاختبار-إعادة الاختبار من HCP (42 موضوع، J=92 منطقة شبكة الوضع الافتراضي)

مؤشرات التقييم

  1. المحاكاة: متوسط الخطأ التربيعي (MSE)، متوسط الانحراف المطلق (MAD)
  2. بيانات HCP: معامل الارتباط المتطابق (CCC) لتقييم موثوقية الاختبار-إعادة الاختبار
  3. بيانات الفئران الميتة: تحليل معدل الإيجابيات الكاذبة

طرق المقارنة

  1. ρ̂CA: مقدّر الارتباط بين المتوسطات التقليدي
  2. ρ̂EBLUE: الارتباط القائم على مقدّر أفضل خطي غير متحيز تجريبي
  3. ρ̂ReML: مقدّر ReML الكامل
  4. ρ̂Vecchia: مقدّر تقريب Vecchia

تفاصيل التنفيذ

  • دوال النواة: النواة الغاوسية H(u;τ) = exp(-τ²u²/2)، نواة Matérn-5/2
  • التحسين: طريقة L-BFGS شبه النيوتونية
  • حجم مجموعة شرط Vecchia: |Jᵢ| = 100
  • اختبار الدلالة: إجراء Benjamini-Yekutieli، FDR < 0.2

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

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

أداء التجارب المحاكاة

  1. الدقة: ρ̂ReML له أصغر انحراف معياري في جميع الإعدادات، والوسيط الأقرب للقيمة الحقيقية
  2. تحليل الانحياز:
    • عندما ρ=0.6، يكون CA و EBLUE منحازين بشكل كبير نحو 0
    • عند التباين المشترك المكاني العالي (ψ=0.8)، يزداد انحياز CA و EBLUE
    • يحافظ ρ̂ReML على الاستقرار في جميع الإعدادات

مثال على النتائج الرقمية

عند قوة إشارة متوسطة (δ=0.5):

  • التباين المشترك المكاني المنخفض (ψ=0.2): MSE لـ ρ̂ReML هو 0.008-0.025، أقل بشكل كبير من CA بـ 0.016-0.033
  • التباين المشترك المكاني العالي (ψ=0.8): الفرق أكثر وضوحاً، MSE لـ ρ̂ReML هو 0.012-0.028، CA هو 0.056-0.194

التجارب الاستئصالية

  1. التحقق من تقريب Vecchia: أداء ρ̂Vecchia و ρ̂ReML متطابقة تقريباً، مما يثبت فعالية طريقة التقريب
  2. قوة تحديد النموذج الخاطئ: يظل ρ̂ReML متفوقاً على الطرق التقليدية حتى تحت هياكل التباين المشترك البديلة
  3. مقارنة المقدّر الأوراكل: يكون فقدان الأداء في الإجراء على مرحلتين ضئيلاً

نتائج البيانات الحقيقية

تحليل مسح الفئران الميتة

  • التحكم في الإيجابيات الكاذبة: لا توجد حواف ذات دلالة إحصائية لـ ρ̂Vecchia عند مستوى الدلالة 5%، بينما لا تزال طريقة CA لديها حواف ذات دلالة
  • التحكم في FDR: كلا الطريقتين لا تحتويان على حواف ذات دلالة بعد التعديل BY (q<0.2)، وهو متوقع

تحليل الاختبار-إعادة الاختبار من HCP

  • تحسن الاتساق: في جميع استراتيجيات بناء الرسم البياني، يكون لدى معظم المواضيع CCC أعلى باستخدام طريقة النموذج المختلط
  • نسبة الحواف: من اختيار الحواف بنسبة 1%-20%، تتفوق طريقة النموذج المختلط باستمرار على CA
  • الدلالة الإحصائية: في أفضل 10% من الحواف، يُظهر حوالي 60-80% من المواضيع اتساقاً أعلى في الاختبار-إعادة الاختبار

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

الاتجاهات البحثية الرئيسية

  1. النمذجة على مستوى الفوكسل: نموذج التأثيرات المختلطة المرتبط بالمهمة من Woolrich et al. (2004)
  2. الاتصالية على مستوى المجموعة: النموذج الهرمي البايزي من Bowman et al. (2008)
  3. الطرق في مجال التردد: نموذج التأثيرات المختلطة في مجال التردد من Kang et al. (2012)
  4. النمذجة الزمكانية: طريقة عملية VAR من Castruccio et al. (2018)

مزايا هذه الورقة

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

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

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

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

القيود

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

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

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

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

المزايا

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

أوجه القصور

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

التأثير

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

السيناريوهات المعمول بها

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

المراجع

تستشهد الورقة بـ 63 مرجعاً ذا صلة، تشمل بشكل أساسي:

  • Achard et al. (2023): التحليل النظري لمقدّرات الارتباط بين المناطق
  • Vecchia (1988): طريقة تقريب الاحتمالية للعمليات المكانية
  • Bowman et al. (2008): النمذجة الهرمية البايزية لبيانات fMRI
  • Kang et al. (2012, 2017): نماذج التأثيرات المختلطة الزمكانية
  • Castruccio et al. (2018): النماذج الزمكانية متعددة الدقة

تقدم هذه الورقة من خلال النمذجة الإحصائية الصارمة والابتكار الحسابي مساهمة منهجية مهمة لتحليل الاتصال الوظيفي fMRI، وتتمتع بقيمة أكاديمية وعملية عالية.