Friedman

# o teste de Friedman pode substituir uma análise de variância em delinemaento de blocos 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 conjunto de dados hipotéticos, onde os tratamentos (quatro tipos de vinho) são avaliados por 5 provadores. A variável resposta é o grau de qualidade da bebida sendo: 1=muito ruim, 2=ruim, 3=mediana, 4=boa, 5=muito boa.

vinhos=rep(c("Vinho_A","Vinho_B","Vinho_C","Vinho_D"),each=5)

provadores=rep(c("a","b","c","d","e"),4)

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

dados=data.frame(vinhos, provadores, resposta)

# observe a ordem das colunas, primeiro os tratamentos, depois os blocos (provadores) e depois a variável resposta

dados

# aplicação do teste, utilizando a função ea1() e o argumento design=15, indicando que será feito o teste de Friedman

resultado=ea1(dados, design=15)

# para observar o resultado

resultado

# observe que ocorre diferença significativa entre tratamentos p<0,05 e no pós teste conclui-se existe diferença estatisticamente significativa entre os vinhos A, D (melhores notas) e um grupo formado por B e C (vinhos inferiores)

# para salvar o resultado em arquivo Word

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