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

من أرابيكا، الموسوعة الحرة
اذهب إلى التنقل اذهب إلى البحث
انحدار لوجستي

الانحدار اللوجستي[1] (بالإنجليزية Logistic regression) هو نموذج إحصائي ينتمي لنماذج الانحدار الخطي يمكن من نمذجة متغير ثنائي الحد بدلالة مجموعة من المتغيرات العشوائية المتوقعة، رقمية كانت أو فئوية. يستخدم الانحدار اللوجستي للتنبؤ باحتمالية وقوع حدث ما بمعرفة إضافية لقيم متغيرات يمكن أن تكون مفسرة أو مرتبطة بهذا الحدث.[2][3][4] يستخدم الانحدارُ اللوجستي عدة متغيرات مُتوقَّعة والتي يمكن أن تكون رقمية أو فئوية. يشتهر الانحدار اللوجستي أيضا بتسميات نموذج لوجيت (Logit) أو المصنف العام للأنتروبية. تستعمل هذه النمذجة بشكل واسع في العديد من التطبيقات العلمية والتجارية وهي من طرق النمذجة الأكثر تطبيقا في مجال التعلم الآلي، حيث تصنف ضمن طرق التعلم الآلي المراقب.

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

تاريخ

يعتبر الطبيب والإحصائي جوزيف بيركسون أول من عرف مفهوم الانحدار اللوجستي وطبقه في مجال الإحصاء الحيوي سنة 1944 ونشرت أول ورقة بحثية يذكر فيها الانحدار اللوجستي باسمه في مجلة الجمعية الإحصائية الأمريكية.[5]

مجالات التطبيق

في ما يلي عينة من مجالات تطبيق الانحدار اللوجستي مع الأمثلة:

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

النموذج الرياضي

اصطلاحات

  • نعتبر عينة Ω سعتها n، يشار لكل فرد (ملاحظة إحصائية) في العينة بالرمز ω.
  • Y هو المتغير النوعي المراد تفسيره ويحتمل قيمتين {0,1}. قيمة Y بالنسبة للفرد ω يرمز لها ب Y(w).
  • X=(X1,X2,...,Xp) هي المتغيرات المفسرة ويمكن أن تكون رقمية أو نوعية. قيمة المتغيرات المفسرة بالنسبة للفرد ω يرمز لها ب X(ω)=(X1(ω),X2(ω),...,Xp(ω)).
  • الاحتمالان القبليان لتحقق القيمتين 0 أو 1 يشار إليهما ب P(Y=1) و P(Y=0). وهما يساويان تباعا النسب الملاحظة للقيمتين 0 و 1 في العينة. الاحتمال هنا قبلي بمعنى أنه لا يستلزم معرفة قبلية بقيم X.
  • الاحتمالان الشرطيان لتحقق قيمة معينة ل Y بمعرفة إضافية لقيم X=(X1,X2,...,Xp) يشار إليهما تباعا ب P(Y=1/X) و P(Y=0/X). يشار للاحتمال الشرطي اللازم تقديره بالنسبة للفرد ω ب π(ω).

نموذج لوجيت (Logit)

نموذج الانحدار اللوجستي يحاول تفسير دالة لوجيت المتغير المراد تفسيره بدلالة توليف خطي للمتغيرات المفسرة:

lnP(Y=1/X)1P(Y=1/X)=a0+a1X1+...+apXp، بحيث ai (i من 0 إلى p هي معاملات النموذج التي يجب تقديرها).

وهو ما يمكن من تقدير الاحتمال الشرطي حسب الصيغة التالية:

P(Y=1/X)=ea0+a1X1+...+apXp1+ea0+a1X1+...+apXp

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

تقدير معاملات النموذج يتم حسب طريقة الإمكانية القصوى (Maximum likelihood estimation):

  • المتغير Y موزع حسب توزيع ثنائي الحدين: Yβ(1,π)، مما يمكن من وضع دالة تقدير الاحتمال الأرجح للعينة على الشكل التالي: L=ωπy(1π)(1y)[0,1]. الإمكانية القصوى تشير إلى احتمال تحقق العينة Ω عبر سحب عشوائي.
  • تهدف الطريقة إلى إيجاد متجهة المعاملات A=(a0,a1,...,ap) التي تحقق القيمة القصوى ل L.
  • لتسهيل التعامل مع هذا المؤشر يتم استعمال لوغاريتم تقدير الاحتمال: LL=ωyln(π)+(1y)ln(1π)],0]

المقدر A^ الذي يحقق القيمة القصوى ل L و LL أن يمتلك الخصائص التالية:

في المراجع الرياضية للانحدار اللوجستي، وأيضا في البرامج الإحصائية، يستعمل أيضا ما يعرف بالانحراف الراسب (Residual Deviance): Dm=2LL[0,+[ وباعتبار هذا المؤشر، يكون المقدر A^ هو الذي يحقق القيمة الدنيا ل Dm.

LL دالة محدبة وبذلك تحتمل حلا وحيدا A^. بخلاف نموذج الانحدار الخطي، لا يمكن إيجاد حل للمقدر بطريقة تحليلية: الحلول المتوفرة هي فقط تجريبية (Heuristic) وأهم الخوارزميات المستعملة في إيجاد المقدر A^ هي خوارزمية نيوتن رافسون (Newton-Raphson).

تقدير المعاملات بطريقة نيوتن رافسون

  • يتم تثبيت متجهة المعاملات A في قيمة بدئية A0 (تكون مقتصرة مثلا فقط على قيمة a0=1 للمعامل الثابت و0 لباقي المعاملات).
  • يتم تكرار حساب قيم A حسب الخوارزمية: Ai+1=Ai(2LAA)1×LA بحيث:
    • A=LA هي متجهة الاشتقاق الجزئي من الدرجة الأولى لدالة تقدير الاحتمال، وهي أيضا متجهة تدرج (Gradient) تؤول إلى 0 في جوار القيمة المثلى. تساوي قيم المتجهة: (aj)=ω(yπ)×xj
    • H(a)=2LAA هي متجهة الاشتقاق الجزئي من الدرجة الثانية لدالة تقدير الاحتمال.
    • Ai هي متجهة الحل المؤقت في المرحلة i من الحساب.

هناك ثلاثة معايير لضبط مآل الخوارزمية (وهي وسائط يمكن ضبطها في البرامج الإحصائية):

  • عدد مرات التكرار أو الاستدعاء الذاتي للطريقة.
  • ϵ المستوى الأدنى المقبول لتغير ΔLL، إذا كانت ΔLL<ϵ توقف.
  • ϵ المستوى الأدنى المقبول لتغير ΔA^، إذا كانت ΔA^<ϵ توقف.

H(A) مصفوفة هيسية، حيزها (P+1)×(P+1) وتكمن أهميتها في كون معكوستها تساوي مصفوفة تغاير المعاملات، وهي بذلك حاسمة في مرحلة اختبارات الدلالة الإحصائية للمعاملات. يمكن أيضا كتابتها (باستعمال المصفوفات) على الشكل التالي:

H(A)=(XVX)1 مع H(j1,j2)=ω(xj1×xj2×π×(1π)) و V هي المصفوفة القطرية المشكلة من القيم π(1π).

يمكن التعبير عن خوارزمية تقدير المعاملات عبر الكتابة المصفوفاتية حسب الشكل التالي:

Ai+1=Ai+1+(tXVX)1X(yp).

تقدير النموذج اللوجستي في البرامج الإحصائية

إكسل

يمكن تطبيق تقنية الانحدار اللوجستي في إكسل يدويا (أي بدون استعمال وظائف معرفة قبليا) فقط عبر توليف الأدوات الحسابية المتوفرة واستعمال البرنامج المساعد Solver لإيجاد المقدرات المثلى للنموذج، وذلك عبر اتباع الخطوات التالية:

  • تنظيم جدول البيانات بالتمييز بين المتغير النوعي المراد تفسيره Y والمتغيرات المفسرة X(ω)=(X1(ω),X2(ω),...,Xp(ω)).
  • تثبيت قيم بدئية لمتجهة المعاملات: a0^=(1,0,...,0).
  • بالنسبة لكل فرد إحصائي في جدول البيانات:
    1. حساب المقدار C(x)=a0+a1x1+...+apxp
    2. ثم حساب المقدار π(x)=11+eC(x) الذي يشكل تقديرا للاحتمال الشرطي بدلالة المعلومات المفسرة المتضمنة في المتغيراتX.
    3. حساب المساهمة الفردية في LL (لوغاريتم تقدير الاحتمال): LL(x)=y×ln(π)+(1y)×ln(1π)
  • حساب الانحراف الراسب Dm=2LL عبر جمع المساهمات الفردية.
  • تشغيل خوارزمية استمثال باستعمال البرنامج المساعد Solver لحل البرنامج التالي:{minDmConstrainta^ والذي سيقوم بمحاكاة مجموعة كبيرة من التوليفات لقيم متجهة المعاملات a^ (التي تمثل قيد برنامج الاستمثال) حتى بلوغ أقل قيمة ممكنة لDm.

آر (R)

توفر بيئة آر العديد من المكتبات لتطبيق الانحدار اللوجستي:

  • وظيفة glm المتوفرة في مكتبة stats
  • وظيفة lrm المتوفرة في مكتبة rms
  • وظيفة gm_logistic المتوفرة في مكتبة Rfast

تقييم نموذج الانحدار اللوجستي

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

على العموم، تنقسم طرق التقييم إلى صنفين:

  • طرق التقييم الذاتي أو الداخلي (Internal Evaluation) والمرتبطة بطريقة النموذج في حد ذاتها (الانحدار اللوجستي).
  • طرق التقييم الخارجي (External Evaluation) المميزة لطرق التعلم الآلي المراقب والتي لا علاقة لها بماهية الطريقة المستعملة.

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

التقييم الذاتي

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

  1. حصر المتغيرات المفسرة: هذه المرحلة تفرض إلماما وحدا أدنى من الخبرة في المجال الذي تطبق فيه النمذجة. اقتراح المتغيرات يجب أن يكون مسندا بحد أدنى من التأصيل النظري الذي يسمح بترشيحها كعوامل محتملة لتفسير المتغير Y.
  2. تقدير المعاملات انطلاقا من عينة إحصائية Ω : يشمل التقييم هنا مدى احترام العينة للمنهج العلمي في عملية الاعتيان وهل تم تفادي (أو تقليل) التحيزات المعرفية أثناء عملية السبر، خصوصا في حالة المتغيرات النوعية الواصفة للأراء الشخصية.
  3. تقييم دقة تقدير المعاملات: بوضع مجالات ثقة للمعاملات ومستوى دلالتها الإحصائية.
  4. تقدير قوة النموذج التفسيرية : بمعنى هل يمكن النموذج من توقع مقبول لقيم Y في العينة.
  5. التصديق على وجود حد أدنى من الارتباط بين Y والمتغيرات المفسرة X.
  6. ترتيب المتغيرات المفسرة X=(X1,X2,...,Xp) حسب مساهمتها في تفسير Y.
  7. تقدير تأثير مكونات العينة على نجاعة النموذج : يجب أن يشمل هذا التقييم حساسية النموذج لتغيرات قيم العينة ومدى تأثير القيم القصوى (الغير اعتيادية) على معاملات النموذج وقوته التفسيرية.

التقييم حسب المصنف المرجعي

يندرج الانحدار اللوجستي ضمن طرق التعلم الآلي المراقب، ومن أهم طرق التقييم في هذا المجال أن تتم مقارنة قوة النموذج الناتج مع نموذج مرجعي يسمى المصنف المرجعي أو البديهي (Default Classifier). النموذج المرجعي هو بكل بساطة نموذج بديهي (Trivial) يفسر Y دون معرفة بقيم المتغيرات X=(X1,X2,...,Xp)، ويعتبر بذلك نموذجا بدئيا (Null Model) يقتصر فقط على المعامل الثابت a0.

انطلاقا من مبدأ التقتير، أحد المبادئ الأساسية لعلم الإحصاء، يجب أن يكون أي نموذج انحدار مقترح على الأقل أحسن من النموذج المرجعي. يتم التقييم عبر المراحل التالية:

  1. تقدير قيمة المعمل الثابت a0 للنموذج المرجعي M0: lnP(Y=1/X)1P(Y=1/X)=a0.
  2. حساب الانحراف الراسب D0.
  3. حساب الشبه - معاملات تحديد (Pseudo-R²) والتي تمثل مقاييسا لجودة النموذج (كلما اقتربت قيمتها من 1، كان النموذج ذا قوة تفسيرية معتبرة) وتكون دائما بدلالة الانحرافين الراسبين Dm و D0.

شبه معاملات التحديد يمكن قراءتها أيضا كتقدير لنسبة مساهمة المتغيرات X=(X1,X2,...,Xp) في تفسير المتغير Y. مثلا معامل بقيمة 0.2 يعني أن 80% من الظواهر المفسرة أو المؤثرة في Y هي مضمرة وأن النموذج المقترح لا يفسر إلا 20 بالمائة من الظاهرة المدروسة. من المهم التأكيد على أن هذا المؤشر لا يحتمل قراءات معاملات التحديد في نموذج الانحدار الخطي (التي يفترض فيها تجاوز مستويات قبول معينة، 0.5 مثلا)، في هذا السياق، يكفي أن يكون غير منعدم (مثلا أكبر من 0.1) لكي يكون النموذج المقترح له قابلية التصديق، بمعنى أنه يفسر المتغير Y أحسن من النموذج البديهي. في ما يلي أكثر المعاملات استعمالا مع صيغها الرياضية:

أكثر شبه معاملات التحديد استعمالا في تقييم نماذج الانحدار اللوجستي
شبه معامل التحديد (Pseudo-R²) الصيغة
R² ماك فادن McFadden RMF2=1LLmLL0
R² كوكس وسنيل Cox & Snell RCS2=1(L0Lm)(2/n)
R² ناغيلكرك Nagelkerke RN2=RCS2(1Lm)(2/n)


التقييم الخارجي

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

مصفوفة الإرباك

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

بنية عامة لمصفوفة إرباك
Y^=1 Y^=0 المجموع
Y=1 a b a+b
Y=0 c d c+d
المجموع a+c b+d a+b+c+d

بحيث:

  • a: عدد أفراد العينة الذين صح توقع قيمتهم الحقيقية 1 وفق النموذج. يصطلح عليهم بالقيم الإيجابية الصحيحة.
  • b: عدد أفراد العينة الذين خاب توقع قيمتهم الحقيقية 1 وفق النموذج. يصطلح عليهم بالقيم السلبية الخاطئة.
  • c: عدد أفراد العينة الذين صح توقع قيمتهم الحقيقية 0 وفق النموذج. يصطلح عليهم بالقيم السلبية الصحيحة.
  • d: عدد أفراد العينة الذين خاب توقع قيمتهم الحقيقية 0 وفق النموذج. يصطلح عليهم بالقيم الإيجابية الخاطئة.

نسبة خطأ النموذج هي نسبة التصنيفات الخاطئة بين أفراد العينة: ϵ=b+cn=1a+dn وتقابلها نسبة النجاح: θ=a+dn.

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

  • حساسية النموذج هي نسبة التوقعات الإيجابية الصحيحة: Se=aa+b
  • دقة النموذج هي نسبة التوقعات الإيجابية الصحيحة ضمن التوقعات الإيجابية: Pr=aa+c ويصطلح عليه في بعض المراجع بالقيمة التوقعية الإيجابية.
  • خصوصية النموذج هي نسبة التوقعات السلبية الصحيحة: Sp=dc+d وتقابلها نسبة القيم الإيجابية الخاطئة: 1Sp=cc+d.
  • مقياس فيشر (F-mesure) وهو مؤشر يأخذ بعين الاعتبار الدقة والحساسية في آن واحد، ويمكن عبر وزن متغير β من ترقية معيار الدقة على حساب الحساسية أو العكس، فأحيانا، في سياقات يكون فيها تحقق القيمة الإيجابية ذا تأثير جسيم ومهم على الفرد الإحصائي، يستحسن تقييم النموذج أساسا حسب قدرته على توقع القيم الإيجابية. بالمقابل، عندما تكون للقيمتين الإيجابية والسلبية نفس الأهمية، يفضل تحقيق التوازن بين الدقة والحساسية:
    • مقياس فيشر لوزن محدد β: F.mesureβ=Fβ=(1+β2)SePrβ2Se+Pr
    • قيمة β=1 هي القيمة الموازنة بين معياري الحساسية والدقة.
    • β<1 تمنح امتيازا للدقة على حساب الحساسية.
    • β>1 تمنح امتيازا للحساسية على حساب الدقة.

للإشارة، فمقياس فيشر هو متوسط توافقي للدقة والحساسية، فباعتبار α بحيث β2=1αα : Fβ=1α1Pr+(1α)1Se

معايير النموذج الأمثل

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

  • نسبتا الخطأ ϵ ونسبة القيم الإيجابية الخاطئة 1Sp يجب أن تؤول إلى 0.
  • الحساسية Se والدقة Pr والخصوصية Sp يجب أن تؤول إلى 1.

المعيار الأخير يشكل تحديا خاصا مقارنة بالأول، لأن Se و Sp يتغيران بطريقة متعارضة. إلا أن هذين المؤشرين لديهما ميزة مهمة، من منظور الاستدلال الإحصائي: هما أقل حساسية لطريقة الاعتيان مقارنة بنسبة الخطأ، ولا يفقدان وجاهتهما الإحصائية حتى في حالة العينات الصغيرة أو تلك التي تكون فيها القيم الإيجابية قليلة. بفضل هذه الخاصية، يتستعملان في تعريف مؤشري تقييم مهمين وجامعين هما مؤشر يودن (Youden Index) ومؤشر نسبة الإمكان (Likelihood Ratio):

  • مؤشر يودن: IY=Se+Sp1[0,1] ويؤول ل 1 بالنسبة للنماذج المثلى.
  • مؤشر نسبة الإمكان : L=Se1Sp ويجب أن يكون L>>1 بالنسبة للنماذج الجيدة.


التعريف

الدالة اللوجستية باللون الأحمر حيث z يمثل المحور الأفقي و (ƒ(z يمثل المحور العمودي

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

f(z)=ezez+1=11+ez

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

z=β0+β1x1+β2x2+β3x3++βkxk

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

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

مراجع

  1. ^ Q111421033، ص. 83، QID:Q111421033
  2. ^ "معلومات عن انحدار لوجستي على موقع catalogue.bnf.fr". catalogue.bnf.fr. مؤرشف من الأصل في 2019-05-02.
  3. ^ "معلومات عن انحدار لوجستي على موقع d-nb.info". d-nb.info. مؤرشف من الأصل في 2019-12-13.
  4. ^ "معلومات عن انحدار لوجستي على موقع id.worldcat.org". id.worldcat.org. مؤرشف من الأصل في 2019-12-13.
  5. ^ بيركسون، جوزيف (سبتمبر 1944). "Application of the logistic function to bio-assay". Journal of the American Statistical Association. مؤرشف من الأصل في 2016-03-10.