Sensitivity analysis for causal mediation: bridge score, sharp sensitivity bounds, and calibration

ملخص تقني: تحليل الحساسية للوساطة السببية عبر درجة الجسر (Bridge Score)، والحدود الحادة، والمعايرة

1. بيان المشكلة

يقوم تحليل الوساطة السببية بتفكيك آثار المعالجة الإجمالية إلى آثار مباشرة طبيعية (NDE) وآثار غير مباشرة طبيعية (NIE). ويعتمد تحديد هذه الآثار عادةً على فرضية التجاهل المتسلسل (sequential ignorability)، والتي تتطلب أن يكون الوسيط (mediator) غير قابل للتجاهل في ظل وجود المعالجة والمتغيرات الأساسية. وهذه الفرضية غير قابلة للاختبار؛ إذ إن أي عامل غير مقاس يدفع بالوسيط والنتيجة معاً ينتهك هذه الفرضية.

بينما يوجد تحليل حساسية لتجاهل الوسيط، حدد المؤلفون ثلاثة تحديات إحصائية محددة في الأدبيات المتعلقة بالتحليل النوعي للحساسية على المقياس الجمعي (additive scale) (أي فرق المتوسط)، وهو المقياس الأكثر شيوعاً في البحث العلمي:

1. القابلية للتفسير مقابل التوطين (Localization): المعلمات القياسية الموجودة (مثل الارتباطات المتبقية) تؤدي إلى انهيار الانتهاكات التي تختلف عبر أذرع المعالجة وملفات تعريف المتغيرات. وعلى العكس من ذلك، فإن دوال الحساسية الموطنة بالكامل يصعب استخراجها ومعايرتها في حالات المتغيرات عالية الأبعاد. هناك نقص في وجود اختزال عام منخفض الأبعاد يحافظ على معلومات الوسيط دون اشتراط التوصيف الكامل للمتغيرات.
2. الافتقار إلى الحدود الجمعية الحادة: توجد حدود غير معلمية حادة لمقياس نسبة المخاطر (المناسب للنتائج الثنائية)، ولكن لا يتوفر نظير لها في شكل غلاف جمعي حاد (mean difference). وتعتمد النهج الجمعية الحالية على صيغ تحيز معلمية مقيدة أو حدود "مانسكي" (Manski-style) غير معلوماتية.
3. المعايرة الهشة: تعتمد المعايرة التشغيلية للمعايير المستمرة غالباً على وحدات أو عرض دعم تعسفي. علاوة على ذلك، فإن الملخصات المتماسكة مع المقياس (مثل قيم E-values) تؤدي إلى انهيار بنية المعلمتين في قيمة قياسية واحدة، مما يفقدها القابلية للتفسير.

2. المنهجية

درجة الجسر (The Bridge Score)
يقدم المؤلفون درجة الجسر، وهي متجه ثنائي الأبعاد يتم تعريفه عند قيمة وسيط مشتركة $m$ ومتغير $x$:
$$B(m, x) = (f_0(m | x), f_1(m | x))$$
حيث $f_a(m|x)$ هي كثافة الوسيط تحت ذراع المعالجة $a$.

• خاصية التوازن: بشرط درجة الجسر $B(m)$، تكون المتغيرات الملحوظة $X$ متوازنة فيما يتعلق بحدث الوسيط $M=m$. وهذا يسمح لدرجة الجسر بأن تعمل كدرجة موازنة لمرحلة الوسيط، على غرار درجة الاحتمالية (propensity score) لتخصيص المعالجة.
• اختزال الأبعاد: تقلل درجة الجسر من كائن الحساسية من دالة لمتجه المتغيرات الكامل $X$ إلى دالة لدرجة الجسر منخفضة الأبعاد $B(m)$، مع الاحتفاظ بقوانين الوسيط الخاصة بكل ذراع واللازمة لتحديد العالم المتقاطع (cross-world identification).

الحدود الحادة للحساسية
يعرف المؤلفون دالة حساسية درجة الجسر $\Delta_a(m, b)$ كخلل (confounding) متبقي مستحث عن طريق الشرط على الوسيط الملحوظ ضمن الذراع $a$، بعد الشرط على درجة الجسر.
تحت فرضية وجود مسبب خفي قياسي ($U$)، يشتق المؤلفون غلافاً جمعياً نقطياً حاداً لـ $\Delta_a(m, b)$:
$$|\Delta_a(m, b)| \leq \Xi_a(m, b) = \eta_a(m, b) \frac{\gamma_a(m, b) - 1}{\gamma_a(m, b)}$$
حيث:

• $\gamma_a(m, b)$ هو نسبة الاختيار المتبقية، وتقيس الميل في توزيع المسبب الخفي الناتج عن الشرط على $M=m$.
• $\eta_a(m, b)$ هو معلم مدى النتيجة، ويقيس انتشار متوسط النتيجة مشروطاً بالمسبب الخفي ضمن طبقة درجة الجسر.

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

نهج المعايرة
لجعل الحدود عملية، تم اقتراح استراتيجيتين للمعايرة:

1. معايرة المعيار (Benchmark Calibration): يتم معايرة المسبب الخفي ضد متغير ملحوظ $W$. يقدم المؤلفون نسخة قائمة على الرتب من هذه المعايرة. من خلال تحويل $W$ إلى رتبته الكسرية، يصبح الغلاف الناتج غير متغير تجاه التعبيرات الرتيبة لـ $W$ وقوياً ضد التوزيعات ذات الذيول الثقيلة، مع الحفاظ على البنية ذات المعلمتين ($\eta, \gamma$).
2. معايرة ميزانية المتبقي (Residual Budget Calibration): يقوم هذا النهج بتغيير حجم معلمة مدى النتيجة $\eta_a$ بناءً على التباين المتبقي للنتيجة بعد التعديل للوسيط ودرجة الجسر. يتحكم معلم يحدده المستخدم $k_a \in [0,1]$ في حصة التباين المتبقي المنسوبة للخلط الخفي، بينما يتحكم قيمة الشبكة $g_a$ في نسبة الاختيار.

خوارزمية الاستدلال
يتم تشغيل الإطار عبر خوارزمية حساب g-البيزي (Bayesian g-computation). يعمل هذا النهج على:

• نشر عدم اليقين من نموذج البيانات الملحوظة (كثافات الوسيط وانحدارات النتيجة).
• دمج عدم اليقين من كميات المعايرة (المعيار أو ميزانية المتبقي).
• استخدام اختزال دالي قياسي لطي الدالة اللانهائية الأبعاد للحساسية إلى معلمات قياسية خاصة بكل ذراع ($\bar{\Delta}_0, \bar{\Delta}_1$)، مما يسمح بفرض توزيعات قبلية محدودة الأبعاد على تصحيحات الحساسية.
• توليد عينات خلفية (posterior draws) للـ NDE والـ NIE تأخذ في الاعتبار جميع مصادر عدم اليقين بشكل كامل.

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

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

4. النتاليات والتوضيح

يوضح البحث الإطار باستخدام تجربة التأطيم للهجرة لـ Brader et al. (2008).

• الإعداد: المعالجة ($A$) هي نبرة سلبية/تأطير لاتيني؛ الوسيط ($M$) هو استجابة عاطفية سلبية؛ النتيجة ($Y$) هي الضرر المتصور من الهجرة.
• النتائج:
  • أنتجت معايرة المعيار القائمة على الرتب أغلفة أضيق من نسخة المقياس الخام، مما أظهر استقرار التحويل الرتبي ضد عرض الدعم التجريبي.
  • أظهرت معايرة ميزانية المتبقي أن الغلاف يتسع بشكل رتيب مع زيادة قوة الاختيار ($g$) وميزانية مدى النتيجة ($k$).
  • أظهر التحليل أن استنتاجات المتانة مدفوعة بعوامل تضخيم قابلة للتفسير (على سبيل المثال، "المسبب الخفي أقوى بـ 1.5 مرة من الدخل") بدلاً من معلمات النماذج الغامضة.
  • نجح حساب g-البيزي في نشر عدم اليقين من تقدير درجة الجسر والمعايرة إلى التوزيعات الخلفية للـ NDE والـ NIE.

5. الأهمية والادعاءات

يضع المؤلفون هذا العمل كـ أداة مكملة وليس بديلاً للأساليب الموجودة مثل تحليل $\rho$-المعلمي (Imai et al., 2010) أو حدود نسبة المخاطر (Ding and VanderWeele, 2016).

تكمن أهمية الإطار في قدرته على:

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

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