Comment lire les résultats du test de Dunn '?
On février 1, 2021 by adminComment lire les résultats de Dunn « s test ? Plus précisément, que me disent les valeurs du tableau ci-dessous?
Jai des données non paramétriques dans 4 groupes, et jai dabord fait un Kruskal-Wallis test pour confirmer que les distributions des groupes étaient différentes les unes des autres et de lensemble de données agrégées. Jai 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
Réponse
Le résultat suivant le test de Kruskal-Wallis fournit toutes les comparaisons possibles par paires (six dans le cas de quatre groupes). Ainsi, celui de la première ligne compare le groupe B au groupe A, le premier sur la deuxième ligne compare le groupe C avec le groupe A, etc.).
Le nombre supérieur pour chaque comparaison est la statistique de test par paires de Dunn z . Le nombre inférieur est dans cet exemple la valeur brute p associée au test (cest-à-dire que vous compareriez à $ \ alpha / 2 $, bien que cette valeur p change en fonction sur loption dajustement des comparaisons multiples du taux derreur par famille ou du taux de fausses découvertes. Pour les ajustements de comparaisons multiples par étapes (par exemple, Holm, Benjamini-Hochberg, etc.), les valeurs p ajustées seront suivies dun astérisque à eux si vous rejetez les hypothèses nulles au niveau de signification spécifié (ce qui nest pas nécessairement directement indiqué par les valeurs ajustées de p puisque le rejet dépend de lordre … voir la documentation et les citations qui y figurent pour plus détails.).
Je suis lauteur de ce paquet (menvoyer un e-mail, comme indiqué explicitement dans la documentation, serait probablement le meilleur moyen de me contacter directement).
kruskal.test
diffèrent légèrement des ceux obtenus avec loptionkw=TRUE
dans le package ci-dessus. Est-ce dû à des approximations ou sagit-il de deux choses différentes?kruskal.test
ne tient pas compte des égalités, alors quedunn.test
ajuste les liens. Comme il est indiqué dans la documentation … juste au-dessus où il donne le-mail à contacter pour les questions dassistance.