Kruskal-Wallis

# o teste de Kruskal-Wallis pode substituir uma análise de variância em delinemaento inteiramente ao acaso em situações onde a variável resposta é qualitativa ordinal ou quando a variável resposta é quantitativa e não atende pressuposições da análise de variância

# para aplicar o teste será usado a função ea1() do pacote "easyanova" que deve ser carregado

require(easyanova)

# No exemplo abaixo segue programação para obter um conjunto de dados hipotéticos, onde os tratamentos são avaliados em 3 grupos de 6 animais. A variável resposta é o grau de doença sendo: 1= doença ausente; 2=doença em estágio inicial, 3=doença em estágio intermediário e 4=doença em estágio final.

tratamentos=c(rep("medicamento_A",6),rep( "medicamento_B",6), rep("placebo",6))

resposta=c(1,2,2,1,1,2,1,4,3,3,4,2,2,3,4,4,4,2)

dados=data.frame(tratamentos,resposta)

dados

# aplicação do teste, utilizando o design=14, que informa a função que será realizado o teste de Kruskal-Wallis

# observe que ocorre diferença significativa entre tratamentos p<0,05 e no pós teste conclui-se que o medicamento B e o placebo são estatisticamente iguais e diferentes do medicamento A

# para salvar o resultado em arquivo word

capture.output(resultado, file="Resultado.doc")