الانحدار القوي
في الإحصاء القوي ، يسعى الانحدار القوي إلى التغلب على بعض قيود تحليل الانحدار التقليدي . يُنمذج تحليل الانحدار العلاقة بين متغير مستقل واحد أو أكثر ومتغير تابع . تتمتع أنواع الانحدار القياسية، مثل المربعات الصغرى العادية ، بخصائص جيدة إذا كانت افتراضاتها الأساسية صحيحة، ولكنها قد تُعطي نتائج مُضللة في غير ذلك (أي أنها ليست قوية في مواجهة انتهاكات الافتراضات). صُممت طرق الانحدار القوي للحد من تأثير انتهاكات الافتراضات من قِبل عملية توليد البيانات الأساسية على تقديرات الانحدار.
على سبيل المثال، تتأثر تقديرات المربعات الصغرى لنماذج الانحدار بشدة بالقيم الشاذة : فالقيمة الشاذة التي يبلغ حجم خطأها ضعف حجم خطأ الملاحظة النموذجية تُسهم بأربعة أضعاف (مربعين) في خسارة مربع الخطأ ، وبالتالي يكون لها تأثير أكبر على تقديرات الانحدار. تُعد دالة خسارة هوبر بديلاً قوياً لخسارة مربع الخطأ القياسية، حيث تُقلل من مساهمة القيم الشاذة في خسارة مربع الخطأ، مما يحد من تأثيرها على تقديرات الانحدار.
التطبيقات
أخطاء التباين غير المتجانس
من الحالات التي يُنصح فيها باللجوء إلى التقدير القوي وجود شك قوي في تباين التباين . في نموذج التباين المتجانس ، يُفترض أن تباين حد الخطأ ثابت لجميع قيم x . يسمح تباين التباين بأن يكون التباين تابعًا لـ x ، وهو ما يُعد أكثر دقة في العديد من السيناريوهات الواقعية. على سبيل المثال، غالبًا ما يكون تباين الإنفاق أكبر لدى الأفراد ذوي الدخل المرتفع مقارنةً بالأفراد ذوي الدخل المنخفض. عادةً ما تستخدم حزم البرامج نموذج التباين المتجانس افتراضيًا، على الرغم من أن هذا النموذج قد يكون أقل دقة من نموذج التباين غير المتجانس. يتمثل أحد الأساليب البسيطة ( توفاليس، 2008 ) في تطبيق طريقة المربعات الصغرى على الأخطاء النسبية، حيث يقلل ذلك من تأثير القيم الأكبر للمتغير التابع مقارنةً بطريقة المربعات الصغرى العادية.
وجود القيم المتطرفة
يُستخدم التقدير القوي في حالات أخرى شائعة، لا سيما عند احتواء البيانات على قيم متطرفة. ففي وجود قيم متطرفة لا تنتمي إلى نفس عملية توليد البيانات، يصبح تقدير المربعات الصغرى غير فعال وقد يكون متحيزًا. ونظرًا لأن تنبؤات المربعات الصغرى تميل نحو القيم المتطرفة، ولأن تباين التقديرات يتضخم بشكل مصطنع، فإن النتيجة هي إمكانية إخفاء القيم المتطرفة. (في كثير من الحالات، بما في ذلك بعض مجالات الإحصاء الجغرافي والإحصاء الطبي، تكون القيم المتطرفة هي محل الاهتمام تحديدًا).
على الرغم من الادعاء أحيانًا بأن طريقة المربعات الصغرى (أو الطرق الإحصائية الكلاسيكية عمومًا) تتمتع بمتانة عالية، إلا أن هذه المتانة تقتصر على عدم ارتفاع معدل الخطأ من النوع الأول عند وجود انتهاكات للنموذج. في الواقع، يميل معدل الخطأ من النوع الأول إلى أن يكون أقل من المستوى الاسمي عند وجود قيم متطرفة، وغالبًا ما يرتفع معدل الخطأ من النوع الثاني بشكل ملحوظ . وقد وُصف انخفاض معدل الخطأ من النوع الأول بأنه تحفظ الطرق الكلاسيكية.
تاريخ وعدم شعبية الانحدار القوي
على الرغم من تفوقها على تقدير المربعات الصغرى في العديد من الحالات، لا تزال الطرق القوية للانحدار غير شائعة الاستخدام. قد تُفسر عدة أسباب عدم شيوعها ( هامبل وآخرون، 1986، 2005 ). أحدها وجود العديد من الطرق المنافسة ، وبداية هذا المجال المتعثرة. كما أن التقديرات القوية تتطلب موارد حاسوبية أكبر بكثير من تقدير المربعات الصغرى ؛ إلا أن هذا الاعتراض أصبح أقل أهمية في السنوات الأخيرة مع ازدياد القدرة الحاسوبية بشكل كبير. سبب آخر محتمل هو فشل بعض حزم البرامج الإحصائية الشائعة في تطبيق هذه الطرق ( سترومبرغ، 2004 ). ولعل أهم أسباب عدم شيوع طرق الانحدار القوية هو أنه عندما يكون تباين الخطأ كبيرًا جدًا أو معدومًا، فإن أي تقدير لمعاملات الانحدار، سواء كان قويًا أم لا، سيكون عديم الجدوى عمليًا لأي مجموعة بيانات معينة، ما لم تكن العينة كبيرة جدًا.
على الرغم من بطء تبني الأساليب القوية، فإن كتب الإحصاء الحديثة السائدة غالبًا ما تتضمن مناقشة لهذه الأساليب (على سبيل المثال، كتب سيبر ولي، وفارواي ؛ وللحصول على وصف عام جيد لكيفية تطور أساليب الانحدار القوي المختلفة من بعضها البعض، انظر كتاب أندرسن ). كما تتضمن حزم البرامج الإحصائية الحديثة مثل R و SAS وStatsmodels و Stata و S-PLUS وظائف كبيرة للتقدير القوي (انظر، على سبيل المثال، كتب فينابلز وريبلي، ومارونا وآخرون ).
أساليب الانحدار القوي
بدائل المربعات الصغرى
أبسط طرق تقدير المعلمات في نموذج الانحدار، والتي تكون أقل حساسية للقيم المتطرفة من تقديرات المربعات الصغرى، هي استخدام أقل الانحرافات المطلقة . ومع ذلك، لا تزال القيم المتطرفة الكبيرة قادرة على التأثير بشكل ملحوظ على النموذج، مما يحفز البحث عن مناهج أكثر قوة.
في عام 1964، قدّم هوبر تقدير M للانحدار. يشير الحرف M في تقدير M إلى "نوع الاحتمال الأقصى". تتميز هذه الطريقة بمتانتها في مواجهة القيم المتطرفة في المتغير التابع، ولكن تبيّن أنها غير مقاومة للقيم المتطرفة في المتغيرات المستقلة ( نقاط التأثير ). في الواقع، عند وجود قيم متطرفة في المتغيرات المستقلة، لا تتمتع هذه الطريقة بأي ميزة على طريقة المربعات الصغرى.
في ثمانينيات القرن الماضي، طُرحت عدة بدائل لتقدير M في محاولة للتغلب على ضعف مقاومته. للاطلاع على مراجعة عملية، يُرجى مراجعة كتاب روسيو وليروي . يُعدّ تقدير المربعات الصغرى المُقَلَّمة (LTS) بديلاً قابلاً للتطبيق، وهو الخيار المُفضَّل حاليًا (2007) وفقًا لروسيو وريان (1997، 2008). يتميز مُقدِّر ثيل-سين بنقطة انهيار أقل من LTS، ولكنه فعال إحصائيًا وشائع الاستخدام. من الحلول المقترحة الأخرى تقدير S. تُحدد هذه الطريقة خطًا (مستوى أو مستوى فائق) يُقلل من تقدير قوي لمقياس البواقي (الذي يُشتق منه حرف S في اسم الطريقة). تتميز هذه الطريقة بمقاومة عالية لنقاط التأثير، كما أنها قوية في مواجهة القيم الشاذة في الاستجابة. مع ذلك، وُجد أن هذه الطريقة غير فعالة أيضًا.
تسعى طريقة التقدير MM إلى الحفاظ على متانة ومقاومة طريقة التقدير S، مع اكتساب كفاءة طريقة التقدير M. تعتمد هذه الطريقة على إيجاد تقدير S عالي المتانة والمقاومة يقلل من تقدير M لمقياس البواقي (الحرف M الأول في اسم الطريقة). ثم يُثبَّت المقياس المُقدَّر أثناء البحث عن تقدير M قريب للمعلمات (الحرف M الثاني).
البدائل البارامترية
يتمثل أحد الأساليب الأخرى لتقدير نماذج الانحدار بدقة في استبدال التوزيع الطبيعي بتوزيع ذي ذيل سميك. وقد أشارت الدراسات إلى أن توزيع t بدرجات حرية تتراوح بين 4 و6 يُعد خيارًا مناسبًا في العديد من الحالات العملية. ويعتمد الانحدار البايزي القوي، كونه بارامتريًا بالكامل، اعتمادًا كبيرًا على هذه التوزيعات.
بافتراض أن البواقي موزعة وفقًا لتوزيع t ، فإن التوزيع ينتمي إلى عائلة التوزيع المكاني-المقياسي. أي،تُسمى درجات حرية توزيع t أحيانًا بمعامل التفرطح . ناقش لانج وليتل وتايلور (1989) هذا النموذج بتفصيل من منظور غير بايزي. أما العرض البايزي فيظهر في جيلمان وآخرون (2003).
يتمثل أحد الأساليب البارامترية البديلة في افتراض أن البواقي تتبع مزيجًا من التوزيعات الطبيعية ( دايمي وآخرون، 2019 )؛ وتحديدًا، توزيعًا طبيعيًا ملوثًا حيث تنتمي غالبية المشاهدات إلى توزيع طبيعي محدد، بينما تنتمي نسبة صغيرة منها إلى توزيع طبيعي ذي تباين أعلى بكثير. أي أن للبواقي احتماليةمن التوزيع الطبيعي بتباين، أينصغير، والاحتمالمن التوزيع الطبيعي بتباينبالنسبة للبعض:
عادة،يُطلق على هذا أحيانًا اسمنموذج التلوث.
تتميز الأساليب البارامترية بميزة أن نظرية الاحتمالية توفر منهجًا جاهزًا للاستدلال (على الرغم من أن نماذج الخليط مثل...في نموذج التلوث (حيث قد لا تنطبق شروط الانتظام المعتادة)، من الممكن بناء نماذج محاكاة من خلال المطابقة. ومع ذلك، تفترض هذه النماذج البارامترية أن النموذج الأساسي صحيح حرفيًا. ولذلك، فهي لا تأخذ في الحسبان التوزيعات المتبقية المنحرفة أو دقة الملاحظة المحدودة.
أوزان الوحدة
من الطرق الفعّالة الأخرى استخدام الأوزان الموحدة ( واينر وثيسن، 1976)، وهي طريقة يمكن تطبيقها عند وجود عدة عوامل تنبؤية لنتيجة واحدة. استخدم إرنست بورغيس (1928) الأوزان الموحدة للتنبؤ بنجاح الإفراج المشروط. قيّم 21 عاملاً إيجابياً على أنها موجودة (مثلاً، "عدم وجود سوابق اعتقال" = 1) أو غائبة ("وجود سوابق اعتقال" = 0)، ثم جمعها للحصول على درجة تنبؤية، والتي ثبت أنها مؤشر مفيد لنجاح الإفراج المشروط. أظهر صموئيل س. ويلكس (1938) أن جميع مجموعات أوزان الانحدار تقريباً، بما في ذلك الأوزان الموحدة، تُجمع لتكوين قيم مركبة شديدة الارتباط ببعضها البعض، وهي نتيجة تُعرف أحياناً بنظرية ويلكس (ري، كاريتا، وإيرلز، 1998). درست روبين داوز (1979) عملية صنع القرار في سياقات تطبيقية، موضحةً أن النماذج البسيطة ذات الأوزان الموحدة غالباً ما تتفوق على الخبراء البشريين. قام بوبكو وروث وباستر (2007) بمراجعة الأدبيات المتعلقة بأوزان الوحدة وخلصوا إلى أن عقودًا من الدراسات التجريبية تُظهر أن أوزان الوحدة تؤدي أداءً مشابهًا لأوزان الانحدار العادية في التحقق المتبادل.
مثال: بيانات الكبد من شركة بوبا
تمت دراسة بيانات الكبد الخاصة بمشروع BUPA من قبل العديد من الباحثين، بمن فيهم بريمان (2001). يمكن الاطلاع على البيانات في صفحة مجموعات البيانات الكلاسيكية ، ويتضمن المقال نقاشًا حول تحويل بوكس-كوكس . يظهر أدناه رسم بياني للوغاريتمات ALT مقابل لوغاريتمات γGT. خطا الانحدار هما الخطان المُقدَّران باستخدام طريقة المربعات الصغرى العادية (OLS) وتقدير MM القوي. أُجري التحليل باستخدام برنامج R الذي وفره فينابلز وريبلي (2002).
يبدو خطّا الانحدار متشابهين للغاية (وهذا ليس بالأمر غير المألوف في مجموعة بيانات بهذا الحجم). مع ذلك، تبرز ميزة الأسلوب القوي عند النظر في تقديرات مقياس البواقي. ففي طريقة المربعات الصغرى العادية، يبلغ تقدير المقياس 0.420، مقارنةً بـ 0.373 في الطريقة القوية. وبالتالي، تبلغ الكفاءة النسبية لطريقة المربعات الصغرى العادية مقارنةً بتقدير MM في هذا المثال 1.266. ويؤدي هذا القصور إلى فقدان القدرة الإحصائية في اختبارات الفرضيات، وإلى فترات ثقة واسعة بشكل غير ضروري للمعلمات المُقدَّرة.
اكتشاف القيم الشاذة
من النتائج الأخرى لعدم كفاءة طريقة المربعات الصغرى العادية إخفاء العديد من القيم الشاذة نتيجة لتضخيم تقدير مقياس البواقي؛ إذ تقترب البواقي المُقاسة من الصفر مقارنةً باستخدام تقدير أكثر دقة للمقياس. تظهر أدناه رسوم بيانية للبواقي المُقاسة من النموذجين. يُمثل المتغير على المحور السيني رقم المشاهدة كما ورد في مجموعة البيانات. يحتوي كتاب روسيو وليروي (1986) على العديد من هذه الرسوم البيانية.
تقع خطوط المرجعية الأفقية عند 2 و -2، لذا يمكن اعتبار أي قيمة متبقية مُقاسة مُلاحظة خارج هذه الحدود قيمةً شاذة. من الواضح أن طريقة المربعات الصغرى تؤدي إلى إخفاء العديد من الملاحظات المهمة.
بينما يمكن الكشف عن القيم الشاذة يدويًا في بُعد أو بُعدين باستخدام الطرق التقليدية، إلا أن مشكلة الحجب قد تجعل تحديد العديد من القيم الشاذة مستحيلاً مع مجموعات البيانات الكبيرة وفي الأبعاد العالية. تكتشف الطرق القوية هذه القيم تلقائيًا، مما يوفر ميزة كبيرة على الطرق التقليدية عند وجود قيم شاذة.
انظر أيضاً
مراجع
- ليو، جيه؛ كوسمان، بي سي؛ راو، بي دي (2018). "الانحدار الخطي القوي عبر تنظيم L0" . معاملات IEEE في معالجة الإشارات . 66 (3): 698-713 . doi : 10.1109/TSP.2017.2771720 .
- أندرسن، ر. (2008). الأساليب الحديثة للانحدار القوي . سلسلة أوراق جامعة سيج حول التطبيقات الكمية في العلوم الاجتماعية، 07-152.
- بن غال، آي.، اكتشاف القيم الشاذة، مؤرشف في 15 ديسمبر 2022 على موقع Wayback Machine ، في: ميمون، أو. وروكاش، ل. (محرران)، دليل استخراج البيانات واكتشاف المعرفة: دليل شامل للممارسين والباحثين، دار نشر كلوير الأكاديمية، 2005، رقم ISBN 0-387-24435-2.
- بوبكو، ب.، روث، ب. ل.، وباستر، م. أ. (2007). "جدوى استخدام أوزان الوحدة في إنشاء الدرجات المركبة: مراجعة أدبية، وتطبيق على صلاحية المحتوى، وتحليل تلوي". أساليب البحث التنظيمي ، المجلد 10، الصفحات 689-709. doi : 10.1177/1094428106294734
- دائمي، عاطفة، هاريبراساد كودامانا، وبياو هوانغ. "نمذجة العمليات الغاوسية باستخدام احتمالية خليط غاوسي". مجلة التحكم في العمليات 81 (2019): 209-220. doi : 10.1016/j.jprocont.2019.06.007
- بريمان، ل. (2001). "النمذجة الإحصائية: الثقافتان" . العلوم الإحصائية . 16 (3): 199-231 . doi : 10.1214/ss/1009213725 . JSTOR 2676681 .
- بورغيس، إي دبليو (1928). "العوامل المحددة للنجاح أو الفشل في الإفراج المشروط". في: أ. أ. بروس (محرر)، آليات قانون الأحكام غير المحددة والإفراج المشروط في إلينوي (ص 205-249). سبرينغفيلد، إلينوي: مجلس الإفراج المشروط لولاية إلينوي. كتب جوجل
- داوز، روبين م. (1979). "الجمال القوي للنماذج الخطية غير المناسبة في صنع القرار". عالم النفس الأمريكي ، المجلد 34، الصفحات 571-582. doi : 10.1037/0003-066X.34.7.571 . ملف PDF مؤرشف
- دريبر، ديفيد (1988). "التحليل القوي للنماذج الخطية القائم على الرتب. الجزء الأول: العرض والمراجعة" . العلوم الإحصائية . 3 (2): 239-257 . doi : 10.1214/ss/1177012915 . JSTOR 2245578 .
- Faraway, JJ (2004). النماذج الخطية باستخدام R. Chapman & Hall/CRC.
- فورنالسكي، ك. و. (2015). "تطبيقات تحليل الانحدار البايزي القوي". المجلة الدولية لعلوم أنظمة المجتمع . 7 (4): 314-333 . doi : 10.1504/IJSSS.2015.073223 .
- جيلمان، أ.؛ جيه بي كارلين؛ إتش إس ستيرن؛ دي بي روبين (2003). تحليل البيانات البايزية ( الطبعة الثانية). تشابمان آند هول/سي آر سي.
- هامبل، الأب؛ إم رونشيتي . بي جي روسيو؛ وا ستاهيل (2005) [1986]. إحصائيات قوية: النهج القائم على وظائف التأثير . وايلي.
- لانج، ك. ل.؛ ليتل، ر. ج. أ.؛ تايلور، ج. م. ج. (1989). "النمذجة الإحصائية القوية باستخدام توزيع t " . مجلة الجمعية الإحصائية الأمريكية . 84 (408): 881-896 . doi : 10.2307/2290063 . JSTOR 2290063 .
- ليرمان، جي.؛ مكوي، إم.؛ تروب، جيه إيه؛ تشانغ تي. (2012). "الحساب القوي للنماذج الخطية، أو كيفية العثور على إبرة في كومة قش" مؤرشف في 2013-09-26 في Wayback Machine ، arXiv : 1202.4044 .
- مارونا، ر.؛ د. مارتن؛ ف. يوهاي (2006). الإحصاءات القوية: النظرية والأساليب . وايلي.
- مكين، جوزيف و. (2004). "التحليل القوي للنماذج الخطية" . العلوم الإحصائية . 19 (4): 562-570 . doi : 10.1214/088342304000000549 . JSTOR 4144426 .
- رادشينكو، إس. جي. (2005). الأساليب القوية لتقدير النماذج الإحصائية: دراسة. (باللغة الروسية) . كييف: دار النشر الروسية «سانسباريل». ص 504. ISBN 978-966-96574-0-4.
- ري، إم جيه، كاريتا، تي آر، وإيرلز، جيه إيه (1998). "في القرارات من أعلى إلى أسفل، لا يهم ترجيح المتغيرات: نتيجة لنظرية ويلكس". أساليب البحث التنظيمي ، المجلد 1 (4)، الصفحات 407-420. doi : 10.1177/109442819814003
- روسيو، بي جيه ؛ إيه إم ليروي (2003) [1986]. الانحدار القوي والكشف عن القيم الشاذة . وايلي.
- ريان، تي بي (2008) [1997]. أساليب الانحدار الحديثة . وايلي.
- سيبر، جي إيه إف؛ إيه جيه لي (2003). تحليل الانحدار الخطي ( الطبعة الثانية). وايلي.
- سترومبرغ، أ. ج. (2004). "لماذا نكتب برامج إحصائية؟ حالة الأساليب الإحصائية القوية" . مجلة البرمجيات الإحصائية . 10 (5). doi : 10.18637/jss.v010.i05 .
- ستروتز، ت. (2016). ملاءمة البيانات وعدم اليقين (مقدمة عملية للمربعات الصغرى الموزونة وما بعدها) . سبرينغر فيو. ISBN 978-3-658-11455-8.
- توفاليس، كريس (2008). "انحدار النسبة المئوية للمربعات الصغرى" . مجلة الأساليب الإحصائية التطبيقية الحديثة . 7 : 526-534 . doi : 10.2139/ssrn.1406472 . hdl : 2299/965 . SSRN 1406472 .
- فينابلز، دبليو إن؛ بي دي ريبلي (2002). الإحصاء التطبيقي الحديث باستخدام S. سبرينغر.
- واينر، هـ. ، وثيسن، د. (1976). "ثلاث خطوات نحو الانحدار القوي". مجلة القياس النفسي ، المجلد 41 (1)، الصفحات 9-34. doi : 10.1007/BF02291695
- ويلكس، إس إس (1938). "أنظمة الترجيح للدوال الخطية للمتغيرات المترابطة عندما لا يوجد متغير تابع". مجلة القياس النفسي ، المجلد 3، الصفحات 23-40. doi : 10.1007/BF02287917
روابط خارجية
- ويكي كتب برمجة R
- ملاحظات بريان ريبلي القوية حول دورة الإحصاء.
- تحتوي ملاحظات نيك فيلر الدراسية حول النمذجة الإحصائية والحساب، المؤرشفة في 2016-03-03 في Wayback Machine، على مواد حول الانحدار القوي.
- نظرة عامة على الإحصاءات القوية لأولفا نصرواي، مؤرشفة بتاريخ 1 أبريل 2022 في أرشيف الإنترنت (Wayback Machine).
- نظرة عامة على التجميع القوي من إعداد أولفا نصرواي، مؤرشفة بتاريخ 1 أبريل 2022 في أرشيف الإنترنت (Wayback Machine).
- لماذا نكتب برامج إحصائية؟ حالة الأساليب الإحصائية القوية، بقلم أ. ج. سترومبرغ. مؤرشف بتاريخ 7 فبراير 2012 في أرشيف الإنترنت.
- برنامج مجاني (Fortran 95) لانحدار L1-norm. تقليل الانحرافات المطلقة بدلاً من المربعات الصغرى.
- تطبيق بايثون مجاني مفتوح المصدر للانحدار غير الخطي القوي.
- الانحدار القوي
