[R-br] Testando frequências de distribuição normal usando o pacote MASS
ASANTOS
alexandresantosbr em yahoo.com.br
Segunda Fevereiro 18 12:38:49 BRT 2013
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 em yahoo.com.br
alexandre.santos em cas.ifmt.edu.br
======================================================================
Mais detalhes sobre a lista de discussão R-br