A função do Leonardo não retornava o resultado devido ao NA em turma1, basta acrescentar na.rm=T em median como abaixo.<br><br>
median.test<-function(y1,y2){<br>
z<-c(y1,y2)<br>
g <- rep(1:2, c(length(y1),length(y2)))<br>
m<-median(z,na.rm=T)<br>
  fisher.test(z<m,g)$p.value<br>
 }<br><br>Att, Eliardo.<br><br><div class="gmail_quote">Em 22 de outubro de 2012 17:40, Zé Augusto <span dir="ltr"><<a href="mailto:jabnaz@gmail.com" target="_blank">jabnaz@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Luiz Henrique<br>
<br>
    Seu script retorna o erro<br>
"Erro em median.default(dados, na.rm = TRUE) : need numeric data"<br>
no R2.15.0 64bits.<br>
<br>
#####################################################<br>
<br>
Leonardo<br>
<br>
    Desculpe a minha ignorância, mas seu script, para<br>
      y1 = turma1<br>
      y2 = turma2<br>
não retorna nada.<br>
<br>
Grato por qualquer informação.<br>
José Augusto<br>
<br>
<br>
Em 22 de outubro de 2012 17:05, Leonard de Assis<br>
<<a href="mailto:assis.leonard@gmail.com">assis.leonard@gmail.com</a>> escreveu:<br>
<div class="HOEnZb"><div class="h5">><br>
median.test<-function(y1,y2){<br>
z<-c(y1,y2)<br>
g <- rep(1:2, c(length(y1),length(y2)))<br>
m<-median(z)<br>
  fisher.test(z<m,g)$p.value<br>
 }<br>
><br>
> Fonte: Lista R<br>
><br>
> []s<br>
> Leonard de Assis<br>
> assis <dot> leonard <at> gmail <dot> com<br>
><br>
> Em 22/10/2012 16:36, Luiz Henrique da Conceicao Leal escreveu:<br>
><br>
> # Dados<br>
>   turma1=c(72,79,93,91,70,95,82,80,74,86,NA)<br>
>   turma2=c(78,66,65,84,69,73,71,75,68,90,76)<br>
>   Notas=data.frame(turma1,turma2)<br>
><br>
><br>
> # Teste da mediana<br>
>   teste.da.mediana=function(dados,alfa=0.05){<br>
>   M=median(dados,na.rm=TRUE)<br>
>   Tabela=table(dados[,1]<=M)<br>
>   for(j in 2:ncol(dados)) Tabela=rbind(Tabela,table(dados[,j]<=M))<br>
>   colnames(Tabela)=c("Acima da mediana","Abaixo da mediana")<br>
>   rownames(Tabela)=colnames(dados)<br>
>   if(length(table(Tabela<=5))==1) {Teste=chisq.test(Tabela,correct = FALSE)}<br>
> else Teste=fisher.test(Tabela)<br>
>   Conclusão=if(Teste$p.value<alfa) {"Rejeita-se a hipótese de igualdade das<br>
> medianas"} else "Não se pode rejeitar a hipótese de igualdade das medianas"<br>
>   Medianas=apply(Notas,2,median,na.rm=TRUE)<br>
>   list(Tabela=Tabela,Teste=Teste,Conclusão=Conclusão,Medianas=Medianas)<br>
>   }<br>
><br>
>   teste.da.mediana(Notas,alfa=0.05)<br>
><br>
><br>
> Em 22 de outubro de 2012 15:54, Márcia Varela da Silva<br>
> <<a href="mailto:marcinhavarela@hotmail.com">marcinhavarela@hotmail.com</a>> escreveu:<br>
>><br>
>> Como fazer o "Teste da Mediana" no R?<br>
>><br>
>><br>
>><br>
>> _______________________________________________<br>
>> R-br mailing list<br>
>> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
>> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>> código mínimo reproduzível.<br>
><br>
><br>
><br>
><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br>
> mínimo reproduzível.<br>
><br>
><br>
><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br>
> mínimo reproduzível.<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</div></div></blockquote></div><br>