Jak czytać wyniki testu Dunna?

12

Jak odczytać wyniki testu Dunna ? Co konkretnie mówią mi wartości z poniższej tabeli?

Posiadam dane nieparametryczne w 4 grupach i najpierw wykonałem test Kruskala-Wallisa, aby potwierdzić, że rozkłady grup były odmienne od siebie i zestaw danych zbiorczych. Następnie użyłem testu Dunna, aby zobaczyć, które grupy różnią się od siebie, a które nie.

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
źródło

Odpowiedzi:

17

Dane wyjściowe po teście Kruskala-Wallisa zapewniają wszystkie możliwe porównania parami (sześć w przypadku czterech grup). Tak więc ten w pierwszym rzędzie porównuje grupę B z grupą A, pierwszy w drugim rzędzie porównuje grupę C z grupą A itd.).

Górna liczba jest dla każdego porównania Dunna parami oo statystyka testowa. Niższa liczba jest w tym przykładzie nieprzetworzoną wartością p związaną z testem (tj. Porównałbyś się z , chociaż ta wartość p zmienia się w zależności od rodzinnego poziomu błędu lub opcji dopasowania wielokrotnych porównań częstotliwości wykrywania. W przypadku wielokrotnych korekt porównawczych krok po kroku (np. Holm, Benjamini-Hochberg itp.), Skorygowane wartości p będą miały gwiazdkę obok nich, jeśli odrzucisz hipotezy zerowe na określonym poziomie istotności (co niekoniecznie jest bezpośrednio wskazane przez skorygowane pα/2)-wartości, ponieważ odrzucenie zależy od zamówienia ... więcej informacji można znaleźć w dokumentacji i cytatach).

Jestem autorem tego pakietu (e-maile, jak wyraźnie zaznaczono w dokumentacji, prawdopodobnie byłyby najlepszym sposobem na bezpośredni kontakt ze mną).

Alexis
źródło
Przepraszam, jeśli to zmierza w innym kierunku, ale kiedy dostosujesz test? Jakie kryteria są niezbędne w danych, aby wybrać korektę?
pocketlizard
Widziałem, że wartości Kruskala-Wallisa uzyskane z kruskal.testnieznacznie różnią się od tych uzyskanych z opcją kw=TRUEw powyższym pakiecie. Czy wynika to z pewnych przybliżeń, czy są to dwie różne rzeczy?
gented
@GennaroTedesco kruskal.testnie dostosowuje się do więzi, dunn.testale dostosowuje się do więzi. Jak napisano w dokumentacji ... tuż powyżej, gdzie podaje e-mail do kontaktu w sprawie pytań dotyczących pomocy technicznej.
Alexis,