Como realizar um teste kruskal-wallis no sas

Um teste de Kruskal-Wallis é usado para determinar se existe ou não uma diferença estatisticamente significativa entre as medianas de três ou mais grupos independentes.

É considerado o equivalente não paramétrico da ANOVA unidirecional .

Este tutorial fornece um exemplo passo a passo de como realizar um teste Kruskal-Wallis no SAS.

Passo 1: Insira os dados

Suponha que os pesquisadores queiram saber se três fertilizantes diferentes levam a diferentes níveis de crescimento das plantas.

Eles selecionam aleatoriamente 30 plantas diferentes e as dividem em três grupos de 10, aplicando um fertilizante diferente em cada grupo. Depois de um mês, medem a altura de cada planta.

Inseriremos os seguintes dados no SAS, que mostram o crescimento total (em polegadas) de cada uma das 10 plantas de cada grupo:

 /*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 ;

Passo 2: Realize o teste Kruskal-Wallis

A seguir, usaremos a instrução proc npar1way para realizar um teste de Kruskal-Wallis para comparar o crescimento médio das plantas entre os três grupos de fertilizantes:

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

Etapa 3: interprete os resultados

A primeira tabela do resultado mostra a estatística geral do teste Qui-quadrado e o valor p correspondente para o teste Kruskal-Wallis:

O valor p do teste é 0,0431 . Como este valor é inferior a 0,05, rejeitamos a hipótese nula de que o crescimento médio das plantas é o mesmo para os três fertilizantes.

Isto significa que temos evidências suficientes para concluir que o tipo de fertilizante utilizado provoca diferenças estatisticamente significativas no crescimento das plantas.

A tabela de resultados finais mostra os valores de p para as comparações pareadas entre cada um dos três grupos:

Nessa tabela podemos perceber que o único valor de p abaixo de 0,05 é a comparação entre o fertilizante 2 e o fertilizante 3, que possui valor de p de 0,0390 .

Isto significa que existe uma diferença estatisticamente significativa no crescimento das plantas entre o fertilizante 2 e o fertilizante 3, mas não entre as outras comparações aos pares.

Recursos adicionais

Os tutoriais a seguir explicam como realizar outros testes estatísticos comuns no SAS:

Como realizar um teste t de uma amostra no SAS
Como realizar um teste t de duas amostras no SAS
Como realizar ANOVA unidirecional no SAS
Como realizar ANOVA bidirecional no SAS