Testando frequências de distribuição normal usando o pacote MASS

Bom tarde Pessoal, Estou tentando testar frequências de distribuição normal para um vetor com 100,000 valores com o pacote MASS, o problema é que quando faço o teste de Qui quadrado entre a distribuição empírica o teste não funciona: Erro em chisq.test(x = table(emp), y = table(est)) : 'x' e 'y' devem ter ao menos 2 níveis Quando faço para binomial negativa funciona, mas para distribuição normal não. Alguém teria uma outra solução, segue CRM: require(MASS) res<-rnorm(100000)### Com 100,000 valores não posso usar shapiro.test ( k <- fitdistr(res,"normal")) par <- k$estimate (md <- par[1])#Média (ep <- par[2])#sd (N <- length(res)) est <-N*dnorm(res,mean=md,sd=ep) ## Estimativas ##Teste Chi - Quadrado para valiar a significância da distribuição fecdf <- ecdf(res) ###ecdf- Função de distribuição cumulativa empírica knotsX <- knots(fecdf) emp <- fecdf(c(knotsX,Inf)) chisq.test(x=table(emp),y=table(est))## Teste de Chi quadrado Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================

Por que não fazer o teste Kolmogorov-Smirnov, ks,test()? À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

Obrigado Walmes, Resolve o problema: a<-rnorm(100000,mean=20,sd=5)## Geração de números com dist. normal res.a<-mean(a)-a ## Resíduo de a b<-rnorm(100000,mean=20,sd=5) ## Geração de números com dist. normal res.b<-mean(b)-b ## Resíduo de b ks.test(res.a,res.b) ##Teste Mas queria ter feito no MASS para testar várias distribuições, inclusive a normal, Abraço, Alexandre Em 19/02/2013 14:56, Walmes Zeviani escreveu:
Por que não fazer o teste Kolmogorov-Smirnov, ks,test()?
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br <mailto:walmes@ufpr.br> skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================
_______________________________________________ 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.
-- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================

Na ks.test você testar qualquer par de amostras e também uma amostra com relação à uma distribuição de teórica. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
ASANTOS
-
Walmes Zeviani