الانحدار التكيفي متعدد المتغيرات

في الإحصاء ، تعد الشرائح الانحدارية التكيفية متعددة المتغيرات ( MARS ) شكلاً من أشكال تحليل الانحدار الذي قدمه جيروم هـ. فريدمان في عام 1991. [ 1 ] وهي تقنية انحدار غير معلمية ويمكن اعتبارها امتدادًا للنماذج الخطية التي تقوم تلقائيًا بنمذجة اللاخطية والتفاعلات بين المتغيرات.

يُعدّ مصطلح "MARS" علامة تجارية مسجلة ومرخصة لشركة سالفورد سيستمز. ولتجنب انتهاكات حقوق الملكية الفكرية، تُسمى العديد من تطبيقات MARS مفتوحة المصدر باسم "Earth". [ 2 ] [ 3 ]

الأساسيات

يقدم هذا القسم نموذج MARS باستخدام بعض الأمثلة. نبدأ بمجموعة بيانات: مصفوفة من متغيرات الإدخال x ، ومتجه من الاستجابات المرصودة y ، حيث تمثل كل صف في x استجابة . على سبيل المثال، قد تكون البيانات كالتالي:

xy
10.516.4
10.718.8
10.819.7
......
20.677.0

هنا يوجد متغير مستقل واحد فقط ، لذا فإن مصفوفة x تتكون من عمود واحد فقط. بناءً على هذه القياسات، نرغب في بناء نموذج يتنبأ بالقيمة المتوقعة لـ y لقيمة معينة لـ x .

نموذج خطي

النموذج الخطي للبيانات المذكورة أعلاه هو y^=-37+5.1x{\displaystyle {\widehat {y}}=-37+5.1x} القبعة علىy^{\displaystyle {\widehat {y}}}يشير إلى أنy^{\displaystyle {\widehat {y}}}يتم تقديرها من البيانات. يوضح الشكل على اليمين رسمًا بيانيًا لهذه الدالة: خط يعطي القيمة المتوقعة.y^{\displaystyle {\widehat {y}}}مقابل x ، مع عرض القيم الأصلية لـ y كنقاط حمراء.

تشير البيانات عند القيم القصوى لـ x إلى أن العلاقة بين y و x قد تكون غير خطية (انظر إلى النقاط الحمراء بالنسبة لخط الانحدار عند القيم المنخفضة والعالية لـ x ). لذلك، نلجأ إلى برنامج MARS لبناء نموذج تلقائيًا يأخذ في الاعتبار العلاقات غير الخطية. يقوم برنامج MARS ببناء نموذج من قيم x و y المعطاة كما يلي:

y^= 25+6.1الأعلى(0،x-13)-3.1الأعلى(0،13-x){\displaystyle {\begin{aligned}{\widehat {y}}=&\ 25\\&{}+6.1\max(0,x-13)\\&{}-3.1\max(0,13-x)\end{aligned}}}

نموذج MARS بسيط لنفس البيانات

يوضح الشكل على اليمين رسمًا بيانيًا لهذه الدالة: القيمة المتوقعةy^{\displaystyle {\widehat {y}}}بالمقارنة مع x ، مع عرض القيم الأصلية لـ y مرة أخرى كنقاط حمراء. أصبحت الاستجابة المتوقعة الآن أكثر توافقًا مع قيم y الأصلية .

قام برنامج MARS تلقائيًا بإنشاء نقطة انعطاف في قيمة y المتوقعة لمراعاة اللاخطية. وتنتج نقطة الانعطاف هذه بواسطة دوال المفصل . دوال المفصل هي التعبيرات التي تبدأ بـالأعلى{\displaystyle \max }(أينالأعلى(أ،ب){\displaystyle \max(a,b)}يكونأ{\displaystyle a}لوأ>ب{\displaystyle a>b}، آخرب{\displaystyle b}). سيتم شرح وظائف المفصلات بمزيد من التفصيل أدناه.

في هذا المثال البسيط، نلاحظ بسهولة من الرسم البياني أن العلاقة بين y و x غير خطية (وقد نستنتج أن y تتناسب طرديًا مع مربع x ). مع ذلك، عادةً ما يكون هناك عدة متغيرات مستقلة ، وبالتالي ستكون العلاقة بين y وهذه المتغيرات غير واضحة ويصعب تحديدها بالرسم البياني. يمكننا استخدام برنامج MARS لاكتشاف هذه العلاقة غير الخطية.

مثال على تعبير MARS مع متغيرات متعددة هو

ozoنهـ= 5.2+0.93الأعلى(0،تهـمص-58)-0.64الأعلى(0،تهـمص-68)-0.046الأعلى(0،234-أنابت)-0.016الأعلى(0،wأناند-7)الأعلى(0،200-vأناs){\displaystyle {\begin{align}\mathrm {ozone} =&\ 5.2\\&{}+0.93\max(0,\mathrm {temp} -58)\\&{}-0.64\max(0,\mathrm {temp} -68)\\&{}-0.046\max(0,234-\mathrm {ibt}) )\\&{}-0.016\max(0,\mathrm {wind} -7)\max(0,200-\mathrm {vis} )\end{محاذاة}}}

التفاعل بين المتغيرات في نموذج MARS

يُمثل هذا التعبير تلوث الهواء (مستوى الأوزون) كدالة لدرجة الحرارة وبعض المتغيرات الأخرى. لاحظ أن الحد الأخير في الصيغة (في السطر الأخير) يتضمن تفاعلاً بينwأناند{\displaystyle \mathrm {ريح} }وvأناs{\displaystyle \mathrm {vis} }.

يوضح الشكل الموجود على اليمين التوقعاتozoنهـ{\displaystyle \mathrm {الأوزون} }مثلwأناند{\displaystyle \mathrm {ريح} }وvأناs{\displaystyle \mathrm {vis} }تختلف المتغيرات، بينما تبقى المتغيرات الأخرى ثابتة عند قيمها المتوسطة. يوضح الشكل أن الرياح لا تؤثر على مستوى الأوزون إلا في حالة انخفاض مستوى الرؤية. نلاحظ أن نموذج MARS قادر على بناء أسطح انحدار مرنة للغاية من خلال دمج دوال المفصل.

للحصول على التعبير أعلاه، يقوم إجراء بناء نموذج MARS تلقائيًا بتحديد المتغيرات التي سيتم استخدامها (بعض المتغيرات مهمة، والبعض الآخر ليس كذلك)، ومواضع الانحناءات في وظائف المفصل، وكيفية دمج وظائف المفصل.

نموذج المريخ

يقوم المريخ ببناء نماذج من الشكل

و^(x)=أنا=1كجأنابأنا(x).{\displaystyle {\widehat {f}}(x)=\sum _{i=1}^{k}c_{i}B_{i}(x).}

النموذج عبارة عن مجموع مرجح لدوال أساسية بأنا(x){\displaystyle B_{i}(x)}. كلجأنا{\displaystyle c_{i}}هو معامل ثابت. على سبيل المثال، كل سطر في صيغة الأوزون أعلاه هو دالة أساسية مضروبة في معاملها.

كل دالة أساسيةبأنا(x){\displaystyle B_{i}(x)}يتخذ أحد الأشكال الثلاثة التالية:

  1. ثابت يساوي 1. يوجد حد واحد فقط من هذا النوع، وهو الحد الثابت . في معادلة الأوزون أعلاه، يبلغ الحد الثابت 5.2.
  2. دالة مفصلية . تأخذ الدالة المفصلية الشكل التالي: الأعلى(0،x-ثابت){\displaystyle \max(0,x-{\text{constant}})}أو الأعلى(0،ثابت-x){\displaystyle \max(0,{\text{constant}}-x)}يختار برنامج MARS تلقائيًا المتغيرات وقيمها لعقد دوال المفصل. ويمكن رؤية أمثلة على هذه الدوال الأساسية في الأسطر الثلاثة الوسطى من معادلة الأوزون.
  3. ناتج عن دمج دالتين مفصليتين أو أكثر. يمكن لهذه الدوال الأساسية أن تُنمذج التفاعل بين متغيرين أو أكثر.

ومن الأمثلة على ذلك السطر الأخير من معادلة الأوزون.

وظائف المفصلات

زوج متناظر من دوال المفصل مع عقدة عند x=3.1

تُعدّ الدوال المفصلية التي تأخذ الشكل جزءًا أساسيًا من نماذج MARSالأعلى(0،x-ج){\displaystyle \max(0,xc)} أو الأعلى(0،ج-x){\displaystyle \max(0,cx)} أينج{\displaystyle c}هو ثابت يُسمى العقدة . يوضح الشكل على اليمين زوجًا متناظرًا من دوال المفصل مع عقدة عند 3.1.

تكون دالة المفصلة صفرًا لجزء من نطاقها، لذا يمكن استخدامها لتقسيم البيانات إلى مناطق منفصلة، ​​يمكن معالجة كل منها بشكل مستقل. على سبيل المثال، يمكن استخدام زوج متناظر من دوال المفصلة في التعبير 6.1الأعلى(0،x-13)-3.1الأعلى(0،13-x){\displaystyle 6.1\max(0,x-13)-3.1\max(0,13-x)} يقوم بإنشاء الرسم البياني الخطي المجزأ الموضح لنموذج MARS البسيط في القسم السابق.

قد يفترض المرء أنه لا يمكن تكوين سوى الدوال الخطية القطعية من دوال المفصل، ولكن يمكن ضرب دوال المفصل معًا لتكوين دوال غير خطية.

تُسمى الدوال المفصلية أيضاً بالدوال المنحدرة أو دوال عصا الهوكي أو الدوال المعدلة . بدلاً منالأعلى{\displaystyle \max }باستخدام الترميز المستخدم في هذه المقالة، غالبًا ما يتم تمثيل دوال المفصلة بواسطة[±(xأنا-ج)]+{\displaystyle [\pm (x_{i}-c)]_{+}}أين[]+{\displaystyle [\cdot ]_{+}}يعني أخذ الجانب الإيجابي.

عملية بناء النموذج

يقوم برنامج MARS ببناء نموذج على مرحلتين: المرور الأمامي والمرور العكسي. هذا النهج ذو المرحلتين هو نفسه المستخدم في أشجار التقسيم المتكررة .

التمريرة الأمامية

يبدأ نموذج MARS بنموذج يتكون فقط من الحد الثابت (وهو متوسط ​​قيم الاستجابة).

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

  1. المصطلحات الموجودة (تسمى المصطلحات الأصلية في هذا السياق)
  2. جميع المتغيرات (لاختيار متغير واحد لدالة الأساس الجديدة)
  3. جميع قيم كل متغير (لعقدة دالة المفصلة الجديدة).

لحساب معامل كل حد، يطبق برنامج MARS الانحدار الخطي على الحدود.

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

عادةً ما يتم البحث في كل خطوة بطريقة البحث الشامل ، ولكن من أهم جوانب خوارزمية MARS أنه نظرًا لطبيعة دوال المفصل، يمكن إجراء البحث بسرعة باستخدام تقنية تحديث المربعات الصغرى السريعة. ويمكن تسريع البحث الشامل باستخدام طريقة استدلالية تقلل عدد الحدود الأصلية التي يتم النظر فيها في كل خطوة ("Fast MARS" [ 4 ] ).

التمريرة الخلفية

عادةً ما يؤدي التمرير الأمامي إلى فرط التخصيص في النموذج. ولإنشاء نموذج ذي قدرة تعميم أفضل، يقوم التمرير العكسي بتقليص النموذج، وحذف الحد الأقل فعالية في كل خطوة حتى يتم العثور على أفضل نموذج فرعي. تتم مقارنة مجموعات النماذج الفرعية باستخدام معيار التحقق المتقاطع المعمم (GCV) الموضح أدناه.

تتمتع عملية المرور العكسي بميزة على عملية المرور الأمامي: ففي أي خطوة يمكنها اختيار أي مصطلح لحذفه، بينما لا يمكن لعملية المرور الأمامي في كل خطوة رؤية سوى الزوج التالي من المصطلحات.

تضيف عملية التمرير الأمامي الحدود في أزواج، بينما تتجاهل عملية التمرير العكسي عادةً أحد طرفي الزوج، ولذلك غالبًا ما لا تُرى الحدود في أزواج في النموذج النهائي. ويمكن ملاحظة مفصل مزدوج في معادلةy^{\displaystyle {\widehat {y}}}في مثال المريخ الأول أعلاه؛ لا توجد أزواج كاملة محفوظة في مثال الأوزون.

التحقق المتقاطع المعمم

تقارن عملية المرور العكسي أداء النماذج المختلفة باستخدام التحقق المتقاطع المعمم (GCV)، وهو شكل معدل من معيار معلومات أكايكي ، يُقارب درجة التحقق المتقاطع بحذف عنصر واحد في الحالة الخاصة التي تكون فيها الأخطاء غاوسية، أو عند استخدام دالة خسارة مربع الخطأ. وقد طُوّر التحقق المتقاطع المعمم بواسطة كرافن ووهبا، ثم وسّعه فريدمان لنظام MARS؛ وتشير القيم المنخفضة للتحقق المتقاطع المعمم إلى نماذج أفضل. صيغة التحقق المتقاطع المعمم هي:

GCV = RSS / ( N · (1 − (العدد الفعلي للمعلمات) / N ) 2 )

حيث RSS هو مجموع مربعات البواقي المقاسة على بيانات التدريب و N هو عدد الملاحظات (عدد الصفوف في مصفوفة x ).

يُعرَّف العدد الفعال للمعاملات على النحو التالي:

(العدد الفعلي للمعاملات) = (عدد حدود مارس) + (العقوبة) × ((عدد حدود مارس) - 1) / 2

حيث تكون العقوبة عادةً 2 (مما يعطي نتائج تعادل معيار معلومات أكايكي ) ولكن يمكن زيادتها من قبل المستخدم إذا رغب في ذلك.

لاحظ أن

(عدد مصطلحات المريخ - 1) / 2

يمثل عدد عقد دالة المفصل، لذا فإن الصيغة تُعاقب إضافة العقد. وبالتالي، تُعدّل صيغة GCV (أي تزيد) مجموع مربعات الاحتمالية للتدريب لمعاقبة النماذج الأكثر تعقيدًا. نُعاقب المرونة لأن النماذج شديدة المرونة ستُحاكي التجسيد المحدد للضوضاء في البيانات بدلًا من مجرد البنية المنهجية للبيانات.

قيود

تم ذكر أحد القيود بالفعل: يمكن للمستخدم تحديد الحد الأقصى لعدد المصطلحات في عملية التمرير الأمامي.

يمكن فرض قيد إضافي على عملية التمرير الأمامي بتحديد الحد الأقصى المسموح به لدرجة التفاعل. عادةً ما يُسمح بدرجة أو درجتين فقط من التفاعل، ولكن يمكن استخدام درجات أعلى عندما تتطلب البيانات ذلك. الحد الأقصى لدرجة التفاعل في مثال MARS الأول أعلاه هو واحد (أي لا تفاعلات أو نموذج إضافي )؛ وفي مثال الأوزون هو اثنان.

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

الإيجابيات والسلبيات

  • تتميز نماذج MARS بسهولة فهمها وتفسيرها. [ 5 ]
  • يمكن لبرنامج MARS التعامل مع البيانات المستمرة والبيانات الفئوية على حد سواء . [ 6 ] [ 7 ]
  • تعتمد طريقة MARS (مثل التقسيم المتكرر) على الاختيار التلقائي للمتغيرات (أي أنها تُدرج المتغيرات المهمة في النموذج وتستبعد المتغيرات غير المهمة). ومع ذلك، قد يكون هناك بعض العشوائية في عملية الاختيار، خاصةً عند وجود متغيرات تنبؤية مترابطة، وهذا قد يؤثر على قابلية التفسير. [ 5 ]
  • غالباً ما يتطلب بناء نماذج MARS القليل من إعداد البيانات أو لا يتطلب أي إعداد على الإطلاق. [ 5 ]
  • الكود من كتاب الأساليب البايزية للتصنيف والانحدار غير الخطي [ 8 ] لـ Bayesian MARS.
  • يمكن دمج النماذج الخطية المعممة (GLMs) في نماذج MARS عن طريق تطبيق دالة الربط بعد بناء نموذج MARS. وبالتالي، على سبيل المثال، يمكن لنماذج MARS دمج الانحدار اللوجستي للتنبؤ بالاحتمالات.
  • يُستخدم الانحدار غير الخطي عندما يكون الشكل الأساسي للدالة معروفًا، ويُستخدم الانحدار فقط لتقدير معلمات تلك الدالة. أما نموذج MARS، فيُقدّر الدوال نفسها، وإن كان ذلك مع قيود صارمة على طبيعة هذه الدوال. (هذه القيود ضرورية لأن اكتشاف نموذج من البيانات يُعدّ مسألة عكسية غير محددة جيدًا بدون قيود على النموذج).
  • التقسيم المتكرر (المعروف باسم CART). يمكن اعتبار MARS تعميمًا للتقسيم المتكرر يسمح بنماذج متصلة، مما يوفر ملاءمة أفضل للبيانات الرقمية.
  • النماذج الإضافية المعممة . على عكس نماذج الانحدار المتعدد الآلي (MARS)، تستخدم نماذج GAM دوال لويس ناعمة أو دوال متعددة الحدود بدلاً من دوال المفصل، ولا تُنمذج تفاعلات المتغيرات تلقائيًا. يقلل التوافق الأكثر سلاسة وعدم وجود حدود انحدار من التباين مقارنةً بنماذج MARS، لكن تجاهل تفاعلات المتغيرات قد يزيد من التحيز.
  • TSMARS . يُستخدم مصطلح "السلاسل الزمنية للمارس" عند تطبيق نماذج MARS في سياق السلاسل الزمنية . في هذا السياق، عادةً ما تكون المتغيرات التنبؤية هي قيم السلاسل الزمنية المتأخرة، مما ينتج عنه نماذج الانحدار الذاتي التكعيبي . هذه النماذج وامتداداتها، بما في ذلك نماذج الانحدار الذاتي التكعيبي للمتوسط ​​المتحرك، موصوفة في " نمذجة السلاسل الزمنية أحادية المتغير والتنبؤ بها باستخدام TSMARS: دراسة لنماذج الانحدار الذاتي التكعيبي للسلاسل الزمنية، والنماذج الموسمية، ونماذج المتوسط ​​المتحرك باستخدام TSMARS".
  • يستخدم نموذج MARS البايزي (BMARS) نفس شكل النموذج، ولكنه يبني النموذج باستخدام منهج بايزي. وقد ينتج عنه نماذج MARS مثلى مختلفة نظرًا لاختلاف منهجية بناء النموذج. وعادةً ما تكون نتيجة BMARS عبارة عن مجموعة من عينات التوزيع الاحتمالي اللاحق لنماذج MARS، مما يسمح بالتنبؤ الاحتمالي. [ 9 ]

انظر أيضاً

مراجع

  1. فريدمان، جيه إتش (1991). "الانحدار التكيفي متعدد المتغيرات باستخدام الدوال التكعيبية". حوليات الإحصاء . 19 ( 1): 1-67 . CiteSeerX 10.1.1.382.970 . doi : 10.1214/aos/1176347963 . JSTOR 2241837. MR 1091842. Zbl 0765.62064 .    
  2. حزمة CRAN الأرضية
  3. الأرض – دوال الانحدار التكيفي متعدد المتغيرات باللون البرتقالي (مكتبة بايثون للتعلم الآلي)
  4. فريدمان، جيه إتش (1993) برنامج فاست مارس ، قسم الإحصاء بجامعة ستانفورد، التقرير الفني رقم 110
  5. 1 2 3 كون، ماكس؛ جونسون، كيل (2013). النمذجة التنبؤية التطبيقية . نيويورك، نيويورك: سبرينغر نيويورك. doi : 10.1007/978-1-4614-6849-3 . ISBN 9781461468486.
  6. فريدمان، جيروم هـ. (1993). "تقدير الدوال للمتغيرات الترتيبية والفئوية المختلطة باستخدام الدوال التكعيبية التكيفية". في ستيفان مورغنثالر؛ إلفيزيو رونشيتي؛ فيرنر ستاهل (محررون). اتجاهات جديدة في تحليل البيانات الإحصائية والمتانة . بيركهاوزر.
  7. فريدمان، جيروم هـ. (1991-06-01). "تقدير الدوال للمتغيرات الترتيبية والفئوية المختلطة باستخدام الدوال التكعيبية التكيفية" . DTIC . مؤرشف من الأصل في 11 أبريل 2022. تم الاسترجاع في 11 أبريل 2022 .
  8. دينيسون، دي جي تي؛ هولمز، سي سي؛ ماليك، بي كيه؛ سميث، إيه إف إم (2002). الأساليب البايزية للتصنيف والانحدار غير الخطي . تشيتشستر، إنجلترا: وايلي. ISBN 978-0-471-49036-4.
  9. دينيسون، دي جي تي؛ ماليك، بي كيه؛ سميث، إيه إف إم (1 ديسمبر 1998). "Bayesian MARS" (ملف PDF) . الإحصاء والحوسبة . 8 (4): 337-346 . doi : 10.1023/A:1008824606259 . ISSN 1573-1375 . S2CID 12570055 .  

للمزيد من القراءة