22 Out
2012
22 Out
'12
20:08
Esqueci de postar na lista
> median.test<-function(y1,y2, na.rm = TRUE){
+ z <- c(y1,y2)
+ g <- rep(1:2, c(length(y1), length(y2)))
+ m <- median(z, na.rm)
+ fisher.test(z<m,g)$p.value
+ }
> turma1=c(72,79,93,91,70,95,82,80,74,86, NA)
> turma2=c(78,66,65,84,69,73,71,75,68,90,76)
> median.test(turma1,turma2)
[1] 0.1983806
[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com
Em 22/10/2012 17:40, Zé Augusto escreveu:
> Luiz Henrique
>
> Seu script retorna o erro
> "Erro em median.default(dados, na.rm = TRUE) : need numeric data"
> no R2.15.0 64bits.
>
> #####################################################
>
> Leonardo
>
> Desculpe a minha ignorância, mas seu script, para
> y1 = turma1
> y2 = turma2
> não retorna nada.
>
> Grato por qualquer informação.
> José Augusto
>
>
> Em 22 de outubro de 2012 17:05, Leonard de Assis
> <assis.leonard@gmail.com> escreveu:
> median.test<-function(y1,y2){
> z<-c(y1,y2)
> g <- rep(1:2, c(length(y1),length(y2)))
> m<-median(z)
> fisher.test(z<m,g)$p.value
> }
>> Fonte: Lista R
>>
>> []s
>> Leonard de Assis
>> assis <dot> leonard <at> gmail <dot> com
>>
>> Em 22/10/2012 16:36, Luiz Henrique da Conceicao Leal escreveu:
>>
>> # Dados
>> turma1=c(72,79,93,91,70,95,82,80,74,86,NA)
>> turma2=c(78,66,65,84,69,73,71,75,68,90,76)
>> Notas=data.frame(turma1,turma2)
>>
>>
>> # Teste da mediana
>> teste.da.mediana=function(dados,alfa=0.05){
>> M=median(dados,na.rm=TRUE)
>> Tabela=table(dados[,1]<=M)
>> for(j in 2:ncol(dados)) Tabela=rbind(Tabela,table(dados[,j]<=M))
>> colnames(Tabela)=c("Acima da mediana","Abaixo da mediana")
>> rownames(Tabela)=colnames(dados)
>> if(length(table(Tabela<=5))==1) {Teste=chisq.test(Tabela,correct = FALSE)}
>> else Teste=fisher.test(Tabela)
>> Conclusão=if(Teste$p.value<alfa) {"Rejeita-se a hipótese de igualdade das
>> medianas"} else "Não se pode rejeitar a hipótese de igualdade das medianas"
>> Medianas=apply(Notas,2,median,na.rm=TRUE)
>> list(Tabela=Tabela,Teste=Teste,Conclusão=Conclusão,Medianas=Medianas)
>> }
>>
>> teste.da.mediana(Notas,alfa=0.05)
>>
>>
>> Em 22 de outubro de 2012 15:54, Márcia Varela da Silva
>> <marcinhavarela@hotmail.com> escreveu:
>>> Como fazer o "Teste da Mediana" no R?
>>>
>>>
>>>
>>> _______________________________________________
>>> 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.
>>
>>
>>
>> _______________________________________________
>> 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.
>>
>>
>>
>> _______________________________________________
>> 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.