Jak wykonać test mediany nastroju w r

Test mediany nastroju służy do porównywania median dwóch lub więcej niezależnych grup.

Do wykonania tego testu w języku R można użyć funkcji median_test z biblioteki monet , która wykorzystuje następującą składnię:

median_test(odpowiedź~grupa, dane)

Złoto:

• odpowiedź: wektor wartości odpowiedzi
• grupa: wektor wartości grupujących
• dane: ramka danych zawierająca wektory odpowiedzi i grup

Poniższy przykład ilustruje sposób wykorzystania tej funkcji do przeprowadzenia testu mediany nastroju w R.

Przykład: Test mediany nastroju w R

Załóżmy, że nauczyciel chce wiedzieć, czy dwie różne metody nauki dają różne wyniki w testach wśród uczniów w jego klasie. Aby to przetestować, losowo prosi 10 uczniów, aby zastosowali jedną metodę nauki, a kolejnych 10 uczniów o inną. Po dwóch tygodniach każdy student przystępuje do tego samego egzaminu.

Postanawia zastosować test mediany Mooda, aby ustalić, czy mediana wyników egzaminu różni się w obu grupach.

Krok 1: Utwórz ramkę danych.

 #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

Krok 2: Wykonaj test średniego nastroju.

 #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

Wartość p testu wynosi 0,6613 . Ponieważ wartość ta jest nie mniejsza niż 0,05, nie możemy odrzucić hipotezy zerowej. Nie mamy wystarczających dowodów, aby stwierdzić, że istnieje statystycznie istotna różnica w medianie wyników egzaminów pomiędzy obiema grupami.

Domyślnie funkcja ta przypisuje wynik 0 obserwacjom dokładnie równym medianie. Można jednak określić tę wartość na 0,5 lub 1, używając argumentu mid.score .

Na przykład poniższy kod przeprowadza dokładnie ten sam test mediany nastroju, ale przypisuje wartość 0,5 obserwacjom równym medianie:

 #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

Testowa wartość p okazuje się wynosić 0,6459 , czyli jest nieco niższa od poprzedniej wartości p wynoszącej 0,6613 . Jednak wniosek z testu jest nadal taki sam: nie mamy wystarczających dowodów, aby stwierdzić, że mediana wyników egzaminów pomiędzy obiema grupami znacząco się różni.