Bayesian inference for partial orders from random linear extensions: power relations from 12th Century Royal Acta

ملخص تقني: الاستدلال البايزي للترتيبات الجزئية من الامتدادات الخطية العشوائية

بيان المشكلة
يتناول البحث تحدياً إحصائياً يتمثل في استنتاج تسلسل هرمي اجتماعي متطور زمنياً من بيانات ترتيب رتبي، حيث لا يكون الهيكل الأساسي بالضرورة ترتيباً كاملاً (كلياً). قام المؤلفون بتحليل 371 قائمة شهود من وثائق "أكتا" (Acta) الملكية الإنجليزية والنورماندية من القرن الثاني عشر، والتي تسجل الترتيب الذي أُدرجت به 67 أسقفاً. يفترض المؤرخون أن هذا الترتيب يعكس تسلسلاً هرمياً اجتماعياً قائماً على المكانة، والأقدمية، والنفوذ السياسي. ومع ذلك، فإن البيانات مشوبة بالضجيج، وغير مكتملة (تظهر مجموعات فرعية من الأساقفة في كل قائمة)، كما أن الهيكل الأساسي قد يحتوي على أزواج غير قابلة للمقارنة بين الفاعلين (ترتيب جزئي)، بدلاً من كونه ترتيباً خطياً واحداً. النماذج الموجودة لبيانات الرتب، مثل نماذج "مالوز" (Mallows) و"بلاكت-لوس" (Plackett-Luce)، تفترض عادةً ترتيباً كلياً أو احتمالات مستقلة عن السياق، وهو ما قد لا يكون مناسباً لهذا السياق التاريخي حيث يُتوقع وجود التعدي (transitivity) ولكن الترتيب يكون جزئياً.

المنهجية
يقترح المؤلفون نموذج ماركوف الخفي (HMM) بايزي، حيث يكون المتغير الكامن هو ترتيب جزئي (poset) يتطور زمنياً، والبيانات المرصودة هي امتدادات خطية عشوائية لهذا الترتب الجزئي.

1. توصيف المتغير الكامن:

• يتم توليد الترتب الجزئي عند الزمن $t$، ويرمز له بـ $h(t)$، من مصفوفة حالة كامنة $Z(t)$ ذات $K$ بُعد.
• لكل أسقف $j$، يتطور متجه سلطة كامن $U^{(t)}_j$ عبر الزمن كنموذج ناقل ذاتي الانحدار (VAR) من الدرجة الأولى.
• يُضاف تأثير متغير مصاحب $\beta_{s_{t,j}}$ بناءً على الأقدمية إلى متجه السلطة لتشكيل متجه الحالة: $Z^{(t)}_j = U^{(t)}_j + \mathbf{1}_K \beta_{s_{t,j}}$.
• يتم تعريف الترتب الجزئي $h(t)$ من خلال تقاطع $K$ من الترتيبات الكلية المستمدة من أعمدة $Z(t)$. يسبق الأسقف $i$ الأسقف $j$ في الترتب الجزئي إذا وفقط إذا كان $Z^{(t)}_{i,k} > Z^{(t)}_{j,k}$ لجميع قيم $k=1,\dots,K$.
• يتم اختيار البعد $K$ ليكون $\lfloor D/2 \rfloor$ (حيث $D$ هو الحد الأقصى لعدد الأساقفة النشطين)، مما يضمن قدرة النموذج على تمثيل أي ترتيب جزئي ممكن (خاصية التمثيل الشامل).

2. نموذج الملاحظة (قفز الدور):

• تُنمذج القوائم المرصودة كامتدادات خطية "مشوشة" للترتيب الجزئي الحقيقي.
• يقدم المؤلفون آلية "قفز الدور": مع احتمال $p$، يتم اختيار الفاعل التالي في القائمة عشوائياً (مع تجاهل قيود الترتب الجزئي)؛ وإلا، يتم اختيار الفاعل وفقاً لامتداد خطي عشوائي للترتيب الفرعي المتبقي.
• يسمح هذا بوجود "ضجيج" حيث يتعارض الترتب المرصود مع الهيكل الأساسي، مع الحفاظ على القيد بأن الامتدادات الخطية الصالحة تكون أكثر احتمالية.
• يتعامل النموذج مع القوائم غير المكتملة (مجموعات فرعية من الأساقفة) من خلال الشرط على الترتب الفرعي المستحث من الفاعلين الحاضرين.

3. النماذج المسبقة والاتساق:

• تم بناء النموذج المسبق لعملية الترتب الجزئي لتكون متسقة هامشياً، مما يعني أن توزيع الترتب الفرعي (المستحث من مجموعة فرعية من الفاعلين) متسق مع توزيع الترتب الكامل؛ ويتحقق ذلك باستخدام نموذج مسبق محدد لمتجهات $U$ ومعامل الارتباط $\rho$.
• يتحكم المعلم الفائق $\rho$ في "عمق" الترتب الجزئي (طول أطول سلسلة). تم تصميم النموذج المسبق ليكون غير متحيز فيما يتعلق بالعمق، لتجنب التركيز على ترتيبات جزئية بعمق 3 كما يحدث في النماذج المسبقة الموحدة.
• تم تقييد تأثيرات الأقدمية $\beta$ لتكون مرتبة ($\beta_1 > \beta_2 > \dots$) لتعكس التوقعات التاريخية بأن الأقدمية ترتبط بالمكانة، رغم أنه تم أيضاً اختبار نماذج غير مقيدة.

4. الاستدلال:

• يتم إجراء الاستدلال البايزي باستخدام طريقة "ماركوف تشين مونت كارلو" (MCMC) مع تحديثات "ميتروبوليس-هاستينغز".
• يتضمن حساب الاحتمالية عملية عد الامتدادات الخطية، وهي مسألة من فئة #P-complete. يستخدم المؤلفون خوارزميات عد تكرارية (تحديداً حزمة lecount) والتي تكون قابلة للتنفيذ للعدد الصغير من الفاعلين ($m \le 22$) في هذه المجموعة من البيانات، لكنها تحد من القدرة على التوسع لمشكلات أكبر.

المساهمات الرئيسية

• فئة نموذج جديدة: يقدم البحث إطاراً عاماً لنمذجة بيانات الرتب المتسلسلة زمنياً حيث يكون الهيكل الأساسي ترتيباً جزئياً، مما يعمم النماذج الموجودة مثل "مالوز" و"بلاكت-لوس".
• الاتساق الهامشي: اشتق المؤلفون وأثبتوا أن نموذجهم للمتغيرات الكامنة للترتيبات الجزئية هو نموذج متسق هامشياً، وهي خاصية غالباً ما تنتهكها النماذج المسبقة المرتجلة للترتيبات الجزئية.
• التمثيل الشامل: أثبتوا أنه مع بُعد كامن $K \ge \lfloor D/2 \rfloor$، يمكن للنموذج تمثيل أي سلسلة زمنية لترتيب جزئي، مما يضمن عدم وجود تحيز هيكلي ضد الهياكل المعقدة.
• التطبيق التاريخي: يعد هذا أول تحليل إحصائي لقوائم شهود "الأكتا" الملكية من القرن الثاني عشر لنمذجة التسلسل الهرمي الاجتماعي صراحةً كترتيب جزئي، مع الفصل بين تأثيرات الأقدمية والسلطة الشخصية.

النتائج

• ملاءمة النموذج: يتفوق النموذج القائم على الترتب الجزئي على البدائل الأبسط (مثل ترتيبات VSP أو ترتيبات الدلاء - Bucket Orders) وعلى نموذج "بلاكت-لوس" المختلط للسلاسل الزمنية من حيث جودة الملاءمة (مقاسة عبر التحقق المتقاطع "Leave-One-Out" وعوامل بايز).
• الأقدمية مقابل السلطة: يؤكد التحليل وجود اعتماد قوي للمكانة على الأقدمية، مع تناقص تأثير الأقدمية مع زيادة رقم الرتبة. ومع ذلك، نجح النموذج في عزل "السلطة" (الانحرافات عن الأقدمية).
• الرؤى التاريخية:
  • الأساقفة النورمان: يرصد النموذج انخفاضاً في مكانة الأساقفة النورمان (مثل أساقفة أفرانش، باييو، وإيفريه) مقارنة بالأساقفة الإنجليز بدءاً من أوائل القرن الثاني عشر، وهو اتجاه لم يتم قياسه كمياً من قبل.
  • النفوذ السياسي: تتذبذب مكانة أساقفة محددين، مثل هنري دي بلوا (شقيق الملك ستيفن) ونايجل من إيلي، بطرق ترتبط بالرضا أو الغضب السياسي، بشكل مستقل عن أقدميتهم.
  • هيكل التسلسل الهرمي: تبين أن ترتيبات الإجماع المستعادة ضحلة نسبياً، وتتكون عادةً من مجموعات قليلة من الأساقفة متوسطي الرتبة مع عدد قليل من الفاعلين فوقهم أو تحتهم، بدلاً من كونها تسلسلاً هرمياً عميقاً وصارماً.

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

يتسم المؤلفون بالتواضع بشأن القدرة على التوسع الحسابي لطريقتهم، حيث يقرون بأن الطبيعة المعقدة (#P-complete) لعد الامتدادات الخطية تحد من النموذج لبيانات تحتوي على عدد صغير من الفاعلين (حوالي 20). ويقترحون أنه بالنسبة لمجموعات البيانات الأكبر، قد تكون الفئات المقيدة من الترتبات الجزئية (مثل VSP أو ترتيبات الدلاء) ضرورية كحلول وسط. كما يشيرون إلى أنه بينما يتناسب نموذجهم جيداً مع البيانات، فإن نموذج "الضجيج" يفترض وقوع الأخطاء في اتجاه واحد (قفز الدور)، وأن العمل المستقبلي قد يستكشف نماذج ضجيج ثنائية الاتجاه.

في النهاية، يجادل البحث بأن افتراض وجود ترتيب كلي (كما في نماذج الترتيب القياسية) هو أمر غير مبرر لهذا النوع من البيانات الاجتماعية، وأن نماذج الترتب الجزئي تقدم تمثيلاً أكثر دقة وقابلية للتفسير للهياكل الاجتماعية التاريخية.