teste=matrix(c(3640,4890,4800,4460,
4200,4550,5320,5500,
4700,6020,5250,5580,
5300,5900,5150,5560) , nrow = 4, byrow =
TRUE,
dimnames = list(1 : 4,c("Trat1", "Trat2", "Trat3","Trat4")));teste
friedman.test(teste)
pvalor=friedman.test(teste)$p.value
alfa=0.05
resposta<-function(pvalor) {
if(pvalor<alfa) { cat("pvalor é ",pvalor, " logo é pequeno >>> rejeita Ho ","\n") }
else if(pvalor>=alfa) {
cat("pvalor é ",pvalor, " logo é grande >>> NÃO rejeita Ho ","\n") }
}
resposta(pvalor)
#COMPARAÇÕES MULTIPLAS
resposta=c(3640,4890,4800,4460,4200,4550,5320,5500,4700,6020,5250,5580,5300,5900,5150,5560)
tratamento = rep(1:4, c(4,4,4,4))
bloco=rep(1:4,4)
require(agricolae)
comparison<-friedman(bloco,tratamento,resposta,alpha=0.05, group=TRUE,main="Teste de Friedman")
Alexandrov