اختبارات حسن التوفيق غير المقاربة واختيار النموذج في نماذج الكتل العشوائية ذات القيم

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

• معرّف الورقة: 2510.13636
• العنوان: اختبارات حسن التوفيق غير المقاربة واختيار النموذج في نماذج الكتل العشوائية ذات القيم
• المؤلفون: Félix Almendra-Hernández, Miles Bakenhus, Vishesh Karwa, Mitsunori Ogawa, Sonja Petrović
• التصنيفات: stat.ME cs.SI math.ST stat.TH
• تاريخ النشر: 16 أكتوبر 2025
• رابط الورقة: https://arxiv.org/abs/2510.13636

الملخص

تدرس هذه الورقة مشكلة اختبارات حسن التوفيق لنماذج الكتل العشوائية ذات القيم (valued stochastic blockmodel, SBM). تمثل نماذج SBM ذات القيم طريقة عامة لنمذجة بيانات الشبكات، حيث يتم تجميع العقد في كتل، وتُقاس الاتصالات بين العقد بواسطة أعداد أو تسميات. تسمح عائلة النموذج هذه بأنظمة أخذ عينات dyad مختلفة، بما في ذلك SBM الكلاسيكي وSBM بواسون والـ SBM المُسمى، بالإضافة إلى الحالات التي تكون فيها بعض ملاحظات الحافة مرقوبة. تركز الورقة على الاختبارات ذات العينات المحدودة لنماذج SBM غير برنولي، وتشتق الحركات الصريحة للقاعدة ماركوفية اللازمة لأخذ عينات من توزيع مرجعي، وتحدد إحصائيات حسن التوفيق لتحديد ملاءمة النموذج. بالنسبة لـ SBM المُسمى (بما في ذلك نماذج الحافة المرقوبة)، يتم دراسة السلوك المقارب لهذه الإحصائيات.

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

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

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

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

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

دافع البحث

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

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

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

شرح الطريقة

تعريف المهمة

بالنظر إلى شبكة ذات قيم $G = (G_{uv})_{1≤u<v≤n}$، حيث يمثل $G_{uv}$ قوة الاتصال (العد أو التسمية) بين زوج العقد $\{u,v\}$، الهدف هو:

1. اختبار ما إذا كانت الشبكة تتوافق مع نموذج SBM ذي قيم معطى
2. إجراء اختبار ملاءمة النموذج عندما يكون التخصيص للكتل غير معروف
3. اختيار عدد الكتل المناسب $k$

معمارية النموذج

1. تعريف SBM ذي القيم

بالنسبة لـ $n$ عقدة و $k$ كتلة، يفترض SBM ذو القيم:

• الاستقلال الشرطي: $G_{uv} \perp\!\!\!\perp G_{u'v'} | Z$ لأي dyad اثنين
• صيغة العائلة الأسية: $P_θ(G = g | Z = z) = \prod_{1≤u<v≤n} \frac{h(g_{uv})\exp⟨T_z(g_{uv}), θ_{z_uz_v}⟩}{ψ(θ_{z_uz_v})}$

حيث $h$ هي القياس الأساسي، $T_z$ هي الإحصائية الكافية، و $θ$ هو متجه المعاملات.

2. الحالات الخاصة

• SBM الكلاسيكي: $G_{uv} | Z = z \sim \text{Bernoulli}(θ_{z_uz_v})$
• SBM بواسون: $G_{uv} | Z = z \sim \text{Poisson}(θ_{z_uz_v})$
• SBM المُسمى: $G_{uv} | Z = z \sim \text{Multinomial}(N, (θ^{(ℓ)}_{z_uz_v})^L_{ℓ=1})$

3. بناء القاعدة ماركوفية

القاعدة ماركوفية لـ SBM بواسون: $B = \{ε_{uv} - ε_{u'v'} : z_u = z_{u'}, z_v = z_{v'}\}$

القاعدة ماركوفية لـ SBM المُسمى: $B = \{ε^{(ℓ)}_{uv} + ε^{(ℓ')}_{u'v'} - ε^{(ℓ')}_{uv} - ε^{(ℓ)}_{u'v'} : ℓ,ℓ' ∈ [L], z_u = z_{u'}, z_v = z_{v'}\}$

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

1. إطار الاختبار الشرطي

• تعريف الألياف: $F_{z,t} := \{g ∈ G : T_z(g) = t\}$
• التوزيع الشرطي: $P(G = g | T_z(g) = t) = \frac{h(g)}{\sum_{g' ∈ F_{z,t}} h(g')}$
• قيمة p الدقيقة: $p(z,g) = P(\text{GoF}_z(G) ≥ \text{GoF}_z(g) | T_z(g))$

2. الطريقة البايزية الجزئية

بالنسبة لتخصيص الكتل غير المعروف، يتم تعريف قيمة p بايزية جزئية: $p_b(g) = \sum_{z ∈ Z_{n,k}} p(z,g)P(Z = z | g)$

تتعامل هذه الطريقة مع عدم اليقين في تخصيص الكتل من خلال حساب المتوسط على التوزيع اللاحق.

3. إحصائيات حسن التوفيق

SBM بواسون: $\text{GoF}_z(g) = \sum_{u=1}^n \sum_{i=1}^k \frac{(m_{ui} - n_iθ̂_{z_ui})^2}{n_iθ̂_{z_ui}}$

SBM المُسمى: $\text{GoF}_z(g) = χ^2_{BC}(g,z) = \sum_{u=1}^n \sum_{i=1}^k \sum_{ℓ=1}^{L-1} \frac{(m^{(ℓ)}_{ui} - n_iθ̂^{(ℓ)}_{z_ui})^2}{n_iθ̂^{(ℓ)}_{z_ui}}$

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

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

1. البيانات المحاكاة:
   • عدد العقد: $n = 50, 100$
   • 4 مصفوفات اتصال مختلفة $θ^{(1)}, θ^{(2)}, θ^{(3)}, θ^{(4)}$
   • إنشاء 100 رسم بياني لكل إعداد
2. البيانات الحقيقية:
   • شبكة الفطريات الطفيلية (154 عقدة)
   • شبكة أنواع الأشجار (51 عقدة)
   • أوزان الحافة تمثل عدد الأنواع المضيفة/الطفيلية المشتركة

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

1. معدل الخطأ من النوع الأول: معدل رفض الفرضية الصفرية عند مستوى الدلالة 0.05
2. قوة الاختبار: معدل رفض النموذج عند أعداد كتل مختلفة
3. دقة اختيار النموذج: المقارنة مع معيار ICL

طرق المقارنة

• معيار ICL (الاحتمالية المتكاملة الكاملة)
• خوارزمية EM المتغيرة لتقدير تخصيص الكتل
• حزمة sbm في R

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

• طول سلسلة MCMC: يتحكم فيه معامل numGraphs
• تقدير تخصيص الكتل: استخدام خوارزمية EM المتغيرة
• عتبة الاحتمالية اللاحقة: $ε = 1/m$، حيث $m$ هو حجم المجموعة الداعمة

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

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

1. التحقق من القوة ومعدل الخطأ من النوع الأول

النتائج عند $n=50$:

النتائج عند $n=100$:

2. تطبيق البيانات الحقيقية

شبكة أنواع الأشجار:

• الكتل 3-7: قيمة p = 0
• الكتل 8-9: قيمة p = 0.01
• الكتلة 10: قيمة p = 0.19
• الكتل 11-15: قيمة p ≥ 0.68

شبكة الفطريات:

• الكتل 3-17: قيمة p = 0
• الكتل 18-21: قيمة p = 0.01
• الكتلة 22: قيمة p = 0.07

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

1. فعالية الطريقة: معدل الرفض قريب من 1 عند استخدام كتلتين أو 3 كتل، وقريب من 0 عند استخدام 4 أو 5 كتل، وهو ما يتوافق مع التوقعات
2. تأثير حجم العينة: حجم عينة أكبر ($n=100$) يوفر قوة إحصائية أقوى
3. الاختلاف عن الطرق الموجودة:
   • تختار هذه الطريقة 10 كتل لشبكة الأشجار و22 كتلة لشبكة الفطريات
   • يختار معيار ICL 7 كتل لشبكة الأشجار و9 كتل لشبكة الفطريات
   • قد يكون الاختلاف بسبب تحفظ الطريقة والمتطلبات الصارمة لملاءمة النموذج

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

اختبارات حسن التوفيق للشبكات

1. الطرق الطيفية: اختبار حسن التوفيق الطيفي لـ Lei (2016)
2. طرق ERGM: طريقة مقارنة التوزيع المرجعي لـ Hunter وآخرين (2008)
3. الاختبارات المحسّنة: طريقة Hu وآخرين (2021) التي تعالج مباشرة تكاليف الحساب والضمانات النظرية

نماذج الكتل العشوائية

1. SBM الكلاسيكي: تمديد الكتل الكامنة لـ Holland وآخرين (1983)
2. الشبكات ذات القيم: تمديد ERGM للشبكات ذات الأعداد لـ Krivitsky (2012)
3. اختيار النموذج: اختيار النموذج القائم على الاحتمالية لـ Wang و Bickel (2017)

الإحصائيات الجبرية

1. القاعدة ماركوفية: النظرية الأساسية لـ Diaconis و Sturmfels (1998)
2. تطبيقات الشبكات: اختبارات العينات المحدودة لـ SBM برنولي لـ Karwa وآخرين (2023)
3. البناء الديناميكي: طريقة القاعدة ماركوفية الديناميكية لـ Gross وآخرين (2014)

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

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

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

القيود

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

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

1. سلاسل ماركوف المركبة: تطوير طرق قادرة على أخذ عينات من ألياف متعددة
2. تحسين الحساب: تحسين التقارب والكفاءة الحسابية لـ MCMC
3. التطبيقات الموسعة: الدمج مع الشبكات الديناميكية والشبكات متعددة الطبقات

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

المزايا

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

أوجه القصور

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

التأثير

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

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

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

المراجع

تتضمن المراجع الرئيسية:

• Diaconis, P. and Sturmfels, B. (1998). خوارزميات جبرية لأخذ عينات من التوزيعات الشرطية
• Holland, P. W., Laskey, K. B., and Leinhardt, S. (1983). نماذج الكتل العشوائية: الخطوات الأولى
• Karwa, V. et al. (2023). اختبارات حسن التوفيق مونت كارلو لنماذج الكتل العشوائية المصححة بالدرجة والنماذج ذات الصلة
• Wang, Y. X. R. and Bickel, P. J. (2017). اختيار النموذج القائم على الاحتمالية لنماذج الكتل العشوائية