Teste de Kruskal-Wallis

O teste de Kruskal-Wallis por postos, teste H de Kruskal-Wallis (que recebe este nome em homenagem a William Kruskal e W. Allen Wallis) ou análise de variância de um fator em postos é um método não paramétrico para testar se amostras se originam da mesma distribuição. É usado para comparar duas ou mais amostras independentes de tamanhos iguais ou diferentes. Ele estende o teste U de Mann-Whitney quando há mais de dois grupos. O equivalente paramétrico do teste de Kruskal-Wallis é o teste F usado na análise de variância de um fator. Um teste de Kruskal-Wallis significante indica que ao menos uma amostra domina estocasticamente uma outra amostra. O teste não identifica onde esta dominância estocástica ocorre ou para quantos pares de grupos se obtém dominância estocástica. O teste de Dunn ou o mais poderoso, mas menos conhecido teste de Conover-Iman ajudariam a analisar os pares específicos de amostras para dominância estocástica em testes post hoc.

Por ser um método não paramétrico, o teste de Kruskal-Wallis não assume uma distribuição normal dos resíduos, diferentemente da análoga análise de variância de um fator. Se o pesquisador puder assumir os pressupostos menos rigorosos de uma distribuição com forma e escala idênticas para todos os grupos, exceto para qualquer diferença nas medianas, então, a hipótese nula é de que as medianas de todos os grupos são iguais e a hipótese alternativa é de que ao menos a mediana de população de um grupo é diferente da mediana de população de ao menos um outro grupo.

Método

1. Classifique todos os dados de todos os grupos juntos, isto é, classifique os dados de 1 a $N$, ignorando a que grupo cada dado pertence. Atribua a quaisquer valores repetidos a média dos postos que eles teriam recebido se não fossem repetidos.
2. A estatística do teste é dada por:

   $H=(N-1){\frac {\sum _{i=1}^{g}n_{i}({\bar {r}}_{i\cdot }-{\bar {r}})^{2}}{\sum _{i=1}^{g}\sum _{j=1}^{n_{i}}(r_{ij}-{\bar {r}})^{2}}},$ em que,

   • $n_{i}$ é o número de observações do grupo $i$;
   • $r_{ij}$ é a classificação (entre todas as observações) de observação $j$ do grupo $i$;
   • $N$ é o número total de observações em todos os grupos;
   • ${\bar {r}}_{i\cdot }={\frac {\sum _{j=1}^{n_{i}}{r_{ij}}}{n_{i}}}$ é a classificação média de todas as observações no grupo $i$;
   • ${\bar {r}}={\tfrac {1}{2}}(N+1)$ é a média de todos os $r_{ij}$.
3. Se os dados não contêm valores repetidos, o denominador da expressão para $H$ é exatamente $(N-1)N(N+1)/12$ e ${\bar {r}}={\frac {N+1}{2}}$. Assim

   ${\begin{aligned}H&={\frac {12}{N(N+1)}}\sum _{i=1}^{g}n_{i}\left({\bar {r}}_{i\cdot }-{\frac {N+1}{2}}\right)^{2}\\&={\frac {12}{N(N+1)}}\sum _{i=1}^{g}n_{i}{\bar {r}}_{i\cdot }^{2}-\ 3(N+1).\end{aligned}}$
   A última fórmula contém apenas os quadrados das classificações médias.

4. Uma correção para valores repetidos, usando-se a fórmula mais curta descrita no ponto anterior, pode ser feita ao dividir $H$ por $1-{\frac {\sum _{i=1}^{G}(t_{i}^{3}-t_{i})}{N^{3}-N}}$, em que $G$ é o número de agrupamento de diferentes postos empatados e $t_{i}$ é o número de observações repetidas dentro do grupo $i$ que são iguais a um valor particular. Esta correção geralmente faz pouca diferença no valor de $H$ a não ser que haja um grande número de valores repetidos.
5. Finalmente, o valor-p é aproximado por $\Pr(\chi _{g-1}^{2}\geq H)$. Se alguns valores $n_{i}$ forem pequenos (isto é, menores que 5), a distribuição de probabilidade de $H$ pode ser bem diferente da distribuição qui-quadrado. Se uma tabela de distribuição qui-quadrado estiver disponível, o valor crítico $\chi _{\alpha :g-1}^{2}$ pode ser encontrado na tabela a $g-1$ graus de liberdade, procurando abaixo da significância desejada ou nível alfa.
6. Se a estatística não for significante, então, não há evidência de dominância estocástica entre as amostras. Entretanto, se o teste for significante , então, ao menos uma amostra domina estocasticamente outra amostra. Por isso, um pesquisador pode usar contrastes de amostras entre pares de amostras individuais ou fazer testes post hoc usando um teste de Dunn, que (1) emprega apropriadamente as mesmas classificações do teste de Kruskal-Wallis e (2) emprega apropriadamente a variância combinada implicada pela hipótese nula do teste de Kruskal-Wallis a fim de determinar quais dos pares de amostras são significantemente diferentes. Quando se realizam múltiplos contrastes de amostras ou testes, a taxa de erro de tipo I tende a se tornar inflada, levantando preocupações quanto às múltiplas comparações.

Tabelas de probabilidades exatas

Uma grande quantidade de recursos computacionais é exigida para computar as probabilidades exatas para o teste de Kruskal-Wallis. Os softwares existentes apenas oferecem probabilidades exatas para amostras com menos de 30 integrantes. Estes programas dependem de aproximação assintótica para amostras de tamanhos maiores. Valores de probabilidades exatas para amostras de tamanhos maiores estão disponíveis. J. D. Spurrier publicou as tabelas de probabilidades exatas para amostras com até 45 integrantes. J. Patrick Meyer e Michael A. Seaman produziram as distribuições de probabilidades exatas para amostras com até 105 participantes.

Ver também

• Análise de variância
• Teste U de Mann-Whitney
• Teste de Friedman

Ligações externas

• Uma versão on-line do teste (em inglês)