# 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.
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 <mailto:marcinhavarela@hotmail.com>> escreveu:
Como fazer o "Teste da Mediana" no R?
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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.
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.
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@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@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.
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.
participantes (5)
-
Eliardo Costa -
Leonard de Assis -
Luiz Henrique da Conceicao Leal -
Márcia Varela da Silva -
Zé Augusto