Hoe de mediane stemmingstest uit te voeren in r
De Mood mediaan test wordt gebruikt om de medianen van twee of meer onafhankelijke groepen te vergelijken.
De median_test- functie van de muntbibliotheek kan worden gebruikt om deze test in R uit te voeren, waarbij de volgende syntaxis wordt gebruikt:
mediaan_test(reactie~groep, gegevens)
Goud:
- respons: een vector van responswaarden
- groep: een vector van groeperingswaarden
- data: een dataframe dat de respons- en groepsvectoren bevat
Het volgende voorbeeld illustreert hoe u deze functie kunt gebruiken om de mediaan stemmingstest in R uit te voeren.
Voorbeeld: stemmingsmediaantest in R
Stel dat een leraar wil weten of twee verschillende studiemethoden verschillende toetsscores opleveren bij zijn leerlingen in zijn klas. Om dit te testen vraagt ze willekeurig 10 studenten om de ene studiemethode te gebruiken en nog eens 10 studenten om een andere te gebruiken. Na twee weken legt elke student hetzelfde examen af.
Ze besluit de mediaantest van Mood te gebruiken om te bepalen of de mediane examenscore tussen de twee groepen verschilt.
Stap 1: Maak het dataframe.
#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
Stap 2: Voer de mediane stemmingstest uit.
#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
De p-waarde van de test is 0,6613 . Omdat deze waarde niet minder dan 0,05 bedraagt, slagen we er niet in de nulhypothese te verwerpen. We hebben niet voldoende bewijs om te zeggen dat er een statistisch significant verschil bestaat in de gemiddelde examenscores tussen de twee groepen.
Standaard kent deze functie een score van 0 toe aan waarnemingen die exact gelijk zijn aan de mediaan. U kunt deze waarde echter opgeven als 0,5 of 1 met behulp van het argument mid.score .
De volgende code voert bijvoorbeeld precies dezelfde gemiddelde stemmingstest uit, maar kent een waarde van 0,5 toe aan waarnemingen die gelijk zijn aan de mediaan:
#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
De test-p-waarde blijkt 0,6459 te zijn, wat iets lager is dan de vorige p-waarde van 0,6613 . De conclusie van de test is echter nog steeds dezelfde: we hebben niet voldoende bewijs om te zeggen dat de mediane examenscores tussen de twee groepen significant verschillend zijn.