Comment lire les résultats du test de Dunn?

12

Comment lire les résultats du test de Dunn ? Plus précisément, que m'indiquent les valeurs du tableau ci-dessous?

J'ai des données non paramétriques dans 4 groupes, et j'ai d'abord fait un test de Kruskal-Wallis pour confirmer que les distributions des groupes étaient différentes les unes des autres et de l'ensemble de données agrégé. J'ai ensuite utilisé le test de Dunn pour voir quels groupes étaient différents les uns des autres et lesquels ne l'étaient pas.

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
la source

Réponses:

17

La sortie suivant le test de Kruskal-Wallis fournit toutes les comparaisons par paires possibles (six dans le cas de quatre groupes). Ainsi, celui de la première ligne compare le groupe B au groupe A, le premier de la deuxième ligne compare le groupe C au groupe A, etc.).

α/2-les valeurs puisque le rejet dépend de la commande ... voir la documentation et les citations qui s'y trouvent pour plus de détails.).

Je suis l'auteur de ce package (m'envoyer un e-mail, comme indiqué explicitement dans la documentation, serait probablement le meilleur moyen de me contacter directement).

Alexis
la source
Désolé si cela prend la question dans une direction différente, mais quand ajusteriez-vous le test? Quels critères sont nécessaires dans vos données pour que vous souhaitiez choisir un ajustement?
pocketlizard
1
@pocketlizard Voir le taux d'erreur par famille et le taux de fausses découvertes .
Alexis
J'ai vu que les valeurs de Kruskal-Wallis obtenues avec kruskal.testdiffèrent légèrement de celles obtenues avec l'option kw=TRUEdans le package ci-dessus. Est-ce dû à des approximations ou s'agit-il de deux choses différentes?
gented
@GennaroTedesco kruskal.testne s'ajuste pas pour les cravates, alors dunn.testqu'il s'ajuste pour les cravates. Comme il est dit dans la documentation ... juste au-dessus de l'endroit où il donne l'e-mail à contacter pour des questions d'assistance.
Alexis