Estimating conditional Mann-Whitney effects using pseudo-observation-based regression

فيما يلي ملخص تقني مفصل لورقة البحث بعنوان "تقدير تأثيرات مان-ويتني الشرطية باستخدام الانحدار القائم على الملاحظات الزائفة" (Estimating conditional Mann-Whitney effects using pseudo-observation-based regression) من إعداد دوبلر، شينك، وشميد.

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

تتناول الورقة تحدي تحليل تأثير مان-ويتني (المعروف أيضاً باسم المؤشر الاحتمالي أو التأثير النسبي)، والذي يُعرف بـ $\theta = P(T_1 > T_2)$، حيث $T_1$ و $T_2$ هما أوقات وقوع الأحداث من مجموعتين مستقلتين. وبينما يعد تأثير مان-ويتني غير الشرطي مقياساً معيارياً للمقارنات بين عينتين، فإن الطرق الحالية غالباً ما تفتقر إلى القدرة على:

1. التعديل وفقاً للمتغيرات المصاحبة (Covariates): دمج الخصائص الخاصة بالمريض ($Z$) لتقدير تأثير شرطي $\theta(z_1, z_2) = P(T_1 > T_2 | Z_1=z_1, Z_2=z_2)$.
2. التعامل مع البيانات المراقبة (Censored data): إدارة بيانات "وقت الحدث" التي تخضع للرقابة اليمينية العشوائية بشكل قوي دون الاعتماد على افتراضات معلمية مقيدة (مثل فرضية المخاطر النسبية المتناسبة).
3. القابلية للتفسير: تقديم النتائج على مقياس الاحتمال بدلاً من مقياس نسبة الخطر (log-hazard scale)، وهو ما يصعب تفسيره سريرياً في كثير من الأحيان.

يهدف المؤلفون إلى تطوير إطار عمل انحداري لا يعتمد على التوزيع (distribution-free) يربط تأثير مان-ويتني بمجموعة خطية من المتغيرات المصاحبة، وهو قابل للتطبيق على كل من البيانات المرصودة بالكامل والبيانات المراقبة يمينياً.

2. المنهجية

يجمع النهج المقترح بين الملاحظات الزائفة (pseudo-observations) ومعادلات التقدير المعممة (GEE).

أ. النموذج

يحدد المؤلفون نموذج انحدار لتأثير مان-وتني الشرطي:
$$\theta(z_1, z_2) = \mu(\beta_0 + \beta_1^\top z_1 + \beta_2^\top z_2)$$
حيث أن:

• $\mu(\cdot)$ هي دالة ربط عكسية رتيبة تماماً (مثل الدالة المحايدة أو اللوجستية).
• $\beta$ هو متجه المعلمات.
• $z_1$ و $z_2$ هما متجهات المتغيرات المصاحبة للأفراد في المجموعة 1 والمجموعة 2 على التوالي.
• يسمح النموذج بمجموعات مختلفة من المتغيرات المصاحبة في كل مجموعة، وإن كان الحالة الخاصة $z_1 = z_2 = z$ تسمح بالتفسير السببي تحت الافتراضات القياسية.

ب. بناء الملاحظات الزائفة

لتطبيق النموذج، قام المؤلفون بتكييف تقنية جاكنيف لعينتين (two-sample jackknife) لتوليد الملاحظات الزائفة ($\tilde{\theta}_{i_1 i_2}$).

• حالة البيانات المرصودة بالكامل: يتم اشتقاق الملاحظة الزائفة للزوج $(i_1, i_2)$ عن طريق إزالة الملاحظة رقم $i_1$ من المجموعة 1 والملاحظة رقم $i_2$ من المجموعة 2، ثم إعادة حساب المقدر، ودمج هذه التقديرات الناتجة عن استبعاد عنصر واحد. في حالة عدم وجود رقابة، يتبسط هذا إلى دالة المؤشر $1\{T_{1i_1} > T_{2i_2}\}$.
• حالة البيانات المراقبة (Censored): يتم استبدال المقدر القياسي $\hat{\theta}$ بتكامل يتضمن مقدرات كابلان-ماير للبقاء ($\hat{S}_1, \hat{S}_2$). وتُبنى الملاحظات الزائفة باستخدام نفس منطق "جاكنيف" ولكن بناءً على مقدرات كابلان-ماير المستمدة من البيانات المراقبة.
• تصحيح التساوي (Tie Correction): يتم تطبيق تصحيح الاستمرارية للتعامل مع أوقات الأحداث المتساوية، وذلك بإضافة $0.5 \cdot P(T_1 = T_2)$ إلى تقدير الاحتمال.

ج. تقدير المعلمات

يتم تقدير معلمات الانحدار $\beta$ عن طريق حل معادلة التقدير المعممة (GEE):
$$U(\beta) = \sum_{i_1} \sum_{i_2} A_{i_1 i_2} (\tilde{\theta}_{i_1 i_2} - \mu_{i_1 i_2}) = 0$$
حيث $A_{i_1 i_2}$ هو مشتق دالة الربط بالنسبة لـ $\beta$، و $\mu_{i_1 i_2}$ هي القيمة المتوقعة للنموذج.

د. النظرية التقاربية والاستدلال

• الاتساق والاعتدالية: أثبت المؤلفون وجود نظرية حد مركزي (CLT) تُظهر أن المقدر $\hat{\beta}$ متسق ويتبع التوزيع الطبيعي تقاربياً تحت شروط انتظام طفيفة (الافتراض 1)، حتى في ظل وجود الرقابة اليمينية.
• اختبار الفرضيات: نظراً لأن هيكل التباين التقاربي معقد (يتضمن المشتقات الثانية للدالة)، اقترح المؤلفون اختبارات فرضية قائمة على البوتستراب (bootstrap) (تحديداً باستخدام "بوتستراب السرعة الفائقة" للكفاءة) لاستنتاج دلالة تأثيرات المتغيرات المصاحبة. تم تقييم أربعة إحصاءات اختبار:
  1. الانحراف المعياري التجريبي.
  2. المدى الربيعي (IQR).
  3. الانحراف المطلق عن الوسيط (MAD).
  4. طريقة المئويات (Percentile method).

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

1. إطار عمل انحداري مبتكر: هذه هي الورقة الأولى التي تقترح نموذج انحدار لتأثير مان-وتني يمكنه التعامل مع البيانات المراقبة يمينياً والبيانات ذات الأوقات المتساوية باستخدام نهج لا يعتمد على التوزيع.
2. توسيع الملاحظات الزائفة: قامت الورقة بتوسيع منهجية الملاحظات الزائفة لـ "جاكنيف" من مشكلات العينة الواحدة إلى إطار العينتين، مما وفر أساساً نظرياً صارماً للخصائص التقاربية للمقدر.
3. التعامل المرن مع المتغيرات المصاحبة: يسمح النموذج بمجموعات مختلفة من المتغيرات المصاحبة في المجموعتين، مما يجعله مناسباً لمقارنة البيانات من مصادر مختلفة (مثل الضوابط التاريخية مقابل التجارب الجديدة).
4. القابلية للتفسير: بخلاف نماذج كوكس التي تقدر نسب الخطر اللوغاريتمية، يقدر هذا الأسلوب التأثيرات على مقياس الاحتمال ($P(T_1 > T_2)$)، مما يوفر تفسيراً سريرياً مباشراً.
5. التفسير السببي: في الحالة الخاصة التي تكون فيها المتغيرات المصاحبة متطابقة ($z_1=z_2$)، يحدد النموذج ضمنياً مقدراً سببياً لاحتمالية الاستفادة من العلاج.

4. النتائج

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

• الإعداد: أُجريت محاكاة واسعة النطاق باستخدام أوقات أحداث تتبع توزيع "ويبول" (Weibull)، مع تغيير أحجام العينات، ونسب الرقابة (تصل إلى 75%)، وأبعاد المتغيرات المصاحبة.
• خطأ النوع الأول (Type-I Error): أظهرت اختبارات البوتستراب القائمة على الانحراف المعياري التجريبي والمئويات التجريبية قدرة جيدة على ضبط معدل خطأ النوع الأول (بالقرب من 5% الاسمية)، حتى في ظل الرقابة العالية. بينما وُجد أن الاختبارات القائمة على IQR و MAD كانت "متساهلة" قليلاً (تؤدي لرفض الفرضية الصفرية أكثر من اللازم).
• القوة (Power): أظهرت الاختبارات المقترحة قوة مرضية، مشابهة لاختبار z من نموذج كوكس، بشرط ألا تتجاوز الرقابة 50%.
• المقارنة مع كوكس: كان اختبار z في نموذج كوكس غالباً ما يكون متساهلاً للغاية في العينات الصغيرة. قدمت الطريقة المقترحة بديلاً قوياً دون افتراض التناسب في المخاطر.

تطبيق على بيانات حقيقية (تجربة SUCCESS-A)

طُبق الأسلوب على تجربة SUCCESS-A لسرطان الثدي (3,754 مريضاً) لتحليل البقاء على قيد الحياة بدون مرض (DFS).

• النتائج: بينما وجدت التحليلات الأصلية عدم وجود فائدة عامة من العلاج، حدد الانحدار المقترح تبايناً في تأثيرات العلاج:
  • الفائدة من التدخل (Gemcitabine): المرضى ذوو مراحل الأورام المتقدمة (pT2–pT4)، والتمايز الضعيف (G3)، والأنواع الجزيئية "ثلاثية السلبية" أظهروا احتمالية أعلى للاستفادة من التدخل.
  • الفائدة من مجموعة الضبط: المرضى ذوو الأورام في مراحل مبكرة، والحالة الهرمونية الإيجابية (ER+/PR+)، وعدم تأثر العقد اللمفاوية، أظهروا احتمالية أعلى للاستفادة من العلاج الكيميائي القياسي وحده.
• المقارنة مع كوكس: حدد النموذج المقترح اتجاهات تأثير مختلفة لبعض المتغيرات (مثل مرحلة الورم) مقارنة بنموذج كوكس. ويعزو المؤلفون ذلك إلى اختلاف المقدرات: يستخدم النموذج المقترح تأثيراً مبتوراً (truncated) عبر أفق زمني محدد ($\tau = 5.5$ سنة)، بينما يقدر نموذج كوكس تأثيراً غير مبتور بافتراض تناسب المخاطر.
• الفائدة السريرية: نجح النموذج في تصنيف المرضى الأفراد إلى مجموعات "مستفيدة من التدخل" و"مستفيدة من الضبط" مع فترات ثقة 95%، مما يبرز إمكانية اتخاذ قرارات علاجية شخصية.

5. الأهمية والخاتمة

توفر هذه الورقة بديلاً قوياً لا يعتمد على التوزيع لنموذج "كوكس للمخاطر النسبية" لتحليل بيانات وقت الحدث. وتكمن أهميتها في:

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

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