Two-stage Estimation for Causal Inference Involving a Semi-continuous Exposure

ملخص تقني لورقة بحثية بعنوان: "التقدير ثنائي المراحل للاستدلال السببي المتعلق بالتعرضات شبه المستمرة"

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

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

• قصور الطرق الحالية: طُورت طرق درجات الميل (Propensity Score - PS) القياسية بشكل جيد للتعرضات الثنائية أو المستمرة تمامًا. ومع ذلك، فإن تطبيق طرق درجات الميل المعممة (الوزن العكسي للكثافة) على البيانات شبه المستمرة غالبًا ما يؤدي إلى مقدرات ذات تباين عالٍ بسبب الأوزان المتطرفة. علاوة على ذلك، تفتقر نماذج "الجزءين" (two-part modeling) الحالية غالبًا إلى إطار سببي موحد يحدد بوضوح ويعرف مقادير سببية منفصلة لـ حالة التعرض (الشرب مقابل عدم الشرب) و جرعة التعرض (الشدة بين الشاربين).
• مثال توضيحي: ركز المؤلفون على التعرض للكحول أثناء الحمل (PAE) وتأثيره على الإدراك لدى الطفل. تُظهر البيانات كتلة نقطية عند الصفر (الأمهات الممتنعات) وتوزيعًا مستمرًا للشاربات. الهدف البحثي هو التمييز بين الأثر السببي لـ أي شرب مقابل عدم الشرب، وبين أثر الجرعة بين الشاربات.

2. المنهجية

يقترح المؤلفون إطارًا سببيًا عامًا يستخدم استراتيجية تقدير ثنائية المراحل قائمة على هيكل درجة ميل مكونة من جزئين (two-part PS).

أ. الإطار السببي والمقادير (Estimands)
يتم تفكيك التعرض $T$ إلى مؤشر ثنائي $A = I(T > 0)$ وجرعة مستمرة $D = \log(T)$ (عندما يكون $A=1$). يحدد الإطار مقدارين سببيين متميزين ضمن نموذج بنيوي هامشي (MSM):

1. $\Delta_D$ (أثر الجرعة): الأثر السببي لزيادة الجرعة بمقدار وحدة واحدة بين مجتمع المتعرضين ($A=1$).
   $$\Delta_D = E\{Y(1, d+1) - Y(1, d) \mid A=1\}$$
2. $\Delta_0(c)$ (أثر الحالة): الأثر السببي لكون الشخص معرضًا عند جرعة مرجعية محددة $c$ مقابل عدم التعرض ($A=0$).
   $$\Delta_0(c) = E\{Y(1, c) - Y(0, 0)\}$$

ب. درجة ميل الجزءين (Two-Part Propensity Score)
لضبط المتغيرات المربكة، قدم المؤلفون ناقل درجة ميل مكون من جزئين $S(X; \alpha) = (S_1, S_2)'$:

• $S_1$ (درجة ميل الجرعة): التوقع الشرطي للجرعة بشرط المعطيات بين المتعرضين: $S_1 = E(D \mid A=1, Z_1)$.
• $S_2$ (درجة ميل الحالة): الاحتمال الشرطي للتعرض بشرط المعطيات: $S_2 = P(A=1 \mid Z_2)$.
• ملاحظة: يمكن أن تختلف مجموعات المتغيرات المربكة $Z_1$ و $Z_2$، مما يسمح بآليات متميزة تقود حالة التعرض مقابل شدة التعرض.

ج. إجراء التقدير ثنائي المراحل
التقدير تسلسلي لفصل تحديد أثر الجرعة عن أثر الحالة:

• المرحلة الأولى: تقدير أثر الجرعة ($\Delta_D$)

  • المستهدف: الأفراد الذين لديهم $A=1$.
  • الطريقة: تعديل الانحدار باستخدام درجة الميل (PS Regression Adjustment).
  • النموذج: انحدار النتيجة $Y$ على الجرعة المرصودة $D$ ودرجة ميل الجرعة $S_1$.
  • المقدر: توفر معامل $D$ تقديرًا متسقًا لأثر الجرعة السببي ($\psi_{11}$)، بشرط أن يكون نموذج $S_1$ محددًا بشكل صحيح.

• المرحلة الثانية: تقدير أثر الحالة ($\Delta_0(c)$)

  • المستهدف: المقارنة بين $A=1$ (عند الجرعة المرجعية $c$) و $A=0$.
  • الاستراتيجية: يتم التعامل مع أثر الجرعة المقدر من المرحلة الأولى ($\hat{\gamma}_{11}$) كـ إزاحة (offset) ثابتة في نموذج المرحلة الثانية. هذا يزيل التباين في النتيجة الناتج عن الانحرافات عن الجرعة المرجعية، مما يعزل أثر "الحالة".
  • الطرق: اقترح المؤلفون ثلاثة نهج للمرحلة الثانية، تعتمد جميعها على درجة ميل الحالة ($S_2$):
    1. تعديل انحدار درجة الميل (PS Regression Adjustment): انحدار $Y$ على $A$ و $S_2$ (مع الإزاحة).
    2. الوزن العكسي للاحتمال (IPW): وزن الملاحظات بـ $1/S_2$ (أو $1/(1-S_2)$).
    3. الوزن العكسي للاحتمال المعزز (AIPW): يجمع بين IPW ونموذج تعويض النتيجة. هذا المقدر ثنائي المتانة (doubly robust): حيث يظل متسقًا إذا تم تحديد نموذج درجة ميل الحالة ($S_2$) بشكل صحيح أو إذا تم تحديد نموذج تعويض النتيجة بشكل صحيح.

د. الخصائص النظرية

• الاتساق والتوزيع الطبيعي التقاربي: أثبت المؤلفون أن المقدرات الناتلة متسقة وتتبع التوزيع الطبيعي التقاربي تحت الشروط المنتظمة القياسية وتحديد النموذج الصحيح.
• تحليل الخطأ في التحديد (Misspecification Analysis): اشتق المؤلفون القيم الحدية للمقدرات عند الخطأ في تحديد نماذج درجة الميل.
  • الخطأ في تحديد المرحلة الأولى ($S_1$) يؤدي إلى تحيز في تقدير الجرعة وينتقل إلى تحيز في تقدير حالة المرحلة الثانية.
  • الخطأ في تحديد المرحلة الثانية ($S_2$) يؤدي إلى تحيز في تقدير الحالة، ما لم يتم استخدام نهج AIPW ويكون نموذج تعويض النتيجة صحيحًا.

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

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

4. النتائج

الدراسات المحاكية

• الأداء: في ظل التحديد الصحيح للنماذج، أظهرت جميع الطرق المقترحة (انحدار PS، IPW، AIPW) تحيزًا ضئيلًا واحتمالات تغطية قريبة من المستوى الاسمي 95%.
• المتانة:
  • الخطأ في تحديد المرحلة الأولى ($S_1$): إذا تم تحديد درجة ميل الجرعة ($S_1$) بشكل خاطئ، فإن كلا التقديرين في المرحلة الأولى والثانية يصبحان متحيزين، مما يسلط الضوء على الأهمية البالغة للنمذجة الصحيحة لمكون التعرض المستمر.
  • الخطأ في تحديد المرحلة الثانية ($S_2$): أظهر مقدر AIPW متانة مزدوجة. عندما تم تحديد درجة ميل الحالة ($S_2$) بشكل خاطئ، ظل المقدر متسقًا بشرط أن يكون نموذج تعويض النتيجة صحيحًا. وعلى العكس من ذلك، إذا كان نموذج التعويض خاطئًا، فقد تم الحفاظ على الاتساق إذا كانت $S_2$ صحيحة.
• التباين: أنتجت الطرق القائمة على الأوزان (IPW/AIPW) عمومًا أخطاء معيارية أكبر مقارنة بتعديل الانحدار، وهو مقايضة معروفة مقابل المتانة.

التطبيق: التعرض للكحول أثناء الحمل (مجموعة ديترويت)

• البيانات: 377 طفلًا من مجموعة ديترويت الطولية، مع ناتج ذكاء كلي (WISC-III).
• النتائج:
  • التعديل التقليدي للمتغيرات المربكة: اقترح وجود أثر إيجابي صغير لحالة الشرب (من المحتمل أن يكون بسبب بقايا التضليل/الارتباك).
  • طرق درجة الميل (PS): أنتجت جميع طرق درجة الميل (الجزءين، درجتين، IPW، AIPW) تقديرات سالبة لكل من أثر الجرعة وأثر الحالة، مما يشير إلى تأثير ضار على الإدراك.
  • الدلالة: بينما أشارت التقديرات النقطية إلى الضرر، لم يصل أي منها إلى الدلالة الإحصائية عند مستوى 5% في هذه العينة المحددة. ومع ذلك، كانت تقديرات PS أكثر اتساقًا مع الأدلة الوبائية السابقة مقارنة بالتعديل الساذج للمتغيرات.
  • التعامل مع البيانات المفقودة: طبق المؤلفون بنجاح التضمين المتعدد (MICE) مع إطار العمل ثنائي المراحل، حيث تم تقدير درجات الميل داخل كل مجموعة بيانات مضمنة للحفاظ على الصلاحية.

5. الأهمية

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

• الإجابة على أسئلة علمية متميزة (مثل: "هل أي شرب ضار؟" مقابل "هل زيادة الشرب تسبب ضررًا أكبر؟").
• تجنب عيوب التقسيم العشوائي أو عدم استقرار الوزن العكسي للكثافة.
• الاستفادة من التقدير ثنائي المتانة (doubly robust) للحماية من الخطأ في تحديد النموذج في مكون التعرض الثنائي.

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