أساسيات تعلم الآلة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 سيئة الاشتراط: بعض تنكمش قرب الصفر. الآن تخيّل أنك تقسم على عدد قريب جداً من الصفر — النتيجة تنفجر. هذا بالضبط ما يحدث: المربعات الصغرى تُضخّم في الاتجاهات الهشّة، فتُخرج معاملات قد تكون غير منحازة نظرياً لكنها تقفز بعنف من عيّنة إلى أخرى — أي أنها عديمة الفائدة عملياً.
القيم الذاتية: مُضخِّمات الضوضاء
لكي نفهم لماذا يحدث هذا الانفجار، ننظر إلى X'X من زاوية القيم الذاتية. كل قيمة ذاتية تُعبّر عن مقدار «الإشارة» المتوفّرة في اتجاه معيّن من فضاء المتنبّئات. حين تعكس المربعات الصغرى المصفوفة X'X، فهي فعلياً تقسم على كل قيمة ذاتية. الاتجاهات القوية — ذات الكبيرة — لا مشكلة فيها. أمّا الاتجاهات ذات الضئيلة فتتضخّم تضخّماً هائلاً: .
هناك مقياس عملي لرصد هذه المشكلة اسمه مُعامل تضخّم التباين (VIF)، ويُحسب كالتالي: ، حيث يقيس مدى إمكانية التنبؤ بالمتغيّر من بقية المتغيّرات. إذا تجاوز VIF القيمة 10، فهذا تحذير واضح: تباين ذلك المعامل مُضخَّم عشر مرات على الأقل.
الحل: حَرْف صغير على القطر يُغيّر كل شيء
الحل أنيق في بساطته: قبل أن تعكس المصفوفة X'X، أضِف ثابتاً موجباً صغيراً إلى كل عنصر على القطر الرئيسي. الأثر العملي أشبه برفع أرضية تحت كل القيم الذاتية — حتى الضئيلة منها أصبحت لا تقلّ عن ، فلم يعد عكسها يُنتج أرقاماً منفجرة.
من الزاوية الهندسية، المربعات الصغرى تبحث عن β الذي يُقلّل . انحدار الحَرْف يُضيف شرطاً إضافياً: قلّل . الحدّ الثاني يفرض على المعاملات ألّا تكبر أكثر من اللازم — يسحبها نحو الصفر بمقدار يتحكّم فيه . هذا هو أقدم شكل لما نسمّيه اليوم .
تخيّل الأمر كلعبة شدّ حبل بين قوّتين: الأولى تسحب β نحو حلّ المربعات الصغرى لكي يُطابق البيانات، والثانية تسحبه نحو الصفر لكي يبقى صغيراً ومنضبطاً. المُعامل يُحدّد أيّ الحبلين أقوى. حين تسيطر المربعات الصغرى وحدها. ومع زيادة ، تنكمش المعاملات أكثر فيزداد وينخفض التباين — وبين الطرفين توجد نقطة مثلى يكون عندها الخطأ الكلي (MSE) في أدنى مستوياته.
المبرهنة: لماذا قد يكون المُقدِّر المنحاز أفضل؟
مبرهنة غاوس-ماركوف تقول إنّ المربعات الصغرى هي أفضل مُقدِّر خطي غير منحاز (BLUE). لكن تأمّل العبارة جيداً: «أفضل بين غير المنحازين» لا يعني «أفضل على الإطلاق». الخطأ الكلي للتنبؤ يتحلّل إلى مركّبتين: . المربعات الصغرى تُحقّق انحيازاً صفرياً، لكنّ تباينها قد يكون ضخماً حين تصغر القيم الذاتية. انحدار الحَرْف يقبل بقدر من الانحياز — يسحب المعاملات منهجياً نحو الصفر — لكنّ المكسب في خفض التباين أكبر بكثير، فينخفض الخطأ الكلي.
النتيجة المحورية لهورل وكينارد تقول: مهما كانت قيمة β الحقيقية، يوجد حتماً يجعل . بعبارة أخرى، قبول قدر من الانحياز يُعطيك نتائج أفضل دائماً حين تكون المتنبّئات مترابطة.
أثر الحَرْف: كيف تختار k بعينيك
اقترح هورل وكينارد أداة بصرية تُسمّى أثر الحَرْف: ترسم كل معامل كدالة في . عند — أي المربعات الصغرى — قد تجد المعاملات كبيرة جداً أو متقلّبة أو حتى بإشارة معكوسة. لكن مع زيادة تبدأ المعاملات بالهدوء والاستقرار، وغالباً تستقرّ عند قيم تتوافق مع المنطق الفيزيائي للمسألة. القاعدة العملية: اختر أصغر تستقرّ عنده المعاملات.
في الممارسة الحديثة، حلّ محلّ أثر الحَرْف في كثير من الحالات — قسّم البيانات إلى أجزاء، جرّب قيماً مختلفة لـ ، واختر القيمة التي تُعطي أقل خطأ تنبؤ على البيانات المحجوبة. لكنّ أثر الحَرْف يظل أداة ثمينة لفهم أيّ المعاملات هشّة ولماذا.
المنظور الهندسي: الدائرة تلتقي بالقطع الناقص
يمكن فهم انحدار الحَرْف بصورة هندسية أنيقة. المربعات الصغرى تُقلّل ، وإذا رسمت خطوط تساوي هذه الخسارة في فضاء المعاملات β، تحصل على قطوع ناقصة مركزها حلّ المربعات الصغرى. عقوبة الحَرْف تُضيف شرطاً: «لا تبتعد كثيراً عن نقطة الأصل» — أي ابقَ داخل كرة نصف قطرها . حلّ الحَرْف هو النقطة التي يلمس فيها أصغر قطع ناقص للخسارة سطحَ كرة القيد.
قارن ذلك مع لاسو: قيد L1 شكله كالمُعيَّن — زواياه على المحاور. ولأنّ نقطة التماس تقع غالباً عند إحدى هذه الزوايا، يُصفّر لاسو بعض المعاملات تماماً — وهذا ما يجعله أداة لانتقاء المتغيّرات. أمّا قيد L2 فشكله دائري بلا زوايا، ولذلك يُقلّص انحدار الحَرْف المعاملات بسلاسة دون أن يُصفّر أيّاً منها.
الفكرة ذاتها في شيفرة برمجية
مبسَّط لإظهار الفكرة — ليس التنفيذ الحقيقي.
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)) # مستقرالصورة الأوسع: من بايز إلى التعلم العميق
انحدار الحَرْف ليس مجرد حيلة جبرية — إنه يتقاطع مع أفكار كبيرة في عدة مجالات. من المنظور البايزي، هو مكافئ تماماً لوضع غاوسية على المعاملات ثم حساب منوال . المُعامل يُعبّر ببساطة عن اعتقاد مسبق بأنّ المعاملات لن تكون كبيرة جداً.
من منظور معالجة الإشارات، انحدار الحَرْف هو مُرشّح وينر (Wiener Filter) الذي يوازن بين كبت الضوضاء والحفاظ على الإشارة الأصلية.
أمّا في ، فإنّ اضمحلال الأوزان — أي إضافة إلى الخسارة — هو بالضبط انحدار حَرْف مُطبَّق على معاملات . حين تكتب weight_decay=0.01 في ، فأنت تستخدم الفكرة ذاتها التي نشرها هورل وكينارد عام 1970.
الأثر الأوسع: كيف وُلد مجال التنظيم
أهمية انحدار الحَرْف تتجاوز المسألة التي حلّها. المبدأ الذي أرساه — أنّ تقييد حرية النموذج يُحسّن قدرته على التنبؤ — أعاد تشكيل الإحصاء و بأكمله، وصار الأساس الذي بُنيت عليه كل أدوات التنظيم الحديثة.
1970
انحدار الحَرْف
هورل وكينارد يُثبتان أنّ إضافة kI إلى X'X تُعطي دائماً خطأ تنبؤ أقل من المربعات الصغرى حين تكون المتنبّئات مترابطة. هنا وُلد التنظيم L2.
1996
لاسو (تيبشيراني)
الفكرة: استبدل عقوبة L2 بعقوبة L1. النتيجة: بعض المعاملات تصل إلى الصفر تماماً — وبذلك أصبح التنظيم أداةً لانتقاء المتغيّرات تلقائياً.
2005
الشبكة المرنة
زو وهاستي يمزجان عقوبتَي L1 وL2 في صيغة واحدة، فيجمعان بين انتقاء المتغيّرات والتعامل الجيّد مع المتغيّرات المترابطة.
2012
التعلم العميق + اضمحلال الأوزان
مع AlexNet وما تلاها، أصبح اضمحلال الأوزان (تنظيم L2 على كل المعاملات) ممارسة قياسية في الشبكات العميقة. فكرة الحَرْف انتقلت من بضعة معاملات إلى مليارات.
2017
AdamW — اضمحلال أوزان مفكوك
لوشتشيلوف وهوتر يكتشفان أنّ اضمحلال الأوزان في Adam يعمل بشكل صحيح فقط حين يُفصل عن خطوة تحديث التدرّج. تصحيح عصري لمبدأ عمره عقود.
الحَرْف أم لاسو: كيف تختار بينهما؟
لاسو (عقوبة L1) وانحدار الحَرْف (عقوبة L2) كلاهما يُقلّص المعاملات، لكنّ آلية التقليص مختلفة تماماً:
-
انحدار الحَرْف يُقلّص كل المعاملات بنِسَب متقاربة دون أن يُلغي أيّاً منها. الخيار المناسب حين تعتقد أنّ جميع المتنبّئات تحمل معلومات مفيدة وأنّ المشكلة هي الارتباط بينها — كما في التحليل الطيفي أو بيانات الجينوم ذات المسارات المتداخلة.
-
لاسو يدفع بعض المعاملات إلى الصفر بالكامل، فيعمل كأداة انتقاء للمتغيّرات. الخيار المناسب حين تتوقّع أنّ عدداً قليلاً فقط من المتنبّئات يؤثّر فعلاً وتريد نموذجاً مختصراً.
-
الشبكة المرنة تمزج بين الاثنين: . مفيدة حين تريد انتقاء المتغيّرات مع الحفاظ على التعامل الجيّد مع المتغيّرات المترابطة.
باختصار: الحَرْف يقول «أبقِ الجميع لكن اخفض أصواتهم»، ولاسو يقول «اطرد من لا فائدة منه». عملياً، جرّب الاثنين ودع التحقق المتقاطع يُقرّر.
المرجعHoerl, Kennard. Ridge Regression: Biased Estimation for Nonorthogonal Problems. Technometrics, 1970.
مصطلحات هذه الورقة
- الانحدار الإحصائيRegression
- الضبط الهيكليRegularization
- الانكماشShrinkage
- الموازنة بين الانحياز والتباعدBias-Variance Tradeoff
- تعدُّد الخطية المشتركMulticollinearity
- المربعات الصغرىLeast Squares
- مُعامل الحَرْفRidge Parameter
- متوسط مربع الخطأMean Squared Error
- القيمة الذاتيةEigenvalue
- مُعامل تضخيم التباينVariance Inflation Factor