Bagaimana cara membaca hasil tes Dunn?

12

Bagaimana cara membaca hasil dari tes Dunn ? Secara khusus, apa yang nilai-nilai dalam tabel di bawah ini katakan padaku?

Saya memiliki data non-parametrik dalam 4 kelompok, dan pertama-tama saya melakukan tes Kruskal-Wallis untuk memastikan bahwa distribusi grup berbeda satu sama lain dan kumpulan data agregat. Saya kemudian menggunakan tes Dunn untuk melihat kelompok mana yang berbeda satu sama lain, dan yang tidak.

library(dunn.test)
dunn.test(data, g=area, kw=TRUE)
Kruskal-Wallis rank sum test

data: x and area
Kruskal-Wallis chi-squared = 1730.4401, df = 3, p-value = 0


                        Comparison of x by area                            
                            (No adjustment)                                
Row Mean-|
Col Mean |          A          B          C
---------+---------------------------------
       B |   20.62310
         |     0.0000
         |
       C |   26.66519  -0.087499
         |     0.0000     0.4651
         |
       D |   39.09084   7.401256   9.469204
         |     0.0000     0.0000     0.0000
bubbalouie
sumber

Jawaban:

17

Output yang mengikuti tes Kruskal-Wallis memberikan semua perbandingan berpasangan (enam dalam kasus empat kelompok). Jadi yang di baris pertama membandingkan grup B dengan grup A, yang pertama di baris kedua membandingkan grup C dengan grup A, dll.).

Angka atas untuk setiap perbandingan adalah statistik uji z berpasangan Dunn . Angka yang lebih rendah dalam contoh ini adalah nilai p- mentah yang terkait dengan tes (yaitu Anda akan membandingkan dengan , meskipun nilai p- ini berubah tergantung pada tingkat kesalahan berdasarkan keluarga atau tingkat penemuan salah beberapa opsi penyesuaian perbandingan perbandingan. Untuk beberapa penyesuaian perbandingan bertahap (mis. Holm, Benjamini-Hochberg, dll.), P -values ​​yang disesuaikan akan memiliki tanda bintang di sebelahnya jika Anda akan menolak hipotesis nol pada tingkat signifikansi yang ditentukan (yang tidak harus secara langsung ditunjukkan oleh disesuaikan pα/2-nilai karena penolakan tergantung pada pemesanan ... lihat dokumentasi dan kutipan di dalamnya untuk lebih jelasnya.).

Saya penulis paket ini (mengirim email kepada saya, seperti yang ditunjukkan secara eksplisit dalam dokumentasi, kemungkinan akan menjadi cara terbaik untuk menghubungi saya secara langsung).

Alexis
sumber
Maaf jika ini membawa pertanyaan ke arah yang berbeda, tetapi kapan Anda akan menyesuaikan tes? Kriteria apa yang diperlukan dalam data Anda agar Anda ingin memilih penyesuaian?
pocketlizard
Saya telah melihat bahwa nilai Kruskal-Wallis diperoleh dengan kruskal.testsedikit berbeda dari yang diperoleh dengan opsi kw=TRUEdalam paket di atas. Apakah ini karena beberapa perkiraan atau apakah mereka dua hal yang berbeda?
gented
@GennaroTedesco kruskal.testtidak menyesuaikan untuk ikatan, sedangkan dunn.testmenyesuaikan untuk ikatan. Seperti yang tertulis di dokumentasi ... tepat di atas di mana ia memberikan email untuk dihubungi untuk pertanyaan dukungan.
Alexis