Segue abaixo um CMR, espero que ajude.
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)
##########CASO FOR NECESSÁRIO##############
#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")
De: Fernanda Balbino <ferbalbino@yahoo.com.br>
Para: Lista R <r-br@listas.c3sl.ufpr.br>
Enviadas: Terça-feira, 22 de Novembro de 2011 15:57
Assunto: [R-br] Teste de Friedman
Olá, boa tarde!
Estou tentando fazer o teste de friedman e utilizei duas ferramentas que o R oferece. Porém me dá resultados diferentes. Utilizei no pacote pgirmess o friedmanmc para comparação multipa e o R me deu o resultado em que em alguns casos não há diferença significativa. E utilizei do pacote agricolae o friedman que me dá como resposta que somente existem
diferenças significativas.
Alguém já utilizou algumas destas funções ou conhece para me explicar por que da diferença? A função friedmanmcdo pacote pgirmessme parece ser mais comleta.
Muito obrigada.
Fernanda Balbino.
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-brLeia o guia de postagem (
http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.