[R-br] Teste no R

[Eliardo Costa]

A função do Leonardo não retornava o resultado devido ao NA em turma1,
basta acrescentar na.rm=T em median como abaixo.

median.test<-function(y1,y2){
z<-c(y1,y2)
g <- rep(1:2, c(length(y1),length(y2)))
m<-median(z,na.rm=T)
  fisher.test(z<m,g)$p.value
 }

Att, Eliardo.

Em 22 de outubro de 2012 17:40, Zé Augusto <jabnaz em gmail.com> 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 em 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 em hotmail.com> escreveu:
> >>
> >> Como fazer o "Teste da Mediana" no R?
> >>
> >>
> >>
> >> _______________________________________________
> >> 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.
> >
> >
> >
> >
> > _______________________________________________
> > 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.
> >
> >
> >
> > _______________________________________________
> > 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.
> _______________________________________________
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121022/e76dc9f9/attachment.html>