أساسيات تعلم الآلة1970تأسيسي9 دقيقة قراءة

انحدار الحَرْف: التقدير المنحاز للمسائل غير المتعامدة

Ridge Regression: Biased Estimation for Nonorthogonal Problems

Hoerl, A. E. · Kennard, R. W. — Technometrics

المشكلة

العادية (OLS) هي الأداة الأساسية في تحليل الانحدار، لكن المشكلة تظهر حين تكون المتغيّرات التنبؤية مترابطة فيما بينها — وهو ما يُعرف بتعدّد الخطية المشترك. في هذه الحالة، أيّ تغيير طفيف في البيانات يدفع المعاملات المُقدَّرة إلى قفزات عنيفة لا يمكن الوثوق بها. السبب الرياضي واضح: المصفوفة X'X تقترب من حالة الانفراد، فقيمها الذاتية الصغيرة تُضخّم β̂ إلى حدود غير عملية. والنتيجة قد يُلائم بيانات التدريب بإتقان، لكنه يفشل عند أول مواجهة مع بيانات جديدة.

الإسهام

الفكرة بسيطة: أضِف ثابتاً موجباً صغيراً k إلى القطر الرئيسي للمصفوفة X'X قبل أن تعكسها: β̂_ridge = (X'X + kI)⁻¹ X'y. هذا «الحَرْف» المُضاف إلى القطر يمنع المصفوفة من الانفجار عند العكس، ويسحب كل معامل باتجاه الصفر، فيُدخل قدراً مضبوطاً من الانحياز. لكنّ المبرهنة المركزية في الورقة تُثبت أنّ هذه المقايضة تستحق دائماً: يوجد حتماً k > 0 يكون عنده الخطأ الكلي (MSE) لـ β̂_ridge أقل من خطأ OLS، لأنّ ما تكسبه من استقرار يفوق ما تخسره من انحياز.

الأثر

غيّر انحدار الحَرْف قناعة راسخة في الإحصاء: أثبت أنّ مُقدِّراً منحازاً قد يتفوّق على المُقدِّر غير المنحاز، وبذلك فتح الباب أمام مجال التنظيم بأكمله. من هذه الفكرة وُلد لاسو (عقوبة L1) والشبكة المرنة وكل نموذج منظَّم نستخدمه اليوم. أمّا مفاضلة الانحياز والتباين التي صاغتها الورقة فأصبحت ركيزة أساسية تُدرَّس في كل مقرّر تعلّم آلي.

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

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

انحدار الحَرْف يُعالج ذلك بسحب كل حَكَم بلطف نحو الدرجات المعتدلة. تأثير الثنائي المترابط يتراجع، والدرجة النهائية قد تنحرف قليلاً عن المركز — لكنها أكثر استقراراً بمراحل. والاستقرار هو ما يهمّك حين تريد التنبؤ بالغطسة القادمة.

المشكلة: ماذا يحدث حين تتحرّك المتنبِّئات معاً؟

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

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

افتح في المختبر
حرّك شريط الارتباط وراقب ما يحدث: كلّما زاد الارتباط بين المتنبّئات، تضطربت معاملات OLS بينما تبقى معاملات الحَرْف مستقرّة.
تستيقظ التجربة عند وصولك…

القيم الذاتية: مُضخِّمات الضوضاء

لكي نفهم لماذا يحدث هذا الانفجار، ننظر إلى X'X من زاوية القيم الذاتية. كل قيمة ذاتية λj\lambda_j تُعبّر عن مقدار «الإشارة» المتوفّرة في اتجاه معيّن من فضاء المتنبّئات. حين تعكس المربعات الصغرى المصفوفة X'X، فهي فعلياً تقسم على كل قيمة ذاتية. الاتجاهات القوية — ذات λj\lambda_j الكبيرة — لا مشكلة فيها. أمّا الاتجاهات ذات λj\lambda_j الضئيلة فتتضخّم تضخّماً هائلاً: Var(β^j)1/λj\text{Var}(\hat{\beta}_j) \propto 1/\lambda_j.

هناك مقياس عملي لرصد هذه المشكلة اسمه مُعامل تضخّم التباين (VIF)، ويُحسب كالتالي: VIFj=1/(1Rj2)\text{VIF}_j = 1/(1 - R_j^2)، حيث Rj2R_j^2 يقيس مدى إمكانية التنبؤ بالمتغيّر jj من بقية المتغيّرات. إذا تجاوز VIF القيمة 10، فهذا تحذير واضح: تباين ذلك المعامل مُضخَّم عشر مرات على الأقل.

افتح في المختبر
راقب ما يحدث حين تنكمش القيم الذاتية لـ X'X: تباين OLS ينفجر. إضافة الحَرْف (kI) ترفع أرضية كل قيمة ذاتية فتمنع هذا الانفجار.
تستيقظ التجربة عند وصولك…

الحل: حَرْف صغير على القطر يُغيّر كل شيء

الحل أنيق في بساطته: قبل أن تعكس المصفوفة X'X، أضِف ثابتاً موجباً صغيراً kk إلى كل عنصر على القطر الرئيسي. الأثر العملي أشبه برفع أرضية تحت كل القيم الذاتية — حتى الضئيلة منها أصبحت لا تقلّ عن kk، فلم يعد عكسها يُنتج أرقاماً منفجرة.

من الزاوية الهندسية، المربعات الصغرى تبحث عن β الذي يُقلّل yXβ2\|y - X\beta\|^2. انحدار الحَرْف يُضيف شرطاً إضافياً: قلّل yXβ2+kβ2\|y - X\beta\|^2 + k\|\beta\|^2. الحدّ الثاني kβ2k\|\beta\|^2 يفرض على المعاملات ألّا تكبر أكثر من اللازم — يسحبها نحو الصفر بمقدار يتحكّم فيه kk. هذا هو أقدم شكل لما نسمّيه اليوم .

β^ridge=(XTX+kI)1XTy\hat{\beta}_{\text{ridge}} = (X^T X + kI)^{-1} X^T y
مُقدِّر الحَرْف — المعادلة المركزية في الورقةX'X = مصفوفة المعادلات السويّة · kI = الحَرْف المُضاف إلى القطر · النتيجة: معاملات منحازة قليلاً نحو الصفر لكنها أكثر استقراراً بكثير من المربعات الصغرى

تخيّل الأمر كلعبة شدّ حبل بين قوّتين: الأولى تسحب β نحو حلّ المربعات الصغرى لكي يُطابق البيانات، والثانية تسحبه نحو الصفر لكي يبقى صغيراً ومنضبطاً. المُعامل kk يُحدّد أيّ الحبلين أقوى. حين k=0k = 0 تسيطر المربعات الصغرى وحدها. ومع زيادة kk، تنكمش المعاملات أكثر فيزداد وينخفض التباين — وبين الطرفين توجد نقطة مثلى يكون عندها الخطأ الكلي (MSE) في أدنى مستوياته.

افتح في المختبر
حرّك k بدءاً من الصفر وراقب المعاملات وهي تنكمش تدريجياً. الانحياز يرتفع والتباين ينخفض — ومنحنى MSE يكشف لك أين تقع النقطة المثلى.
تستيقظ التجربة عند وصولك…

المبرهنة: لماذا قد يكون المُقدِّر المنحاز أفضل؟

مبرهنة غاوس-ماركوف تقول إنّ المربعات الصغرى هي أفضل مُقدِّر خطي غير منحاز (BLUE). لكن تأمّل العبارة جيداً: «أفضل بين غير المنحازين» لا يعني «أفضل على الإطلاق». الخطأ الكلي للتنبؤ يتحلّل إلى مركّبتين: MSE=Bias2+Variance\text{MSE} = \text{Bias}^2 + \text{Variance}. المربعات الصغرى تُحقّق انحيازاً صفرياً، لكنّ تباينها قد يكون ضخماً حين تصغر القيم الذاتية. انحدار الحَرْف يقبل بقدر من الانحياز — يسحب المعاملات منهجياً نحو الصفر — لكنّ المكسب في خفض التباين أكبر بكثير، فينخفض الخطأ الكلي.

النتيجة المحورية لهورل وكينارد تقول: مهما كانت قيمة β الحقيقية، يوجد حتماً k>0k > 0 يجعل MSE(β^ridge)<MSE(β^OLS)\text{MSE}(\hat{\beta}_{\text{ridge}}) < \text{MSE}(\hat{\beta}_{\text{OLS}}). بعبارة أخرى، قبول قدر من الانحياز يُعطيك نتائج أفضل دائماً حين تكون المتنبّئات مترابطة.

MSE(β^ridge)=σ2j=1pλj(λj+k)2+k2j=1pαj2(λj+k)2\text{MSE}(\hat{\beta}_{\text{ridge}}) = \sigma^2 \sum_{j=1}^{p} \frac{\lambda_j}{(\lambda_j + k)^2} + k^2 \sum_{j=1}^{p} \frac{\alpha_j^2}{(\lambda_j + k)^2}
تفكيك MSE — التباين (المجموع الأول) + الانحياز² (المجموع الثاني)λⱼ = القيم الذاتية لـ X'X · αⱼ = المعاملات الحقيقية في أساس المتجهات الذاتية · k يُقلّص حدّ التباين (من خلال (λⱼ+k)²) لكنه يُنمّي حدّ الانحياز (من خلال k²) — منحنى MSE على شكل U وله نقطة دنيا عند k الأمثل
افتح في المختبر
راقب منحنيَي الانحياز² والتباين أثناء تغيير k. مجموعهما — أي MSE — يرسم شكل حرف U، وقاعه هو قيمة الحَرْف المثلى.
تستيقظ التجربة عند وصولك…

أثر الحَرْف: كيف تختار k بعينيك

اقترح هورل وكينارد أداة بصرية تُسمّى أثر الحَرْف: ترسم كل معامل β^j\hat{\beta}_j كدالة في kk. عند k=0k = 0 — أي المربعات الصغرى — قد تجد المعاملات كبيرة جداً أو متقلّبة أو حتى بإشارة معكوسة. لكن مع زيادة kk تبدأ المعاملات بالهدوء والاستقرار، وغالباً تستقرّ عند قيم تتوافق مع المنطق الفيزيائي للمسألة. القاعدة العملية: اختر أصغر kk تستقرّ عنده المعاملات.

في الممارسة الحديثة، حلّ محلّ أثر الحَرْف في كثير من الحالات — قسّم البيانات إلى أجزاء، جرّب قيماً مختلفة لـ kk، واختر القيمة التي تُعطي أقل خطأ تنبؤ على البيانات المحجوبة. لكنّ أثر الحَرْف يظل أداة ثمينة لفهم أيّ المعاملات هشّة ولماذا.

افتح في المختبر
راقب كيف تهدأ المعاملات تدريجياً مع زيادة k. الخط المتقطّع العمودي يُشير إلى k الأمثل الذي حدّده التحقق المتقاطع.
تستيقظ التجربة عند وصولك…

المنظور الهندسي: الدائرة تلتقي بالقطع الناقص

يمكن فهم انحدار الحَرْف بصورة هندسية أنيقة. المربعات الصغرى تُقلّل yXβ2\|y - X\beta\|^2، وإذا رسمت خطوط تساوي هذه الخسارة في فضاء المعاملات β، تحصل على قطوع ناقصة مركزها حلّ المربعات الصغرى. عقوبة الحَرْف kβ2k\|\beta\|^2 تُضيف شرطاً: «لا تبتعد كثيراً عن نقطة الأصل» — أي ابقَ داخل كرة نصف قطرها tt. حلّ الحَرْف هو النقطة التي يلمس فيها أصغر قطع ناقص للخسارة سطحَ كرة القيد.

قارن ذلك مع لاسو: قيد L1 شكله كالمُعيَّن — زواياه على المحاور. ولأنّ نقطة التماس تقع غالباً عند إحدى هذه الزوايا، يُصفّر لاسو بعض المعاملات تماماً — وهذا ما يجعله أداة لانتقاء المتغيّرات. أمّا قيد L2 فشكله دائري بلا زوايا، ولذلك يُقلّص انحدار الحَرْف المعاملات بسلاسة دون أن يُصفّر أيّاً منها.

افتح في المختبر
القطع الناقص يُمثّل خطوط تساوي خسارة OLS، الدائرة هي قيد الحَرْف، والمعيَّن هو قيد لاسو. حرّك k وراقب أين يلمس القطع الناقص كل شكل.
تستيقظ التجربة عند وصولك…

الفكرة ذاتها في شيفرة برمجية

انحدار الحَرْف من الصفرpython

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

import numpy as np

def ridge(X, y, k):
    """انحدار الحَرْف: (X'X + kI)^{-1} X'y"""
    p = X.shape[1]
    return np.linalg.solve(X.T @ X + k * np.eye(p), X.T @ y)

def ols(X, y):
    """المربعات الصغرى العادية — انحدار حَرْف بـ k=0"""
    return ridge(X, y, k=0)

def ridge_mse(X, y, beta_true, k):
    """احسب MSE = الانحياز² + التباين لقيمة k معطاة"""
    p = X.shape[1]
    W = np.linalg.inv(X.T @ X + k * np.eye(p)) @ X.T
    beta_hat = W @ y
    bias_sq = np.sum((beta_hat - beta_true) ** 2)
    # أثر التباين: σ² * trace(W @ W')
    return bias_sq  # مُبسَّط — النسخة الكاملة تُضيف حدّ التباين

# جرّبها: متنبّئات مترابطة
rng = np.random.default_rng(42)
n, p = 100, 5
X = rng.standard_normal((n, p))
X[:, 1] = X[:, 0] + 0.01 * rng.standard_normal(n)  # شبه خطية
beta_true = np.array([1, 2, 0.5, -1, 0.3])
y = X @ beta_true + 0.5 * rng.standard_normal(n)

print("OLS:  ", ols(X, y).round(2))        # غير مستقر
print("Ridge:", ridge(X, y, k=1.0).round(2))  # مستقر

الصورة الأوسع: من بايز إلى التعلم العميق

انحدار الحَرْف ليس مجرد حيلة جبرية — إنه يتقاطع مع أفكار كبيرة في عدة مجالات. من المنظور البايزي، هو مكافئ تماماً لوضع غاوسية βN(0,σ2/kI)\beta \sim N(0, \sigma^2/k \cdot I) على المعاملات ثم حساب منوال . المُعامل kk يُعبّر ببساطة عن اعتقاد مسبق بأنّ المعاملات لن تكون كبيرة جداً.

من منظور معالجة الإشارات، انحدار الحَرْف هو مُرشّح وينر (Wiener Filter) الذي يوازن بين كبت الضوضاء والحفاظ على الإشارة الأصلية.

أمّا في ، فإنّ اضمحلال الأوزان — أي إضافة λθ2\lambda \|\theta\|^2 إلى الخسارة — هو بالضبط انحدار حَرْف مُطبَّق على معاملات . حين تكتب weight_decay=0.01 في ، فأنت تستخدم الفكرة ذاتها التي نشرها هورل وكينارد عام 1970.

الأثر الأوسع: كيف وُلد مجال التنظيم

أهمية انحدار الحَرْف تتجاوز المسألة التي حلّها. المبدأ الذي أرساه — أنّ تقييد حرية النموذج يُحسّن قدرته على التنبؤ — أعاد تشكيل الإحصاء و بأكمله، وصار الأساس الذي بُنيت عليه كل أدوات التنظيم الحديثة.

  1. 1970

    انحدار الحَرْف

    هورل وكينارد يُثبتان أنّ إضافة kI إلى X'X تُعطي دائماً خطأ تنبؤ أقل من المربعات الصغرى حين تكون المتنبّئات مترابطة. هنا وُلد التنظيم L2.

  2. 1996

    لاسو (تيبشيراني)

    الفكرة: استبدل عقوبة L2 بعقوبة L1. النتيجة: بعض المعاملات تصل إلى الصفر تماماً — وبذلك أصبح التنظيم أداةً لانتقاء المتغيّرات تلقائياً.

  3. 2005

    الشبكة المرنة

    زو وهاستي يمزجان عقوبتَي L1 وL2 في صيغة واحدة، فيجمعان بين انتقاء المتغيّرات والتعامل الجيّد مع المتغيّرات المترابطة.

  4. 2012

    التعلم العميق + اضمحلال الأوزان

    مع AlexNet وما تلاها، أصبح اضمحلال الأوزان (تنظيم L2 على كل المعاملات) ممارسة قياسية في الشبكات العميقة. فكرة الحَرْف انتقلت من بضعة معاملات إلى مليارات.

  5. 2017

    AdamW — اضمحلال أوزان مفكوك

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

الحَرْف أم لاسو: كيف تختار بينهما؟

لاسو (عقوبة L1) وانحدار الحَرْف (عقوبة L2) كلاهما يُقلّص المعاملات، لكنّ آلية التقليص مختلفة تماماً:

  • انحدار الحَرْف يُقلّص كل المعاملات بنِسَب متقاربة دون أن يُلغي أيّاً منها. الخيار المناسب حين تعتقد أنّ جميع المتنبّئات تحمل معلومات مفيدة وأنّ المشكلة هي الارتباط بينها — كما في التحليل الطيفي أو بيانات الجينوم ذات المسارات المتداخلة.

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

  • الشبكة المرنة تمزج بين الاثنين: αβ1+(1α)β22\alpha \|\beta\|_1 + (1-\alpha)\|\beta\|_2^2. مفيدة حين تريد انتقاء المتغيّرات مع الحفاظ على التعامل الجيّد مع المتغيّرات المترابطة.

باختصار: الحَرْف يقول «أبقِ الجميع لكن اخفض أصواتهم»، ولاسو يقول «اطرد من لا فائدة منه». عملياً، جرّب الاثنين ودع التحقق المتقاطع يُقرّر.

افتح في المختبر
قارن مسارات المعاملات في الحَرْف (L2) ولاسو (L1) جنباً إلى جنب. لاحظ أنّ لاسو يُوصل بعض المعاملات إلى الصفر تماماً، بينما الحَرْف يُقرّبها من الصفر دون أن يُلغيها.
تستيقظ التجربة عند وصولك…

المرجعHoerl, Kennard. Ridge Regression: Biased Estimation for Nonorthogonal Problems. Technometrics, 1970.

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