Come leggere i risultati del test di Dunn '?
Su Febbraio 1, 2021 da adminCome leggo i risultati di Dunn “s test ? Nello specifico, cosa mi dicono i valori nella tabella sottostante?
Ho dati non parametrici in 4 gruppi e ho prima fatto un Kruskal-Wallis test per confermare che le distribuzioni dei gruppi erano dissimili luna dallaltra e dal set di dati aggregato. Ho quindi utilizzato il test di Dunn per vedere quali gruppi erano dissimili tra loro e quali no.
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
Risposta
Loutput che segue il test di Kruskal-Wallis fornisce tutti i possibili confronti a coppie (sei nel caso di quattro gruppi). Quindi quello sulla prima riga confronta il gruppo B con il gruppo A, il primo sulla seconda riga confronta il gruppo C con il gruppo A, ecc.).
Il numero superiore per ogni confronto è la statistica del test z a coppie di Dunn. Il numero più basso è in questo esempio il valore grezzo p associato al test (cioè si confronterebbe con $ \ alpha / 2 $, sebbene questo valore p cambi a seconda sullopzione di regolazione del confronto multiplo del tasso di errore familiare o del tasso di false scoperte. Per gli aggiustamenti graduali del confronto multiplo (ad esempio Holm, Benjamini-Hochberg, ecc.), i valori di p aggiustati avranno un asterisco accanto a loro se rifiutassi le ipotesi nulle al livello di significatività specificato (che non è necessariamente indicato direttamente dai valori p aggiustati poiché il rifiuto dipende dallordinamento … vedi la documentazione e le citazioni ivi contenute per ulteriori informazioni dettagli.).
Sono lautore di questo pacchetto (inviarmi une-mail, come esplicitamente indicato nella documentazione, sarebbe probabilmente il modo migliore per contattarmi direttamente).
kruskal.test
differiscono leggermente dal quelli ottenuti con lopzionekw=TRUE
nel pacchetto precedente. Ciò è dovuto ad alcune approssimazioni o sono due cose diverse?kruskal.test
non si adatta ai legami, mentredunn.test
si regola per i pareggi. Come dice nella documentazione … appena sopra dove fornisce lemail da contattare per domande di supporto.