[R-br] Como calcular o rank dos dados para uma análise não paramétrica

Daniel Marcelino dmsilva.br em gmail.com
Quarta Outubro 7 00:20:57 BRT 2015


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 em 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-for.html
>
>
>
>
> Em Segunda-feira, 5 de Outubro de 2015 11:46, Simone D. Sartorio
> <sisartorio em 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 em 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 em 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.


Mais detalhes sobre a lista de discussão R-br