الانحدار اللوجستي

رسم بياني توضيحي لمنحنى الانحدار اللوجستي المُطابق للبيانات. يُظهر المنحنى الاحتمالية المُقدَّرة لاجتياز الامتحان (متغير تابع ثنائي) مقابل ساعات الدراسة (متغير مستقل عددي). انظر قسم "  مثال" للاطلاع على تفاصيل الحل.

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

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

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

تُقدَّر معلمات الانحدار اللوجستي عادةً باستخدام طريقة تقدير الاحتمال الأقصى (MLE). لا يوجد لهذه الطريقة صيغة مغلقة، على عكس طريقة المربعات الصغرى الخطية ؛ انظر قسم "  ملاءمة النموذج" . يؤدي الانحدار اللوجستي باستخدام طريقة تقدير الاحتمال الأقصى دورًا أساسيًا مشابهًا للاستجابات الثنائية أو الفئوية، كما يؤديه الانحدار الخطي باستخدام طريقة المربعات الصغرى العادية (OLS) للاستجابات العددية : فهو نموذج أساسي بسيط ومُحلَّل جيدًا؛ انظر قسم "  المقارنة مع الانحدار الخطي" للمزيد من التفاصيل. طُوِّر الانحدار اللوجستي كنموذج إحصائي عام، وانتشر استخدامه بشكل أساسي على يد جوزيف بيركسون ، [ 5 ] بدءًا من كتابه (1944) ، حيث صاغ مصطلح "logit"؛ انظر قسم "  التاريخ" .

التطبيقات

عام

يُستخدم الانحدار اللوجستي في مجالات متنوعة، تشمل التعلم الآلي، ومعظم المجالات الطبية، والعلوم الاجتماعية. فعلى سبيل المثال، طُوِّر مقياس شدة الصدمة والإصابة ( TRISS )، الذي يُستخدم على نطاق واسع للتنبؤ بمعدل الوفيات لدى المرضى المصابين، في الأصل بواسطة بويد وآخرون باستخدام الانحدار اللوجستي. [ 6 ] كما طُوِّرت العديد من المقاييس الطبية الأخرى المستخدمة لتقييم شدة حالة المريض باستخدام الانحدار اللوجستي. [ 7 ] [ 8 ] [ 9 ] [ 10 ] ويمكن استخدام الانحدار اللوجستي للتنبؤ بخطر الإصابة بمرض معين (مثل داء السكري ، أو مرض الشريان التاجي )، بناءً على الخصائص الملحوظة للمريض (العمر، والجنس، ومؤشر كتلة الجسم ، ونتائج فحوصات الدم المختلفة ، وما إلى ذلك). [ 11 ] [ 12 ] مثال آخر هو التنبؤ بما إذا كان الناخب النيبالي سيصوت لحزب المؤتمر النيبالي أو الحزب الشيوعي النيبالي أو لأي حزب آخر، بناءً على العمر، والدخل، والجنس، والعرق، ومكان الإقامة، والتصويت في الانتخابات السابقة، وما إلى ذلك. [ 13 ] يمكن استخدام هذه التقنية أيضًا في الهندسة ، وخاصة للتنبؤ باحتمالية فشل عملية أو نظام أو منتج معين. [ 14 ] [ 15 ] كما تُستخدم في تطبيقات التسويق ، مثل التنبؤ بميل العميل لشراء منتج أو إلغاء اشتراك، وما إلى ذلك. [ 16 ] في الاقتصاد ، يمكن استخدامها للتنبؤ باحتمالية انضمام شخص ما إلى سوق العمل، ومن تطبيقاتها التجارية التنبؤ باحتمالية تخلف مالك منزل عن سداد قرض عقاري . تُستخدم الحقول العشوائية الشرطية ، وهي امتداد للانحدار اللوجستي للبيانات المتسلسلة، في معالجة اللغة الطبيعية . يعتمد مخططو الكوارث والمهندسون على هذه النماذج للتنبؤ بالقرارات التي يتخذها أصحاب المنازل أو شاغلو المباني في عمليات الإخلاء الصغيرة والكبيرة، مثل حرائق المباني وحرائق الغابات والأعاصير وغيرها. [ 17 ] [ 18 ] [ 19 ] تساعد هذه النماذج في تطوير خطط موثوقة لإدارة الكوارث وتصميم بيئة مبنية أكثر أمانًا .

التعلم الآلي الخاضع للإشراف

الانحدار اللوجستي هو خوارزمية تعلم آلي خاضعة للإشراف، تُستخدم على نطاق واسع في مهام التصنيف الثنائي ، مثل تحديد ما إذا كانت رسالة البريد الإلكتروني بريدًا عشوائيًا أم لا، وتشخيص الأمراض من خلال تقييم وجود أو غياب حالات مرضية محددة بناءً على نتائج فحوصات المرضى. يعتمد هذا النهج على الدالة اللوجستية (أو الدالة السينية) لتحويل توليفة خطية من خصائص الإدخال إلى قيمة احتمالية تتراوح بين 0 و1. تشير هذه الاحتمالية إلى مدى ترجيح أن يتوافق مُدخل مُعين مع إحدى فئتين مُحددتين مسبقًا. تستند الآلية الأساسية للانحدار اللوجستي إلى قدرة الدالة اللوجستية على نمذجة احتمالية النتائج الثنائية بدقة. بفضل منحنى الدالة اللوجستية المميز على شكل حرف S، فإنها تُحول أي عدد حقيقي إلى قيمة ضمن النطاق من 0 إلى 1. هذه الميزة تجعلها مناسبة بشكل خاص لمهام التصنيف الثنائي، مثل فرز رسائل البريد الإلكتروني إلى "بريد عشوائي" أو "غير عشوائي". من خلال حساب احتمالية تصنيف المتغير التابع ضمن مجموعة محددة، يوفر الانحدار اللوجستي إطارًا احتماليًا يدعم اتخاذ القرارات المستنيرة. [ 20 ]

مثال

مشكلة

كمثال بسيط، يمكننا استخدام الانحدار اللوجستي مع متغير تفسيري واحد وفئتين للإجابة على السؤال التالي:

تقضي مجموعة من 20 طالباً ما بين صفر وست ساعات في الدراسة لامتحان. كيف يؤثر عدد ساعات الدراسة على احتمالية نجاح الطالب في الامتحان؟

يكمن سبب استخدام الانحدار اللوجستي في هذه المسألة في أن قيم المتغير التابع، النجاح والرسوب، على الرغم من تمثيلها بالرقمين "1" و"0"، ليست أعدادًا أصلية . لو تم تغيير المسألة بحيث استُبدل النجاح/الرسوب بالدرجات من 0 إلى 100 (أعداد أصلية)، لأمكن استخدام تحليل الانحدار البسيط.

يوضح الجدول عدد الساعات التي قضاها كل طالب في الدراسة، وما إذا كان قد نجح (1) أو رسب (0).

ساعات ( س ك )0.500.751.001.251.501.751.752.002.252.502.753.003.253.504.004.254.504.755.005.50
مرر ( y k )00000010101010111111

نرغب في تطبيق دالة لوجستية على البيانات التي تتكون من ساعات الدراسة ( xk ) ونتيجة الاختبار ( yk = 1  للنجاح، 0 للرسوب). يتم فهرسة نقاط البيانات بواسطة الرمز السفلي k الذي يبدأ منك=1{\displaystyle k=1}لك=ك=20{\displaystyle k=K=20}. يُطلق على المتغير x اسم " المتغير التفسيري "، ويُطلق على المتغير y اسم " المتغير الفئوي " الذي يتكون من فئتين: "ناجح" أو "راسب" بما يتوافق مع القيم الفئوية 1 و 0 على التوالي.

نموذج

رسم بياني لمنحنى الانحدار اللوجستي المُطابق لبيانات ( x m , y m ). يوضح المنحنى احتمالية اجتياز الامتحان مقابل ساعات الدراسة.

تأخذ الدالة اللوجستية الشكل التالي:

ص(x)=11+هـ-(x-μ)/s{\displaystyle p(x)={\frac {1}{1+e^{-(x-\mu )/s}}}}

حيث μ هو معامل الموقع (نقطة منتصف المنحنى، حيثص(μ)=1/2{\displaystyle p(\mu )=1/2}) و s هو معامل قياس . يمكن إعادة كتابة هذا التعبير على النحو التالي:

ص(x)=11+هـ-(β0+β1x){\displaystyle p(x)={\frac {1}{1+e^{-(\beta _{0}+\beta _{1}x)}}}}

أينβ0=-μ/s{\displaystyle \beta _{0}=-\mu /s}ويُعرف باسم نقطة التقاطع (وهي نقطة التقاطع الرأسية أو نقطة التقاطع مع المحور الصادي للخط).y=β0+β1x{\displaystyle y=\beta _{0}+\beta _{1}x})، وβ1=1/s{\displaystyle \beta _{1}=1/s}(معامل المقياس العكسي أو معامل المعدل ): وهما نقطة تقاطع المحور الصادي وميل منحنى اللوغاريتم النسبي كدالة لـ x . وبالعكس،μ=-β0/β1{\displaystyle \mu =-\beta _{0}/\beta _{1}}وs=1/β1{\displaystyle s=1/\beta _{1}}.

لاحظ أن هذا النموذج هو في الواقع تبسيط مفرط، لأنه يفترض أن كل طالب سينجح إذا درس إلى أجل غير مسمى (الحد = 1).

ملائم

يُستخدم عادةً مقياس جودة المطابقة للانحدار اللوجستي باستخدام دالة الخسارة اللوجستية (أو دالة الخسارة اللوغاريتمية )، وهي دالة الاحتمال اللوغاريتمي السالبة . بالنسبة لقيم معينة لـ x<sub> k</sub> و y<sub> k </sub>، اكتبصك=ص(xك){\displaystyle p_{k}=p(x_{k})}الـصك{\displaystyle p_{k}} هي احتمالات أنyك{\displaystyle y_{k}} ستساوي واحدًا، و1-صك{\displaystyle 1-p_{k}}تمثل احتمالات أن تكون هذه القيم صفرًا (انظر توزيع برنولي ). نريد إيجاد قيم هذه الاحتمالات .β0{\displaystyle \beta _{0}}وβ1{\displaystyle \beta _{1}}والتي تُعطي "أفضل تطابق" مع البيانات. وللمقارنة بين "أفضل تطابق" مع البيانات، انظر حالة الانحدار الخطي . هناك، يُعتبر مجموع مربعات انحرافات التطابق عن نقاط البيانات ( yk ) ، أي خسارة الخطأ التربيعي ، مقياسًا لجودة التطابق، ويتم الحصول على أفضل تطابق عندما تكون هذه الخسارة في أدنى حد .

الخسارة اللوغاريتمية للنقطة kك{\displaystyle \ell _{k}}هو :

ك={-lnصك لو yك=1،-ln(1-صك) لو yك=0.{\displaystyle \ell _{k}={\begin{cases}-\ln p_{k}&{\text{ if }}y_{k}=1,\\-\ln(1-p_{k})&{\text{ if }}y_{k}=0.\end{cases}}}

يمكن تفسير خسارة اللوغاريتم على أنها " مفاجأة " النتيجة الفعلية .yك{\displaystyle y_{k}}مقارنة بالتوقعصك{\displaystyle p_{k}}وهو مقياس لمحتوى المعلومات . يكون فقدان اللوغاريتم دائمًا أكبر من أو يساوي صفرًا، ويساوي صفرًا فقط في حالة التنبؤ المثالي (أي عندماصك=1{\displaystyle p_{k}=1}وyك=1{\displaystyle y_{k}=1}، أوصك=0{\displaystyle p_{k}=0}وyك=0{\displaystyle y_{k}=0}وتقترب من اللانهاية كلما ساءت التوقعات (أي عندماyك=1{\displaystyle y_{k}=1}وصك0{\displaystyle p_{k}\to 0}أوyك=0{\displaystyle y_{k}=0}وصك1{\displaystyle p_{k}\to 1}مما يعني أن النتيجة الفعلية "أكثر إثارة للدهشة". وبما أن قيمة الدالة اللوجستية تقع دائمًا بين الصفر والواحد، فإن الخسارة اللوغاريتمية تكون دائمًا أكبر من الصفر وأقل من اللانهاية. على عكس الانحدار الخطي، حيث يمكن أن تكون خسارة النموذج صفرًا عند نقطة ما بالمرور عبر نقطة بيانات (وخسارة صفرية إجمالًا إذا كانت جميع النقاط على خط مستقيم)، فإنه في الانحدار اللوجستي لا يمكن أن تكون الخسارة صفرًا عند أي نقطة، لأنyك{\displaystyle y_{k}}إما أن تكون القيمة 0 أو 1، ولكن0<صك<1{\displaystyle 0<p_{k}<1} .

يمكن دمج هذه العناصر في تعبير واحد:

ك=-yكlnصك-(1-yك)ln(1-صك).{\displaystyle \ell _{k}=-y_{k}\ln p_{k}-(1-y_{k})\ln(1-p_{k}).}

يُعرف هذا التعبير بشكل أكثر رسمية باسم الإنتروبيا المتقاطعة للتوزيع المتوقع(صك،(1-صك)){\displaystyle {\big (}p_{k},(1-p_{k}){\big )}}من التوزيع الفعلي(yك،(1-yك)){\displaystyle {\big (}y_{k},(1-y_{k}){\big )}}، كتوزيعات احتمالية على فضاء العنصرين (ناجح، فاشل).

مجموع هذه الخسائر، أي الخسارة الكلية، هو الاحتمالية اللوغاريتمية السالبة الإجمالية-{\displaystyle -\ell }ويتم الحصول على أفضل تطابق لتلك الخيارات منβ0{\displaystyle \beta _{0}}وβ1{\displaystyle \beta _{1}} الذي-{\displaystyle -\ell }يتم تقليلها إلى الحد الأدنى .

بدلاً من ذلك، بدلاً من تقليل الخسارة، يمكن للمرء أن يزيد من معكوسها، وهو دالة الاحتمال اللوغاريتمي (الموجبة):

=ك:yك=1ln(صك)+ك:yك=0ln(1-صك)=ك=1ك(yكln(صك)+(1-yك)ln(1-صك)){\displaystyle \ell =\sum _{k:y_{k}=1}\ln(p_{k})+\sum _{k:y_{k}=0}\ln(1-p_{k})=\sum _{k=1}^{K}\left(\,y_{k}\ln(p_{k})+(1-y_{k})\ln(1-p_{k})\right)}

أو بمعنى آخر، تعظيم دالة الاحتمال نفسها، وهي احتمال أن تكون مجموعة البيانات المعطاة ناتجة عن دالة لوجستية معينة:

ل=ك:yك=1صكك:yك=0(1-صك){\displaystyle L=\prod _{k:y_{k}=1}p_{k}\,\prod _{k:y_{k}=0}(1-p_{k})}

تُعرف هذه الطريقة باسم تقدير الاحتمال الأقصى .

تقدير المعلمات

بما أن غير خطية في β0{\displaystyle \beta _{0}}وβ1{\displaystyle \beta _{1}}يتطلب تحديد قيمها المثلى استخدام طرق عددية. إحدى طرق تعظيم هي حساب مشتقات بالنسبة إلىβ0{\displaystyle \beta _{0}}وβ1{\displaystyle \beta _{1}}أن يكون صفرًا:

0=β0=ك=1ك(yك-صك){\displaystyle 0={\frac {\partial \ell }{\partial \beta _{0}}}=\sum _{k=1}^{K}(y_{k}-p_{k})}
0=β1=ك=1ك(yك-صك)xك{\displaystyle 0={\frac {\partial \ell }{\partial \beta _{1}}}=\sum _{k=1}^{K}(y_{k}-p_{k})x_{k}}

ويمكن إنجاز عملية التعظيم عن طريق حل المعادلتين أعلاه لإيجاد قيمة β0{\displaystyle \beta _{0}}وβ1{\displaystyle \beta _{1}}وهذا ، مرة أخرى، سيتطلب عمومًا استخدام الأساليب العددية.

قيم β0{\displaystyle \beta _{0}}وβ1{\displaystyle \beta _{1}}تم التوصل إلى القيم التي تزيد من قيمتي و L باستخدام البيانات المذكورة أعلاه على النحو التالي:

β0-4.1{\displaystyle \beta _{0}\approx -4.1}
β11.5{\displaystyle \beta _{1}\approx 1.5}

مما ينتج عنه قيمة لـ μ و s كالتالي:

μ=-β0/β12.7{\displaystyle \mu =-\beta _{0}/\beta _{1}\approx 2.7}
s=1/β10.67{\displaystyle s=1/\beta _{1}\approx 0.67}

التوقعات

الـβ0{\displaystyle \beta _{0}}وβ1{\displaystyle \beta _{1}}يمكن إدخال المعاملات في معادلة الانحدار اللوجستي لتقدير احتمالية اجتياز الامتحان.

على سبيل المثال، بالنسبة للطالب الذي يدرس ساعتين، أدخل القيمةx=2{\displaystyle x=2}بإدخال هذه القيمة في المعادلة، نحصل على الاحتمالية التقديرية لاجتياز الامتحان وهي 0.25:

ت=β0+2β1-4.1+21.5=-1.1{\displaystyle t=\beta _{0}+2\beta _{1}\approx -4.1+2\cdot 1.5=-1.1}
ص=11+هـ-ت0.25=احتمالية اجتياز الامتحان{\displaystyle p={\frac {1}{1+e^{-t}}}\approx 0.25={\text{Probability of passing exam}}}

وبالمثل، بالنسبة للطالب الذي يدرس 4 ساعات، فإن الاحتمالية المقدرة لاجتياز الامتحان هي 0.87:

ت=β0+4β1-4.1+41.5=1.9{\displaystyle t=\beta _{0}+4\beta _{1}\approx -4.1+4\cdot 1.5=1.9}
ص=11+هـ-ت0.87=احتمالية اجتياز الامتحان{\displaystyle p={\frac {1}{1+e^{-t}}}\approx 0.87={\text{Probability of passing exam}}}

يوضح هذا الجدول الاحتمالية المقدرة لاجتياز الامتحان لعدة قيم من ساعات الدراسة.

ساعات الدراسة ( س )اجتياز الامتحان
احتمالات اللوغاريتم ( t )الاحتمالات ( هـ ت )الاحتمال ( p )
1-2.570.076 ≈ 1:13.10.07
2-1.070.34 ≈ 1:2.910.26
μ2.7{\displaystyle \mu \approx 2.7}01١ / ٢ = ٠ ٫ ٥٠
30.441.550.61
41.946.960.87
53.4531.40.97

تقييم النموذج

يُعطي تحليل الانحدار اللوجستي النتائج التالية.

المعاملالخطأ المعياريقيمة zقيمة p (والد)
نقطة التقاطع ( β 0 )-4.11.8-2.30.021
الساعات ( β 1 )1.50.9 1.70.017

بحسب اختبار والد ، تشير النتائج إلى أن عدد ساعات الدراسة يرتبط ارتباطًا وثيقًا باحتمالية اجتياز الامتحان (ص=0.017{\displaystyle p=0.017}بدلاً من طريقة والد، فإن الطريقة الموصى بها [ 21 ] لحساب قيمة p للانحدار اللوجستي هي اختبار نسبة الاحتمال (LRT)، والذي يعطي لهذه البياناتص0.00064{\displaystyle p\approx 0.00064}(انظر §  اختبارات الانحراف ونسبة الاحتمال أدناه).

التعميمات

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

خلفية

الشكل 1. الدالة اللوجستية القياسيةσ(ت){\displaystyle \sigma (t)}؛σ(ت)(0،1){\displaystyle \sigma (t)\in (0,1)}للجميعت{\displaystyle t}.

تعريف الدالة اللوجستية

يمكن البدء بشرح الانحدار اللوجستي من خلال شرح الدالة اللوجستية القياسية . الدالة اللوجستية هي دالة سيجمويد ، تأخذ أي مدخل حقيقي .ت{\displaystyle t}وتُخرج قيمة بين الصفر والواحد. [ 2 ] بالنسبة للوغيت، يُفسَّر هذا على أنه يأخذ مُدخلات لوغاريتمية للنسبة الاحتمالية ويُخرج احتمالًا . دالة اللوجستي القياسيةσ:R(0،1){\displaystyle \sigma :\mathbb {R} \rightarrow (0,1)} تُعرَّف على النحو التالي:

σ(ت)=هـتهـت+1=11+هـ-ت{\displaystyle \sigma (t)={\frac {e^{t}}{e^{t}+1}}={\frac {1}{1+e^{-t}}}}

يظهر في الشكل 1 رسم بياني للدالة اللوجستية على الفترة الزمنية t (−6,6).

لنفترض أنت{\displaystyle t}هي دالة خطية لمتغير تفسيري واحدx{\displaystyle x}(الحالة التيت{\displaystyle t}(يُعامل التركيب الخطي لمتغيرات تفسيرية متعددة بشكل مماثل). يمكننا بعد ذلك التعبير عنت{\displaystyle t}على النحو التالي:

ت=β0+β1x{\displaystyle t=\beta _{0}+\beta _{1}x}

والدالة اللوجستية العامةص:R(0،1){\displaystyle p:\mathbb {R} \rightarrow (0,1)}يمكن الآن كتابتها على النحو التالي:

ص(x)=σ(ت)=11+هـ-(β0+β1x){\displaystyle p(x)=\sigma (t)={\frac {1}{1+e^{-(\beta _{0}+\beta _{1}x)}}}}

في النموذج اللوجستي،ص(x){\displaystyle p(x)}يُفسر على أنه احتمال المتغير التابعY{\displaystyle Y}يُعادل ذلك نجاحاً/حالةً بدلاً من فشل/عدم وجود حالة. من الواضح أن متغيرات الاستجابةYأنا{\displaystyle Y_{i}}لا تتوزع بشكل متطابق:P(Yأنا=1|X){\displaystyle P(Y_{i}=1\mid X)}يختلف عن نقطة بيانات واحدةXأنا{\displaystyle X_{i}}بالنسبة لبعضها البعض، على الرغم من استقلالها بالنظر إلى مصفوفة التصميم.X{\displaystyle X}والمعايير المشتركةβ{\displaystyle \beta }[ 11 ]

تعريف معكوس الدالة اللوجستية

يمكننا الآن تعريف دالة اللوجيت (اللوغاريتم الاحتمالي) على أنها الدالة العكسيةز=σ-1{\displaystyle g=\sigma ^{-1}}من الدالة اللوجستية القياسية. من السهل ملاحظة أنها تحقق ما يلي:

ز(ص(x))=σ-1(ص(x))=لوجيتص(x)=ln(ص(x)1-ص(x))=β0+β1x،{\displaystyle g(p(x))=\sigma ^{-1}(p(x))=\operatorname {logit} p(x)=\ln \left({\frac {p(x)}{1-p(x)}}\right)=\beta _{0}+\beta _{1}x,}

وبالمثل، بعد رفع كلا الطرفين إلى الأس، نحصل على الاحتمالات:

ص(x)1-ص(x)=هـβ0+β1x.{\displaystyle {\frac {p(x)}{1-p(x)}}=e^{\beta _{0}+\beta _{1}x}.}

تفسير هذه المصطلحات

في المعادلات أعلاه، تكون الحدود كما يلي:

  • ز{\displaystyle g}هي دالة اللوجيت. معادلة لـز(ص(x)){\displaystyle g(p(x))}يوضح ذلك أن اللوجيت (أي اللوغاريتم الطبيعي للنسبة الاحتمالية) يعادل تعبير الانحدار الخطي.
  • ln{\displaystyle \ln }يرمز إلى اللوغاريتم الطبيعي .
  • ص(x){\displaystyle p(x)}هي احتمالية أن يساوي المتغير التابع حالة معينة، بالنظر إلى توليفة خطية معينة من المتغيرات المستقلة. صيغة لـص(x){\displaystyle p(x)}يوضح هذا أن احتمال تساوي المتغير التابع لحالة معينة يساوي قيمة الدالة اللوجستية لمعادلة الانحدار الخطي. وهذا مهم لأنه يُبين أن قيمة معادلة الانحدار الخطي يمكن أن تتراوح من سالب ما لا نهاية إلى موجب ما لا نهاية، ومع ذلك، بعد التحويل، فإن الصيغة الناتجة للاحتمال تبقى ثابتة.ص(x){\displaystyle p(x)}تتراوح قيمتها بين 0 و 1.
  • β0{\displaystyle \beta _{0}}يمثل هذا التقاطع من معادلة الانحدار الخطي (قيمة المعيار عندما يكون المتنبئ مساوياً للصفر).
  • β1x{\displaystyle \beta _{1}x}هو معامل الانحدار مضروبًا في قيمة معينة للمتغير التنبؤي.
  • قاعدةهـ{\displaystyle e}يرمز إلى الدالة الأسية.

تعريف الاحتمالات

احتمالات أن يساوي المتغير التابع حالة معينة (بالنظر إلى تركيبة خطية معينة)x{\displaystyle x}(من المتغيرات التنبؤية) يُعادل الدالة الأسية لمعادلة الانحدار الخطي. يوضح هذا كيف تعمل دالة اللوجيت كدالة ربط بين الاحتمالية ومعادلة الانحدار الخطي. ونظرًا لأن قيمة اللوجيت تتراوح بين سالب ما لا نهاية وموجب ما لا نهاية، فإنها توفر معيارًا مناسبًا لإجراء الانحدار الخطي، ويمكن تحويل قيمة اللوجيت بسهولة إلى نسبة الاحتمالات. [ 2 ]

لذا نُعرّف احتمالات تساوي المتغير التابع مع حالة معينة (بالنظر إلى تركيبة خطية معينة).x{\displaystyle x}(من المتغيرات التنبؤية) على النحو التالي:

احتمال=هـβ0+β1x.{\displaystyle {\text{odds}}=e^{\beta _{0}+\beta _{1}x}.}

نسبة الاحتمالات

بالنسبة للمتغير المستقل المستمر، يمكن تعريف نسبة الاحتمالات على النحو التالي:

تمثل الصورة مخططًا لما تبدو عليه نسبة الاحتمالات كتابةً، من خلال نموذج بالإضافة إلى مثال درجة الاختبار في قسم "أمثلة" من المحتويات. بعبارة أخرى، إذا افترضنا أن نسبة الاحتمالات هي 2 إلى 1، فيمكننا القول: "مقابل كل زيادة بمقدار وحدة واحدة في ساعات الدراسة، تكون احتمالات النجاح (المجموعة 1) أو الرسوب (المجموعة 0) (كما هو متوقع) 2 إلى 1" (دينيس، 2019).
ياR=احتمال(x+1)احتمال(x)=(ص(x+1)1-ص(x+1))(ص(x)1-ص(x))=هـβ0+β1(x+1)هـβ0+β1x=هـβ1{\displaystyle \mathrm {OR} ={\frac {\operatorname {odds} (x+1)}{\operatorname {odds} (x)}}={\frac {\left({\frac {p(x+1)}{1-p(x+1)}}\right)}{\left({\frac {p(x)}{1-p(x)}}\right)}}={\frac {e^{\beta _{0}+\beta _{1}(x+1)}}{e^{\beta _{0}+\beta _{1}x}}}=e^{\beta _{1}}}

توفر هذه العلاقة الأسية تفسيراً لـβ1{\displaystyle \beta _{1}}تتضاعف الاحتمالات بـهـβ1{\displaystyle e^{\beta _{1}}}لكل زيادة بمقدار وحدة واحدة في x. [ 22 ]

بالنسبة لمتغير مستقل ثنائي، تُعرَّف نسبة الاحتمالات على النحو التالي:أدبج{\displaystyle {\frac {ad}{bc}}}حيث تمثل a و b و c و d خلايا في جدول توافق 2×2 . [ 23 ]

متغيرات تفسيرية متعددة

في حالة وجود متغيرات تفسيرية متعددة، فإن التعبير أعلاهβ0+β1x{\displaystyle \beta _{0}+\beta _{1}x}يمكن تعديلها إلىβ0+β1x1+β2x2++βمxم=β0+أنا=1مβأناxأنا{\displaystyle \beta _{0}+\beta _{1}x_{1}+\beta _{2}x_{2}+\cdots +\beta _{m}x_{m}=\beta _{0}+\sum _{i=1}^{m}\beta _{i}x_{i}}ثم عند استخدام هذا في المعادلة التي تربط اللوغاريتم الطبيعي لاحتمالية النجاح بقيم المتغيرات التنبؤية، سيصبح الانحدار الخطي انحدارًا متعددًا مع m من المتغيرات التفسيرية؛ المعاملاتβأنا{\displaystyle \beta _{i}}للجميعأنا=0،1،2،...،م{\displaystyle i=0,1,2,\dots ,m}جميعها تقديرات.

ومرة أخرى، المعادلات الأكثر تقليدية هي:

سجلص1-ص=β0+β1x1+β2x2++βمxم{\displaystyle \log {\frac {p}{1-p}}=\beta _{0}+\beta _{1}x_{1}+\beta _{2}x_{2}+\cdots +\beta _{m}x_{m}}

و

ص=11+ب-(β0+β1x1+β2x2++βمxم){\displaystyle p={\frac {1}{1+b^{-(\beta _{0}+\beta _{1}x_{1}+\beta _{2}x_{2}+\cdots +\beta _{m}x_{m})}}}}

حيث عادةب=هـ{\displaystyle b=e}.

تعريف

تحتوي مجموعة البيانات على N نقطة. تتكون كل نقطة i من مجموعة من m متغيرات إدخال x1 , i ... xm , i (وتُسمى أيضًا المتغيرات المستقلة ، أو المتغيرات التفسيرية، أو المتغيرات التنبؤية، أو السمات، أو الخصائص)، ومتغير نتيجة ثنائي Yi ( يُعرف أيضًا بالمتغير التابع ، أو متغير الاستجابة، أو متغير الإخراج، أو الفئة)، أي أنه لا يمكن أن يأخذ إلا القيمتين المحتملتين 0 (والتي غالبًا ما تعني "لا" أو "فشل") أو 1 (والتي غالبًا ما تعني "نعم" أو "نجاح"). يهدف الانحدار اللوجستي إلى استخدام مجموعة البيانات لإنشاء نموذج تنبؤي لمتغير النتيجة.

كما هو الحال في الانحدار الخطي، يُفترض أن تعتمد متغيرات النتيجة Y i على المتغيرات التفسيرية x i ... x m،i .

المتغيرات التفسيرية

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

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

المتغيرات الناتجة

بصورة رسمية، تُوصف النتائج Yᵢ بأنها بيانات موزعة وفقًا لتوزيع برنولي ، حيث تُحدد كل نتيجة باحتمال غير مُلاحظ pᵢ خاص بالنتيجة قيد الدراسة، ولكنه مرتبط بالمتغيرات التفسيرية. ويمكن التعبير عن ذلك بأي من الصيغ المكافئة التالية :

Yأنا|x1،أنا،...،xم،أنا برنولي(صأنا)هـ[Yأنا|x1،أنا،...،xم،أنا]=صأنابرو(Yأنا=y|x1،أنا،...،xم،أنا)={صأنالو y=11-صأنالو y=0برو(Yأنا=y|x1،أنا،...،xم،أنا)=صأناy(1-صأنا)(1-y){\displaystyle {\begin{aligned}Y_{i}\mid x_{1,i},\ldots ,x_{m,i}\ &\sim \operatorname {Bernoulli} (p_{i})\\[5pt]\operatorname {\mathbb {E} } [Y_{i}\mid x_{1,i},\ldots ,x_{m,i}]&=p_{i}\\[5pt]\Pr(Y_{i}=y\mid x_{1,i},\ldots ,x_{m,i})&={\begin{cases}p_{i}&{\text{if }}y=1\\1-p_{i}&{\text{if }}y=0\end{cases}}\\[5pt]\Pr(Y_{i}=y\mid x_{1,i},\ldots ,x_{m,i})&=p_{i}^{y}(1-p_{i})^{(1-y)}\end{aligned}}}

معاني هذه الأسطر الأربعة هي:

  1. يُعبّر السطر الأول عن التوزيع الاحتمالي لكل قيمة من قيم Yᵢ : بافتراض ثبات المتغيرات التفسيرية، فإنه يتبع توزيع برنولي بمعاملات pᵢ ، وهي احتمال الحصول على النتيجة 1 في المحاولة i . وكما ذُكر سابقًا، لكل محاولة احتمال نجاحها الخاص، تمامًا كما لكل محاولة متغيراتها التفسيرية الخاصة. لا يُلاحظ احتمال النجاح pᵢ ، وإنما تُلاحظ فقط نتيجة محاولة برنولي الفردية باستخدام هذا الاحتمال. 
  2. يُعبّر السطر الثاني عن حقيقة أن القيمة المتوقعة لكل Yᵢ تساوي احتمال النجاح pᵢ ، وهي خاصية عامة لتوزيع برنولي. بعبارة أخرى، إذا أجرينا عددًا كبيرًا من تجارب برنولي باستخدام نفس احتمال النجاح pᵢ ، ثم أخذنا متوسط ​​جميع النتائج 1 و0، فسيكون الناتج قريبًا من pᵢ . وذلك لأن حساب المتوسط ​​بهذه الطريقة يُحسب ببساطة نسبة النجاحات المُشاهدة، والتي نتوقع أن تتقارب مع احتمال النجاح الأصلي.
  3. يكتب السطر الثالث دالة الكتلة الاحتمالية لتوزيع برنولي، ويحدد احتمال رؤية كل من النتيجتين المحتملتين.
  4. السطر الرابع هو طريقة أخرى لكتابة دالة الكتلة الاحتمالية، تتجنب كتابة الحالات المنفصلة وتُعدّ أكثر ملاءمة لأنواع معينة من الحسابات. يعتمد هذا على حقيقة أن Yᵢ لا يمكن أن تأخذ إلا القيمة 0 أو 1. في كلتا الحالتين، سيكون أحد الأسين 1، ما يعني "اختيار" القيمة التي تحته، بينما يكون الآخر 0، ما يعني "إلغاء" القيمة التي تحته. وبالتالي، تكون النتيجة إما pᵢ أو 1 − pᵢ ، كما في السطر السابق.  
دالة التنبؤ الخطي

تعتمد الفكرة الأساسية للانحدار اللوجستي على استخدام الآلية المطورة بالفعل للانحدار الخطي من خلال نمذجة الاحتمالية pᵢ باستخدام دالة تنبؤ خطية ، أي توليفة خطية من المتغيرات التفسيرية ومجموعة من معاملات الانحدار الخاصة بالنموذج قيد الدراسة ولكنها نفسها لجميع التجارب.و(أنا){\displaystyle f(i)}بالنسبة لنقطة بيانات معينة i ، يتم كتابتها على النحو التالي:

و(أنا)=β0+β1x1،أنا++βمxم،أنا،{\displaystyle f(i)=\beta _{0}+\beta _{1}x_{1,i}+\cdots +\beta _{m}x_{m,i},}

أينβ0،...،βم{\displaystyle \beta _{0},\ldots ,\beta _{m}}هي معاملات الانحدار التي تشير إلى التأثير النسبي لمتغير تفسيري معين على النتيجة.

عادةً ما يتم وضع النموذج في شكل أكثر إيجازًا على النحو التالي:

  • يتم تجميع معاملات الانحدار β 0 ، β 1 ، ... ، β m في متجه واحد β بحجم m  +  1.
  • لكل نقطة بيانات i ، يتم إضافة متغير زائف تفسيري إضافي x 0, i ، بقيمة ثابتة 1، تتوافق مع معامل التقاطع β 0 .
  • ثم يتم تجميع المتغيرات التفسيرية الناتجة x 0, i , x 1, i , ..., x m,i في متجه واحد X i بحجم m  +  1.

وهذا يجعل من الممكن كتابة دالة التنبؤ الخطي على النحو التالي:

و(أنا)=βXأنا،{\displaystyle f(i)={\boldsymbol {\beta }}\cdot \mathbf {X} _{i},}

باستخدام الترميز الخاص بالضرب النقطي بين متجهين.

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

العديد من المتغيرات التفسيرية، فئتان

يمكن تعميم المثال المذكور أعلاه للانحدار اللوجستي الثنائي على متغير تفسيري واحد ليشمل الانحدار اللوجستي الثنائي على أي عدد من المتغيرات التفسيرية x1 ، x2 ، ... وأي عدد من القيم الفئوية .y=0،1،2،...{\displaystyle y=0,1,2,\dots }.

بدايةً، يمكننا النظر في نموذج لوجستي ذي M متغيرًا تفسيريًا، x1 ، x2 ، ... ، xM ، وكما في المثال أعلاه، قيمتين فئويتين ( y = 0 و1). بالنسبة لنموذج الانحدار اللوجستي الثنائي البسيط، افترضنا وجود علاقة خطية بين المتغير التنبؤي واللوغاريتم الطبيعي لنسبة احتمالية وقوع الحدث (المعروف أيضًا باسم اللوجيت ).y=1{\displaystyle y=1}يمكن تعميم هذه العلاقة الخطية على حالة وجود M من المتغيرات التفسيرية:

ت=سجلبص1-ص=β0+β1x1+β2x2++βمxم{\displaystyle t=\log _{b}{\frac {p}{1-p}}=\beta _{0}+\beta _{1}x_{1}+\beta _{2}x_{2}+\cdots +\beta _{M}x_{M}}

حيث t هو اللوغاريتم الطبيعي للنسبة الاحتمالية وβأنا{\displaystyle \beta _{i}}هي معلمات النموذج. وقد تم إدخال تعميم إضافي لا يقتصر فيه أساس النموذج ( b ) على عدد أويلر e . في معظم التطبيقات، يكون الأساسب{\displaystyle b}عادةً ما يتم اعتبار اللوغاريتم هو e . ومع ذلك، في بعض الحالات قد يكون من الأسهل توصيل النتائج بالعمل في الأساس 2 أو الأساس 10.

لتبسيط الترميز، سنحدد المتغيرات التفسيرية ومعاملات β على النحو التالي :(م+1){\displaystyle (M+1)}المتجهات ذات الأبعاد :

x={x0،x1،x2،...،xم}{\displaystyle {\boldsymbol {x}}=\{x_{0},x_{1},x_{2},\dots ,x_{M}\}}
β={β0،β1،β2،...،βم}{\displaystyle {\boldsymbol {\beta }}=\{\beta _{0},\beta _{1},\beta _{2},\dots ,\beta _{M}\}}

مع إضافة متغير تفسيري x 0 =1. يمكن الآن كتابة اللوجيت على النحو التالي:

ت=م=0مβمxم=βx{\displaystyle t=\sum _{m=0}^{M}\beta _{m}x_{m}={\boldsymbol {\beta }}\cdot x}

إيجاد الاحتمال p أنy=1{\displaystyle y=1}ينتج عنه:

ص(x)=بβx1+بβx=11+ب-βx=Sب(ت){\displaystyle p({\boldsymbol {x}})={\frac {b^{{\boldsymbol {\beta }}\cdot {\boldsymbol {x}}}}{1+b^{{\boldsymbol {\beta }}\cdot {\boldsymbol {x}}}}}={\frac {1}{1+b^{-{\boldsymbol {\beta }}\cdot {\boldsymbol {x}}}}}=S_{b}(t)}،

أينSب{\displaystyle S_{b}}هي دالة سيجمويد ذات أساسب{\displaystyle b}تُظهر الصيغة أعلاه أنه بمجردβم{\displaystyle \beta _{m}}إذا كانت المتغيرات ثابتة، فيمكننا بسهولة حساب احتمالات اللوغاريتم التيy=1{\displaystyle y=1}بالنسبة لملاحظة معينة، أو احتمال أنy=1{\displaystyle y=1}بالنسبة لملاحظة معينة. تتمثل حالة الاستخدام الرئيسية للنموذج اللوجستي في إعطاء ملاحظةx{\displaystyle {\boldsymbol {x}}}وتقدير الاحتماليةص(x){\displaystyle p({\boldsymbol {x}})}الذي - التيy=1{\displaystyle y=1}يمكن إيجاد معاملات بيتا المثلى مرة أخرى عن طريق تعظيم دالة الاحتمال اللوغاريتمي. بالنسبة لقياسات K ، يتم تعريفxك{\displaystyle {\boldsymbol {x}}_{k}}باعتباره المتجه التفسيري للقياس رقم k ، وyك{\displaystyle y_{k}}باعتبارها النتيجة التصنيفية لهذا القياس، يمكن كتابة دالة الاحتمال اللوغاريتمي في شكل مشابه جدًا للشكل البسيطم=1{\displaystyle M=1}الحالة المذكورة أعلاه:

=ك=1كyكسجلب(ص(xك))+ك=1ك(1-yك)سجلب(1-ص(xك)){\displaystyle \ell =\sum _{k=1}^{K}y_{k}\log _{b}(p({\boldsymbol {x_{k}}}))+\sum _{k=1}^{K}(1-y_{k})\log _{b}(1-p({\boldsymbol {x_{k}}}))}

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

βم=0=ك=1كyكxمك-ك=1كص(xك)xمك{\displaystyle {\frac {\partial \ell }{\partial \beta _{m}}}=0=\sum _{k=1}^{K}y_{k}x_{mk}-\sum _{k=1}^{K}p({\boldsymbol {x}}_{k})x_{mk}}

حيث x mk هي قيمة المتغير التفسيري x m من القياس رقم k .

لنأخذ مثالاً معم=2{\displaystyle M=2}المتغيرات التفسيرية،ب=10{\displaystyle b=10}، والمعاملاتβ0=-3{\displaystyle \beta _{0}=-3}،β1=1{\displaystyle \beta _{1}=1}، وβ2=2{\displaystyle \beta _{2}=2}والتي تم تحديدها بالطريقة المذكورة أعلاه. وبشكل أكثر تحديداً، فإن النموذج هو:

ت=سجل10ص1-ص=-3+x1+2x2{\displaystyle t=\log _{10}{\frac {p}{1-p}}=-3+x_{1}+2x_{2}}
ص=بβx1+بβx=بβ0+β1x1+β2x21+بβ0+β1x1+β2x2=11+ب-(β0+β1x1+β2x2){\displaystyle p={\frac {b^{{\boldsymbol {\beta }}\cdot {\boldsymbol {x}}}}{1+b^{{\boldsymbol {\beta }}\cdot x}}}={\frac {b^{\beta _{0}+\beta _{1}x_{1}+\beta _{2}x_{2}}}{1+b^{\beta _{0}+\beta _{1}x_{1}+\beta _{2}x_{2}}}}={\frac {1}{1+b^{-(\beta _{0}+\beta _{1}x_{1}+\beta _{2}x_{2})}}}}،

حيث p هو احتمال وقوع الحدث الذيy=1{\displaystyle y=1}ويمكن تفسير ذلك على النحو التالي:

  • β0=-3{\displaystyle \beta _{0}=-3}يمثل نقطة تقاطع المحور الصادي . وهو اللوغاريتم الطبيعي لاحتمالية وقوع الحدث الذيy=1{\displaystyle y=1}عندما تكون المتغيرات التنبؤيةx1=x2=0{\displaystyle x_{1}=x_{2}=0}من خلال الرفع إلى الأس، يمكننا أن نرى أنه عندماx1=x2=0{\displaystyle x_{1}=x_{2}=0}احتمالات وقوع الحدث الذيy=1{\displaystyle y=1}هي من 1 إلى 1000، أو10-3{\displaystyle 10^{-3}}وبالمثل، فإن احتمال وقوع الحدث الذيy=1{\displaystyle y=1}متىx1=x2=0{\displaystyle x_{1}=x_{2}=0}يمكن حسابها على النحو التالي1/(1000+1)=1/1001.{\displaystyle 1/(1000+1)=1/1001.}
  • β1=1{\displaystyle \beta _{1}=1}زيادةx1{\displaystyle x_{1}}يزيد الرقم بمقدار 1 من احتمالات اللوغاريتم بمقدار1{\displaystyle 1}لذا إذاx1{\displaystyle x_{1}}تزداد احتمالات أنy=1{\displaystyle y=1}زيادة بمعامل101{\displaystyle 10^{1}}احتماليةy=1{\displaystyle y=1}وقد ازدادت أيضاً، لكنها لم تزد بنفس القدر الذي ازدادت به الاحتمالات.
  • β2=2{\displaystyle \beta _{2}=2}زيادةx2{\displaystyle x_{2}}يزيد الرقم بمقدار 1 من احتمالات اللوغاريتم بمقدار2{\displaystyle 2}لذا إذاx2{\displaystyle x_{2}}تزداد احتمالات أنy=1{\displaystyle y=1}زيادة بمعامل102.{\displaystyle 10^{2}.}لاحظ كيف أثرx2{\displaystyle x_{2}}يكون تأثير اللوغاريتم الاحتمالي ضعف تأثيرx1{\displaystyle x_{1}}لكن التأثير على الاحتمالات أكبر بعشر مرات. لكن التأثير على احتماليةy=1{\displaystyle y=1}ليس أكبر بعشر مرات، إنما التأثير على الاحتمالات هو الأكبر بعشر مرات فقط.

الانحدار اللوجستي متعدد الحدود: العديد من المتغيرات التفسيرية والعديد من الفئات

في الحالات المذكورة أعلاه لفئتين (الانحدار اللوجستي الثنائي)، تم ترقيم الفئات بـ "0" و"1"، وكان لدينا احتمالان: احتمال أن تكون النتيجة في الفئة 1 يُعطى بواسطةص(x){\displaystyle p({\boldsymbol {x}})}واحتمالية أن تكون النتيجة في الفئة 0 تُعطى بواسطة1-ص(x){\displaystyle 1-p({\boldsymbol {x}})}مجموع هذه الاحتمالات يساوي 1، وهو أمر يجب أن يكون صحيحًا، لأن "0" و "1" هما الفئتان الوحيدتان الممكنتان في هذا الإعداد.

بشكل عام ، إذا كان لدينام+1{\displaystyle M+1}المتغيرات التفسيرية (بما في ذلك x 0 ) وشمال+1{\displaystyle N+1}سنحتاج إلى تصنيفاتشمال+1{\displaystyle N+1}احتمالات منفصلة ، ​​احتمال واحد لكل فئة، مُفهرسة بالرمز n ، تصف احتمال أنتكون النتيجة الفئوية y في الفئة y=n ، بشرط وجود متجه المتغيرات المساعدة x . يجب أن يساوي مجموع هذه الاحتمالات على جميع الفئات 1. باستخدام الأساس الرياضي e ، تكون هذه الاحتمالات كما يلي:

صن(x)=هـβنx1+u=1شمالهـβux{\displaystyle p_{n}({\boldsymbol {x}})={\frac {e^{{\boldsymbol {\beta }}_{n}\cdot {\boldsymbol {x}}}}{1+\sum _{u=1}^{N}e^{{\boldsymbol {\beta }}_{u}\cdot {\boldsymbol {x}}}}}}لن=1،2،...،شمال{\displaystyle n=1,2,\dots ,N}
ص0(x)=1-ن=1شمالصن(x)=11+u=1شمالهـβux{\displaystyle p_{0}({\boldsymbol {x}})=1-\sum _{n=1}^{N}p_{n}({\boldsymbol {x}})={\frac {1}{1+\sum _{u=1}^{N}e^{{\boldsymbol {\beta }}_{u}\cdot {\boldsymbol {x}}}}}}

كل الاحتمالات باستثناءص0(x){\displaystyle p_{0}({\boldsymbol {x}})}سيكون لكل منها مجموعة خاصة بها من معاملات الانحدارβن{\displaystyle {\boldsymbol {\beta }}_{n}}يتضح أنه، كما هو مطلوب، مجموعصن(x){\displaystyle p_{n}({\boldsymbol {x}})}في جميع الفئات، يكون عدد الفئات n هو 1. اختيارص0(x){\displaystyle p_{0}({\boldsymbol {x}})}إن تعريفها بدلالة الاحتمالات الأخرى أمرٌ مصطنع. كان من الممكن اختيار أيٍّ من الاحتمالات لتعريفها بهذه الطريقة. تُسمى هذه القيمة الخاصة لـ n "مؤشر المحور"، ويتم التعبير عن اللوغاريتم الطبيعي للنسبة ( t n ) بدلالة احتمال المحور، ويتم التعبير عنها بدورها كتركيبة خطية للمتغيرات التفسيرية.

تن=ln(صن(x)ص0(x))=βنx{\displaystyle t_{n}=\ln \left({\frac {p_{n}({\boldsymbol {x}})}{p_{0}({\boldsymbol {x}})}}\right)={\boldsymbol {\beta }}_{n}\cdot {\boldsymbol {x}}}

لاحظ أيضًا أنه في الحالة البسيطة لـشمال=1{\displaystyle N=1}تم استرداد القضية المصنفة ضمن فئتين، معص(x)=ص1(x){\displaystyle p({\boldsymbol {x}})=p_{1}({\boldsymbol {x}})}وص0(x)=1-ص1(x){\displaystyle p_{0}({\boldsymbol {x}})=1-p_{1}({\boldsymbol {x}})}.

يمكن الآن حساب احتمالية توليد مجموعة معينة من K قياسات أو نقاط بيانات بواسطة الاحتمالات المذكورة أعلاه. وبفهرسة كل قياس بالرمز k ، فلنرمز إلى المجموعة k من المتغيرات التفسيرية المقاسة بالرمز .xك{\displaystyle {\boldsymbol {x}}_{k}}ويمكن الإشارة إلى نتائجها التصنيفية بـyك{\displaystyle y_{k}}والتي يمكن أن تساوي أي عدد صحيح في الفترة [0، N]. وبالتالي، فإن دالة الاحتمال اللوغاريتمي هي:

=ك=1كن=0شمالΔ(ن،yك)ln(صن(xك)){\displaystyle \ell =\sum _{k=1}^{K}\sum _{n=0}^{N}\Delta (n,y_{k})\,\ln(p_{n}({\boldsymbol {x}}_{k}))}

أينΔ(ن،yك){\displaystyle \Delta (n,y_{k})}هي دالة مؤشر تساوي 1 إذا كان y <sub>k </sub> = n، وتساوي صفرًا فيما عدا ذلك. في حالة وجود متغيرين تفسيريين، عُرِّفت هذه الدالة على أنها y<sub> k</sub> عندما n = 1، و 1 - y<sub> k</sub> عندما n = 0. كان هذا مناسبًا، ولكنه ليس ضروريًا. [ 24 ] مرة أخرى، يمكن إيجاد معاملات بيتا المثلى عن طريق تعظيم دالة الاحتمال اللوغاريتمي باستخدام الطرق العددية عمومًا. إحدى طرق الحل الممكنة هي جعل مشتقات دالة الاحتمال اللوغاريتمي بالنسبة لكل معامل بيتا تساوي صفرًا، ثم حل المعادلة لإيجاد معاملات بيتا.

βنم=0=ك=1كΔ(ن،yك)xمك-ك=1كصن(xك)xمك{\displaystyle {\frac {\partial \ell }{\partial \beta _{nm}}}=0=\sum _{k=1}^{K}\Delta (n,y_{k})x_{mk}-\sum _{k=1}^{K}p_{n}({\boldsymbol {x}}_{k})x_{mk}}

أينβنم{\displaystyle \beta _{nm}}هو المعامل رقم m منβن{\displaystyle {\boldsymbol {\beta }}_{n}}متجه وxمك{\displaystyle x_{mk}}يمثل المتغير التفسيري رقم m للقياس رقم k . بمجرد تقدير معاملات بيتا من البيانات، سنتمكن من تقدير احتمالية أن تؤدي أي مجموعة لاحقة من المتغيرات التفسيرية إلى أي من فئات النتائج الممكنة.

التفسيرات

توجد مواصفات وتفسيرات مكافئة مختلفة للانحدار اللوجستي، والتي تتناسب مع أنواع مختلفة من النماذج الأكثر عمومية، وتسمح بتعميمات مختلفة.

كنموذج خطي معمّم

إن النموذج المحدد المستخدم في الانحدار اللوجستي، والذي يميزه عن الانحدار الخطي القياسي وعن أنواع تحليل الانحدار الأخرى المستخدمة للنتائج ذات القيم الثنائية ، هو الطريقة التي يتم بها ربط احتمالية حدوث نتيجة معينة بدالة التنبؤ الخطية:

لوجيت(هـ[Yأنا|x1،أنا،...،xم،أنا])=لوجيت(صأنا)=ln(صأنا1-صأنا)=β0+β1x1،أنا++βمxم،أنا{\displaystyle \operatorname {logit} (\operatorname {\mathbb {E} } [Y_{i}\mid x_{1,i},\ldots ,x_{m,i}])=\operatorname {logit} (p_{i})=\ln \left({\frac {p_{i}}{1-p_{i}}}\right)=\beta _{0}+\beta _{1}x_{1,i}+\cdots +\beta _{m}x_{m,i}}

مكتوبة باستخدام الصيغة المختصرة الموضحة أعلاه، وهي كالتالي:

لوجيت(هـ[Yأنا|Xأنا])=لوجيت(صأنا)=ln(صأنا1-صأنا)=βXأنا{\displaystyle \operatorname {logit} (\operatorname {\mathbb {E} } [Y_{i}\mid \mathbf {X} _{i}])=\operatorname {logit} (p_{i})=\ln \left({\frac {p_{i}}{1-p_{i}}}\right)={\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}

تعبر هذه الصيغة عن الانحدار اللوجستي كنوع من النموذج الخطي المعمم ، والذي يتنبأ بالمتغيرات بأنواع مختلفة من التوزيعات الاحتمالية عن طريق ملاءمة دالة تنبؤ خطية بالشكل المذكور أعلاه لنوع من التحويل التعسفي للقيمة المتوقعة للمتغير.

تم شرح فكرة التحويل باستخدام دالة اللوجيت (اللوغاريتم الطبيعي للنسبة الاحتمالية) أعلاه . ولها أيضًا تأثير عملي يتمثل في تحويل الاحتمالية (التي تتراوح قيمتها بين 0 و1) إلى متغير يتراوح بين 0 و1.(-،+){\displaystyle (-\infty ,+\infty )}— وبالتالي مطابقة النطاق المحتمل لدالة التنبؤ الخطي على الجانب الأيمن من المعادلة.

كل من الاحتمالات pᵢ ومعاملات الانحدار غير قابلة للملاحظة، ولا تُعدّ طريقة تحديدها جزءًا من النموذج نفسه. عادةً ما يتم تحديدها من خلال نوع من إجراءات التحسين، مثل تقدير الاحتمال الأقصى ، الذي يجد القيم التي تُطابق البيانات المرصودة على أفضل وجه (أي التي تُعطي أدق التنبؤات للبيانات المرصودة بالفعل)، وعادةً ما تخضع لشروط تنظيمية تسعى إلى استبعاد القيم غير المحتملة، مثل القيم الكبيرة جدًا لأي من معاملات الانحدار. يُعادل استخدام شرط تنظيمي إجراء تقدير الاحتمال اللاحق الأقصى (MAP)، وهو امتداد لتقدير الاحتمال الأقصى. (يتم التنظيم عادةً باستخدام دالة تنظيمية تربيعية ، وهي تُعادل وضع توزيع غاوسي مسبق بمتوسط ​​صفري على المعاملات، ولكن توجد منظمات أخرى ممكنة أيضًا). سواءً تم استخدام التنظيم أم لا، فإنه عادةً ما يكون من غير الممكن إيجاد حل مغلق. بدلاً من ذلك، يجب استخدام طريقة عددية تكرارية، مثل المربعات الصغرى الموزونة بشكل تكراري (IRLS) أو، بشكل أكثر شيوعًا هذه الأيام، طريقة شبه نيوتن مثل طريقة L-BFGS . [ 25 ]

يُفسَّر تقدير المعامل β j على أنه التأثير التراكمي على لوغاريتم احتمالات حدوث تغيير بمقدار وحدة واحدة في المتغير التفسيري j . في حالة المتغير التفسيري الثنائي، على سبيل المثال، الجنسهـβ{\displaystyle e^{\beta }}هو تقدير احتمالات حدوث النتيجة، على سبيل المثال، للذكور مقارنة بالإناث.

تستخدم الصيغة المكافئة معكوس دالة اللوجيت، وهي الدالة اللوجستية ، أي:

هـ[Yأنا|Xأنا]=صأنا=لوجيت-1(βXأنا)=11+هـ-βXأنا{\displaystyle \operatorname {\mathbb {E} } [Y_{i}\mid \mathbf {X} _{i}]=p_{i}=\operatorname {logit} ^{-1}({\boldsymbol {\beta }}\cdot \mathbf {X} _{i})={\frac {1}{1+e^{-{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}}}}

يمكن أيضًا كتابة الصيغة على شكل توزيع احتمالي (وتحديدًا باستخدام دالة الكتلة الاحتمالية ):

برو(Yأنا=y|Xأنا)=صأناy(1-صأنا)1-y=(هـβXأنا1+هـβXأنا)y(1-هـβXأنا1+هـβXأنا)1-y=هـβXأناy1+هـβXأنا{\displaystyle \Pr(Y_{i}=y\mid \mathbf {X} _{i})={p_{i}}^{y}(1-p_{i})^{1-y}=\left({\frac {e^{{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}}{1+e^{{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}}}\right)^{y}\left(1-{\frac {e^{{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}}{1+e^{{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}}}\right)^{1-y}={\frac {e^{{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}\cdot y}}{1+e^{{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}}}}

كنموذج متغير كامن

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

تخيل أنه لكل تجربة i ، يوجد متغير كامن مستمر Y i * (أي متغير عشوائي غير ملاحظ ) يتم توزيعه على النحو التالي:

Yأنا*=βXأنا+εأنا{\displaystyle Y_{i}^{\ast }={\boldsymbol {\beta }}\cdot \mathbf {X} _{i}+\varepsilon _{i}\,}

أين

εأناالخدمات اللوجستية(0،1){\displaystyle \varepsilon _{i}\sim \operatorname {Logistic} (0,1)\,}

أي أن المتغير الكامن يمكن كتابته مباشرة بدلالة دالة التنبؤ الخطي ومتغير الخطأ العشوائي الإضافي الذي يتم توزيعه وفقًا للتوزيع اللوجستي القياسي .

ويمكن اعتبار Y i مؤشراً على ما إذا كان هذا المتغير الكامن موجباً أم لا:

Yأنا={1لو Yأنا*>0  أي -εأنا<βXأنا،0خلاف ذلك.{\displaystyle Y_{i}={\begin{cases}1&{\text{if }}Y_{i}^{\ast }>0\ {\text{ i.e. }}{-\varepsilon _{i}}<{\boldsymbol {\beta }}\cdot \mathbf {X} _{i},\\0&{\text{otherwise.}}\end{cases}}}

قد يبدو اختيار نمذجة متغير الخطأ باستخدام التوزيع اللوجستي القياسي، بدلاً من التوزيع اللوجستي العام مع تحديد الموقع والمقياس بقيم عشوائية، مقيدًا، ولكنه في الواقع ليس كذلك. يجب أن نضع في اعتبارنا أنه بإمكاننا اختيار معاملات الانحدار بأنفسنا، وفي كثير من الأحيان يمكننا استخدامها لتعويض التغيرات في معلمات توزيع متغير الخطأ. على سبيل المثال، يُكافئ التوزيع اللوجستي لمتغير الخطأ بمعامل موقع غير صفري μ (الذي يحدد المتوسط) توزيعًا بمعامل موقع صفري، حيث تمت إضافة μ إلى معامل التقاطع. تُنتج كلتا الحالتين القيمة نفسها لـ Y i * بغض النظر عن إعدادات المتغيرات التفسيرية. وبالمثل، يُكافئ معامل المقياس العشوائي s ضبط معامل المقياس على 1 ثم قسمة جميع معاملات الانحدار على s . في الحالة الأخيرة، ستكون القيمة الناتجة لـ Y i * أصغر بمعامل s مقارنة بالحالة الأولى، لجميع مجموعات المتغيرات التفسيرية - ولكن الأهم من ذلك، أنها ستبقى دائمًا على نفس جانب 0، وبالتالي ستؤدي إلى نفس اختيار Y i .

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

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

برو(εأنا<x)=لوجيت-1(x){\displaystyle \Pr(\varepsilon _{i}<x)=\operatorname {logit} ^{-1}(x)}

ثم:

برو(Yأنا=1|Xأنا)=برو(Yأنا*>0|Xأنا)=برو(βXأنا+εأنا>0)=برو(εأنا>-βXأنا)=برو(εأنا<βXأنا)(لأن التوزيع اللوجستي متماثل)=لوجيت-1(βXأنا)=صأنا(انظر أعلاه){\displaystyle {\begin{aligned}\Pr(Y_{i}=1\mid \mathbf {X} _{i})&=\Pr(Y_{i}^{\ast }>0\mid \mathbf {X} _{i})\\[5pt]&=\Pr({\boldsymbol {\beta }}\cdot \mathbf {X} _{i}+\varepsilon _{i}>0)\\[5pt]&=\Pr(\varepsilon _{i}>-{\boldsymbol {\beta }}\cdot \mathbf {X} _{i})\\[5pt]&=\Pr(\varepsilon _{i}<{\boldsymbol {\beta }}\cdot \mathbf {X} _{i})&&{\text{(because the logistic distribution is symmetric)}}\\[5pt]&=\operatorname {logit} ^{-1}({\boldsymbol {\beta }}\cdot \mathbf {X} _{i})&\\[5pt]&=p_{i}&&{\text{(see above)}}\end{aligned}}}

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

نموذج المتغيرات الكامنة ثنائي الاتجاه

وهناك صيغة أخرى تستخدم متغيرين كامنين منفصلين:

Yأنا0*=β0Xأنا+ε0Yأنا1*=β1Xأنا+ε1{\displaystyle {\begin{aligned}Y_{i}^{0\ast }&={\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}+\varepsilon _{0}\,\\Y_{i}^{1\ast }&={\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}+\varepsilon _{1}\,\end{aligned}}}

أين

ε0سيارة كهربائية1(0،1)ε1سيارة كهربائية1(0،1){\displaystyle {\begin{aligned}\varepsilon _{0}&\sim \operatorname {EV} _{1}(0,1)\\\varepsilon _{1}&\sim \operatorname {EV} _{1}(0,1)\end{aligned}}}

حيث EV 1 (0,1) هو توزيع القيمة القصوى القياسي من النوع الأول : أي

برو(ε0=x)=برو(ε1=x)=هـ-xهـ-هـ-x{\displaystyle \Pr(\varepsilon _{0}=x)=\Pr(\varepsilon _{1}=x)=e^{-x}e^{-e^{-x}}}

ثم

Yأنا={1لو Yأنا1*>Yأنا0*،0خلاف ذلك.{\displaystyle Y_{i}={\begin{cases}1&{\text{if }}Y_{i}^{1\ast }>Y_{i}^{0\ast },\\0&{\text{otherwise.}}\end{cases}}}

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

يبدو اختيار توزيع القيمة القصوى من النوع الأول تعسفيًا إلى حد ما، ولكنه يجعل الرياضيات تعمل، وقد يكون من الممكن تبرير استخدامه من خلال نظرية الاختيار العقلاني .

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

β=β1-β0{\displaystyle {\boldsymbol {\beta }}={\boldsymbol {\beta }}_{1}-{\boldsymbol {\beta }}_{0}}
ε=ε1-ε0{\displaystyle \varepsilon =\varepsilon _{1}-\varepsilon _{0}}

ينبع هذا الحدس من حقيقة أنه بما أننا نختار بناءً على القيمة القصوى بين قيمتين، فإن الفرق بينهما هو المهم فقط، وليس القيم نفسها، وهذا يُلغي فعليًا درجة واحدة من درجات الحرية . ومن الحقائق المهمة الأخرى أن الفرق بين متغيرين من النوع الأول يتبع توزيع القيم المتطرفة يتبع التوزيع اللوجستي، أيε=ε1-ε0الخدمات اللوجستية(0،1).{\displaystyle \varepsilon =\varepsilon _{1}-\varepsilon _{0}\sim \operatorname {Logistic} (0,1).}يمكننا توضيح المكافئ على النحو التالي:

برو(Yأنا=1|Xأنا)=برو(Yأنا1*>Yأنا0*|Xأنا)=برو(Yأنا1*-Yأنا0*>0|Xأنا)=برو(β1Xأنا+ε1-(β0Xأنا+ε0)>0)=برو((β1Xأنا-β0Xأنا)+(ε1-ε0)>0)=برو((β1-β0)Xأنا+(ε1-ε0)>0)=برو((β1-β0)Xأنا+ε>0)(بديل ε كما سبق)=برو(βXأنا+ε>0)(بديل β كما سبق)=برو(ε>-βXأنا)(الآن، نفس النموذج المذكور أعلاه)=برو(ε<βXأنا)=لوجيت-1(βXأنا)=صأنا{\displaystyle {\begin{aligned}\Pr(Y_{i}=1\mid \mathbf {X} _{i})={}&\Pr \left(Y_{i}^{1\ast }>Y_{i}^{0\ast }\mid \mathbf {X} _{i}\right)&\\[5pt]={}&\Pr \left(Y_{i}^{1\ast }-Y_{i}^{0\ast }>0\mid \mathbf {X} _{i}\right)&\\[5pt]={}&\Pr \left({\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}+\varepsilon _{1}-\left({\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}+\varepsilon _{0}\right)>0\right)&\\[5pt]={}&\Pr \left(({\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}-{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i})+(\varepsilon _{1}-\varepsilon _{0})>0\right)&\\[5pt]={}&\Pr(({\boldsymbol {\beta }}_{1}-{\boldsymbol {\beta }}_{0})\cdot \mathbf {X} _{i}+(\varepsilon _{1}-\varepsilon _{0})>0)&\\[5pt]={}&\Pr(({\boldsymbol {\beta }}_{1}-{\boldsymbol {\beta }}_{0})\cdot \mathbf {X} _{i}+\varepsilon >0)&&{\text{(substitute }}\varepsilon {\text{ as above)}}\\[5pt]={}&\Pr({\boldsymbol {\beta }}\cdot \mathbf {X} _{i}+\varepsilon >0)&&{\text{(substitute }}{\boldsymbol {\beta }}{\text{ as above)}}\\[5pt]={}&\Pr(\varepsilon >-{\boldsymbol {\beta }}\cdot \mathbf {X} _{i})&&{\text{(now, same as above model)}}\\[5pt]={}&\Pr(\varepsilon <{\boldsymbol {\beta }}\cdot \mathbf {X} _{i})&\\[5pt]={}&\operatorname {logit} ^{-1}({\boldsymbol {\beta }}\cdot \mathbf {X} _{i})\\[5pt]={}&p_{i}\end{aligned}}}

كنموذج "لوغاريتمي خطي"

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

هنا، بدلاً من كتابة لوغاريتم الاحتمالات p i كمتنبئ خطي، نقوم بفصل المتنبئ الخطي إلى اثنين، واحد لكل من النتيجتين:

lnبرو(Yأنا=0)=β0Xأنا-lnZlnبرو(Yأنا=1)=β1Xأنا-lnZ{\displaystyle {\begin{aligned}\ln \Pr(Y_{i}=0)&={\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}-\ln Z\\\ln \Pr(Y_{i}=1)&={\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}-\ln Z\end{aligned}}}

تم إدخال مجموعتين منفصلتين من معاملات الانحدار، تمامًا كما هو الحال في نموذج المتغير الكامن ثنائي الاتجاه، وتظهر المعادلتان في شكل يكتب لوغاريتم الاحتمال المرتبط كمتنبئ خطي، مع حد إضافي.-lnZ{\displaystyle -\ln Z}في النهاية. هذا الحد، كما اتضح، يعمل كعامل توحيد يضمن أن تكون النتيجة توزيعًا. ويمكن ملاحظة ذلك برفع كلا الطرفين إلى الأس:

برو(Yأنا=0)=1Zهـβ0Xأنابرو(Yأنا=1)=1Zهـβ1Xأنا{\displaystyle {\begin{aligned}\Pr(Y_{i}=0)&={\frac {1}{Z}}e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}\\[5pt]\Pr(Y_{i}=1)&={\frac {1}{Z}}e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}\end{aligned}}}

في هذا الشكل ، يتضح أن الغرض من Z هو ضمان أن يكون التوزيع الناتج على Yᵢ توزيعًا احتماليًا بالفعل ، أي أن مجموعها يساوي 1. وهذا يعني أن Z هو ببساطة مجموع جميع الاحتمالات غير المعيارية، وبقسمة كل احتمال على Z ، تصبح الاحتمالات " معيارية ". أي:

Z=هـβ0Xأنا+هـβ1Xأنا{\displaystyle Z=e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}+e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}

والمعادلات الناتجة هي

برو(Yأنا=0)=هـβ0Xأناهـβ0Xأنا+هـβ1Xأنابرو(Yأنا=1)=هـβ1Xأناهـβ0Xأنا+هـβ1Xأنا.{\displaystyle {\begin{aligned}\Pr(Y_{i}=0)&={\frac {e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}}{e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}+e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}}\\[5pt]\Pr(Y_{i}=1)&={\frac {e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}{e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}+e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}}.\end{aligned}}}

أو بشكل عام:

برو(Yأنا=ج)=هـβجXأناحهـβحXأنا{\displaystyle \Pr(Y_{i}=c)={\frac {e^{{\boldsymbol {\beta }}_{c}\cdot \mathbf {X} _{i}}}{\sum _{h}e^{{\boldsymbol {\beta }}_{h}\cdot \mathbf {X} _{i}}}}}

يوضح هذا بوضوح كيفية تعميم هذه الصيغة لتشمل أكثر من نتيجتين، كما هو الحال في نموذج اللوجيت متعدد الحدود . هذه الصيغة العامة هي بالضبط دالة softmax كما في

برو(Yأنا=ج)=سوفتماكس(ج،β0Xأنا،β1Xأنا،...).{\displaystyle \Pr(Y_{i}=c)=\operatorname {softmax} (c,{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i},{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i},\dots ).}

لإثبات أن هذا النموذج مكافئ للنموذج السابق، نبدأ بإدراك أن النموذج المذكور أعلاه مُفرط التحديد، بمعنى أنبرو(Yأنا=0){\displaystyle \Pr(Y_{i}=0)}وبرو(Yأنا=1){\displaystyle \Pr(Y_{i}=1)}لا يمكن تحديدها بشكل مستقل: بلبرو(Yأنا=0)+برو(Yأنا=1)=1{\displaystyle \Pr(Y_{i}=0)+\Pr(Y_{i}=1)=1}لذا فإن معرفة أحدهما تحدد الآخر تلقائيًا. ونتيجة لذلك، يكون النموذج غير قابل للتحديد ، حيث توجد تركيبات متعددة منβ0{\displaystyle {\boldsymbol {\beta }}_{0}}وβ1{\displaystyle {\boldsymbol {\beta }}_{1}}سينتج عن ذلك نفس الاحتمالات لجميع المتغيرات التفسيرية الممكنة. في الواقع، يمكن ملاحظة أن إضافة أي متجه ثابت إلى كليهما سينتج عنه نفس الاحتمالات.

برو(Yأنا=1)=هـ(β1+ج)Xأناهـ(β0+ج)Xأنا+هـ(β1+ج)Xأنا=هـβ1XأناهـجXأناهـβ0XأناهـجXأنا+هـβ1XأناهـجXأنا=هـجXأناهـβ1XأناهـجXأنا(هـβ0Xأنا+هـβ1Xأنا)=هـβ1Xأناهـβ0Xأنا+هـβ1Xأنا.{\displaystyle {\begin{aligned}\Pr(Y_{i}=1)&={\frac {e^{({\boldsymbol {\beta }}_{1}+\mathbf {C} )\cdot \mathbf {X} _{i}}}{e^{({\boldsymbol {\beta }}_{0}+\mathbf {C} )\cdot \mathbf {X} _{i}}+e^{({\boldsymbol {\beta }}_{1}+\mathbf {C} )\cdot \mathbf {X} _{i}}}}\\[5pt]&={\frac {e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}e^{\mathbf {C} \cdot \mathbf {X} _{i}}}{e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}e^{\mathbf {C} \cdot \mathbf {X} _{i}}+e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}e^{\mathbf {C} \cdot \mathbf {X} _{i}}}}\\[5pt]&={\frac {e^{\mathbf {C} \cdot \mathbf {X} _{i}}e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}{e^{\mathbf {C} \cdot \mathbf {X} _{i}}(e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}+e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}})}}\\[5pt]&={\frac {e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}{e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}+e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}}.\end{aligned}}}

ونتيجة لذلك، يمكننا تبسيط الأمور واستعادة إمكانية التحديد عن طريق اختيار قيمة عشوائية لأحد المتجهين. نختار أن نضبطβ0=0.{\displaystyle {\boldsymbol {\beta }}_{0}=\mathbf {0} .} ثم،

هـβ0Xأنا=هـ0Xأنا=1{\displaystyle e^{{\boldsymbol {\beta }}_{0}\cdot \mathbf {X} _{i}}=e^{\mathbf {0} \cdot \mathbf {X} _{i}}=1}

وهكذا

برو(Yأنا=1)=هـβ1Xأنا1+هـβ1Xأنا=11+هـ-β1Xأنا=صأنا{\displaystyle \Pr(Y_{i}=1)={\frac {e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}{1+e^{{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}}={\frac {1}{1+e^{-{\boldsymbol {\beta }}_{1}\cdot \mathbf {X} _{i}}}}=p_{i}}

مما يدل على أن هذه الصيغة مكافئة بالفعل للصيغة السابقة. (كما هو الحال في صيغة المتغير الكامن ثنائي الاتجاه، أي إعدادات حيثβ=β1-β0{\displaystyle {\boldsymbol {\beta }}={\boldsymbol {\beta }}_{1}-{\boldsymbol {\beta }}_{0}}(سيؤدي ذلك إلى نتائج مماثلة.)

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

كشبكة عصبية أحادية الطبقة

للنموذج صيغة مكافئة

صأنا=11+هـ-(β0+β1x1،أنا++βكxك،أنا).{\displaystyle p_{i}={\frac {1}{1+e^{-(\beta _{0}+\beta _{1}x_{1,i}+\cdots +\beta _{k}x_{k,i})}}}.\,}

يُطلق على هذا الشكل الوظيفي عادةً اسم الشبكة العصبية الاصطناعية أحادية الطبقة أو الشبكة العصبية أحادية الطبقة . تُحسب الشبكة العصبية أحادية الطبقة مخرجًا متصلًا بدلًا من دالة متدرجة . تُحسب مشتقة pᵢ بالنسبة إلى X = ( x₁ , ... , xₖ ) من الشكل العام التالي:  

y=11+هـ-و(X){\displaystyle y={\frac {1}{1+e^{-f(X)}}}}

حيث f ( X ) دالة تحليلية في X. بهذا الاختيار، تصبح الشبكة العصبية أحادية الطبقة مطابقة لنموذج الانحدار اللوجستي. تتميز هذه الدالة بمشتقة متصلة، مما يسمح باستخدامها في خوارزمية الانتشار العكسي . كما تُفضّل هذه الدالة لسهولة حساب مشتقتها.

دyدX=y(1-y)دودX.{\displaystyle {\frac {\mathrm {d} y}{\mathrm {d} X}}=y(1-y){\frac {\mathrm {d} f}{\mathrm {d} X}}.\,}

فيما يتعلق بالبيانات ذات الحدين

يفترض نموذج وثيق الصلة أن كل i لا يرتبط بتجربة برنولي واحدة ولكن بـ n i تجارب مستقلة موزعة بشكل متطابق ، حيث تكون الملاحظة Y i هي عدد النجاحات الملاحظة (مجموع المتغيرات العشوائية الفردية الموزعة وفقًا لتوزيع برنولي)، وبالتالي تتبع توزيعًا ذا حدين :

Yأناسلة المهملات(نأنا،صأنا)، ل أنا=1،...،ن{\displaystyle Y_{i}\,\sim \operatorname {Bin} (n_{i},p_{i}),{\text{ for }}i=1,\dots ,n}

ومن الأمثلة على هذا التوزيع نسبة البذور ( pi ) التي تنبت بعد زراعة ni .

من حيث القيم المتوقعة ، يتم التعبير عن هذا النموذج على النحو التالي:

صأنا=هـ[Yأنانأنا|Xأنا]،{\displaystyle p_{i}=\operatorname {\mathbb {E} } \left[\left.{\frac {Y_{i}}{n_{i}}}\,\right|\,\mathbf {X} _{i}\right]\,,}

لهذا السبب.

لوجيت(هـ[Yأنانأنا|Xأنا])=لوجيت(صأنا)=ln(صأنا1-صأنا)=βXأنا،{\displaystyle \operatorname {logit} \left(\operatorname {\mathbb {E} } \left[\left.{\frac {Y_{i}}{n_{i}}}\,\right|\,\mathbf {X} _{i}\right]\right)=\operatorname {logit} (p_{i})=\ln \left({\frac {p_{i}}{1-p_{i}}}\right)={\boldsymbol {\beta }}\cdot \mathbf {X} _{i}\,,}

أو ما يعادل ذلك:

برو(Yأنا=y|Xأنا)=(نأناy)صأناy(1-صأنا)نأنا-y=(نأناy)(11+هـ-βXأنا)y(1-11+هـ-βXأنا)نأنا-y.{\displaystyle \Pr(Y_{i}=y\mid \mathbf {X} _{i})={n_{i} \choose y}p_{i}^{y}(1-p_{i})^{n_{i}-y}={n_{i} \choose y}\left({\frac {1}{1+e^{-{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}}}\right)^{y}\left(1-{\frac {1}{1+e^{-{\boldsymbol {\beta }}\cdot \mathbf {X} _{i}}}}\right)^{n_{i}-y}\,.}

يمكن ملاءمة هذا النموذج باستخدام نفس أنواع الأساليب المستخدمة في النموذج الأساسي المذكور أعلاه.

ملاءمة النموذج

تقدير الاحتمال الأقصى (MLE)

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

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

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

المربعات الصغرى المعاد ترجيحها بشكل متكرر (IRLS)

الانحدار اللوجستي الثنائي (y=0{\displaystyle y=0}أوy=1{\displaystyle y=1}يمكن، على سبيل المثال، حسابها باستخدام طريقة المربعات الصغرى الموزونة بشكل متكرر (IRLS)، والتي تعادل تعظيم دالة الاحتمال اللوغاريتمي لعملية موزعة وفقًا لتوزيع برنولي باستخدام طريقة نيوتن . إذا كُتبت المسألة في شكل مصفوفة متجهة، مع معلماتwتي=[β0،β1،β2،...]{\displaystyle \mathbf {w} ^{T}=[\beta _{0},\beta _{1},\beta _{2},\ldots ]}المتغيرات التفسيريةx(أنا)=[1،x1(أنا)،x2(أنا)،...]تي{\displaystyle \mathbf {x} (i)=[1,x_{1}(i),x_{2}(i),\ldots ]^{T}}والقيمة المتوقعة لتوزيع برنوليμ(أنا)=11+هـ-wتيx(أنا){\displaystyle \mu (i)={\frac {1}{1+e^{-\mathbf {w} ^{T}\mathbf {x} (i)}}}}، المعاييرw{\displaystyle \mathbf {w} }يمكن إيجادها باستخدام الخوارزمية التكرارية التالية:

wك+1=(XتيSكX)-1Xتي(SكXwك+y-μك){\displaystyle \mathbf {w} _{k+1}=\left(\mathbf {X} ^{T}\mathbf {S} _{k}\mathbf {X} \right)^{-1}\mathbf {X} ^{T}\left(\mathbf {S} _{k}\mathbf {X} \mathbf {w} _{k}+\mathbf {y} -\mathbf {\boldsymbol {\mu }} _{k}\right)}

أينS=التشخيص(μ(أنا)(1-μ(أنا))){\displaystyle \mathbf {S} =\operatorname {diag} (\mu (i)(1-\mu (i)))}هي مصفوفة ترجيح قطرية،μ=[μ(1)،μ(2)،...]{\displaystyle {\boldsymbol {\mu }}=[\mu (1),\mu (2),\ldots ]}متجه القيم المتوقعة،

X=[1x1(1)x2(1)...1x1(2)x2(2)...]{\displaystyle \mathbf {X} ={\begin{bmatrix}1&x_{1}(1)&x_{2}(1)&\ldots \\1&x_{1}(2)&x_{2}(2)&\ldots \\\vdots &\vdots &\vdots \end{bmatrix}}}

مصفوفة الانحدار وy(أنا)=[y(1)،y(2)،...]تي{\displaystyle \mathbf {y} (i)=[y(1),y(2),\ldots ]^{T}}متجه متغيرات الاستجابة. يمكن الاطلاع على مزيد من التفاصيل في المراجع. [ 29 ]

بايزي

مقارنة الدالة اللوجستية بدالة بروبيت معكوسة مُقاسة (أي دالة التوزيع التراكمي للتوزيع الطبيعي )، مقارنةσ(x){\displaystyle \sigma (x)}مقابلΦ(π8x){\textstyle \Phi ({\sqrt {\frac {\pi }{8}}}x)}مما يجعل المنحدرات متساوية عند نقطة الأصل. وهذا يُظهر الذيول الأثقل للتوزيع اللوجستي.

في سياق الإحصاء البايزي ، تُوضع التوزيعات الاحتمالية المسبقة عادةً على معاملات الانحدار، على سبيل المثال في صورة توزيعات غاوسية . لا يوجد توزيع احتمالي مسبق مترافق لدالة الاحتمال في الانحدار اللوجستي. عندما كان الاستدلال البايزي يُجرى تحليليًا، كان هذا يجعل حساب التوزيع الاحتمالي اللاحق صعبًا إلا في الأبعاد المنخفضة جدًا. أما الآن، فتتيح البرامج الآلية مثل OpenBUGS و JAGS و PyMC و Stan و Turing.jl حساب هذه التوزيعات اللاحقة باستخدام المحاكاة، لذا فإن عدم وجود التوزيع الاحتمالي المترافق لا يُشكل مشكلة. مع ذلك، عندما يكون حجم العينة أو عدد المعلمات كبيرًا، قد تكون المحاكاة البايزية الكاملة بطيئة، ولذا يلجأ الباحثون غالبًا إلى طرق تقريبية مثل الطرق البايزية التباينية ونشر التوقع .

"قاعدة العشرة"

تُعرف قاعدة " واحد من كل عشرة " على نطاق واسع، وتنص على أن نماذج الانحدار اللوجستي تُعطي قيمًا ثابتة للمتغيرات التفسيرية إذا استندت إلى حد أدنى يبلغ حوالي 10 أحداث لكل متغير تفسيري (EPV)؛ حيث يشير الحدث إلى الحالات التي تنتمي إلى الفئة الأقل تكرارًا في المتغير التابع. وبالتالي، فإن الدراسة المصممة لاستخدامك{\displaystyle k}المتغيرات التفسيرية لحدث (مثل احتشاء عضلة القلب ) المتوقع حدوثه بنسبة معينةص{\displaystyle p}سيحتاج عدد المشاركين في الدراسة إلى ما مجموعه10ك/ص{\displaystyle 10k/p}المشاركون. ومع ذلك، ثمة جدل كبير حول موثوقية هذه القاعدة، التي تستند إلى دراسات محاكاة وتفتقر إلى أساس نظري متين. [ 30 ] ووفقًا لبعض المؤلفين [ 31 فإن القاعدة متحفظة للغاية في بعض الحالات، حيث ذكروا: "إذا اعتبرنا (بشكل شخصي إلى حد ما) أن تغطية فاصل الثقة أقل من 93%، أو أن الخطأ من النوع الأول أكبر من 7%، أو أن التحيز النسبي أكبر من 15%، يمثل مشكلة، فإن نتائجنا تشير إلى أن المشاكل شائعة إلى حد ما مع 2-4 EPV، وغير شائعة مع 5-9 EPV، ولا تزال تُلاحظ مع 10-16 EPV. لم تكن أسوأ حالات كل مشكلة حادة مع 5-9 EPV، وعادةً ما تكون مماثلة لتلك التي تحدث مع 10-16 EPV". [ 32 ]

وقد توصل آخرون إلى نتائج لا تتفق مع ما سبق، باستخدام معايير مختلفة. ومن المعايير المفيدة ما إذا كان من المتوقع أن يحقق النموذج المُطابق نفس القدرة على التمييز التنبؤي في عينة جديدة كما حققها في عينة تطوير النموذج. ولتحقيق هذا المعيار، قد يلزم 20 حدثًا لكل متغير مرشح. [ 33 ] كما يمكن القول إن 96 مشاهدة كافية لتقدير قيمة التقاطع للنموذج بدقة كافية بحيث يكون هامش الخطأ في الاحتمالات المتوقعة ±0.1 عند مستوى ثقة 0.95. [ 13 ]

خطأ وأهمية المطابقة

اختبار الانحراف ونسبة الاحتمال - حالة بسيطة

في أي عملية ملاءمة، تُحسّن إضافة مُعامل ملاءمة آخر إلى النموذج (مثل مُعاملات بيتا في نموذج الانحدار اللوجستي) قدرة النموذج على التنبؤ بالنتائج المقاسة في أغلب الأحيان. ويصدق هذا حتى لو لم يكن للمُعامل الإضافي أي قيمة تنبؤية، لأن النموذج سيُفرط في ملاءمة البيانات. والسؤال المطروح هو: هل التحسن الناتج عن إضافة مُعامل ملاءمة آخر كافٍ للتوصية بإدراجه، أم أنه مجرد تحسن متوقع من فرط الملاءمة؟

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

تتشابه نماذج الانحدار الخطي والانحدار اللوجستي في العديد من الجوانب. فعلى سبيل المثال، في الانحدار الخطي البسيط، يتم مطابقة مجموعة من K نقطة بيانات ( x<sub> k</sub> , y<sub> k</sub> ) مع دالة نموذجية مقترحة على النحو التالي:y=ب0+ب1x{\displaystyle y=b_{0}+b_{1}x}يتم الحصول على المطابقة عن طريق اختيار معلمات b التي تقلل من مجموع مربعات البواقي (حد الخطأ التربيعي) لكل نقطة بيانات:

ε2=ك=1ك(ب0+ب1xك-yك)2.{\displaystyle \varepsilon ^{2}=\sum _{k=1}^{K}(b_{0}+b_{1}x_{k}-y_{k})^{2}.}

سيتم الإشارة إلى القيمة الدنيا التي تمثل المطابقة بـε^2{\displaystyle {\hat {\varepsilon }}^{2}}

يمكن تقديم فكرة النموذج الصفري ، حيث يُفترض أن المتغير x ليس له فائدة في التنبؤ بنتائج y k : يتم مطابقة نقاط البيانات مع دالة نموذج صفري على شكل y  = b 0 مع حد خطأ تربيعي: 

ε2=ك=1ك(ب0-yك)2.{\displaystyle \varepsilon ^{2}=\sum _{k=1}^{K}(b_{0}-y_{k})^{2}.}

تتضمن عملية التوفيق اختيار قيمة b 0 التي تقلل منε2{\displaystyle \varepsilon ^{2}}من مدى ملاءمة النموذج الصفري، المشار إليه بـ εφ2{\displaystyle \varepsilon _{\varphi }^{2}}حيثφ{\displaystyle \varphi }يشير الرمز السفلي إلى النموذج الصفري. ويُلاحظ أن النموذج الصفري مُحسَّن بواسطةب0=y¯{\displaystyle b_{0}={\overline {y}}}أينy¯{\displaystyle {\overline {y}}}هو متوسط ​​قيم y k ، والقيمة المُحسَّنةεφ2{\displaystyle \varepsilon _{\varphi }^{2}}يكون:

ε^φ2=ك=1ك(y¯-yك)2{\displaystyle {\hat {\varepsilon }}_{\varphi }^{2}=\sum _{k=1}^{K}({\overline {y}}-y_{k})^{2}}

وهو يتناسب مع مربع الانحراف المعياري للعينة (غير المصحح) لنقاط بيانات y k .

يمكننا تخيل حالة يتم فيها تعيين نقاط البيانات y k عشوائيًا لمختلف قيم x k ، ثم يتم تركيبها باستخدام النموذج المقترح. على وجه التحديد، يمكننا النظر في تركيب النموذج المقترح لكل تبديل لنتائج y k . يمكن إثبات أن الخطأ الأمثل لأي من هذه التركيبات لن يكون أبدًا أقل من الخطأ الأمثل للنموذج الصفري، وأن الفرق بين هذين الخطأين الأدنى سيتبع توزيع كاي تربيع ، بدرجات حرية تساوي درجات حرية النموذج المقترح مطروحًا منها درجات حرية النموذج الصفري، والتي ستكون في هذه الحالة2-1=1{\displaystyle 2-1=1}باستخدام اختبار مربع كاي ، يمكننا بعد ذلك تقدير عدد مجموعات y k المبدلة التي ستؤدي إلى خطأ أدنى أقل من أو يساوي الخطأ الأدنى باستخدام y k الأصلي ، وبالتالي يمكننا تقدير مدى أهمية التحسن الناتج عن تضمين المتغير x في النموذج المقترح.

في الانحدار اللوجستي، يُعد مقياس جودة المطابقة دالة الاحتمال L ، أو لوغاريتمها، دالة الاحتمال اللوغاريتمي . دالة الاحتمال L مماثلة لـε2{\displaystyle \varepsilon ^{2}}في حالة الانحدار الخطي، باستثناء أن دالة الاحتمال تُعظّم بدلاً من أن تُصغّر. نرمز إلى دالة الاحتمال اللوغاريتمي المُعظّمة للنموذج المقترح بـ^{\displaystyle {\hat {\ell }}}.

في حالة الانحدار اللوجستي الثنائي البسيط، يتم تركيب مجموعة نقاط البيانات K بشكل احتمالي على دالة بالشكل التالي:

ص(x)=11+هـ-ت{\displaystyle p(x)={\frac {1}{1+e^{-t}}}}

أينص(x){\displaystyle p(x)} هو احتمال أنy=1{\displaystyle y=1}يتم حساب اللوغاريتمات النسبية باستخدام الصيغة التالية:

ت=β0+β1x{\displaystyle t=\beta _{0}+\beta _{1}x}

ودالة الاحتمال اللوغاريتمي هي:

=ك=1ك(yكln(ص(xك))+(1-yك)ln(1-ص(xك))){\displaystyle \ell =\sum _{k=1}^{K}\left(y_{k}\ln(p(x_{k}))+(1-y_{k})\ln(1-p(x_{k}))\right)}

بالنسبة للنموذج الصفري، فإن احتمال أنy=1{\displaystyle y=1}يُعطى بواسطة:

صφ(x)=11+هـ-تφ{\displaystyle p_{\varphi }(x)={\frac {1}{1+e^{-t_{\varphi }}}}}

يتم إعطاء اللوغاريتمات الاحتمالية للنموذج الصفري بواسطة:

تφ=β0{\displaystyle t_{\varphi }=\beta _{0}}

ودالة الاحتمال اللوغاريتمي هي:

φ=ك=1ك(yكln(صφ)+(1-yك)ln(1-صφ)){\displaystyle \ell _{\varphi }=\sum _{k=1}^{K}\left(y_{k}\ln(p_{\varphi })+(1-y_{k})\ln(1-p_{\varphi })\right)}

بما أنناصφ=y¯{\displaystyle p_{\varphi }={\overline {y}}}عند الحد الأقصى لـ L ، تكون قيمة الاحتمالية اللوغاريتمية القصوى للنموذج الصفري هي

^φ=ك(y¯ln(y¯)+(1-y¯)ln(1-y¯)){\displaystyle {\hat {\ell }}_{\varphi }=K(\,{\overline {y}}\ln({\overline {y}})+(1-{\overline {y}})\ln(1-{\overline {y}}))}

الأمثلβ0{\displaystyle \beta _{0}}يكون:

β0=ln(y¯1-y¯){\displaystyle \beta _{0}=\ln \left({\frac {\overline {y}}{1-{\overline {y}}}}\right)}

أينy¯{\displaystyle {\overline {y}}}وهو مرة أخرى متوسط ​​قيم y k . ومرة ​​أخرى، يمكننا أن ننظر نظرياً في مدى ملاءمة النموذج المقترح لكل تبديل لـ y ويمكن إثبات أن الحد الأقصى لقيمة الاحتمال اللوغاريتمي لهذه التبديلات لن يكون أبداً أصغر من قيمة النموذج الصفري.

^^φ{\displaystyle {\hat {\ell }}\geq {\hat {\ell }}_{\varphi }}

كذلك، على غرار خطأ حالة الانحدار الخطي، يمكننا تعريف انحراف نموذج الانحدار اللوجستي على النحو التالي:

د=ln(ل^2ل^φ2)=2(^-^φ){\displaystyle D=\ln \left({\frac {{\hat {L}}^{2}}{{\hat {L}}_{\varphi }^{2}}}\right)=2({\hat {\ell }}-{\hat {\ell }}_{\varphi })}

والتي ستكون دائمًا موجبة أو صفرًا. والسبب في هذا الاختيار هو أن الانحراف ليس فقط مقياسًا جيدًا لجودة المطابقة، بل إنه أيضًا يتبع تقريبًا توزيع كاي تربيع، مع تحسن التقريب بزيادة عدد نقاط البيانات ( K )، ليصبح توزيع كاي تربيع تمامًا في حالة عدد لا نهائي من نقاط البيانات. وكما هو الحال في الانحدار الخطي، يمكننا استخدام هذه الحقيقة لتقدير احتمالية أن تعطي مجموعة عشوائية من نقاط البيانات مطابقة أفضل من المطابقة التي تم الحصول عليها بواسطة النموذج المقترح، وبالتالي الحصول على تقدير لمدى تحسن النموذج بشكل ملحوظ من خلال تضمين xk نقطة بيانات في النموذج المقترح.

بالنسبة للنموذج البسيط لدرجات اختبار الطلاب الموصوف أعلاه، فإن القيمة القصوى للوغاريتم الاحتمالية للنموذج الصفري هي^φ=-13.8629...{\displaystyle {\hat {\ell }}_{\varphi }=-13.8629\ldots }القيمة القصوى للوغاريتم الاحتمالية للنموذج البسيط هي^=-8.02988...{\displaystyle {\hat {\ell }}=-8.02988\ldots }بحيث يكون الانحرافد=2(^-^φ)=11.6661...{\displaystyle D=2({\hat {\ell }}-{\hat {\ell }}_{\varphi })=11.6661\ldots }

باستخدام اختبار مربع كاي للدلالة، فإن تكامل توزيع مربع كاي بدرجة حرية واحدة من 11.6661... إلى ما لا نهاية يساوي 0.00063649...

هذا يعني فعليًا أنه من المتوقع أن يكون حوالي 6 من كل 10000 عملية مطابقة لقيمة y k عشوائية أكثر ملاءمة (انحرافًا أصغر) من قيمة y k المعطاة ، وبالتالي يمكننا استنتاج أن تضمين المتغير x والبيانات في النموذج المقترح يُعد تحسينًا كبيرًا جدًا مقارنةً بالنموذج الصفري. بعبارة أخرى، نرفض الفرضية الصفرية .1-د99.94%{\displaystyle 1-D\approx 99.94\%}ثقة.

ملخص مدى ملاءمة المنتج

يُقاس مدى جودة المطابقة في نماذج الانحدار الخطي عمومًا باستخدام . ولأن هذا لا يوجد له نظير مباشر في الانحدار اللوجستي، يمكن استخدام طرق مختلفة [ 34 ] : الفصل 21، بما في ذلك ما يلي.

اختبارات الانحراف ونسبة الاحتمال

في تحليل الانحدار الخطي، يُعنى بتقسيم التباين عبر حساب مجموع المربعات ، حيث يُقسم التباين في المعيار أساسًا إلى تباين تُفسره المتغيرات التنبؤية وتباين متبقٍ. أما في تحليل الانحدار اللوجستي، فيُستخدم الانحراف بدلًا من حساب مجموع المربعات. [ 35 ] يُشابه الانحراف حساب مجموع المربعات في الانحدار الخطي [ 2 ] ، وهو مقياس لمدى عدم تطابق نموذج الانحدار اللوجستي مع البيانات. [ 35 ] عند توفر نموذج "مُشبع" (نموذج ذو تطابق مثالي نظريًا)، يُحسب الانحراف بمقارنة النموذج المُعطى بالنموذج المُشبع. [ 2 ] تُعطي هذه العملية الحسابية اختبار نسبة الاحتمال : [ 2 ]

د=-2lnاحتمالية النموذج المُطابقاحتمالية النموذج المشبع.{\displaystyle D=-2\ln {\frac {\text{likelihood of the fitted model}}{\text{likelihood of the saturated model}}}.}

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

عندما لا يكون النموذج المشبع متاحًا (وهي حالة شائعة)، يتم حساب الانحراف ببساطة على أنه −2·(log احتمالية النموذج الملائم)، ويمكن إزالة الإشارة إلى log احتمالية النموذج المشبع من كل ما يلي دون ضرر.

يُعدّ مقياسان للانحراف في الانحدار اللوجستي مهمين للغاية: الانحراف الصفري وانحراف النموذج. يُمثّل الانحراف الصفري الفرق بين نموذج يحتوي على الحد الثابت فقط (أي "بدون مُتنبئات") والنموذج المُشبع. بينما يُمثّل انحراف النموذج الفرق بين نموذج يحتوي على مُتنبئ واحد على الأقل والنموذج المُشبع. [ 35 ] في هذا السياق، يُوفّر النموذج الصفري أساسًا لمقارنة نماذج المُتنبئات. وبما أن الانحراف هو مقياس للفرق بين نموذج مُعيّن والنموذج المُشبع، فإن القيم الأصغر تُشير إلى مُلاءمة أفضل. لذا، لتقييم مُساهمة مُتنبئ أو مجموعة من المُتنبئات، يُمكن طرح انحراف النموذج من الانحراف الصفري وتقييم الفرق بناءً علىχs-ص2،{\displaystyle \chi _{s-p}^{2},} توزيع مربع كاي بدرجات حرية [ 2 ] تساوي الفرق في عدد المعلمات المقدرة.

يترك

دباطل=-2lnاحتمالية النموذج الصفرياحتمالية النموذج المشبعدمُجهز=-2lnاحتمالية النموذج المُلائماحتمالية النموذج المشبع.{\displaystyle {\begin{aligned}D_{\text{null}}&=-2\ln {\frac {\text{likelihood of null model}}{\text{likelihood of the saturated model}}}\\[6pt]D_{\text{fitted}}&=-2\ln {\frac {\text{likelihood of fitted model}}{\text{likelihood of the saturated model}}}.\end{aligned}}}

إذن، الفرق بينهما هو:

دباطل-دمُجهز=-2(lnاحتمالية النموذج الصفرياحتمالية النموذج المشبع-lnاحتمالية النموذج المُلائماحتمالية النموذج المشبع)=-2ln(احتمالية النموذج الصفرياحتمالية النموذج المشبع)(احتمالية النموذج المُلائماحتمالية النموذج المشبع)=-2lnاحتمالية النموذج الصفرياحتمالية النموذج المُلائم.{\displaystyle {\begin{aligned}D_{\text{null}}-D_{\text{fitted}}&=-2\left(\ln {\frac {\text{likelihood of null model}}{\text{likelihood of the saturated model}}}-\ln {\frac {\text{likelihood of fitted model}}{\text{likelihood of the saturated model}}}\right)\\[6pt]&=-2\ln {\frac {\left({\dfrac {\text{likelihood of null model}}{\text{likelihood of the saturated model}}}\right)}{\left({\dfrac {\text{likelihood of fitted model}}{\text{likelihood of the saturated model}}}\right)}}\\[6pt]&=-2\ln {\frac {\text{likelihood of the null model}}{\text{likelihood of fitted model}}}.\end{aligned}}}

إذا كان انحراف النموذج أصغر بكثير من الانحراف الصفري، فيمكن الاستنتاج أن المتغير التنبؤي أو مجموعة المتغيرات التنبؤية تُحسّن بشكل ملحوظ من ملاءمة النموذج. وهذا يُشابه اختبار F المستخدم في تحليل الانحدار الخطي لتقييم دلالة التنبؤ. [ 35 ]

معامل التحديد الزائف (R²)

في الانحدار الخطي، يُستخدم معامل الارتباط المتعدد التربيعي (R² ) لتقييم جودة المطابقة، إذ يُمثل نسبة التباين في المعيار الذي تُفسره المتغيرات التنبؤية. [ 35 ] أما في تحليل الانحدار اللوجستي، فلا يوجد مقياس مماثل مُتفق عليه، ولكن توجد عدة مقاييس مُتنافسة، لكل منها قيودها. [ 35 ] [ 36 ]

تتناول هذه الصفحة أربعة من أكثر المؤشرات استخداماً ومؤشراً واحداً أقل استخداماً:

  • نسبة الاحتمال R 2 L
  • كوكس وسنيل آر 2 سي إس
  • ناغلكيرك آر 2 شمال
  • مكفادين آر 2 مكف
  • تيور آر 2 تي

اختبار هوسمر-ليميشو

يستخدم اختبار هوسمر -ليميشو إحصائية اختبار تتبع بشكل تقاربيχ2{\displaystyle \chi ^{2}}يُستخدم التوزيع لتقييم ما إذا كانت معدلات الأحداث المرصودة تتطابق مع معدلات الأحداث المتوقعة في مجموعات فرعية من مجتمع النموذج. ويعتبر بعض الإحصائيين هذا الاختبار قديمًا نظرًا لاعتماده على تصنيف عشوائي للاحتمالات المتوقعة وقوته النسبية المنخفضة. [ 37 ]

أهمية المعامل

بعد ملاءمة النموذج، من المرجح أن يرغب الباحثون في دراسة مساهمة كل متغير تنبؤي على حدة. وللقيام بذلك، سيرغبون في فحص معاملات الانحدار. في الانحدار الخطي، تمثل معاملات الانحدار التغير في المعيار لكل وحدة تغير في المتغير التنبؤي. [ 35 ] أما في الانحدار اللوجستي، فتمثل معاملات الانحدار التغير في اللوجيت لكل وحدة تغير في المتغير التنبؤي. ونظرًا لأن اللوجيت ليس بديهيًا، فمن المرجح أن يركز الباحثون على تأثير المتغير التنبؤي على الدالة الأسية لمعامل الانحدار - نسبة الاحتمالات (انظر التعريف ). في الانحدار الخطي، يتم تقييم دلالة معامل الانحدار عن طريق حساب اختبار t . أما في الانحدار اللوجستي، فهناك عدة اختبارات مختلفة مصممة لتقييم دلالة كل متغير تنبؤي على حدة، وأبرزها اختبار نسبة الاحتمال وإحصائية والد.

اختبار نسبة الاحتمال

يُعد اختبار نسبة الاحتمال المذكور أعلاه لتقييم مدى ملاءمة النموذج الإجراء الموصى به أيضًا لتقييم مساهمة كل "متغير تنبؤي" في نموذج معين. [ 2 ] [ 26 ] [ 35 ] في حالة نموذج متغير تنبؤي واحد، تتم مقارنة انحراف نموذج المتغير التنبؤي بانحراف النموذج الصفري على توزيع كاي تربيع بدرجة حرية واحدة. إذا كان انحراف نموذج المتغير التنبؤي أصغر بشكل ملحوظ (مقارنةً بتوزيع كاي تربيع باستخدام الفرق في درجات حرية النموذجين)، فيمكن الاستنتاج بوجود ارتباط ذي دلالة إحصائية بين المتغير التنبؤي والنتيجة. على الرغم من أن بعض الحزم الإحصائية الشائعة (مثل SPSS) توفر إحصائيات اختبار نسبة الاحتمال، إلا أنه بدون هذا الاختبار الذي يتطلب موارد حاسوبية كبيرة، سيكون من الصعب تقييم مساهمة كل متغير تنبؤي في حالة الانحدار اللوجستي المتعدد. لتقييم مساهمة كل متغير تنبؤي، يمكن إدخال المتغيرات التنبؤية بشكل هرمي، ومقارنة كل نموذج جديد بالنموذج السابق لتحديد مساهمة كل متغير. [ 35 ] ثمة جدل بين الإحصائيين حول مدى ملاءمة ما يُسمى بالإجراءات "التدريجية". ويكمن الخوف في أنها قد لا تحافظ على الخصائص الإحصائية الاسمية، وقد تُصبح مُضللة. [ 38 ]

إحصائية والد

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

دبليوج=βج2Sهـβج2{\displaystyle W_{j}={\frac {\beta _{j}^{2}}{SE_{\beta _{j}}^{2}}}}

على الرغم من أن العديد من البرامج الإحصائية (مثل SPSS وSAS) تُبلغ عن إحصائية والد لتقييم مساهمة كل مُتنبئ على حدة، إلا أن هذه الإحصائية لها بعض القيود. فعندما يكون معامل الانحدار كبيرًا، يميل الخطأ المعياري لهذا المعامل إلى أن يكون أكبر أيضًا، مما يزيد من احتمالية الخطأ من النوع الثاني . كما تميل إحصائية والد إلى التحيز عندما تكون البيانات قليلة. [ 35 ]

أخذ العينات باستخدام أسلوب دراسة الحالات والشواهد

لنفترض أن الحالات نادرة. عندها قد نرغب في أخذ عينات منها بوتيرة أعلى من معدل انتشارها في المجتمع. على سبيل المثال، لنفترض وجود مرض يصيب شخصًا واحدًا من بين كل 10,000، ولجمع بياناتنا نحتاج إلى إجراء فحص بدني شامل. قد يكون إجراء آلاف الفحوصات البدنية لأشخاص أصحاء للحصول على بيانات لعدد قليل فقط من المصابين مكلفًا للغاية. لذا، قد نقوم بتقييم عدد أكبر من المصابين، وربما جميع الحالات النادرة. يُعرف هذا أيضًا بأخذ العينات بأثر رجعي، أو ما يُسمى بالبيانات غير المتوازنة. وكقاعدة عامة، فإن أخذ عينات من الضوابط بمعدل خمسة أضعاف عدد الحالات سيوفر بيانات كافية للضوابط. [ 39 ]

يتميز الانحدار اللوجستي بكونه فريدًا من نوعه، إذ يمكن تقديره على بيانات غير متوازنة، بدلًا من بيانات مأخوذة عشوائيًا، ومع ذلك يُعطي تقديرات صحيحة لمعاملات تأثير كل متغير مستقل على النتيجة. بمعنى آخر، إذا قمنا بتكوين نموذج لوجستي من هذه البيانات، وإذا كان النموذج صحيحًا في عموم السكان، فإنβج{\displaystyle \beta _{j}}جميع المعاملات صحيحة باستثناءβ0{\displaystyle \beta _{0}}يمكننا تصحيح ذلكβ0{\displaystyle \beta _{0}}إذا عرفنا الانتشار الحقيقي على النحو التالي: [ 39 ]

β^0*=β^0+سجلπ1-π-سجلπ~1-π~{\displaystyle {\widehat {\beta }}_{0}^{*}={\widehat {\beta }}_{0}+\log {\frac {\pi }{1-\pi }}-\log {{\tilde {\pi }} \over {1-{\tilde {\pi }}}}}

أينπ{\displaystyle \pi }هو الانتشار الحقيقي وπ~{\displaystyle {\tilde {\pi }}}يمثل هذا الانتشار في العينة.

مناقشة

على غرار أنواع تحليل الانحدار الأخرى ، يستخدم الانحدار اللوجستي متغيرًا تنبؤيًا واحدًا أو أكثر، قد يكون متصلًا أو فئويًا. ولكن على عكس الانحدار الخطي العادي، يُستخدم الانحدار اللوجستي للتنبؤ بمتغيرات تابعة تنتمي إلى عدد محدود من الفئات (حيث يُعامل المتغير التابع في حالة الانحدار اللوجستي الثنائي كنتيجة لتجربة برنولي )، بدلًا من كونه نتيجة متصلة. ونظرًا لهذا الاختلاف، تُنتهك افتراضات الانحدار الخطي. وعلى وجه الخصوص، لا يمكن أن تتبع البواقي التوزيع الطبيعي. إضافةً إلى ذلك، قد يُعطي الانحدار الخطي تنبؤات غير منطقية لمتغير تابع ثنائي. لذا، نحتاج إلى طريقة لتحويل المتغير الثنائي إلى متغير متصل يمكن أن يأخذ أي قيمة حقيقية (سالبة أو موجبة). ولتحقيق ذلك، يحسب الانحدار اللوجستي الثنائي أولًا احتمالية وقوع الحدث لمستويات مختلفة من كل متغير مستقل، ثم يأخذ لوغاريتم هذه الاحتمالية لإنشاء معيار متصل كنسخة مُحوّلة من المتغير التابع. لوغاريتم الاحتمالات هو لوجيت الاحتمال، ويُعرَّف اللوجيت على النحو التالي:لوجيتص=lnص1-صل 0<ص<1.{\displaystyle \operatorname {logit} p=\ln {\frac {p}{1-p}}\quad {\text{for }}0<p<1\,.}

على الرغم من أن المتغير التابع في الانحدار اللوجستي هو برنولي، فإن دالة اللوجيت تكون على مقياس غير مقيد. [ 2 ] دالة اللوجيت هي دالة الربط في هذا النوع من النماذج الخطية المعممة، أي لوجيتهـ(Y)=β0+β1x{\displaystyle \operatorname {logit} \operatorname {\mathcal {E}} (Y)=\beta _{0}+\beta _{1}x}

Y هو متغير الاستجابة الموزع وفقًا لتوزيع برنولي و x هو متغير التنبؤ؛ قيم β هي المعلمات الخطية.

ثم يُطبَّق لوغاريتم احتمالية النجاح على المتغيرات التنبؤية. تُحوَّل القيمة المتوقعة للوغاريتم إلى احتمالات متوقعة، باستخدام معكوس اللوغاريتم الطبيعي - الدالة الأسية . وبالتالي، على الرغم من أن المتغير التابع المرصود في الانحدار اللوجستي الثنائي هو متغير ثنائي (0 أو 1)، فإن الانحدار اللوجستي يُقدِّر احتمالات كون المتغير التابع "ناجحًا"، كمتغير مستمر. في بعض التطبيقات، تكون الاحتمالات كافية. وفي تطبيقات أخرى، يلزم تنبؤ محدد بنعم أو لا لتحديد ما إذا كان المتغير التابع "ناجحًا" أم لا؛ ويمكن أن يستند هذا التنبؤ الفئوي إلى احتمالات النجاح المحسوبة، حيث تُترجم الاحتمالات المتوقعة التي تتجاوز قيمة حدية مختارة إلى تنبؤ بالنجاح.

التعلم الآلي ودالة خسارة الإنتروبيا المتقاطعة

في تطبيقات التعلم الآلي حيث يتم استخدام الانحدار اللوجستي للتصنيف الثنائي، فإن MLE يقلل من دالة خسارة الإنتروبيا المتقاطعة .

يُعدّ الانحدار اللوجستي خوارزمية مهمة في مجال التعلّم الآلي . والهدف منه هو نمذجة احتمالية متغير عشوائي.Y{\displaystyle Y}تكون القيمة 0 أو 1 بناءً على البيانات التجريبية. [ 40 ]

لنفترض دالة نموذج خطي معمّم مُعَلمة بواسطةθ{\displaystyle \theta }،

حθ(X)=11+هـ-θتيX=برو(Y=1|X؛θ){\displaystyle h_{\theta }(X)={\frac {1}{1+e^{-\theta ^{T}X}}}=\Pr(Y=1\mid X;\theta )}

لذلك،

برو(Y=0|X؛θ)=1-حθ(X){\displaystyle \Pr(Y=0\mid X;\theta )=1-h_{\theta }(X)}

ومنذ ذلك الحينY{0،1}{\displaystyle Y\in \{0,1\}}، نرى ذلكبرو(y|X؛θ){\displaystyle \Pr(y\mid X;\theta )}يُعطى بواسطةبرو(y|X؛θ)=حθ(X)y(1-حθ(X))(1-y).{\displaystyle \Pr(y\mid X;\theta )=h_{\theta }(X)^{y}(1-h_{\theta }(X))^{(1-y)}.}نقوم الآن بحساب دالة الاحتمال بافتراض أن جميع المشاهدات في العينة موزعة توزيع برنولي بشكل مستقل.

ل(θ|y؛x)=برو(Y|X؛θ)=أنابرو(yأنا|xأنا؛θ)=أناحθ(xأنا)yأنا(1-حθ(xأنا))(1-yأنا){\displaystyle {\begin{aligned}L(\theta \mid y;x)&=\Pr(Y\mid X;\theta )\\&=\prod _{i}\Pr(y_{i}\mid x_{i};\theta )\\&=\prod _{i}h_{\theta }(x_{i})^{y_{i}}(1-h_{\theta }(x_{i}))^{(1-y_{i})}\end{aligned}}}

عادةً، يتم تعظيم دالة الاحتمال اللوغاريتمي.

شمال-1سجلل(θ|y؛x)=شمال-1أنا=1شمالسجلبرو(yأنا|xأنا؛θ){\displaystyle N^{-1}\log L(\theta \mid y;x)=N^{-1}\sum _{i=1}^{N}\log \Pr(y_{i}\mid x_{i};\theta )}

والتي يتم تحقيق أقصى قدر منها باستخدام تقنيات التحسين مثل انحدار التدرج .

بافتراض(x،y){\displaystyle (x,y)}يتم سحب الأزواج بشكل منتظم من التوزيع الأساسي، ثم في حالة N الكبيرة ، 

ليمشمال+شمال-1أنا=1شمالسجلبرو(yأنا|xأنا؛θ)=xXyYبرو(X=x،Y=y)سجلبرو(Y=y|X=x؛θ)=xXyYبرو(X=x،Y=y)(-سجلبرو(Y=y|X=x)برو(Y=y|X=x؛θ)+سجلبرو(Y=y|X=x))=-دكوالالمبور(YYθ)-ح(Y|X){\displaystyle {\begin{aligned}&\lim \limits _{N\rightarrow +\infty }N^{-1}\sum _{i=1}^{N}\log \Pr(y_{i}\mid x_{i};\theta )=\sum _{x\in {\mathcal {X}}}\sum _{y\in {\mathcal {Y}}}\Pr(X=x,Y=y)\log \Pr(Y=y\mid X=x;\theta )\\[6pt]={}&\sum _{x\in {\mathcal {X}}}\sum _{y\in {\mathcal {Y}}}\Pr(X=x,Y=y)\left(-\log {\frac {\Pr(Y=y\mid X=x)}{\Pr(Y=y\mid X=x;\theta )}}+\log \Pr(Y=y\mid X=x)\right)\\[6pt]={}&-D_{\text{KL}}(Y\parallel Y_{\theta })-H(Y\mid X)\end{aligned}}}

أينح(Y|X){\displaystyle H(Y\mid X)}هي الإنتروبيا الشرطية ودكوالالمبور{\displaystyle D_{\text{KL}}}هو تباعد كولباك-لايبير . وهذا يقودنا إلى استنتاج بديهي مفاده أنه من خلال تعظيم دالة الاحتمال اللوغاريتمي لنموذج ما، فإنك تقلل من تباعد كولباك-لايبير لهذا النموذج عن توزيع الإنتروبيا القصوى. وبالتالي، نبحث بشكل بديهي عن النموذج الذي يعتمد على أقل عدد من الافتراضات في معاييره.

مقارنة بالانحدار الخطي

يمكن اعتبار الانحدار اللوجستي حالة خاصة من النموذج الخطي المعمم، وبالتالي فهو مماثل للانحدار الخطي . مع ذلك، يعتمد نموذج الانحدار اللوجستي على افتراضات مختلفة تمامًا (حول العلاقة بين المتغيرات التابعة والمستقلة) عن تلك الخاصة بالانحدار الخطي. على وجه الخصوص، يمكن ملاحظة الاختلافات الرئيسية بين هذين النموذجين في السمتين التاليتين للانحدار اللوجستي. أولًا، التوزيع الشرطيy|x{\displaystyle y\mid x}يتبع التوزيع توزيع برنولي وليس التوزيع الغاوسي ، لأن المتغير التابع ثنائي. ثانيًا، القيم المتوقعة هي احتمالات، وبالتالي فهي محصورة بين (0،1) من خلال دالة التوزيع اللوجستي، لأن الانحدار اللوجستي يتنبأ باحتمالية نتائج معينة وليس بالنتائج نفسها.

البدائل

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

يُعدّ الانحدار اللوجستي بديلاً لطريقة فيشر لعام 1936، وهي تحليل التمييز الخطي . [ 42 ] إذا تحققت افتراضات تحليل التمييز الخطي، يُمكن عكس عملية التكييف لإنتاج الانحدار اللوجستي. إلا أن العكس غير صحيح، لأن الانحدار اللوجستي لا يتطلب افتراض التوزيع الطبيعي متعدد المتغيرات الخاص بتحليل التمييز. [ 43 ]

يمكن بسهولة تخفيف افتراض التأثيرات الخطية للمتنبئ باستخدام تقنيات مثل دوال التجزئة . [ 13 ]

تاريخ

يُقدّم كريمر (2002) تاريخًا مُفصّلًا للانحدار اللوجستي . طُوّرت الدالة اللوجستية كنموذج لنمو السكان ، وأطلق عليها بيير فرانسوا فيرهولست اسم "لوجستي" في ثلاثينيات وأربعينيات القرن التاسع عشر، بتوجيه من أدولف كيتليه ؛ انظر قسم تاريخ الدالة اللوجستية  لمزيد من التفاصيل. [ 44 ] في ورقته البحثية الأولى (1838)، لم يُحدّد فيرهولست كيفية مُلاءمة المنحنيات للبيانات. [ 45 ] [ 46 ] في ورقته البحثية الأكثر تفصيلًا (1845)، حدّد فيرهولست مُعاملات النموذج الثلاثة بجعل المنحنى يمرّ بثلاث نقاط مُرصودة، مما أسفر عن تنبؤات ضعيفة. [ 47 ] [ 48 ]

طُوِّرت الدالة اللوجستية بشكل مستقل في الكيمياء كنموذج للتفاعل التحفيزي الذاتي ( فيلهلم أوستوالد ، 1883). [ 49 ] التفاعل التحفيزي الذاتي هو تفاعل يكون فيه أحد النواتج عاملًا محفزًا للتفاعل نفسه، بينما يكون إمداد أحد المتفاعلات ثابتًا. وهذا يُؤدي بطبيعة الحال إلى ظهور المعادلة اللوجستية لنفس سبب نمو السكان: فالتفاعل يُعزز نفسه بنفسه ولكنه مُقيد.

أُعيد اكتشاف الدالة اللوجستية بشكل مستقل كنموذج لنمو السكان عام 1920 على يد ريموند بيرل ولويل ريد ، ونُشرت دراستهما تحت عنوان " بيرل وريد (1920)" ، مما أدى إلى استخدامها في الإحصاءات الحديثة. لم يكونا على دراية في البداية بعمل فيرهولست، ويُفترض أنهما تعرّفا عليه من ل. غوستاف دو باسكييه ، لكنهما لم يُنصفاه ولم يتبنيا مصطلحاته. [ 50 ] تم الاعتراف بأسبقية فيرهولست، وأعاد أودني يول إحياء مصطلح "لوجستي" عام 1925، والذي لا يزال يُستخدم منذ ذلك الحين. [ 51 ] طبّق بيرل وريد النموذج لأول مرة على سكان الولايات المتحدة، وقاما أيضًا في البداية بمطابقة المنحنى بجعله يمر بثلاث نقاط؛ وكما هو الحال مع فيرهولست، أسفرت هذه الطريقة أيضًا عن نتائج ضعيفة. [ 52 ]

في ثلاثينيات القرن العشرين، طُوِّر نموذج البروبيت ووُضِّح بشكل منهجي على يد تشيستر إيتنر بليس ، الذي صاغ مصطلح "بروبيت" في كتابه (1934) ، وجون غادوم في كتابه (1933) . ثم قام رونالد أ. فيشر بتطبيق النموذج باستخدام تقدير الاحتمال الأقصى في كتابه (1935) ، كإضافة إلى عمل بليس. استُخدم نموذج البروبيت بشكل أساسي في المقايسات الحيوية ، وسبقته أعمال سابقة تعود إلى عام 1860؛ انظر نموذج البروبيت § التاريخ . أثَّر نموذج البروبيت على التطوير اللاحق لنموذج اللوجيت، وتنافس هذان النموذجان فيما بينهما. [ 53 ] 

يُرجّح أن يكون إدوين بيدويل ويلسون وطالبته جين وورسيستر قد استخدما النموذج اللوجستي لأول مرة كبديل لنموذج البروبيت في التحليل البيولوجي في كتاب ويلسون وورسيستر (1943) . [ 54 ] مع ذلك، يُعزى تطوير النموذج اللوجستي كبديل عام لنموذج البروبيت بشكل أساسي إلى عمل جوزيف بيركسون على مدى عقود، بدءًا من كتابه (1944) حيث صاغ مصطلح "لوجيت" قياسًا على "بروبيت"، واستمر ذلك في كتابه (1951) والسنوات اللاحقة. [ 55 ] في البداية، رُفض نموذج اللوجيت باعتباره أقل كفاءة من نموذج البروبيت، لكنه "حقق تدريجيًا مكانة مساوية لنموذج البروبيت"، [ 56 ] لا سيما بين عامي 1960 و1970. وبحلول عام 1970، أصبح نموذج اللوجيت مكافئًا لنموذج البروبيت المستخدم في المجلات الإحصائية، ثم تفوق عليه بعد ذلك. يعود هذا الانتشار النسبي إلى اعتماد نموذج اللوجيت خارج نطاق التحليل البيولوجي، بدلاً من استبدال نموذج البروبيت داخل التحليل البيولوجي، واستخدامه غير الرسمي في الممارسة العملية؛ ويعزى انتشار نموذج اللوجيت إلى بساطته الحسابية وخصائصه الرياضية وعموميته، مما يسمح باستخدامه في مجالات متنوعة. [ 3 ]

وقد طرأت تحسينات مختلفة خلال تلك الفترة، ولا سيما من قبل ديفيد كوكس ، كما في كوكس (1958) . [ 4 ]

طُرح نموذج اللوجيت متعدد الحدود بشكل مستقل في دراسات كوكس (1966) وثيل (1969) ، مما وسّع نطاق تطبيقه وزاد من شعبيته بشكل كبير. [ 57 ] وفي عام 1973، ربط دانيال مكفادين نموذج اللوجيت متعدد الحدود بنظرية الاختيار المنفصل ، وتحديدًا بديهية لوس للاختيار ، موضحًا أن نموذج اللوجيت متعدد الحدود ينبع من افتراض استقلال البدائل غير ذات الصلة ، ومفسرًا احتمالات البدائل على أنها تفضيلات نسبية؛ [ 58 ] وقد شكّل هذا أساسًا نظريًا للانحدار اللوجستي. [ 57 ]

الإضافات

هناك عدد كبير من الامتدادات:

انظر أيضاً

مراجع

  1. تولز، جوليانا؛ ميورر، ويليام جيه (2016). "الانحدار اللوجستي الذي يربط خصائص المريض بالنتائج". مجلة الجمعية الطبية الأمريكية . 316 (5): 533-534 . doi : 10.1001/jama.2016.7653 . ISSN 0098-7484 . OCLC 6823603312. PMID 27483067 .   
  2. 1 2 3 4 5 6 7 8 9 10 11 هوسمر، ديفيد دبليو؛ ليميشو، ستانلي (2000). الانحدار اللوجستي التطبيقي ( الطبعة الثانية). وايلي. ISBN  978-0-471-35632-5.
  3. 1 2 Cramer 2002 ، ص. 10–11.
  4. 1 2 ووكر، إس إتش؛ دنكان، دي بي (1967). "تقدير احتمالية وقوع حدث ما كدالة لعدة متغيرات مستقلة". Biometrika . 54 (1/2): 167–178 . doi : 10.2307/2333860 . JSTOR 2333860 . 
  5. كريمر 2002 ، ص 8.
  6. بويد، سي آر؛ تولسون، إم إيه؛ كوبس، دبليو إس (1987). "تقييم رعاية الإصابات: طريقة TRISS. درجة الصدمة ودرجة شدة الإصابة" . مجلة الصدمات . 27 (4): 370-378 . doi : 10.1097/00005373-198704000-00005 . PMID 3106646 . 
  7. كولوغلو، م.؛ إلكر، د.؛ ألتون، هـ.؛ سايك، إ. (2001). "التحقق من صحة مؤشر التروية العضلية (MPI) ومؤشر التهاب الصفاق المناعي (PIA II) في مجموعتين مختلفتين من المرضى المصابين بالتهاب الصفاق الثانوي". مجلة أمراض الكبد والجهاز الهضمي . 48 (37): 147-151 . PMID 11268952 . 
  8. ^ بيوندو، إس. راموس، إي. ديروس، م.؛ راغي، جي إم؛ دي أوكا، J .؛ مورينو، ب. فران، ل.؛ جوريتا، إي. (2000). “العوامل النذير للوفيات في التهاب الصفاق القولوني الأيسر: نظام تسجيل جديد”. مجلة الكلية الأمريكية للجراحين . 191 (6): 635– 42. دوى : 10.1016/S1072-7515(00)00758-4 . بميد 11129812 . 
  9. مارشال، جيه سي؛ كوك، دي جيه؛ كريستو، إن في؛ برنارد، جي آر؛ سبرونغ، سي إل؛ سيبالد، دبليو جيه (1995). "مؤشر اختلال وظائف الأعضاء المتعددة: وصف موثوق لنتائج سريرية معقدة". طب العناية المركزة . 23 (10): 1638-1652 . doi : 10.1097/00003246-199510000-00007 . PMID 7587228 . 
  10. لو غال، الابن؛ ليميشو، س.؛ سولنييه، ف. (1993). "مؤشر جديد مبسط لتقييم الحالة الفسيولوجية الحادة (SAPS II) بناءً على دراسة متعددة المراكز أوروبية/أمريكية شمالية". مجلة الجمعية الطبية الأمريكية . 270 (24): 2957-2963 . Bibcode : 1993JAMA..270.2957L . doi : 10.1001/jama.1993.03510240069035 . PMID 8254858 . 
  11. 1 2 ديفيد أ. فريدمان (2009). النماذج الإحصائية: النظرية والتطبيق . مطبعة جامعة كامبريدج . ص 128. 
  12. ترويت، ج؛ كورنفيلد، ج؛ كانيل، و (1967). "تحليل متعدد المتغيرات لمخاطر الإصابة بأمراض القلب التاجية في فرامنغهام". مجلة الأمراض المزمنة . 20 (7): 511-24 . doi : 10.1016/0021-9681(67)90082-3 . PMID 6028270 . 
  13. 1 2 3 هاريل، فرانك إي. (2015). استراتيجيات نمذجة الانحدار . سلسلة سبرينغر في الإحصاء ( الطبعة الثانية). نيويورك؛ سبرينغر. doi : 10.1007/978-3-319-19425-7 . ISBN  978-3-319-19424-0.
  14. م. سترانو؛ ب. م. كولوسيمو (2006). "تحليل الانحدار اللوجستي لتحديد مخططات حد التشكيل تجريبياً". المجلة الدولية لأدوات وآلات التصنيع . 46 (6): 673-682 . doi : 10.1016/j.ijmachtools.2005.07.005 .
  15. بالي، إس كيه؛ داس، إس كيه (2009). "نموذج الانحدار اللوجستي للتنبؤ بمخاطر انهيار السقف في أعمال الحفر والأعمدة في مناجم الفحم: منهجية". علوم السلامة . 47 : 88-96 . doi : 10.1016/j.ssci.2008.01.002 .
  16. بيري، مايكل جيه إيه (1997). تقنيات استخراج البيانات للتسويق والمبيعات ودعم العملاء . وايلي. ص 10. 
  17. ^ ميسا أرانجو، رودريجو؛ حسن، صامويل؛ أوكوسوري، ساتيش ف.؛ موراي تويت ، باميلا (فبراير 2013). "نموذج على مستوى الأسرة لاختيار نوع وجهة الإخلاء أثناء الإعصار باستخدام بيانات إعصار إيفان" . مراجعة المخاطر الطبيعية . 14 (1): 11– 20. بيب كود : 2013NHRev..14...11M . دوى : 10.1061/(ASCE)NH.1527-6996.0000083 . ردمك 1527-6988 . 
  18. ويبنماير، ماثيو جيه؛ هاند، مايكل إس؛ كالكين، ديفيد إي؛ فين، تايرون جيه؛ طومسون، ماثيو بي (يونيو 2013). "تفضيلات المخاطرة في اتخاذ القرارات الاستراتيجية المتعلقة بحرائق الغابات: تجربة اختيار مع مديري حرائق الغابات في الولايات المتحدة" . تحليل المخاطر . 33 (6): 1021-1037 . Bibcode : 2013RiskA..33.1021W . doi : 10.1111/j.1539-6924.2012.01894.x . ISSN 0272-4332 . PMID 23078036. S2CID 45282555 .   
  19. لوفريجليو، روجيرو؛ بوري، دينو؛ ديلوليو، لويجي؛ إيباس، أنجيل (2014-02-01). "نموذج اختيار منفصل قائم على المنافع العشوائية لاختيار المخرج في عمليات الإخلاء الطارئة" . علوم السلامة . 62 : 418-426 . doi : 10.1016/j.ssci.2013.10.004 . ISSN 0925-7535 . 
  20. "الانحدار اللوجستي" . CORP-MIDS1 (MDS) . تم الاسترجاع في 16-03-2024 .
  21. نيمان، جبيرسون، إي إس (1933)، "حول مشكلة الاختبارات الأكثر كفاءة للفرضيات الإحصائية" (ملف PDF) ، المعاملات الفلسفية للجمعية الملكية في لندن أ ، 231 ( 694-706 ): 289-337 ، Bibcode : 1933RSPTA.231..289N ، doi : 10.1098/rsta.1933.0009 ، JSTOR 91247 
  22. "كيفية تفسير نسبة الاحتمالات في الانحدار اللوجستي؟" . معهد البحوث والتعليم الرقمي.
  23. إيفريت، برايان (1998). قاموس كامبريدج للإحصاء . كامبريدج، المملكة المتحدة: نيويورك: مطبعة جامعة كامبريدج. ISBN 978-0-521-59346-5.
  24. على سبيل المثال، يمكن تعريف دالة المؤشر في هذه الحالة على النحو التاليΔ(ن،y)=1-(y-ن)2{\displaystyle \Delta (n,y)=1-(y-n)^{2}}
  25. مالوف، روبرت (2002). "مقارنة بين خوارزميات تقدير معلمات الإنتروبيا القصوى" . وقائع المؤتمر السادس حول تعلم اللغة الطبيعية (CoNLL-2002) . الصفحات 49-55 . doi : 10.3115/1118853.1118871 . 
  26. 1 2 3 4 5 6 7 مينارد، سكوت و. (2002). الانحدار اللوجستي التطبيقي ( الطبعة الثانية). سيج. ISBN  978-0-7619-2208-7.
  27. غورييرو، كريستيان؛ مونفور، آلان (1981). "الخصائص التقاربية لمُقدِّر الاحتمال الأقصى في نماذج اللوجيت الثنائية". مجلة الاقتصاد القياسي . 17 (1): 83-97 . doi : 10.1016/0304-4076(81)90060-9 .
  28. بارك، بيونغ يو؛ سيمار، ليوبولد؛ زيلينيوك، فالنتين (2017). "التقدير غير البارامتري لنماذج الاختيار المنفصل الديناميكي لبيانات السلاسل الزمنية" (ملف PDF) . الإحصاءات الحاسوبية وتحليل البيانات . 108 : 97-120 . doi : 10.1016/j.csda.2016.10.024 .
  29. مورفي، كيفن ب. (2012). التعلم الآلي - منظور احتمالي . مطبعة معهد ماساتشوستس للتكنولوجيا. ص 245. ISBN  978-0-262-01802-9.
  30. فان سميدن، م.؛ دي غروت، ج. أ.؛ مونز، ك. ج.؛ كولينز، ج. س.؛ ألتمان، د. ج.؛ إيكمانز، م. ج.؛ ريتسما، ج. ب. (2016). "لا يوجد أساس منطقي لمعيار متغير واحد لكل 10 أحداث في تحليل الانحدار اللوجستي الثنائي" . مجلة BMC لمنهجية البحوث الطبية . 16 (1): 163. doi : 10.1186/s12874-016-0267-3 . PMC 5122171. PMID 27881078 .  
  31. بيدوزي، ب؛ كونكاتو، ج؛ كيمبر، إ؛ هولفورد، ت ر؛ فاينشتاين، أ ر (ديسمبر 1996). "دراسة محاكاة لعدد الأحداث لكل متغير في تحليل الانحدار اللوجستي" . مجلة علم الأوبئة السريرية . 49 (12): 1373-1379 . doi : 10.1016/s0895-4356(96)00236-3 . PMID 8970487 . 
  32. فيتينغهوف، إي.؛ ماكولوتش، سي إي (12 يناير 2007). "تخفيف قاعدة عشرة أحداث لكل متغير في الانحدار اللوجستي وانحدار كوكس" . المجلة الأمريكية لعلم الأوبئة . 165 (6): 710-718 . doi : 10.1093/aje/kwk052 . PMID 17182981 . 
  33. فان دير بلوغ، تيرد؛ أوستن، بيتر سي؛ ستايربيرغ، إيووت دبليو (2014). "تقنيات النمذجة الحديثة تتطلب كميات هائلة من البيانات: دراسة محاكاة للتنبؤ بالنتائج الثنائية" . مجلة BMC لمنهجية البحوث الطبية . 14 137. doi : 10.1186/1471-2288-14-137 . PMC 4289553. PMID 25532820 .  
  34. غرين، ويليام ن. (2003). التحليل الاقتصادي القياسي ( الطبعة الخامسة). برنتيس هول. ISBN  978-0-13-066189-0.
  35. 1 2 3 4 5 6 7 8 9 10 كوهين، جاكوب؛ كوهين، باتريشيا؛ ويست، ستيفن ج.؛ أيكن، ليونا س. (2002). تحليل الانحدار/الارتباط المتعدد التطبيقي للعلوم السلوكية ( الطبعة الثالثة). روتليدج. ISBN  978-0-8058-2223-6.
  36. أليسون، بول د. "مقاييس ملاءمة الانحدار اللوجستي" (ملف PDF) . شركة ستاتيستيكال هورايزونز وجامعة بنسلفانيا.
  37. هوسمر، د. و. (1997). "مقارنة اختبارات جودة المطابقة لنموذج الانحدار اللوجستي". الإحصاء الطبي . 16 (9): 965-980 . doi : 10.1002/(sici)1097-0258(19970515)16:9 < 965::aid-sim509 > 3.3.co ; 2-f . PMID 9160492 . 
  38. هاريل، فرانك إي. (2010). استراتيجيات نمذجة الانحدار: مع تطبيقات على النماذج الخطية، والانحدار اللوجستي، وتحليل البقاء . نيويورك: سبرينغر. ISBN 978-1-4419-2918-1.
  39. 1 2 https://class.stanford.edu/c4x/HumanitiesScience/StatLearning/asset/classification.pdf الشريحة 16
  40. نغ، أندرو (2000). "ملاحظات محاضرة CS229" (ملف PDF) . ملاحظات محاضرة CS229 : 16-19 .
  41. رودريغيز، ج. (2007). ملاحظات المحاضرات حول النماذج الخطية المعممة . الفصل 3، الصفحة 45. 
  42. غاريث جيمس؛ دانييلا ويتن؛ تريفور هاستي؛ روبرت تيبشيراني (2013). مقدمة في التعلم الإحصائي . سبرينغر. ص 6. 
  43. ^ بوهار، ماجا؛ بلاس، ماتيا؛ ترك ، ساندرا (2004). "مقارنة الانحدار اللوجستي والتحليل التمييزي الخطي: دراسة محاكاة" . ميتودولوسكي زفيزكي . 1 (1).
  44. كريمر 2002 ، ص 3-5.
  45. ^ فيرهولست، بيير فرانسوا (1838). "لاحظ sur la loi que la السكان يصبون في تزايدهم" (PDF) . المراسلات الرياضيات والفيزياء . 10 : 113 – 121 . تم الاسترجاع 3 ديسمبر 2014 .
  46. كريمر 2002 ، ص 4، "لم يذكر كيف قام بتركيب المنحنيات". 
  47. ^ فيرهلست، بيير فرانسوا (1845). "Recherches mathématiques sur la loi d'accroissement de la السكان" [ أبحاث رياضية في قانون زيادة النمو السكاني ] . مذكرات جديدة من الأكاديمية الملكية للعلوم والآداب الجميلة في بروكسل . 18 . تم الاسترجاع 2013/02/18 .
  48. كريمر 2002 ، ص 4.
  49. كريمر 2002 ، ص 7.
  50. كريمر 2002 ، ص. 6.
  51. Cramer 2002 ، ص 6-7.
  52. كريمر 2002 ، ص 5.
  53. Cramer 2002 ، ص 7-9.
  54. كريمر 2002 ، ص 9.
  55. كريمر 2002 ، ص 8، "على حد علمي، فإن إدخال اللوجستيات كبديل لدالة الاحتمال الطبيعي هو عمل شخص واحد، جوزيف بيركسون (1899-1982)، ..." 
  56. كريمر 2002 ، ص 11.
  57. 1 2 Cramer 2002 ، ص. 13.
  58. ماكفادين، دانيال (1973). "تحليل اللوجيت الشرطي لسلوك الاختيار النوعي" (ملف PDF) . في: زاريمبكا، ب. (محرر). آفاق في الاقتصاد القياسي . نيويورك: أكاديميك برس. ص 105-142 . مؤرشف من الأصل (ملف PDF) بتاريخ 27 نوفمبر 2018. تم الاطلاع عليه بتاريخ 20 أبريل 2019 . 

مصادر

  • أغريستي، آلان. (2002). تحليل البيانات الفئوية . نيويورك: وايلي-إنترساينس. ISBN 978-0-471-36093-3.
  • أميميا، تاكيشي (1985). "نماذج الاستجابة النوعية" . الاقتصاد القياسي المتقدم . أكسفورد: باسل بلاكويل. ص 267-359 . ISBN  978-0-631-13345-2.
  • بالاكريشنان، ن. (1991). دليل التوزيع اللوجستي . دار مارسيل ديكر للنشر. رقم ISBN 978-0-8247-8587-1.
  • غوريرو، كريستيان (2000). "الثنائية البسيطة" . الاقتصاد القياسي للمتغيرات التابعة النوعية . نيويورك: مطبعة جامعة كامبريدج. ص 6-37 . ISBN  978-0-521-58985-7.
  • غرين، ويليام هـ. (2003). التحليل الاقتصادي القياسي، الطبعة الخامسة . برنتيس هول. ISBN 978-0-13-066189-0.
  • هيلبي، جوزيف م. (2009). نماذج الانحدار اللوجستي . تشابمان آند هول/سي آر سي برس. رقم ISBN 978-1-4200-7575-5.
  • هوسمر، ديفيد (2013). الانحدار اللوجستي التطبيقي . هوبوكين، نيو جيرسي: وايلي. ISBN 978-0-470-58247-3.
  • هاول، ديفيد سي. (2010). الأساليب الإحصائية لعلم النفس، الطبعة السابعة . بلمونت، كاليفورنيا: تومسون وادزورث. ISBN 978-0-495-59786-5.
  • بيدوزي، ب.؛ ج. كونكاتو؛ إ. كيمبر؛ ت. ر. هولفورد؛ أ. ر. فاينشتاين (1996). "دراسة محاكاة لعدد الأحداث لكل متغير في تحليل الانحدار اللوجستي" . مجلة علم الأوبئة السريرية . 49 (12): 1373-1379 . doi : 10.1016/s0895-4356(96)00236-3 . PMID 8970487 . 
  • بيري، مايكل جيه إيه؛ لينوف، جوردون (1997). تقنيات استخراج البيانات للتسويق والمبيعات ودعم العملاء . وايلي.