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

Daniel Marcelino dmsilva.br em gmail.com
Quinta Outubro 8 16:06:34 BRT 2015


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 em 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 em 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 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