Muito bom, Benilton! Ainda prefiro, pra esses casos, um bom qq-plot, apesar de usar a estatística do K-S como critério de parada em alguns métodos.
Em 16 de novembro de 2011 09:16, Benilton Carvalho < [email protected]> escreveu: > O teste de KS para 1 amostra (qdo vc usa punif) mede a distancia entre > a distribuicao empirica da sua amostra 'x' e a distribuicao teorica da > uniforme (via punif). O teste de KS para 2 amostras (seu segundo > caso), mede a distancia entre as duas distribuicoes empiricas (note a > diferenca). > > Por sinal, note que, ao rejeitar a hipotese nula, vc nao pode dizer > que ela e' falsa... Tudo o que vc pode afirmar e' que, com os dados > que voce possui, vc nao tem evidencias que suportem a hipotese nula. > E' uma diferenca sutil, mas essencial para todos os que usam testes de > hipotese. > > Para ilustracao: gere trocentas amostras com distribuicao normal... > faca o teste de hipotese (usando KS) de normalidade para cada uma das > amostras... Se vc preferir o ponto de corte de 5% para o p-valor e > usa-lo como criterio de rejeicao, vc observara' que vc rejeitara' > cerca de 5% dos conjuntos de dados e ainda assim a hipotese nula eh > verdadeira (afinal todos os dados sao normais).... > > Para facilitar sua vida no "dever de casa": > > prop.table(table(replicate(10000, ks.test(rnorm(10), pnorm)$p.value < > .05))) > > b > > 2011/11/16 julio cesar araujo <[email protected]>: > > Kolmogorov - Smirnov test no R, > > > > Pessoal, quero verificar se um conjunto de dados tem distribuição > uniforme > > (0,1). Daí através da função ks.test, uso, por exemplo, o comando > > ks.test(x, "punif"). Para um conjutno de dados que possuo esse comando > > apresenta valor abaixo de 0.5 o que significa que meus dados não possuem > a > > distribuição uniforme. Só que essa mesma função pode ser utilizada para > > comparar se duas amostras possuem a mesma distribuição, por exemplo > > ks.test(x, y), daí criei uma variável y, através de um runif (uma > uniforme > > randomica) e testei com a X e resultado foi aceito (acima de 0.20, no meu > > caso). Como isso é possível???? seria x então uniforme ou não, como é > > possível quando utilizado "punif" dar que não oe quando comparado a uma > > variável uniforme criada aleatorimente por mim acusar ter a mesma > > distribuição?? > > > > Desde já agradeço e peço desculpas pela "trivialidade" da pergunta. > > > > _______________________________________________ > > R-br mailing list > > [email protected] > > 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 > [email protected] > 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. > -- Fernando A.B. Colugnati
_______________________________________________ R-br mailing list [email protected] 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.
