Como ler os resultados do teste de Dunn?

12

Como leio os resultados do teste de Dunn ? Especificamente, o que os valores na tabela abaixo me dizem?

Eu tenho dados não paramétricos em 4 grupos e fiz primeiro um teste de Kruskal-Wallis para confirmar que as distribuições dos grupos eram diferentes umas das outras e do conjunto de dados agregado. Eu então usei o teste de Dunn para ver quais grupos eram diferentes entre si e quais não eram.

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
fonte

Respostas:

17

A saída após o teste de Kruskal-Wallis fornece todas as comparações possíveis em pares (seis no caso de quatro grupos). Portanto, o da primeira linha compara o grupo B com o grupo A, o primeiro da segunda linha compara o grupo C com o grupo A, etc.).

α/2 , embora esse valor p mude dependendo da taxa de erro familiar ou da opção de ajuste de múltiplas comparações da taxa de descoberta falsa. Se houver vários ajustes comparativos passo a passo (por exemplo, Holm, Benjamini-Hochberg, etc.), os valores p ajustados terão um asterisco ao lado deles se você rejeitar as hipóteses nulas no nível de significância especificado (o que não é necessariamente indicado diretamente pelo ajuste) p-valores desde que a rejeição depende do pedido ... consulte a documentação e as citações para obter mais detalhes.).

Eu sou o autor deste pacote (enviar um e-mail para mim, conforme indicado explicitamente na documentação, provavelmente seria a melhor maneira de entrar em contato diretamente comigo).

Alexis
fonte
Desculpe se isso está levando a pergunta em uma direção diferente, mas quando você ajustaria o teste? Quais critérios são necessários em seus dados para você escolher um ajuste?
pocketlizard
1
@pocketlizard Veja taxa de erro família-wise e taxa de detecção falsa .
Alexis
Vi que os valores de Kruskal-Wallis obtidos kruskal.testdiferem ligeiramente dos obtidos com a opção kw=TRUEno pacote acima. Isso se deve a algumas aproximações ou são duas coisas diferentes?
gented
@GennaroTedesco kruskal.testnão se ajusta a empates, enquanto dunn.testse ajusta a empates. Como diz a documentação ... logo acima, onde é enviado o email para perguntas de suporte.
Alexis