Pessoal, obrigada pela ajuda!O rank parece que vai funcionar...estou em teste aqui. Consegui um livro que descreve como fazer o teste no SPSS, Minitab e Excel file:///C:/Users/Simone/Desktop/Dytham%202011%20statistics.pdf ...mas no R, ainda não encontrei nada pronto...obrigada mais uma vez!abraçosSimone ***---***---***---***---***---***---***---***---***---***---***---***---**** Simone Daniela Sartorio* Professora Adjunta II no CCA/UFSCar, Araras/SP.* Doutora e Mestre em Estatística e Experimentação Agronômica - ESALQ/USP.* Licenciada em Matemática - UNESP/Rio Claro.***---***---***---***---***---***---***---***---***---***---***---***---*** Tenha um bom dia! ;) Em Quarta-feira, 7 de Outubro de 2015 0:20, Daniel Marcelino <dmsilva.br@gmail.com> escreveu: Se você não uma função no R para realizar o teste, isso por si só já pode ser considerado um sinal de que o teste não é muito popular. Veja o exemplo no livro: # Sokal RR, Rohlf FJ, 1995. Biometry, 3rd ed. New York: Freeman. rat <- c(709,679,699,657,594,677,592,538,476,508,505,539) fat <- gl(2,6,labels=c("Fresh","Rancid")) sex <- gl(2,3,12,labels=c("man","woman")) srh(rat,sex,fat) Função para gerar os resultados: srh <-function(r, pf , sf, n=3){ abn = (nlevels(pf)*nlevels(sf)*n) lm_rank <- lm(rank(r)~pf*sf) # the variance must be correct if there are ties aov_lm <- anova(lm_rank) ans <- aov_lm[1:3,1:3] ans[,4] <- ans[,2]/(length(r)*(length(r)+1)/abn) ans[,5] <- (1-pchisq(ans[,4],ans[,1])) colnames(ans)[4:5] <- c("H","p-value") ans } 2015-10-06 9:30 GMT-03:00 Antonio Moita <awmoita@yahoo.com.br>:
Caros,
Dei uma uma olhada no google e encontrei este blog que tem uma forma em seis passos para calcular. Este assunto me interessa, qdo tiver mais tempo eu vou dar "assuntar" melhor.
No R tem o pacote ??? que faz muito mais do que voce pretende realizar, se eu nao me engano é o nparLD
Att, Moita http://blog.excelmasterseries.com/2014/05/scheirer-ray-hope-test-alternative...
Em Segunda-feira, 5 de Outubro de 2015 11:46, Simone D. Sartorio <sisartorio@yahoo.com.br> escreveu:
Pessoal, bom dia!
Eu tentando aplicar um teste não paramétrico, em particular o "teste Scheirer-Ray-Hare". Ele é uma extensão do Kruskal-Wallis para experimento com 2 fatores (em esquema fatorial). Já fucei no R e não encontrei essa função pronta. Logo, preciso obter apenas o rank dos dados para poder aplicar o teste, mas não obtive sucesso ao procurar como obter esse rank com a função kruskal.test...alguém poderia me ajudar?
Aqui vai um exemplo fictício para poderem me ajudar...lembrando que esse exemplo é de um fator apenas. O que eu quero é apenas a obtenção dos ranks.
#---------------------------------------------------------- A<-c(27,14,8,18,7) B<-c(48,18,32,51,22) C<-c(11,0,3,15,8) D<-c(44,72,81,55,39)
orq.dados<-data.frame(Campo<-gl(4,5), Orquideas<-c(A,B,C,D)); orq.dados a<- kruskal.test(Orquideas~Campo, orq.dados) summary(a) names(a) ?kruskal.test #----------------------------------------------------------
desde já, muito obrigada! Simone
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne� c�igo m�imo reproduz�el.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.