Медианный тест настроения используется для сравнения медиан двух или более независимых групп.
Для выполнения этого теста в R можно использовать функцию median_test из библиотеки монет , которая использует следующий синтаксис:
median_test (группа ответов, данные)
куда:
- ответ: вектор значений ответа
- group: вектор значений группировки
- данные: кадр данных, содержащий вектор ответа и группы
В следующем примере показано, как использовать эту функцию для выполнения медианного теста настроения в R.
Пример: Медианный тест настроения в R
Предположим, учитель хочет знать, дают ли два разных метода обучения разные экзаменационные баллы среди его учеников. Чтобы проверить это, она случайным образом поручает 10 ученикам использовать один метод обучения, а другим 10 ученикам — другой. Через две недели каждый студент сдает один и тот же экзамен.
Она решает использовать медианный тест настроения, чтобы определить, различается ли средний балл на экзамене между двумя группами.
Шаг 1: Создайте фрейм данных.
#create data
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)
#view data
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(coin)
#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.Однако вывод теста остается прежним: у нас нет достаточных доказательств, чтобы сказать, что медианы экзаменационных баллов между двумя группами значительно различаются.