Como faço para descobrir o poder de um teste de hipótese
Olá, Estou gerando duas amostras aleatórias de poisson: amostra1<-rpois(500,lambda=2.5) amostra2<-rpois(500,lambda=3.5) e aplicando o teste de wilcoxon e o teste t para dados pareados. wilcox.test(amostra1,amostra2, paired=T) t.test(amostra1,amostra2, paired=T) e gostaria de obter o poder de cada um dos testes acima. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>
Esse post pode ser de grande valia: http://r.789695.n4.nabble.com/How-to-compute-the-power-of-a-wilcoxon-test-td... abs Bruno Barros @brunosoab Em 21 de fevereiro de 2014 23:15, Jodavid Ferreira <jodavid.arts@gmail.com>escreveu:
Olá, Estou gerando duas amostras aleatórias de poisson:
amostra1<-rpois(500,lambda=2.5) amostra2<-rpois(500,lambda=3.5)
e aplicando o teste de wilcoxon e o teste t para dados pareados.
wilcox.test(amostra1,amostra2, paired=T) t.test(amostra1,amostra2, paired=T)
e gostaria de obter o poder de cada um dos testes acima.
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>
_______________________________________________ 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.
Qual é epistemologia de fazer um teste de poder *depois* de ter feito a amostragem? Se você precisa definir um poder para determinar como vai fazer seu experimento (que para a estatística significa determinar o tamanho da amostra), isso deve ser feito a priori. (v. por exemplo, http://homepage.stat.uiowa.edu/~rlenth/Power/ em "Advice" o 1º *bullet*. Se você está curioso quanto ao Erro Tipo II do teste t quando a amostra vem de uma população com distribuição de Poisson você deveria fazer um grande número de vezes os testes tabelando quantas vezes o teste "acerta" ou não. *Caveat emptor*: para que isso seja realmente útil você teria que se certificar que o gerador de números aleatórios não tem nenhum "defeito" que enviese seus resultados, o ideal seria usar um baseado em um HW apropriado para tal ou usar uma tabela de números aleatórios seguramente livre de vícios. HTH -- Cesar Rabak 2014-02-21 23:15 GMT-03:00 Jodavid Ferreira <jodavid.arts@gmail.com>:
Olá, Estou gerando duas amostras aleatórias de poisson:
amostra1<-rpois(500,lambda=2.5) amostra2<-rpois(500,lambda=3.5)
e aplicando o teste de wilcoxon e o teste t para dados pareados.
wilcox.test(amostra1,amostra2, paired=T) t.test(amostra1,amostra2, paired=T)
e gostaria de obter o poder de cada um dos testes acima.
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>
_______________________________________________ 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 (3)
-
Bruno Barros -
Cesar Rabak -
Jodavid Ferreira