Jak wykonać test kruskala-wallisa w sas

Test Kruskala-Wallisa służy do określenia, czy istnieje statystycznie istotna różnica między medianami trzech lub więcej niezależnych grup.

Uważa się, że jest to nieparametryczny odpowiednik jednokierunkowej analizy wariancji ANOVA .

W tym samouczku przedstawiono krok po kroku przykład wykonania testu Kruskala-Wallisa w SAS-ie.

Krok 1: Wprowadź dane

Załóżmy, że badacze chcą wiedzieć, czy trzy różne nawozy prowadzą do różnych poziomów wzrostu roślin.

Losowo wybierają 30 różnych roślin i dzielą je na trzy grupy po 10, stosując do każdej grupy inny nawóz. Po miesiącu mierzą wysokość każdej rośliny.

Wprowadzimy następujące dane do SAS, które pokazują całkowity wzrost (w calach) każdej z 10 roślin w każdej grupie:

 /*create dataset*/
data fertilizer_data;
    input fertilizer $growth;
    datalines ;
fert1 7
fert1 14
fert1 14
fert1 13
fert1 12
fert1 9
fert1 6
fert1 14
fert1 12
fert1 8
fert2 15
fert2 17
fert2 13
fert2 15
fert2 15
fert2 13
fert2 9
fert2 12
fert2 10
fert2 8
fert3 6
fert3 8
fert3 8
fert3 9
fert3 5
fert3 14
fert3 13
fert3 8
fert3 10
fert3 9
;
run ;

Krok 2: Wykonaj test Kruskala-Wallisa

Następnie użyjemy instrukcji proc npar1way do wykonania testu Kruskala-Wallisa w celu porównania średniego wzrostu roślin pomiędzy trzema grupami nawozów:

 /*perform Kruskal-Wallis test*/
proc npar1way data =fertilizer_data wilcoxon dscf ;
    class fertilizer;
    vargrowth ;
run ;

Krok 3: Interpretacja wyników

Pierwsza tabela wyników przedstawia ogólną statystykę testu Chi-kwadrat i odpowiadającą jej wartość p dla testu Kruskala-Wallisa:

Wartość p testu wynosi 0,0431 . Ponieważ wartość ta jest mniejsza niż 0,05, odrzucamy hipotezę zerową, że średni wzrost roślin jest taki sam dla wszystkich trzech nawozów.

Oznacza to, że mamy wystarczające dowody, aby stwierdzić, że rodzaj użytego nawozu powoduje statystycznie istotne różnice we wzroście roślin.

Ostateczna tabela wyników pokazuje wartości p dla porównań parami pomiędzy każdą z trzech grup:

Z tej tabeli widzimy, że jedyną wartością p poniżej 0,05 jest porównanie nawozu 2 z nawozem 3, dla którego wartość p wynosi 0,0390 .

Oznacza to, że istnieje statystycznie istotna różnica we wzroście roślin pomiędzy nawozem 2 i nawozem 3, ale nie pomiędzy pozostałymi porównaniami parami.

Dodatkowe zasoby

Poniższe samouczki wyjaśniają, jak wykonywać inne popularne testy statystyczne w SAS-ie:

Jak wykonać test t dla jednej próby w SAS
Jak wykonać test t dla dwóch próbek w SAS
Jak wykonać jednokierunkową ANOVA w SAS
Jak wykonać dwukierunkową ANOVA w SAS