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

Cesar Rabak cesar.rabak em gmail.com
Qua Mar 6 20:41:09 -03 2019


Na minha tentativa de rodar s/exemplo antes da linha
> res

As seguintes msgs de alerta são impressas:
Warning messages:
1: In densfun(x, parm[1], parm[2], ...) : NaNs produzidos
2: In chisq.test(table(emp), table(est), correct = TRUE) :
  Aproximação do qui-quadrado pode estar incorreta
3: In densfun(x, parm[1], parm[2], ...) : NaNs produzidos
4: In densfun(x, parm[1], parm[2], ...) : NaNs produzidos
5: In chisq.test(table(emp), table(est), correct = TRUE) :
  Aproximação do qui-quadrado pode estar incorreta
6: In chisq.test(table(emp), table(est), correct = TRUE) :
  Aproximação do qui-quadrado pode estar incorreta

CONFERE?



On Tue, Mar 5, 2019 at 5:07 PM ASANTOS por (R-br) <r-br em listas.c3sl.ufpr.br>
wrote:

> 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/
> ======================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190306/a49abf5d/attachment.html>


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