كيفية إجراء اختبار دن في بايثون

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

إذا كانت نتائج اختبار كروسكال واليس ذات دلالة إحصائية، فمن المناسب إجراء اختبار دن لتحديد المجموعات المختلفة بالضبط.

يشرح هذا البرنامج التعليمي كيفية إجراء اختبار Dunn في بايثون.

مثال: اختبار دن في بايثون

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

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

لإجراء اختبار Dunn في بايثون، يمكننا استخدام الدالة posthoc_dunn() ‎ من مكتبة scikit-posthocs.

يوضح الكود التالي كيفية استخدام هذه الوظيفة:

الخطوة 1: تثبيت scikit-posthocs.

نحتاج أولاً إلى تثبيت مكتبة scikit-posthocs:

 pip install scikit-posthocs

الخطوة الثانية: إجراء اختبار دان.

ثم يمكننا إنشاء البيانات وإجراء اختبار Dunn:

 #specify the growth of the 10 plants in each group
group1 = [7, 14, 14, 13, 12, 9, 6, 14, 12, 8]
group2 = [15, 17, 13, 15, 15, 13, 9, 12, 10, 8]
group3 = [6, 8, 8, 9, 5, 14, 13, 8, 10, 9]
data = [group1, group2, group3]

#perform Dunn's test using a Bonferonni correction for the p-values
import scikit_posthocs as sp
sp. posthoc_dunn (data, p_adjust = ' bonferroni ')

               1 2 3
1 1.000000 0.550846 0.718451
2 0.550846 1.000000 0.036633
3 0.718451 0.036633 1.000000

لاحظ أننا اخترنا استخدام تصحيح Bonferroni للقيم p للتحكم في معدل الخطأ العائلي ، لكن الاختيارات المحتملة الأخرى للوسيطة p_adjust تتضمن ما يلي:

• صديق
• هولم سيداك
• سيمز هوشبيرج
• manl
• fdr_bh
• fdr_by
• fdr_tsbh

راجع الوثائق للحصول على مزيد من التفاصيل حول كل طريقة من طرق تعديل القيمة p.

الخطوة 3: تفسير النتائج.

ومن نتائج اختبار دان نلاحظ ما يلي:

• القيمة p المعدلة للفرق بين المجموعة 1 والمجموعة 2 هي 0.550846 .
• القيمة p المعدلة للفرق بين المجموعة 1 والمجموعة 3 هي 0.718451 .
• القيمة p المعدلة للفرق بين المجموعة 2 والمجموعة 3 هي 0.036633 .

وبالتالي، فإن المجموعتين الوحيدتين اللتين تختلفان بشكل كبير إحصائيًا عند α = 0.05 هما المجموعتان 2 و3.

مصادر إضافية

مقدمة لاختبار دن للمقارنات المتعددة
كيفية إجراء اختبار دن في R