[R-br] Teste de Qui-quadrado aplicado para a distribuição binomial

ASANTOS alexandre.santos em cas.ifmt.edu.br
Ter Mar 5 17:06:59 -03 2019


Prezados Membros,

     Estou tentando realizar um teste de aderência (Qui-quadrado) 
aplicado a dados com distribuição binomial negativa sem sucesso. Em meu CRM:

#Pacote MASS
library(MASS)

#Crio um banco de dados simulado
Features<-sort(rep(1:3,100))
Resp1<-rnbinom(200,0.75,0.10) ##60 dados com distribuição binomial negativa
Resp2<-rpois(100,10) ##30 dados com distribuição de Poisson
Resp<-c(Resp1,Resp2)
d<-as.data.frame(cbind(Features,Resp))

#Criando um loop para realizar o teste de Chi quadrado para os três 
conjuntos de dados - Feature 1, 2 e 3
uniq <- unique(unlist(d$Features))
res<-NULL
for (i in 1:length(uniq)){
     data_1 <- subset(d, Features == uniq[i])
     k <- fitdistr(data_1$Resp,"negative binomial") #Extração dos 
parâmetros média e mu da binomial negativa
     par <- k$estimate
     size <- par[1]#Parâmetro k
     mu <- par[2]#Média
     N <- length(data_1$Resp)
     est <-N*dnbinom(data_1$Resp,size=size,mu=mu)  ## Criando valores 
estimados
     fecdf <- ecdf(data_1$Resp) ###ecdf- Cria a distribuição cumulativa 
empírica
     knotsX <- knots(fecdf)
     emp <- fecdf(c(knotsX,Inf))  # Distribuição empírica
     testChi<-chisq.test(table(emp),table(est),correct=TRUE) #Teste de 
Qui-quadrado sobre a contagem dos dados empíricos versus estimados
     result<-ifelse(testChi$p.value>0.05,"Binomial Negative","Other 
distribution") ## Classifica os testes em binomial negativa ou outra 
distribuição
     res<-rbind(res,c(uniq[i],result))
     colnames(res)<-c("Features","Distribution")
}
res

Meu output é: Features Distribution [1,] "1" "Other distribution" [2,] 
"2" "Other distribution" [3,] "3" "Binomial Negative" Ou seja, dá o 
resultado errado, na verdade as Features 1 e 2 é que são binomial 
negativa e a Feature 3 que são dados com distribuição de Poisson que tem 
que dar "Other distribution". Alguém teria alguma sugestão para dar? 
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 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
OrcID: orcid.org/0000-0001-8232-6722   -   ResearcherID: A-5790-2016
Researchgate: www.researchgate.net/profile/Alexandre_Santos10
LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635
Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190305/f8adb555/attachment.html>


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