التراجع متساوي التوتر

مثال على الانحدار المتساوي التوتر (الخط الأحمر المتصل) مقارنةً بالانحدار الخطي على نفس البيانات، وكلاهما مصمم لتقليل متوسط ​​مربع الخطأ . خاصية الشكل الحر للانحدار المتساوي التوتر تعني أن الخط يمكن أن يكون أكثر انحدارًا حيث تكون البيانات أكثر انحدارًا؛ أما قيد التساوي التوتري فيعني أن الخط لا يتناقص.

في الإحصاء والتحليل العددي ، يعتبر الانحدار المتساوي أو الانحدار الرتيب أسلوبًا لتركيب خط حر الشكل على سلسلة من الملاحظات بحيث يكون الخط الملائم غير متناقص (أو غير متزايد) في كل مكان، ويقع بالقرب من الملاحظات قدر الإمكان.

التطبيقات

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

من التطبيقات الأخرى التحجيم متعدد الأبعاد غير المتري [ 1 ] ، حيث يُبحث عن تمثيل منخفض الأبعاد لنقاط البيانات بحيث يتطابق ترتيب المسافات بين النقاط في التمثيل مع ترتيب الاختلاف بينها. ويُستخدم الانحدار المتساوي بشكل تكراري لضبط المسافات المثالية للحفاظ على ترتيب الاختلاف النسبي.

يُستخدم الانحدار المتساوي التوتر أيضًا في التصنيف الاحتمالي لمعايرة الاحتمالات المتوقعة لنماذج التعلم الآلي الخاضعة للإشراف . [ 2 ]

الانحدار المتساوي التوتر للحالة المرتبة ببساطة مع متغير واحدx،y{\displaystyle x,y}وقد طُبقت هذه الطريقة لتقدير علاقات الاستجابة للجرعة المستمرة في مجالات مثل التخدير وعلم السموم. وبشكل أدق، لا يوفر الانحدار متساوي التوتر سوى تقديرات نقطية عند القيم المرصودة لـx.{\displaystyle x.}يتم عادةً تقدير منحنى الاستجابة للجرعة الكامل دون أي افتراضات إضافية عن طريق الاستيفاء الخطي بين تقديرات النقاط. [ 3 ]

تم تطوير برامج لحساب الانحدار المتساوي (الرتيب) للغات R ، [ 4 ] [ 5 ] [ 6 ] و Stata و Python . [ 7 ]

بيان المشكلة والخوارزميات

يترك(x1،y1)،...،(xن،yن){\displaystyle (x_{1},y_{1}),\ldots ,(x_{n},y_{n})}لنفترض مجموعة معينة من الملاحظات، حيثyأناR{\displaystyle y_{i}\in \mathbb {R} }وxأنا{\displaystyle x_{i}}تقع ضمن مجموعة مرتبة جزئياً . وللتعميم، كل ملاحظة(xأنا،yأنا){\displaystyle (x_{i},y_{i})}قد يُعطى وزناًwأنا0{\displaystyle w_{i}\geq 0}، على الرغم من شيوعهاwأنا=1{\displaystyle w_{i}=1}للجميعأنا{\displaystyle i}.

يسعى الانحدار المتساوي التوتر إلى إيجاد أفضل ملاءمة باستخدام طريقة المربعات الصغرى الموزونةy^أناyأنا{\displaystyle {\hat {y}}_{i}\approx y_{i}}للجميعأنا{\displaystyle i}مع مراعاة القيد التالي:y^أناy^ج{\displaystyle {\hat {y}}_{i}\leq {\hat {y}}_{j}}حينماxأناxج{\displaystyle x_{i}\leq x_{j}}وهذا يُعطي البرنامج التربيعي التالي (QP) بدلالة المتغيراتy^1،...،y^ن{\displaystyle {\hat {y}}_{1},\ldots ,{\hat {y}}_{n}}:

مينأنا=1نwأنا(y^أنا-yأنا)2{\displaystyle \min \sum _{i=1}^{n}w_{i}({\hat {y}}_{i}-y_{i})^{2}}رهناً بـy^أناy^ج للجميع (أنا،ج)هـ{\displaystyle {\hat {y}}_{i}\leq {\hat {y}}_{j}{\text{ لجميع }}(i,j)\in E}

أينهـ={(أنا،ج):xأناxج}{\displaystyle E=\{(i,j):x_{i}\leq x_{j}\}}يحدد الترتيب الجزئي للمدخلات المرصودةxأنا{\displaystyle x_{i}}(ويمكن اعتبارها مجموعة حواف رسم بياني موجه غير دوري (dag) ذي رؤوس1،2،...ن{\displaystyle 1,2,\ldots n}). يمكن حل المشكلات من هذا النوع باستخدام تقنيات البرمجة التربيعية العامة.

في السياق المعتاد حيثxأنا{\displaystyle x_{i}}تندرج القيم ضمن مجموعة مرتبة ترتيباً كلياً مثلR{\displaystyle \mathbb {R} }، يمكننا أن نفترض دون فقدان العمومية أن المشاهدات قد تم فرزها بحيثx1x2xن{\displaystyle x_{1}\leq x_{2}\leq \cdots \leq x_{n}}وخذهـ={(أنا،أنا+1):1أنا<ن}{\displaystyle E=\{(i,i+1):1\leq i<n\}}في هذه الحالة، تُعد خوارزمية تجميع المخالفين المتجاورين خوارزمية تكرارية بسيطة لحل البرنامج التربيعي . في المقابل، درس بيست وشاكرافارتي [ 8 ] المشكلة باعتبارها مشكلة تحديد مجموعة فعالة ، واقترحا خوارزمية أولية. يمكن اعتبار هاتين الخوارزميتين متقابلتين، ولكل منهما تعقيد حسابي قدرهيا(ن){\displaystyle O(n)}على البيانات المصنفة مسبقًا. [ 8 ]

لإكمال مهمة الانحدار المتساوي التوتر، يمكننا بعد ذلك اختيار أي دالة غير متناقصةو(x){\displaystyle f(x)}بحيثو(xأنا)=y^أنا{\displaystyle f(x_{i})={\hat {y}}_{i}}لكل i. من الواضح أن أي دالة من هذا القبيل تحل

مينوأنا=1نwأنا(و(xأنا)-yأنا)2{\displaystyle \min _{f}\sum _{i=1}^{n}w_{i}(f(x_{i})-y_{i})^{2}}رهناً بـو{\displaystyle f}كونها غير متناقصة

ويمكن استخدامها للتنبؤ بـy{\displaystyle y}قيم القيم الجديدة لـx{\displaystyle x}خيار شائع عندماxأناR{\displaystyle x_{i}\in \mathbb {R} }يتمثل الحل في إجراء استيفاء خطي بين النقاط(xأنا،y^أنا){\displaystyle (x_{i},{\hat {y}}_{i})}كما هو موضح في الشكل، مما ينتج عنه دالة خطية متصلة ومجزأة:

و(x)={y^1لو xx1y^أنا+x-xأناxأنا+1-xأنا(y^أنا+1-y^أنا)لو xأناxxأنا+1y^نلو xxن{\displaystyle f(x)={\begin{cases}{\hat {y}}_{1}&{\text{إذا كان }}x\leq x_{1}\\{\hat {y}}_{i}+{\frac {x-x_{i}}{x_{i+1}-x_{i}}}({\hat {y}}_{i+1}-{\hat {y}}_{i})&{\text{إذا كان }}x_{i}\leq x\leq x_{i+1}\\{\hat {y}}_{n}&{\text{إذا كان }}x\geq x_{n}\end{cases}}}

الانحدار المتساوي التوتر المركزي

كما يوضح الشكل الأول في هذه المقالة، في حالة وجود انتهاكات للرتابة، فإن المنحنى الناتج عن الاستيفاء سيكون له فترات ثابتة. ومن المعروف عادةً في تطبيقات استجابة الجرعة أنو(x){\displaystyle f(x)}ليست رتيبة فحسب، بل سلسة أيضًا . الفواصل الموسيقية المسطحة لا تتوافق معو(x){\displaystyle f(x)}الشكل المفترض، ويمكن إثبات تحيزه. وقد طوّر أورون وفلورنوي تحسينًا بسيطًا لمثل هذه التطبيقات، يُسمى الانحدار المتساوي التوتر المركزي (CIR)، وأثبت أنه يقلل بشكل كبير من خطأ التقدير في كل من تطبيقات استجابة الجرعة وتحديد الجرعة. [ 9 ] يتم تطبيق كل من CIR والانحدار المتساوي التوتر القياسي للحالة أحادية المتغير والمرتبة ببساطة، في حزمة R المسماة "cir". [ 4 ] توفر هذه الحزمة أيضًا تقديرات فترات الثقة التحليلية.

مراجع

  1. كروسكال، ج. ب. (1964). "القياس متعدد الأبعاد غير المتري: طريقة عددية". مجلة القياس النفسي . 29 (2): 115-129 . doi : 10.1007/BF02289694 . S2CID 11709679 . 
  2. نيكوليسكو-ميزيل، ألكساندرو؛ كاروانا، ريتش (2005). "التنبؤ باحتمالات جيدة باستخدام التعلم الخاضع للإشراف". في: دي رايدت، لوك؛ وروبل، ستيفان (محرران). وقائع المؤتمر الدولي الثاني والعشرين للتعلم الآلي (ICML 2005)، بون، ألمانيا، 7-11 أغسطس 2005. سلسلة وقائع المؤتمرات الدولية لجمعية آلات الحوسبة. المجلد 119. جمعية آلات الحوسبة. الصفحات 625-632 . doi : 10.1145/1102351.1102430 .  
  3. ستيليانو، إم بي؛ فلورنوي، إن (2002). "تحديد الجرعة باستخدام تصميم العملة المتحيزة صعودًا وهبوطًا والانحدار المتساوي التوتر". القياسات الحيوية . 58 (1): 171-177 . doi : 10.1111/j.0006-341x.2002.00171.x . PMID 11890313. S2CID 8743090 .  
  4. 1 2 أورون، عساف. "حزمة 'cir'"" . CRAN . مؤسسة R للحوسبة الإحصائية . تم الاطلاع عليه بتاريخ 26 ديسمبر 2020 .
  5. ليو، جان دي؛ هورنيك، كورت؛ ماير، باتريك (2009). "تحسين التساوي في R: خوارزمية المخالفين المجاورين المجمعين (PAVA) وطرق المجموعة النشطة" . مجلة البرمجيات الإحصائية . 32 (5): 1-24 . doi : 10.18637/jss.v032.i05 . ISSN 1548-7660 . 
  6. ^ شو، تشيبنغ؛ صن، تشينكاي؛ كاروناكاران، أمان. "حزمة UniIsoRegression" (PDF) . كران . مؤسسة R للحوسبة الإحصائية . تم الاسترجاع في 29 أكتوبر 2021 .
  7. بيدريجوسا، فابيان؛ وآخرون (2011). "Scikit-learn: التعلم الآلي في بايثون". مجلة أبحاث التعلم الآلي . 12 : 2825-2830 . arXiv : 1201.0490 . Bibcode : 2011JMLR...12.2825P . 
  8. 1 2 بيست، مايكل جيه؛ تشاكرافارتي، نيلوتبال (1990). "خوارزميات المجموعة النشطة للانحدار المتساوي التوتر؛ إطار موحد" . البرمجة الرياضية . 47 ( 1-3 ): 425-439 . doi : 10.1007/bf01580873 . ISSN 0025-5610 . S2CID 31879613 .  
  9. أورون، أ.ب.؛ فلورنوي، ن. (2017). "الانحدار المتساوي التوتر المركزي: تقدير النقطة والفترة لدراسات استجابة الجرعة". الإحصاء في البحوث الصيدلانية الحيوية . 9 (3): 258-267 . arXiv : 1701.05964 . doi : 10.1080/19466315.2017.1286256 . S2CID 88521189 . 

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

  • روبرتسون، ت.؛ رايت، ف. ت.؛ ديكسترا، ر. ل. (1988). الاستدلال الإحصائي المقيد بالترتيب . نيويورك: وايلي. ISBN 978-0-471-91787-8.
  • بارلو، ر. إي.؛ بارثولوميو، د. ج.؛ بريمنر، ج. م.؛ برونك، هـ. د. (1972). الاستدلال الإحصائي في ظل قيود الترتيب؛ نظرية وتطبيق الانحدار المتساوي . نيويورك: وايلي. ISBN 978-0-471-04970-8.
  • شيفلي، تي إس، ساجر، تي دبليو، ووكر، إس جي (2009). "نهج بايزي لتقدير الدوال الرتيبة غير البارامترية". مجلة الجمعية الإحصائية الملكية، السلسلة ب . 71 (1): 159-175 . CiteSeerX 10.1.1.338.3846 . doi : 10.1111/j.1467-9868.2008.00677.x . S2CID 119761196 .  {{cite journal}}: صيانة CS1: أسماء متعددة: قائمة المؤلفين ( رابط )
  • وو، دبليو بي ؛ وودروف، إم ؛ مينتز، جي (2001). "الانحدار المتساوي التوتر: نظرة أخرى على مشكلة نقطة التغير". بيومتريكا . 88 (3): 793-804 . doi : 10.1093/biomet/88.3.793 .