{"id":1050,"date":"2023-07-27T20:35:01","date_gmt":"2023-07-27T20:35:01","guid":{"rendered":"https:\/\/statorials.org\/nl\/dunns-test-in-r\/"},"modified":"2023-07-27T20:35:01","modified_gmt":"2023-07-27T20:35:01","slug":"dunns-test-in-r","status":"publish","type":"post","link":"https:\/\/statorials.org\/nl\/dunns-test-in-r\/","title":{"rendered":"Hoe de dunn-test uit te voeren in r"},"content":{"rendered":"

<\/p>\n

\n

Een Kruskal-Wallis-test<\/a> wordt gebruikt om te bepalen of er al dan niet een statistisch significant verschil bestaat tussen de medianen van drie of meer onafhankelijke groepen. Het wordt beschouwd als het niet-parametrische equivalent van eenrichtings-ANOVA<\/a> .<\/span><\/p>\n

Als de resultaten van een Kruskal-Wallis-test statistisch significant zijn, is het passend om de Dunn-test<\/a> uit te voeren om precies te bepalen welke groepen verschillend zijn.<\/span><\/p>\n

In deze tutorial wordt uitgelegd hoe u de Dunn-test uitvoert in R.<\/span><\/p>\n

Voorbeeld: Dunn-test in R<\/strong><\/span><\/h3>\n

Een onderzoeker wil weten of drie medicijnen verschillende effecten hebben op rugpijn. Daarom recruteert hij 30 mensen die allemaal aan soortgelijke rugpijn lijden en verdeelt ze willekeurig in drie groepen die medicijn A, medicijn B of medicijn C krijgen. Een maand na inname van het medicijn vraagt de onderzoeker elk individu om hun rugpijn te beoordelen. een schaal van 1 tot 100, waarbij 100 de meest ernstige pijn aangeeft.<\/span><\/p>\n

De onderzoeker voert een Kruskal-Wallis-test uit met een significantieniveau van 0,05 om te bepalen of er een statistisch significant verschil bestaat tussen de gemiddelde rugpijnbeoordelingen in deze drie groepen.<\/span><\/p>\n

De volgende code laat zien hoe u het dataframe in R maakt en een Kruskal-Wallis-test uitvoert:<\/span><\/p>\n

 #make this example reproducible\nset.seed(0)\n<\/span>\n#create data frame\ndata <- data.frame(drug = rep(c(\"A\", \"B\", \"C\"), each = 10),\n                   bread = c(runif(10, 40, 60),\n                            runif(10, 45, 65),\n                            runif(10, 55, 70)))<\/span>\n\n#view first six rows of data frame\nhead(data)\n<\/span>\n# drug pain\n#1 A 57.93394\n#2 A 45.31017\n#3 A 47.44248\n#4 A 51.45707\n#5 A 58.16416\n#6 A 44.03364\n\n#perform Kruskal-Wallis Test<\/span>\nkruskal.test(pain ~ drug, data = data)\n\n\tKruskal-Wallis rank sum test\n\ndata: pain by drug\nKruskal-Wallis chi-squared = 11.105, df = 2, p-value = 0.003879\n<\/span><\/span><\/strong><\/pre>\n
 Aangezien de totale p-waarde ( 0,003879<\/strong> ) kleiner is dan 0,05, betekent dit dat er een statistisch significant verschil is in de gerapporteerde pijnniveaus tussen de drie medicijnen. We kunnen dus de Dunn-test uitvoeren om precies te bepalen welke medicijnen verschillend zijn.<\/span><\/p>\n
 #loadlibrary\nlibrary(FSA)\n\n#perform Dunn's Test with Bonferroni correction for p-values<\/span>\ndunnTest(pain ~ drug,\n         data=data,\n         method=\" bonferroni<\/span> \")\n\nDunn (1964) Kruskal-Wallis multiple comparison\n  p-values adjusted with the Bonferroni method.\n\n  Comparison Z P.unadj P.adj\n1 A - B -0.8890009 0.374002602 1.000000000\n2 A - C -3.2258032 0.001256197 0.003768591\n3 B - C -2.3368023 0.019449464 0.058348393<\/span><\/span><\/strong><\/pre>\n
 Merk op dat we ervoor hebben gekozen om een Bonferroni-correctie te gebruiken voor de p-waarden van meerdere vergelijkingen, maar andere mogelijke opties zijn onder meer:<\/span><\/p>\n
 De volgende code laat zien hoe u de Dunn-test in R uitvoert met behulp van de dunnTest()-<\/a> functie uit de FSA()-<\/strong> bibliotheek:<\/span><\/p>\n