Bagaimana melakukan tes dunn di r

Uji Kruskal-Wallis digunakan untuk menentukan apakah terdapat perbedaan yang signifikan secara statistik antara median tiga atau lebih kelompok independen. Ini dianggap setara nonparametrik dari ANOVA satu arah .

Jika hasil uji Kruskal-Wallis signifikan secara statistik, maka uji Dunn sebaiknya dilakukan untuk menentukan secara pasti kelompok mana yang berbeda.

Tutorial ini menjelaskan cara melakukan tes Dunn di R.

Contoh: Tes Dunn di R

Seorang peneliti ingin mengetahui apakah ketiga obat memiliki efek berbeda terhadap sakit punggung. Jadi dia merekrut 30 orang yang semuanya menderita nyeri punggung serupa dan secara acak membagi mereka menjadi tiga kelompok untuk menerima obat A, obat B, atau obat C. Sebulan setelah meminum obat tersebut, peneliti meminta setiap individu untuk menilai nyeri punggung mereka berdasarkan skala 1 sampai 100, dengan 100 menunjukkan nyeri paling parah.

Peneliti melakukan uji Kruskal-Wallis menggunakan tingkat signifikansi 0,05 untuk menentukan apakah terdapat perbedaan yang signifikan secara statistik antara median peringkat nyeri punggung pada ketiga kelompok ini.

Kode berikut menunjukkan cara membuat bingkai data di R dan melakukan pengujian Kruskal-Wallis:

 #make this example reproducible
set.seed(0)

#create data frame
data <- data.frame(drug = rep(c("A", "B", "C"), each = 10),
                   bread = c(runif(10, 40, 60),
                            runif(10, 45, 65),
                            runif(10, 55, 70)))

#view first six rows of data frame
head(data)

# drug pain
#1 A 57.93394
#2 A 45.31017
#3 A 47.44248
#4 A 51.45707
#5 A 58.16416
#6 A 44.03364

#perform Kruskal-Wallis Test
kruskal.test(pain ~ drug, data = data)

	Kruskal-Wallis rank sum test

data: pain by drug
Kruskal-Wallis chi-squared = 11.105, df = 2, p-value = 0.003879

Karena nilai p keseluruhan ( 0,003879 ) kurang dari 0,05, ini berarti terdapat perbedaan yang signifikan secara statistik dalam tingkat nyeri yang dilaporkan antara ketiga obat. Jadi kita dapat melakukan uji Dunn untuk menentukan obat mana yang berbeda.

Kode berikut menunjukkan cara melakukan pengujian Dunn di R menggunakan fungsi dunnTest() dari pustaka FSA() :

 #loadlibrary
library(FSA)

#perform Dunn's Test with Bonferroni correction for p-values
dunnTest(pain ~ drug,
         data=data,
         method=" bonferroni ")

Dunn (1964) Kruskal-Wallis multiple comparison
  p-values adjusted with the Bonferroni method.

  Comparison Z P.unadj P.adj
1 A - B -0.8890009 0.374002602 1.000000000
2 A - C -3.2258032 0.001256197 0.003768591
3 B - C -2.3368023 0.019449464 0.058348393

Perhatikan bahwa kami memilih untuk menggunakan koreksi Bonferroni untuk nilai p dari beberapa perbandingan, namun opsi lain yang memungkinkan meliputi:

• “sidak” (penyesuaian sidak)
• “holm” (penyesuaian holom)
• “hs” (penyesuaian Holm-Sidak)
• “bs” (penyesuaian Bonferroni-Sidak)
• “oleh” (penyesuaian Benjamin-Yekuteili)
• “bh” ( prosedur Benjamini-Hochberg )

Pada α = 0,05, obat A dan C adalah dua obat yang secara statistik berbeda secara signifikan satu sama lain (nilai p yang disesuaikan = 0,003768 ).