كيفية إجراء اختبار المزاج المتوسط في r

يُستخدم اختبار المزاج المتوسط لمقارنة متوسطات مجموعتين مستقلتين أو أكثر.

يمكن استخدام وظيفة median_test الخاصة بمكتبة العملات المعدنية لإجراء هذا الاختبار في R، والذي يستخدم بناء الجملة التالي:

اختبار متوسط (الاستجابة ~ المجموعة، البيانات)

ذهب:

• الاستجابة: ناقل لقيم الاستجابة
• المجموعة: ناقل لقيم التجميع
• البيانات: إطار بيانات يحتوي على الاستجابة ومتجهات المجموعة

يوضح المثال التالي كيفية استخدام هذه الوظيفة لإجراء اختبار المزاج المتوسط في R.

مثال: اختبار متوسط المزاج في R

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

قررت استخدام اختبار المزاج المتوسط لتحديد ما إذا كانت درجة الاختبار المتوسطة تختلف بين المجموعتين.

الخطوة 1: إنشاء إطار البيانات

 #createdata
method = rep(c('method1', 'method2'), each=10)
score = c(75, 77, 78, 83, 83, 85, 89, 90, 91, 97, 77, 80, 84, 84, 85, 90, 92, 92, 94, 95)
examData = data.frame(method, score)

#viewdata
examData

    method score
1 method1 75
2 method1 77
3 method1 78
4 method1 83
5 method1 83
6 method1 85
7 method1 89
8 method1 90
9 method1 91
10 method1 97
11 method2 77
12 method2 80
13 method2 84
14 method2 84
15 method2 85
16 method2 90
17 method2 92
18 method2 92
19 method2 94
20 method2 95

الخطوة 2: إجراء اختبار المزاج المتوسط.

 #load the coin library
library(corner)

#perform Mood's Median Test
median_test(score~method, data = examData)

#output
	Asymptotic Two-Sample Brown-Mood Median Test

data: score by method (method1, method2)
Z = -0.43809, p-value = 0.6613
alternative hypothesis: true mu is not equal to 0

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

افتراضيًا، تقوم هذه الدالة بتعيين درجة 0 للملاحظات المساوية تمامًا للوسيط. ومع ذلك، يمكنك تحديد هذه القيمة لتكون 0.5 أو 1 باستخدام الوسيطة mid.score .

على سبيل المثال، الكود التالي ينفذ بالضبط نفس اختبار الحالة المزاجية المتوسطة، لكنه يعين قيمة 0.5 للملاحظات المساوية للمتوسط:

 #perform Mood's Median Test
median_test(score~method, mid.score="0.5" , data = examData)

#output
	Asymptotic Two-Sample Brown-Mood Median Test

data: score by method (method1, method2)
Z = -0.45947, p-value = 0.6459
alternative hypothesis: true mu is not equal to 00

تبين أن قيمة الاختبار p هي 0.6459 ، وهي أقل قليلاً من القيمة p السابقة البالغة 0.6613 . ومع ذلك، فإن نتيجة الاختبار لا تزال هي نفسها: ليس لدينا أدلة كافية لنقول أن متوسط درجات الاختبار بين المجموعتين تختلف بشكل كبير.