الاحتمال والإحصاء1958تأسيسي9 دقيقة قراءة

تحليل الانحدار للمتتاليات الثنائية

The Regression Analysis of Binary Sequences

Cox, D. R. — Journal of the Royal Statistical Society: Series B

المشكلة

الانحدار الخطي التقليدي يتنبأ بقيمة مستمرة، لكن كثيراً من المسائل العملية نتائجها ثنائية (0 أو 1). حين تمرر خطاً مستقيماً على بيانات ثنائية، قد يُخرج لك «احتمالات» أعلى من 1 أو أقل من 0، وهذا بلا معنى. كان الباحثون بحاجة إلى إطار إحصائي متين يصف كيف يتغيّر احتمال حدث ثنائي بحسب متغير تفسيري أو أكثر، مع اختبارات صارمة وطرق مضبوطة لتقدير المعلمات.

الإسهام

صاغ كوكس الانحدار اللوجستي صياغةً رسمية: نمذجة لوغاريتم أرجحية النتيجة الثنائية كدالة خطية في المتغيرات التفسيرية. تحوّل هذه القيمة الخطية إلى احتمال يقع بين 0 و1. واشتقّ طريقة للمعاملات، وحدّد الإحصاءات الكافية — وهي الملخصات التي تحتوي كامل معلومات البيانات عن كل معلمة — وبنى استدلالاً شرطياً دقيقاً يتخلّص من المعلمات المُزعجة دون أي تقريب. بذلك أصبح للبيانات الثنائية نظرية انحدار صارمة تُوازي الانحدار الخطي الكلاسيكي.

الأثر

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

تخيّل مفتاح إضاءة تدريجياً يتحكم في مصباح. الخطي يشبه قرصاً دوّاراً بلا حدود: يمكنك لفّه إلى ما بعد «الإضاءة الكاملة» أو ما قبل «الإطفاء الكامل»، فتدخل منطقة لا وجود لها — لأن المصباح لا يمكن أن يكون أسطع من 100% أو أخفت من 0%.

الانحدار اللوجستي يستبدل ذلك القرص بـآلية على شكل حرف S: مهما أدرتَ القرص، تتشبّع الإضاءة بسلاسة عند 0% و100%. موضع يدك (المُدخل) يُترجَم إلى درجة سطوع () لها دائماً معنى فيزيائي.

ورقة كوكس هي المخطط الهندسي لتلك الآلية.

المشكلة: الخطوط المستقيمة لا تصلح للتنبؤ بنعم أو لا

في الانحدار الخطي التقليدي تتنبأ بمخرج مستمر yy انطلاقاً من مُدخل xx: حيث y=α+βxy = \alpha + \beta x. لكن ماذا لو كانت yy لا تأخذ سوى 0 أو 1 — مريض يعيش أو يموت، قطعة سليمة أو معيبة؟

حين تمرر خطاً مستقيماً على أصفار وآحاد، سيتجاوز الخط حتماً القيمة 1 أو يهبط دون 0 عند القيم المتطرفة لـxx. هذه التنبؤات لا معنى لها كاحتمالات. ما تحتاجه هو نموذج مخرجاته دائماً احتمال صحيح — بين 0 و1.

افتح في المختبر
بدّل بين الخط المستقيم والمنحنى اللوجستي لترى لماذا يفشل الانحدار الخطي مع البيانات الثنائية.
تستيقظ التجربة عند وصولك…

الحل: الدالة اللوجستية

استخدم كوكس الدالة اللوجستية (المعروفة أيضاً بدالة سيجمويد) لتحويل أي عدد حقيقي إلى احتمال يقع بين 0 و1. الفكرة الجوهرية تمرّ بثلاث مراحل متتالية:

المرحلة الأولى — القيمة الخطية. نحسب مجموعاً مُرجّحاً z=α+βxz = \alpha + \beta x تماماً كما في الانحدار العادي. هذه القيمة يمكن أن تمتد من -\infty إلى ++\infty.

المرحلة الثانية — ضغط سيجمويد. ندخل zz في دالة سيجمويد ذات الشكل S:‏ σ(z)=1/(1+ez)\sigma(z) = 1 / (1 + e^{-z}). هذه الدالة «تضغط» أي عدد حقيقي ليقع داخل الفترة (0,1)(0, 1).

المرحلة الثالثة — التفسير كاحتمال. المخرج σ(z)\sigma(z) هو تقدير النموذج للقيمة P(Y=1x)P(Y = 1 \mid x) — أي احتمال أن تكون النتيجة 1 بمعلومية المُدخل xx.

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

P(Y=1x)=σ(α+βx)=11+e(α+βx)P(Y=1 \mid x) = \sigma(\alpha + \beta x) = \frac{1}{1 + e^{-(\alpha + \beta x)}}
نموذج الانحدار اللوجستي — الاحتمال كدالة سيجمويد لقيمة خطيةα = نقطة التقاطع (لوغاريتم الأرجحية الأساسي) · β = مقدار انزياح لوغاريتم الأرجحية عند زيادة x بوحدة واحدة · دالة سيجمويد σ تضغط القيمة الخطية داخل المجال [0, 1]
افتح في المختبر
اسحب α وβ لترى كيف ينزاح منحنى سيجمويد ويزداد حدّة.
تستيقظ التجربة عند وصولك…

اللوجيت: قلب الاحتمال رأساً على عقب

دالة سيجمويد تحوّل القيمة الخطية إلى احتمال. ومعكوسها — دالة اللوجيت — يحوّل الاحتمال عائداً إلى قيمة خطية. هذه النظرة المعكوسة هي ما يمنح الانحدار اللوجستي أناقته وقابليته للتفسير.

لنبدأ من أرجحية الحدث: إذا كان الاحتمال pp، فالأرجحية هي p/(1p)p / (1 - p). مثلاً احتمال 0.75 يعني أرجحية 3:1 — أي أن الحدث مرجّح ثلاث مرات أكثر من عدم حدوثه. بعد ذلك نأخذ اللوغاريتم الطبيعي للأرجحية فنحصل على لوغاريتم الأرجحية: logit(p)=ln(p/(1p))\text{logit}(p) = \ln(p / (1 - p)).

الملاحظة الجوهرية عند كوكس: لوجيت الاحتمال هو دالة خطية في المتغيرات التفسيرية. هذه أبسط علاقة يمكن تصوّرها بين المدخلات والمخرج المُحوَّل — وهذا تحديداً سبب تسمية النموذج بالانحدار اللوجستي (من log-odds).

عند هذه النقطة يصبح تفسير المعاملات بديهياً: β\beta هو مقدار التغيّر في لوغاريتم الأرجحية لكل وحدة تغيّر في xx. إذا كانت β\beta موجبة فزيادة xx تجعل الحدث أكثر ترجيحاً، وإذا كانت سالبة فزيادة xx تُضعف احتمال وقوعه.

logit(p)=ln ⁣(p1p)=α+βx\text{logit}(p) = \ln\!\left(\frac{p}{1-p}\right) = \alpha + \beta x
رابطة اللوجيت — لوغاريتم الأرجحية خطي في المتغيرات التفسيريةp/(1−p) = الأرجحية · ln يحوّلها إلى لوغاريتم الأرجحية · α + βx خط مستقيم — الصيغة ذاتها في الانحدار الخطي، لكن في فضاء لوغاريتم الأرجحية بدلاً من فضاء الاحتمال
افتح في المختبر
حرّك قيمة الاحتمال على اليسار وشاهد كيف تنتقل عبر الأرجحية → لوغاريتم الأرجحية على اليمين.
تستيقظ التجربة عند وصولك…

إيجاد أفضل منحنى: تقدير الأرجحية القصوى

في الانحدار العادي نقلّل مربعات الأخطاء. لكن مع النتائج الثنائية، الخطأ التربيعي ليس المقياس المناسب — فالتنبؤ بقيمة 0.99 لنتيجة فعلية تساوي 1 يجب أن يُكافأ أكثر من التنبؤ بقيمة 0.51.

استخدم كوكس تقدير الأرجحية القصوى: أوجد قيم α\alpha وβ\beta التي تجعل البيانات المُلاحظة أكثر ترجيحاً وفق النموذج. لكل نقطة بيانات: إذا كانت yi=1y_i = 1 فإسهام النموذج هو pip_i، وإذا كانت yi=0y_i = 0 فإسهامه هو 1pi1 - p_i. الأرجحية الكلية هي حاصل ضرب كل هذه الإسهامات، والهدف تعظيمها.

عملياً نعظّم لوغاريتم الأرجحية — مجموع بدلاً من حاصل ضرب — لأن الجمع أكثر استقراراً عددياً وأسهل في الاشتقاق. هذا اللوغاريتم هو بالضبط سالب ما يُسمّى في التعلم العميق الحديث الثنائية. كل شبكة عصبية دُرّبت على التصنيف تحلّ عملياً مسألة ذاتها التي صاغها كوكس عام 1958.

(α,β)=i=1n[yilnpi+(1yi)ln(1pi)]\ell(\alpha, \beta) = \sum_{i=1}^{n} \left[ y_i \ln p_i + (1 - y_i) \ln(1 - p_i) \right]
لوغاريتم الأرجحية — دالة الهدف في الانحدار اللوجستيكل حد يكافئ النموذج على منحه احتمالاً عالياً للنتيجة التي وقعت فعلاً · تعظيم هذا المجموع ≡ تقليل الإنتروبيا التقاطعية الثنائية
افتح في المختبر
اسحب المنحنى لمطابقة البيانات — شاهد لوغاريتم الأرجحية يرتفع كلما تحسّنت المطابقة.
تستيقظ التجربة عند وصولك…

الإحصاءات الكافية والاستدلال الشرطي

من أبرز إسهامات كوكس تحديده الإحصاءات الكافية للانحدار اللوجستي. هو ملخّص للبيانات يلتقط كامل المعلومات التي تحملها عن معينة — بمعنى أنك تستطيع الاستغناء عن البيانات الخام تماماً دون أن تفقد أي معلومة.

بيّن كوكس أن النموذج اللوجستي يملك الإحصاءات الكافية ذاتها الموجودة في النموذج الخطي الكلاسيكي: Y=yiY = \sum y_i (العدد الإجمالي لقيم 1) وX=yixiX = \sum y_i x_i (مجموع مُرجّح لقيم المتغير التفسيري). YY يختزل كل المعلومات عن نقطة التقاطع α\alpha، وXX يختزل كل المعلومات عن الميل β\beta.

هذا فتح الباب لتقنية فعّالة: . لتقدير β\beta دون الحاجة إلى التعامل مع المعلمة المُزعجة α\alpha، يكفي أن نشترط على القيمة الملاحظة لـYY. بذلك نتخلّص من α\alpha كلياً دون أي تقريب. هذا الأسلوب هو النظير الثنائي لاختبار فيشر الدقيق في جداول التوافق.

متغيرات تفسيرية متعددة: التوسّع إلى أبعاد أعلى

يمتد النموذج بشكل طبيعي ليشمل عدة متغيرات تفسيرية. مع pp من المتغيرات x1,,xpx_1, \ldots, x_p تصبح القيمة الخطية:

z=α+β1x1+β2x2++βpxpz = \alpha + \beta_1 x_1 + \beta_2 x_2 + \cdots + \beta_p x_p

كل βj\beta_j يقيس تأثير المتغير jj على لوغاريتم الأرجحية مع تثبيت بقية المتغيرات. المنطق هو ذاته في الانحدار الخطي المتعدد، لكنه يعمل في فضاء لوغاريتم الأرجحية بدلاً من فضاء القيم المستمرة.

بالترميز الحديث نكتب z=wx+bz = \mathbf{w}^\top \mathbf{x} + b — وهو مع ، أي اللبنة الأساسية في كل من طبقات الشبكة العصبية. طبقة المخرجات في أي مصنّف ثنائي ضمن التعلم العميق هي هذه المعادلة متبوعة بدالة سيجمويد — أي حرفياً نموذج الانحدار اللوجستي.

P(Y=1x)=11+e(wx+b)P(Y=1 \mid \mathbf{x}) = \frac{1}{1 + e^{-(\mathbf{w}^\top \mathbf{x} + b)}}
الانحدار اللوجستي متعدد المتغيرات بالترميز المتّجهي الحديثw = متّجه الأوزان · x = متّجه السمات · b = الانحياز · هذا هو عصبون المخرجات في كل مصنّف ثنائي
افتح في المختبر
حرّك متّجه الأوزان والانحياز لترى كيف تتحرّك حدود القرار في فضاء ثنائي الأبعاد.
تستيقظ التجربة عند وصولك…

الفكرة نفسها بالشيفرة

الانحدار اللوجستي من الصفرpython

مبسَّط لإظهار الفكرة — ليس التنفيذ الحقيقي.

import numpy as np

def sigmoid(z):
    """الدالة اللوجستية — تحوّل أي عدد حقيقي إلى (0, 1)."""
    return 1 / (1 + np.exp(-z))

def log_likelihood(X, y, w):
    """لوغاريتم الأرجحية: ما مدى جودة تفسير الأوزان للبيانات؟"""
    p = sigmoid(X @ w)
    return np.sum(y * np.log(p + 1e-12) + (1 - y) * np.log(1 - p + 1e-12))

def fit_logistic(X, y, lr=0.01, steps=1000):
    """الصعود التدريجي على لوغاريتم الأرجحية."""
    w = np.zeros(X.shape[1])
    for _ in range(steps):
        p = sigmoid(X @ w)
        gradient = X.T @ (y - p)      # اتجاه أشدّ صعود
        w += lr * gradient             # ادفع الأوزان لزيادة الأرجحية
    return w

# هذا كل شيء — سيجمويد + لوغاريتم الأرجحية + الصعود التدريجي هو
# كامل الخوارزمية. المكتبات الحديثة تُضيف تنظيماً ومُحسِّنات
# أذكى، لكن الجوهر لم يتغيّر منذ ورقة كوكس (1958).

السياق التاريخي: من التحليل الحيوي إلى التعلم العميق

الدالة اللوجستية نفسها تعود إلى فيرهولست (1845) الذي استخدمها لنمذجة النمو السكاني. ثم طبّقها بيركسون (1944) في التحليل الحيوي لتقدير الجرعة الدوائية التي يستجيب عندها 50% من الأفراد. لكن ورقة كوكس عام 1958 هي التي جمعت هذه الخيوط في إطار انحدار عام مع استدلال صارم، فحوّلت النموذج اللوجستي إلى أداة صالحة لـأي مسألة ذات نتائج ثنائية.

  1. 1845

    منحنى فيرهولست اللوجستي

    قدّم بيير-فرانسوا فيرهولست الدالة اللوجستية لنمذجة النمو السكاني بسقف استيعابي — أول ظهور للمنحنى على شكل حرف S.

  2. 1944

    نموذج بيركسون اللوجستي

    طبّق جوزيف بيركسون الدالة اللوجستية على منحنيات الاستجابة للجرعة في التحليل الحيوي، وصاغ مصطلح «لوجيت» ودافع عن تفوّقه على نموذج البروبيت.

  3. 1958

    كوكس: تحليل الانحدار للمتتاليات الثنائية

    صاغ الانحدار اللوجستي كإطار إحصائي عام مع تقدير الأرجحية القصوى والإحصاءات الكافية والاستدلال الشرطي.

  4. 1972

    النماذج الخطية المعمّمة

    وحّد نيلدر وودربورن الانحدار اللوجستي وانحدار بواسون وغيرهما في إطار واحد — النماذج الخطية المعمّمة — مع خوارزمية المربعات الصغرى المُعاد ترجيحها تكرارياً.

  5. 1986

    الانتشار العكسي وعصبون سيجمويد

    استخدم روملهارت وهينتون وويليامز دالة سيجمويد (اللوجستية) كدالة تنشيط في الشبكات العصبية المُدرَّبة بالانتشار العكسي — وارثين دالة كوكس مباشرةً.

  6. 2012

    ثورة التعلم العميق

    فازت شبكة AlexNet بمسابقة ImageNet مستخدمةً طبقات مخرجات سيجمويد/softmax مع الإنتروبيا التقاطعية — وهي التعميم متعدد الفئات للوغاريتم الأرجحية الثنائي الذي وضعه كوكس.

  7. 2026

    لا يزال النموذج المرجعي الافتراضي

    لا يزال الانحدار اللوجستي المصنّف المعتمد كنقطة انطلاق. وتُظهر الدراسات باستمرار أنه يُضاهي النماذج المعقدة أو ينافسها في كثير من السيناريوهات الواقعية.

لماذا لا تزال هذه الورقة مهمة

المرجعCox, D. R.. The Regression Analysis of Binary Sequences. Journal of the Royal Statistical Society: Series B, 1958.

مصطلحات هذه الورقة