Como calcular o rank dos dados para uma análise não paramétrica
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.dadosa<- kruskal.test(Orquideas~Campo, orq.dados)summary(a)names(a)?kruskal.test #---------------------------------------------------------- desde já, muito obrigada!Simone
Não é só usar a função rank()? W.
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.dadosa<- 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.
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.
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.
Simone, você cehgou a testar a função que eu postei? Para min ela funciona melhor do que a solução apresentada no livro que você mencionou: dados = read.table(text="Male,LongDay,FoodIntake 1,1,69.50 1,0,70.00 1,1,71.10 1,1,72.10 1,0,72.30 1,0,72.90 1,1,73.20 1,0,73.30 0,1,75.50 0,1,76.30 0,1,78.10 0,0,80.90 0,1,81.20 0,0,82.40 0,0,83.00 0,0,88.20", sep=",", header=TRUE)
with(dados, srh(FoodIntake, as.factor(LongDay), as.factor(Male))) Df Sum Sq Mean Sq H p-value pf 1 25 25 1.1029 0.29362 sf 1 256 256 11.2941 0.00078 pf:sf 1 4 4 0.1765 0.67442
2015-10-07 9:59 GMT-03:00 Simone D. Sartorio <sisartorio@yahoo.com.br>:
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ços Simone
***---***---***---***---***---***---***---***---***---***---***---***---*** * 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.
participantes (4)
-
Antonio Moita -
Daniel Marcelino -
Simone D. Sartorio -
Walmes Zeviani