تراجع الحافة

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

طُرحت هذه النظرية لأول مرة من قِبل هورل وكينارد عام 1970 في ورقتيهما البحثيتين في مجلة تكنومتركس بعنوان "انحدارات ريدج: التقدير المتحيز للمسائل غير المتعامدة" و"انحدارات ريدج: تطبيقات في المسائل غير المتعامدة". [ 5 ] [ 6 ] [ 1 ]

طُوِّرَ انحدار ريدج كحلٍّ مُحتملٍ لعدم دقة مُقدِّرات المربعات الصغرى عندما تحتوي نماذج الانحدار الخطي على بعض المتغيرات المستقلة متعددة الارتباط (شديدة الارتباط) - وذلك عن طريق إنشاء مُقدِّر انحدار ريدج (RR). يُوفِّر هذا تقديرًا أكثر دقةً لمعاملات ريدج، حيث يكون تباينه ومُقدِّر متوسط ​​مربعاته غالبًا أصغر من مُقدِّرات المربعات الصغرى المُشتقة سابقًا. [ 7 ] [ 2 ]

ملخص

في حل المربعات الصغرى العادية لـ

Y=Xβ+ε،{\displaystyle \mathbf {Y} =\mathbf {X} {\boldsymbol {\beta }}+{\boldsymbol {\varepsilon }},\,}

مشكلة مصفوفة العزم شبه المفردةXتيX{\displaystyle \mathbf {X} ^{\mathsf {T}}\mathbf {X} }يتم تخفيف هذه المشكلة بإضافة عناصر موجبة إلى الأقطار ، مما يقلل من رقم شرطها . بالمقارنة مع مُقدِّر المربعات الصغرى العادي، يحتوي مُقدِّر ريدج البسيط على حد إضافي.λأنا{\displaystyle \lambda \mathbf {I} }في المقام: β^λ=(XتيX+λأنا)-1XتيY{\displaystyle {\hat {\boldsymbol {\beta}}}_{\lambda }=\left(\mathbf {X} ^{\mathsf {T}}\mathbf {X} +\lambda \mathbf {I} \right)^{-1}\mathbf {X} ^{\mathsf {T}}\mathbf {Y} } أينY{\displaystyle \mathbf {Y} }هو المتغير التابع أو متجه الاستجابة،X{\displaystyle \mathbf {X} }هي مصفوفة التصميم ،أنا{\displaystyle \mathbf {I} }هي مصفوفة الوحدة ، ومعامل تنظيم ريدج (أو تيكهونوف)λ0{\displaystyle \lambda \geq 0}يُستخدم كثابت لإزاحة أقطار مصفوفة العزوم. [ 8 ] يمكن إثبات أن هذا المُقدِّر هو حل لمسألة المربعات الصغرى مع مراعاة القيد.βتيβ=ج{\displaystyle {\boldsymbol {\beta }}^{\mathsf {T}}{\boldsymbol {\beta }}=c}، والتي يمكن التعبير عنها كعملية تصغير لاغرانجية: أرجينينβY-Xβ2+λ(βتيβ-ج){\displaystyle {\text{argmin}}_{\boldsymbol {\beta }}\,\|\mathbf {Y} -\mathbf {X} {\boldsymbol {\beta }}\|^{2}+\lambda \left({\boldsymbol {\beta }}^{\mathsf {T}}{\boldsymbol {\beta }}-c\right)} مما يدل على أنλ{\displaystyle \lambda }ليس سوى مُضاعِف لاغرانج للقيد. [ 9 ] في الواقع، توجد علاقة تناظرية بينج{\displaystyle c}وλ{\displaystyle \lambda }وبما أننا، عملياً، لا نعرفج{\displaystyle c}، نحن نحددλ{\displaystyle \lambda }أو إيجادها من خلال استراتيجيات إضافية لملاءمة البيانات، انظر تحديد معلمة تيخونوف أدناه.

لاحظ ذلك كماλ0{\displaystyle \lambda \downarrow 0}يصبح القيد في النهاية غير ملزم ، ويتقارب مقدر ريدج إلى مقدر المربعات الصغرى العادية ذي المعيار الأدنى ، والذي يُشار إليه هنا بـ β^=β^0{\displaystyle {\hat {\boldsymbol {\beta }}}={\hat {\boldsymbol {\beta }}}_{0}}:

ليمλ0β^λ=X+Y=β^0،{\displaystyle \lim _{\lambda \downarrow 0}{\hat {\boldsymbol {\beta }}}_{\lambda }=\mathbf {X} ^{+}\mathbf {Y} ={\hat {\boldsymbol {\beta }}}_{0},} مع X+{\displaystyle \mathbf {X} ^{+}}يشير إلى المعكوس الزائف لـX{\displaystyle \mathbf {X} }.

تحديد معامل تيخونوف

معامل التنظيم الأمثلλ{\displaystyle \lambda }عادةً ما يكون غير معروف، وفي الواقع العملي يحتاج إلى تقدير. عادةً ما يتم اختيار معامل تنظيم تيخونوف بناءً على البيانات.λ{\displaystyle \lambda } يتم إنجاز ذلك إما عن طريق التحقق المتبادل ، أو عن طريق إجراء إضافي، كما يلي.

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

خيار شائع يعتمد على البيانات لـλ{\displaystyle \lambda }هو المُقلِّل لخسارة التحقق المتقاطع أو تعميماتها. على سبيل المثال، أثبتت غريس وهبة أن المعامل الأمثل، بمعنى التحقق المتقاطع المعمم، يُقلِّل [ 10 ] [ 11 ].

جي=RSSτ2=Xβ^-Y2[tr(أنا-X(XتيX+λ2أنا)-1Xتي)]2،{\displaystyle G={\frac {\operatorname {RSS} }{\tau ^{2}}}={\frac {\left\|\mathbf {X} {\hat {\boldsymbol {\beta }}}-\mathbf {Y} \right\|^{2}}{\left[\operatorname {tr} \left(\mathbf {I} -\mathbf {X} \left(\mathbf {X} ^{\mathsf {T}}\mathbf {X} +\lambda ^{2}\mathbf {I} \right)^{-1}\mathbf {X} ^{\mathsf {T}}\right)\right]^{2}}},} أينRSS{\displaystyle \operatorname {RSS} }هو مجموع مربعات البواقي ، وτ{\displaystyle \tau }يمثل العدد الفعال لدرجات الحرية .

مُقدِّر إضافي

افترض أنX{\displaystyle \mathbf {X} }هو ن×ص{\displaystyle n\times p}المصفوفة وتعريفهاΩ:=(XX/ن)+أوميغا :=(\mathbf {X} ^{\top }\mathbf {X} /n)^{+}} . ثم، ضع في اعتبارك الخيار التالي لمعامل تنظيم تيخونوف:

λ*:=ς2تر(Ω)βΩβ+3ς2تر(Ω2)/ن،{\displaystyle \lambda ^{*}:={\frac {\varsigma ^{2}\mathrm {tr} (\Omega )}{{\boldsymbol {\beta }}^{\top }\Omega {\boldsymbol {\beta }}+3\varsigma ^{2}\mathrm {tr} (\Omega ^{2})/n}},}

أينς2{\displaystyle \varsigma ^{2}}يمثل تباين الضوضاء ε=Y-Xβ{\displaystyle {\boldsymbol {\varepsilon}}=\mathbf {Y} -\mathbf {X} {\boldsymbol {\beta }}}، إنه، Vأر(ε)=ς2أنا{\displaystyle \mathrm {Var} ({\boldsymbol {\varepsilon }})=\varsigma ^{2}\mathbf {I} }يمكن إثبات [ 12 ] أن مُقدِّر الحافةβ^λ*{\displaystyle {\hat {\boldsymbol {\beta }}}_{\lambda ^{*}}}يتمتع بمخاطر متوقعة أقل داخل العينة مقارنةً بمُقدِّر المربعات الصغرى ذي المعيار الأدنىβ^0=X+Y{\displaystyle {\hat {\boldsymbol {\beta }}}_{0}=\mathbf {X} ^{+}\mathbf {Y} }وبشكل أدق،

هـY-Xβ^02هـY-Xβ^λ*2+ς2نλ*تر(Ω)،{\displaystyle \mathbb {E} \|\mathbf {Y} '-\mathbf {X} {\hat {\boldsymbol {\beta }}}_{0}\|^{2}\geq \mathbb {E} \|\mathbf {Y} '-\mathbf {X} {\hat {\boldsymbol {\beta }}}_{\lambda ^{*}}\|^{2}+{\frac {\varsigma ^{2}}{n}}\lambda ^{*}\mathrm {tr} (\Omega ),}

حيث تتعامل التوقعاتX{\displaystyle \mathbf {X} }كما هو ثابت وY{\displaystyle \mathbf {Y} '}بيانات استجابة الاختبار ، مستقلة عنY{\displaystyle \mathbf {Y} }(وبالتالي مستقل عن المقدرات)β^0{\displaystyle {\hat {\boldsymbol {\beta }}}_{0}}وβ^λ*{\displaystyle {\hat {\boldsymbol {\beta }}}_{\lambda ^{*}}}).

بالطبع، في الواقع العملي، تكون الصيغة لـλ*{\displaystyle \lambda ^{*}} يتم استخدامها عن طريق إدخال المقدرات الإحصائية للمعلمات المجهولة β{\displaystyle {\boldsymbol {\beta }}}وς2{\displaystyle \varsigma ^{2}}. متىن>ص{\displaystyle n>p}، وأكثر المقدرات الطبيعية لهذه المعلمات هي مقدرات المربعات الصغرى المعتادة:

β^=X+Y،ς^2=Y-Xβ^2ن-ص.{\displaystyle {\hat {\boldsymbol {\beta }}}=\mathbf {X} ^{+}\mathbf {Y} ,\qquad {\hat {\varsigma }}^{2}={\frac {\|\mathbf {Y} -\mathbf {X} {\hat {\boldsymbol {\beta }}}\|^{2}}{n-p}}.} استبدال المجهولβ،ς2{\displaystyle {\boldsymbol {\beta }},\varsigma ^{2}}في صيغةλ*{\displaystyle \lambda ^{*}}مع ما يقابلهβ^،ς^2{\displaystyle {\hat {\boldsymbol {\beta }}},{\hat {\varsigma }}^{2}}وهذا يعطي ما يسمى بالمُقدِّر الإضافيλ^*{\displaystyle {\widehat {\lambda }}^{*}}للحصول على الأمثلλ*{\displaystyle \lambda ^{*}}.

تشمل الأساليب البديلة لاختيار معلمة تنظيم تيخونوف القائمة على البيانات مبدأ التباين ، وطريقة منحنى L ، [ 13 ] الحد الأقصى المقيد للاحتمالية .

تاريخ

تم ابتكار تنظيم تيخونوف بشكل مستقل في سياقات متعددة. واكتسب شهرة واسعة من خلال تطبيقه على المعادلات التكاملية في أعمال أندريه تيخونوف [ 14 ] [ 15 ] [ 16 ] [ 17 ] [ 18 ] وديفيد ل. فيليبس [ 19 ] . ويستخدم بعض المؤلفين مصطلح تنظيم تيخونوف-فيليبس . وقد شرح آرثر إي. هورل الحالة ذات الأبعاد المحدودة ، معتمدًا على منهج إحصائي [ 20 ] ، وكذلك مانوس فوستر، الذي فسر هذه الطريقة على أنها مرشح وينر-كولموغوروف (كريغينغ) [ 21 ] . وتبعًا لهورل، يُعرف في الأدبيات الإحصائية باسم انحدار ريدج [ 22 ] ، نسبةً إلى تحليل ريدج ("ريدج" تشير إلى المسار من القيمة القصوى المقيدة) [ 23 ] .

تنظيم تيكهونوف للمعادلات الخطية

لنفترض أنه بالنسبة لمصفوفة حقيقية معروفةأ{\displaystyle A}والمتجهب{\displaystyle \mathbf {b} }نريد إيجاد متجهx{\displaystyle \mathbf {x} }بحيث أx=ب،{\displaystyle A\mathbf {x} =\mathbf {b} ,} أينx{\displaystyle \mathbf {x} }وب{\displaystyle \mathbf {b} }قد تكون بأحجام مختلفة وأ{\displaystyle A}بل قد يكون غير مربع.

الأسلوب القياسي هو الانحدار الخطي باستخدام طريقة المربعات الصغرى العادية . ومع ذلك، إذا لم يكن هناكx{\displaystyle \mathbf {x} }يحقق المعادلة أو أكثر من معادلة واحدةx{\displaystyle \mathbf {x} }إذا لم يكن الحل فريدًا، يُقال إن المسألة غير محددة جيدًا . في مثل هذه الحالات، يؤدي تقدير المربعات الصغرى العادية إلى نظام معادلات زائد التحديد ، أو في أغلب الأحيان نظام معادلات ناقص التحديد . تؤثر معظم الظواهر الواقعية تأثير مرشحات التمرير المنخفض في الاتجاه الأمامي حيثأ{\displaystyle A}خرائطx{\displaystyle \mathbf {x} }لب{\displaystyle \mathbf {b} }لذلك، عند حل المسألة العكسية، تعمل عملية التحويل العكسي كمرشح تمرير عالي ، مما يؤدي إلى تضخيم الضوضاء ( تكون القيم الذاتية /القيم المفردة أكبر ما يمكن في التحويل العكسي، بينما كانت أصغر ما يمكن في التحويل الأمامي). بالإضافة إلى ذلك، فإن طريقة المربعات الصغرى العادية تلغي ضمنيًا كل عنصر من عناصر النسخة المُعاد بناؤها.x{\displaystyle \mathbf {x} }ذلك في الفضاء الصفري لـأ{\displaystyle A}بدلاً من السماح باستخدام نموذج كمعلومات مسبقة لـx{\displaystyle \mathbf {x} }تسعى طريقة المربعات الصغرى العادية إلى تقليل مجموع مربعات البواقي ، والتي يمكن كتابتها باختصار على النحو التالي: أx-ب22،{\displaystyle \left\|A\mathbf {x} -\mathbf {b} \right\|_{2}^{2},} أين2{\displaystyle \|\cdot \|_{2}}هو المعيار الإقليدي .

لإعطاء الأفضلية لحل معين ذي خصائص مرغوبة، يمكن تضمين مصطلح تنظيمي في عملية التصغير هذه: أx-ب22+Γx22=أx-ب22،{\displaystyle \left\|A\mathbf {x} -\mathbf {b} \right\|_{2}^{2}+\left\|\Gamma \mathbf {x} \right\|_{2}^{2}=\left\|{\mathcal {A}}\mathbf {x} -{\mathcal {b}}\right\|_{2}^{2},} أينأ=(أΓ){\displaystyle {\mathcal {A}}={\begin{pmatrix}A\\\Gamma \end{pmatrix}}}وب=(ب0){\displaystyle {\mathcal {b}}={\begin{pmatrix}\mathbf {b} \\{\boldsymbol {0}}\end{pmatrix}}}، بالنسبة لمصفوفة تيخونوف المختارة بشكل مناسبΓ{\displaystyle \Gamma }في كثير من الحالات، يتم اختيار هذه المصفوفة كمضاعف قياسي لمصفوفة الوحدة (Γ=αأنا{\displaystyle \Gamma =\alpha I}يُفضّل هذا الأسلوب الحلول ذات المعايير الأصغر ، ويُعرف باسم تنظيم L2 . [ 24 ] في حالات أخرى، يمكن استخدام مُعاملات التمرير العالي (مثل مُعامل الفرق أو مُعامل فورييه الموزون ) لفرض السلاسة إذا كان يُعتقد أن المتجه الأساسي متصل في الغالب. يُحسّن هذا التنظيم حالة المسألة، مما يُتيح حلًا عدديًا مباشرًا. ويمكن التعامل معها كمسألة مربعات صغرى عادية مع مصفوفات مُوسّعة.أ{\displaystyle {\mathcal {A}}}وب{\displaystyle {\mathcal {b}}}الحل هو x^=(أتيأ)-1أتيب=(أتيأ+ΓتيΓ)-1أتيب.{\displaystyle {\hat {\mathbf {x} }}=({\mathcal {A}}^{\mathsf {T}}{\mathcal {A}})^{-1}{\mathcal {A}}^{\mathsf {T}}\mathbf {\mathcal {b}} =(A^{\mathsf {T}}A+\Gamma ^{\mathsf {T}}\Gamma )^{-1}A^{\mathsf {T}}\mathbf {b} .} قد يختلف تأثير التنظيم باختلاف حجم المصفوفةΓ{\displaystyle \Gamma }. لΓ=0{\displaystyle \Gamma =0}يُختزل هذا إلى حل المربعات الصغرى غير المنتظم، بشرط وجود ( ATA ) −1 . لاحظ أنه في حالة المصفوفة المركبةأ{\displaystyle A}وكالعادة، يتم نقلهاأتي{\displaystyle A^{\mathsf {T}}}يجب استبدالها بالتحويل الهرميتيأح{\displaystyle A^{\mathsf {H}}}.

يُستخدم تنظيم L2 في العديد من السياقات إلى جانب الانحدار الخطي، مثل التصنيف باستخدام الانحدار اللوجستي أو آلات المتجهات الداعمة ، [ 25 ] وتحليل المصفوفات. [ 26 ]

تطبيق على نتائج الملاءمة الحالية

بما أن تنظيم تيكهونوف يضيف ببساطة حدًا تربيعيًا إلى دالة الهدف في مسائل التحسين، فمن الممكن القيام بذلك بعد إجراء التحسين غير المنظم. على سبيل المثال، إذا كانت المسألة المذكورة أعلاه معΓ=0{\displaystyle \Gamma =0}ينتج الحلx^0{\displaystyle {\hat {\mathbf {x} }}_{0}}، الحل في وجودΓ0{\displaystyle \Gamma \neq 0}يمكن التعبير عنها على النحو التالي: x^=بx^0،{\displaystyle {\hat {\mathbf {x} }}=B{\hat {\mathbf {x} }}_{0},} باستخدام "مصفوفة التنظيم"ب=(أتيأ+ΓتيΓ)-1أتيأ{\displaystyle B=\left(A^{\mathsf {T}}A+\Gamma ^{\mathsf {T}}\Gamma \right)^{-1}A^{\mathsf {T}}A}.

إذا كانت عملية ملاءمة المعلمات تأتي مع مصفوفة تباين لعدم اليقين في المعلمات المقدرةV0{\displaystyle V_{0}}إذن، ستكون مصفوفة التنظيم هي ب=(V0-1+ΓتيΓ)-1V0-1،{\displaystyle B=(V_{0}^{-1}+\Gamma ^{\mathsf {T}}\Gamma )^{-1}V_{0}^{-1},} وستكون للنتيجة المنتظمة تباين مشترك جديد V=بV0بتي.{\displaystyle V=BV_{0}B^{\mathsf {T}}.}

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

تنظيم تيخونوف المعمم

بالنسبة للتوزيعات الطبيعية متعددة المتغيرات العامة لـx{\displaystyle \mathbf {x} }وباستخدام خطأ البيانات، يمكن تطبيق تحويل على المتغيرات لتقليلها إلى الحالة المذكورة أعلاه. وبالمثل، يمكن البحث عنx{\displaystyle \mathbf {x} }لتقليل أx-بP2+x-x0سؤال2،{\displaystyle \left\|A\mathbf {x} -\mathbf {b} \right\|_{P}^{2}+\left\|\mathbf {x} -\mathbf {x} _{0}\right\|_{Q}^{2},} حيث استخدمناxسؤال2{\displaystyle \left\|\mathbf {x} \right\|_{Q}^{2}}لترمز إلى مربع المعيار المرجحxتيسؤالx{\displaystyle \mathbf {x} ^{\mathsf {T}}Q\mathbf {x} }(قارن بمسافة ماهالانوبيس ). في التفسير البايزيP{\displaystyle P}هي مصفوفة التغاير العكسية لـب{\displaystyle \mathbf {b} }،x0{\displaystyle \mathbf {x} _{0}}القيمة المتوقعة لـx{\displaystyle \mathbf {x} }، وسؤال{\displaystyle Q}هي مصفوفة التغاير العكسية لـx{\displaystyle \mathbf {x} }.

لا يتم تضمين مصفوفة تيخونوف بشكل صريح لأن مصطلح التنظيم المقابل Γx-x0سؤال2{\displaystyle \left\|\Gamma \mathbf {x} -\mathbf {x} _{0}'\right\|_{Q'}^{2}}يُختزل إلى ما سبق معΓx0=x0{\displaystyle \Gamma \mathbf {x} _{0}=\mathbf {x} _{0}'}وسؤال=ΓتيسؤالΓ{\displaystyle Q=\Gamma ^{T}Q'\Gamma }بالنسبة للتنظيم الطبيعي حيثسؤال=أنا{\displaystyle Q'=I}ثم تظهر مصفوفة تيخونوف في تحليل تشوليسكيسؤال=ΓتيΓ{\displaystyle Q=\Gamma ^{\mathsf {T}}\Gamma }ويعتبر مرشحاً لتبييض البشرة .

لهذه المسألة المعممة حل أمثلx^{\displaystyle {\hat {\mathbf {x} }}}والتي يمكن كتابتها صراحة باستخدام الصيغة x^=(أتيPأ+سؤال)-1(أتيPب+سؤالx0)=x0+(أتيPأ+سؤال)-1(أتيP(ب-أx0)).{\displaystyle \mathbf {\hat {\mathbf {x} }} =\left(A^{\mathsf {T}}PA+Q\right)^{-1}\left(A^{\mathsf {T}}P\mathbf {b} +Q\mathbf {x} _{0}\right)=\mathbf {x} _{0}+\left(A^{\mathsf {T}}PA+Q\right)^{-1}\left(A^{\mathsf {T}}P\left(\mathbf {b} -A\mathbf {x} _{0}\right)\right).}

تنظيم لافرينتييف

في بعض الحالات، يمكن تجنب استخدام النقلأتي{\displaystyle A^{\mathsf {T}}}كما اقترح ميخائيل لافرينتييف . [ 28 ] على سبيل المثال، إذاأ{\displaystyle A}متماثلة موجبة محددة، أيأ=أتي>0{\displaystyle A=A^{\mathsf {T}}>0}وكذلك معكوسهأ-1{\displaystyle A^{-1}}وبالتالي، يمكن استخدام ذلك لإنشاء معيار التربيع المرجح.xP2=xتيأ-1x{\displaystyle \left\|\mathbf {x} \right\|_{P}^{2}=\mathbf {x} ^{\mathsf {T}}A^{-1}\mathbf {x} }في تنظيم تيخونوف المعمم، مما يؤدي إلى تقليل أx-بأ-12+x-x0سؤال2{\displaystyle \left\|A\mathbf {x} -\mathbf {b} \right\|_{A^{-1}}^{2}+\left\|\mathbf {x} -\mathbf {x} _{0}\right\|_{Q}^{2}} أو، بشكل مكافئ حتى حد ثابت ، xتي(أ+سؤال)x-2xتي(ب+سؤالx0).{\displaystyle \mathbf {x} ^{\mathsf {T}}\left(A+Q\right)\mathbf {x} -2\mathbf {x} ^{\mathsf {T}}\left(\mathbf {b} +Q\mathbf {x} _{0}\right).}

لهذه المسألة التصغيرية حل أمثلx*{\displaystyle \mathbf {x} ^{*}}والتي يمكن كتابتها صراحة باستخدام الصيغة x*=(أ+سؤال)-1(ب+سؤالx0)،{\displaystyle \mathbf {x} ^{*}=\left(A+Q\right)^{-1}\left(\mathbf {b} +Q\mathbf {x} _{0}\right),} وهو ليس سوى حل مسألة تيخونوف المعممة حيثأ=أتي=P-1.{\displaystyle A=A^{\mathsf {T}}=P^{-1}.}

يُعدّ تنظيم لافرينتييف، إن أمكن تطبيقه، أفضل من تنظيم تيخونوف الأصلي، لأن مصفوفة لافرينتييفأ+سؤال{\displaystyle A+Q}يمكن تحسين حالتها، أي أن يكون لها رقم حالة أصغر ، مقارنةً بمصفوفة تيخونوف.أتيأ+ΓتيΓ.{\displaystyle A^{\mathsf {T}}A+\Gamma ^{\mathsf {T}}\Gamma .}

التنظيم في فضاء هيلبرت

عادةً ما تنتج مسائل المعادلات التكاملية الخطية المنفصلة سيئة التكييف عن تجزئة المعادلات التكاملية ، ويمكن صياغة تنظيم تيخونوف في السياق الأصلي ذي الأبعاد اللانهائية. في ما سبق، يمكننا تفسيرأ{\displaystyle A}كعامل مضغوط على فضاءات هيلبرت ، وx{\displaystyle x}وب{\displaystyle b}كعناصر في نطاق ومجالأ{\displaystyle A}المشغلأ*أ+ΓتيΓ{\displaystyle A^{*}A+\Gamma ^{\mathsf {T}}\Gamma }إذن، يكون مؤثرًا ذاتيًا مترافقًا ومحدودًا وقابلًا للعكس.

العلاقة بتحليل القيم المفردة ومرشح وينر

معΓ=αأنا{\displaystyle \Gamma =\alpha I}يمكن تحليل حل المربعات الصغرى هذا بطريقة خاصة باستخدام تحليل القيم المفردة . وبالنظر إلى تحليل القيم المفردة أ=يوΣVتي{\displaystyle A=U\Sigma V^{\mathsf {T}}} بقيم مفردةσأنا{\displaystyle \sigma _{i}}يمكن التعبير عن حل تيخونوف المنتظم على النحو التالي: x^=Vديوتيب،{\displaystyle {\hat {x}}=VDU^{\mathsf {T}}b,} أيند{\displaystyle D}يحتوي على قيم قطرية دأناأنا=σأناσأنا2+α2{\displaystyle D_{ii}={\frac {\sigma _{i}}{\sigma _{i}^{2}+\alpha ^{2}}}} وتكون قيمتها صفرًا في أي مكان آخر. وهذا يُبيّن تأثير مُعامل تيخونوف على رقم شرط المسألة المُنتظمة. أما في الحالة المُعمّمة، فيُمكن اشتقاق تمثيل مُشابه باستخدام تحليل القيم المفردة المُعمّم . [ 29 ]

وأخيراً، يتعلق الأمر بمرشح وينر : x^=أنا=1qوأناuأناتيبσأناvأنا،{\displaystyle {\hat {x}}=\sum _{i=1}^{q}f_{i}{\frac {u_{i}^{\mathsf {T}}b}{\sigma _{i}}}v_{i},} حيث توجد أوزان وينروأنا=σأنا2σأنا2+α2{\displaystyle f_{i}={\frac {\sigma _{i}^{2}}{\sigma _{i}^{2}+\alpha ^{2}}}}وq{\displaystyle q}هي رتبةأ{\displaystyle A}.

العلاقة بالصياغة الاحتمالية

يُدخل الصياغة الاحتمالية للمسألة العكسية (عندما تكون جميع حالات عدم اليقين غاوسية) مصفوفة تباين مشتركةجم{\displaystyle C_{M}}تمثل هذه المصفوفة حالات عدم اليقين المسبقة بشأن معلمات النموذج، ومصفوفة التغاير.جد{\displaystyle C_{D}}تمثل هذه القيم عدم اليقين في المعلمات المرصودة. [ 30 ] في الحالة الخاصة عندما تكون هاتان المصفوفتان قطريتين ومتناحيتين،جم=σم2أنا{\displaystyle C_{M}=\sigma _{M}^{2}I}وجد=σد2أنا{\displaystyle C_{D}=\sigma _{D}^{2}I}وفي هذه الحالة، تختزل معادلات النظرية العكسية إلى المعادلات المذكورة أعلاه، معα=σد/σم{\displaystyle \alpha ={\sigma _{D}}/{\sigma _{M}}}[ 31 ] [ 32 ]

التفسير البايزي

على الرغم من أن اختيار الحل لهذه المشكلة المنتظمة قد يبدو مصطنعًا في البداية، وبالفعل المصفوفةΓ{\displaystyle \Gamma }يبدو الأمر تعسفيًا إلى حد ما، لكن يمكن تبرير هذه العملية من وجهة نظر بايزية . [ 33 ] تجدر الإشارة إلى أنه بالنسبة للمسألة غير المحددة جيدًا، يجب بالضرورة إدخال بعض الافتراضات الإضافية للحصول على حل فريد. إحصائيًا، فإن التوزيع الاحتمالي المسبق لـx{\displaystyle x}يُعتبر أحيانًا توزيعًا طبيعيًا متعدد المتغيرات . [ 34 ] ولتبسيط الأمر، تُفترض الافتراضات التالية: المتوسطات تساوي صفرًا؛ مكوناتها مستقلة؛ المكونات لها نفس الانحراف المعياريσx{\displaystyle \sigma _{x}}كما أن البيانات عرضة للأخطاء، والأخطاء فيب{\displaystyle b}يُفترض أيضًا أن تكون مستقلة بمتوسط ​​صفر وانحراف معياريσب{\displaystyle \sigma _{b}}في ظل هذه الافتراضات، يُعدّ الحل المُنتظم وفقًا لنموذج تيكهونوف هو الحل الأكثر احتمالًا بالنظر إلى البيانات والتوزيع المسبق لـx{\displaystyle x}، وفقًا لنظرية بايز . [ 35 ]

إذا تم استبدال افتراض التوزيع الطبيعي بافتراضات تجانس التباين وعدم ارتباط الأخطاء ، وإذا استمر افتراض المتوسط ​​الصفري، فإن نظرية جاوس-ماركوف تستلزم أن الحل هو المقدر الخطي الأدنى غير المتحيز . [ 36 ]

انظر أيضاً

ملحوظات

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

مراجع

  1. 1 2 هيلت، دونالد إي.؛ سيغريست، دونالد دبليو. (1977). ريدج، برنامج حاسوبي لحساب تقديرات انحدار ريدج . doi : 10.5962/bhl.title.68934 .
  2. 1 2 غروبر، مارفن (1998). تحسين الكفاءة عن طريق التقليص: مقدرات انحدار جيمس-شتاين وريدج . مطبعة سي آر سي. ص 2. ISBN  978-0-8247-0156-7.
  3. كينيدي، بيتر (2003). دليل الاقتصاد القياسي ( الطبعة الخامسة). كامبريدج: مطبعة معهد ماساتشوستس للتكنولوجيا. الصفحات 205-206 . ISBN   0-262-61183-X.
  4. غروبر، مارفن (1998). تحسين الكفاءة عن طريق التقليص: مُقدِّرات جيمس-شتاين وانحدار ريدج . بوكا راتون: مطبعة سي آر سي. الصفحات 7-15 . ISBN  0-8247-0156-9.
  5. هورل، آرثر إي.؛ كينارد، روبرت دبليو. (1970). "انحدار ريدج: تقدير متحيز للمسائل غير المتعامدة". تكنومتركس . 12 (1): 55-67 . doi : 10.2307/1267351 . JSTOR 1267351 . 
  6. هورل، آرثر إي.؛ كينارد، روبرت دبليو. (1970). "انحدار ريدج: تطبيقات على المسائل غير المتعامدة". تكنومتركس . 12 (1): 69-82 . doi : 10.2307/1267352 . JSTOR 1267352 . 
  7. جوليف، آي تي ​​(2006). تحليل المكونات الرئيسية . سبرينغر ساينس آند بيزنس ميديا. ص 178. ISBN  978-0-387-22440-4.
  8. لاختيارλ{\displaystyle \lambda }للاطلاع على التفاصيل العملية، انظر: خلف، غضبان؛ شكور، غازي (2005). "اختيار معامل ريدج لمسائل الانحدار". مجلة الاتصالات في الإحصاء - النظرية والأساليب . 34 (5): 1177-1182 . doi : 10.1081/STA-200056836 . S2CID 122983724 . 
  9. ^ فان فيرينجن ، فيسيل (2021-05-31). “ملاحظات محاضرة عن انحدار التلال”. أرخايف : 1509.09169 [ stat.ME ].
  10. وهبة، ج. (1990). "نماذج الدوال التكعيبية للبيانات الرصدية". سلسلة مؤتمرات CBMS-NSF الإقليمية في الرياضيات التطبيقية . جمعية الرياضيات الصناعية والتطبيقية. Bibcode : 1990smod.conf.....W .
  11. غولوب، ج.؛ هيث، م.؛ وهبة، ج. (1979). "التحقق المتقاطع المعمم كطريقة لاختيار معلمة ريدج جيدة" (ملف PDF) . تكنومتركس . 21 (2): 215-223 . doi : 10.1080/00401706.1979.10489751 .
  12. بوتيف، زدرافكو إي.؛ كروس، ديرك ب.؛ تايمر، توماس (2025). علم البيانات والتعلم الآلي: الأساليب الرياضية والإحصائية ( الطبعة الثانية). بوكا راتون ؛ لندن: مطبعة سي آر سي. ص 267-268. ISBN    978-1-032-48868-4.
  13. بي سي هانسن، "منحنى L واستخدامه في المعالجة العددية للمسائل العكسية"،
  14. ^ تيخونوف، أندريه نيكولايفيتش (1943). "حول الثبات في المشكلات العكسية" [ حول استقرار المشكلات العكسية ] . دوكلادي أكاديمي ناوك SSSR . 39 (5): 195– 198. مؤرشفة من الأصلي بتاريخ 27-02-2005.
  15. ^ تيخونوف، أن (1963). "حول القرارات غير الصحيحة بشأن فترة ما بعد الولادة وطريقة التنظيم". دوكلادي أكاديمي ناوك SSSR . 151 : 501 – 504.. مترجم في "حل المسائل المصاغة بشكل غير صحيح وطريقة التنظيم". الرياضيات السوفيتية . 4 : 1035-1038 .
  16. تيخونوف، أ.ن.؛ في.ي. أرسينين (1977). حل المسائل غير المحددة جيدًا . واشنطن: وينستون وأولاده. ISBN 0-470-99124-0.
  17. تيخونوف، أندريه نيكولايفيتش؛ غونشارسكي، أ.؛ ستيبانوف، ف. ف.؛ ياغولا، أناتولي غريغوريفيتش (30 يونيو 1995). الطرق العددية لحل المسائل غير المحددة جيدًا . هولندا: سبرينغر هولندا. ISBN 0-7923-3583-Xتم الاطلاع عليه بتاريخ 9 أغسطس 2018 .
  18. ^ تيخونوف، أندريه نيكولاييفيتش؛ ليونوف، ألكسندر س. ياجولا، أناتوليج جريجوريفيتش (1998). المشكلات غير الخطية المطروحة . لندن: تشابمان وهال. رقم ISBN 0-412-78660-5تم الاطلاع عليه بتاريخ 9 أغسطس 2018 .
  19. فيليبس، د. ل. (1962). "تقنية للحل العددي لبعض المعادلات التكاملية من النوع الأول" . مجلة ACM . 9 : 84-97 . doi : 10.1145/321105.321114 . S2CID 35368397 . 
  20. هورل، آرثر إي. (1962). "تطبيق تحليل ريدج على مسائل الانحدار". التقدم في الهندسة الكيميائية . 58 (3): 54-59 .
  21. فوستر، م. (1961). "تطبيق نظرية وينر-كولموغوروف للتنعيم على معكوس المصفوفة". مجلة جمعية الرياضيات الصناعية والتطبيقية . 9 (3): 387-392 . doi : 10.1137/0109031 .
  22. هورل، أ. إي.؛ كينارد، ر. و. (1970). "انحدار ريدج: تقدير متحيز للمسائل غير المتعامدة". تكنومتركس . 12 (1): 55-67 . doi : 10.1080/00401706.1970.10488634 .
  23. هورل، روجر و. (2020-10-01). "انحدار ريدج: سياق تاريخي" . تكنومتركس . 62 (4): 420-425 . doi : 10.1080/00401706.2020.1742207 . ISSN 0040-1706 . 
  24. نغ، أندرو واي. (2004). اختيار الميزات، والتنظيم L1 مقابل L2، والثبات الدوراني (PDF) . وقائع المؤتمر الدولي للتعلم الآلي .
  25. ر. إ. فان؛ ك. و. تشانغ؛ س. ج. هسيه؛ إكس. ر. وانغ؛ س. ج. لين (2008). "LIBLINEAR: مكتبة للتصنيف الخطي واسع النطاق". مجلة أبحاث تعلم الآلة . 9 : 1871-1874 .
  26. غوان، نايانغ؛ تاو، داتشنغ؛ لو، تشيغانغ؛ يوان، بو (2012). "تحليل المصفوفات غير السالبة عبر الإنترنت باستخدام تقريب عشوائي قوي". معاملات IEEE في الشبكات العصبية وأنظمة التعلم . 23 (7): 1087-1099 . Bibcode : 2012ITNNL..23.1087G . doi : 10.1109/TNNLS.2012.2197827 . PMID 24807135. S2CID 8755408 .  
  27. كوتش، لوكاس (2022). "التنظيم اللاحق لقياسات المقطع العرضي غير المطوي". مجلة الأجهزة . 17 (10) 10021. arXiv : 2207.02125 . Bibcode : 2022JInst..17P0021K . doi : 10.1088/1748-0221/17/10/P10021 .
  28. لافرينتييف، م.م. (1967). بعض المسائل المطروحة بشكل غير صحيح في الفيزياء الرياضية . نيويورك: سبرينغر.
  29. هانسن، بير كريستيان (1 يناير 1998). مسائل الرتبة الناقصة والمسائل غير المحددة جيدًا المنفصلة: الجوانب العددية للانعكاس الخطي ( الطبعة الأولى). فيلادلفيا، الولايات المتحدة الأمريكية: SIAM. ISBN  978-0-89871-403-6.
  30. تارانتولا، ألبرت (2005). نظرية المسائل العكسية وطرق تقدير معلمات النموذج ( الطبعة الأولى). فيلادلفيا: جمعية الرياضيات الصناعية والتطبيقية (SIAM). ISBN  0-89871-792-2تم الاطلاع عليه بتاريخ 9 أغسطس 2018 .
  31. هوانغ، يونفي؛ وآخرون (2019). "مجهر قوة الجر مع تنظيم مُحسَّن واختيار تلقائي لمعاملات بايز لمقارنة الخلايا" . التقارير العلمية . 9 (1) 539: 537. arXiv : 1810.05848 . Bibcode : 2019NatSR...9..539H . doi : 10.1038/ s41598-018-36896 -x . PMC 6345967. PMID 30679578 .   
  32. هوانغ، يونفي؛ غومبر، جيرهارد؛ ساباس، بينيديكت (2020). "طريقة مجهرية لقوة الجر باستخدام نظرية بايز مع إزالة التشويش التلقائية في حزمة برمجية سهلة الاستخدام". مجلة اتصالات الفيزياء الحاسوبية . 256 107313. arXiv : 2005.01377 . Bibcode : 2020CoPhC.25607313H . doi : 10.1016/j.cpc.2020.107313 .
  33. غرينبيرغ، إدوارد؛ ويبستر، تشارلز إي. الابن (1983). الاقتصاد القياسي المتقدم: جسر إلى الأدبيات . نيويورك: جون وايلي وأولاده. ص 207-213 . ISBN  0-471-09077-8.
  34. هوانغ، يونفي؛ وآخرون (2019). "مجهر قوة الجر مع تنظيم مُحسَّن واختيار تلقائي لمعاملات بايز لمقارنة الخلايا" . التقارير العلمية . 9 (1) 539: 537. arXiv : 1810.05848 . Bibcode : 2019NatSR...9..539H . doi : 10.1038/ s41598-018-36896 -x . PMC 6345967. PMID 30679578 .   
  35. فوغل، كورتيس ر. (2002). الأساليب الحسابية للمسائل العكسية . فيلادلفيا: جمعية الرياضيات الصناعية والتطبيقية. ISBN 0-89871-550-4.
  36. أميميا، تاكيشي (1985). الاقتصاد القياسي المتقدم . مطبعة جامعة هارفارد. ص 60-61 . ISBN  0-674-00560-0.

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