[R-br] Teste de distribuição de erros tabajara

ASANTOS alexandresantosbr em yahoo.com.br
Quarta Outubro 28 16:25:38 BRST 2015


Prezados listeiros,

       Estou tentando criar uma função para testar a distribuição de 
erros de normal, poisson e binomial negativa por enquanto, porém não sei 
o que estou fazendo de errado na minha função e solicito a ajuda de 
vocês para conseguir achar o erro, segue CRM:

## Teste de distribuição de erros tabajara

require(MASS)

BD<-rpois(1000,10) ### Valores a serem testados

Dist<-c("negative binomial","normal","Poisson")

for(typeD in 1:Dist){
k <- fitdistr(BD,typeD[Dist])
par <- k$estimate
size <- par[1]#Parametro
mu <- par[2]#Média
SD<-sd(BD)
N <- length(BD)

if (typeD[Dist]=='negative binomial'){       ### Para binomial negativa
est <-N*dnbinom(BD,size=size,mu=mu)  ## Estimativas
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica
knotsX <- knots(fecdf)
emp <- fecdf(c(knotsX,Inf))  # Empírica
chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado 
para valiar a significância da distribuição
}
if (typeD[Dist]=='normal'){          ### Para normal
est <-N*dnorm(BD,mean=mu, sd=SD)  ## Estimativas
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica
knotsX <- knots(fecdf)
emp <- fecdf(c(knotsX,Inf))  # Empírica
chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado 
para valiar a significância da distribuição
}
  if (typeD[Dist]=='Poisson'){        ### Para Poisson
est <-N*dpois(BD,lambda=mu)  ## Estimativas
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica
knotsX <- knots(fecdf)
emp <- fecdf(c(knotsX,Inf))  # Empírica
chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado 
para valiar a significância da distribuição
}
}
#

Obrigado,

-- 
======================================================================
Alexandre dos Santos
Proteção 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
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================



Mais detalhes sobre a lista de discussão R-br