<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p><font size="-1" face="Arial">Prezados Membros,</font></p>
<p><font size="-1" face="Arial"> Estou tentando realizar um teste
de aderência (Qui-quadrado) aplicado a dados com distribuição
binomial negativa sem sucesso. Em meu CRM:<br>
</font></p>
<p><font size="-1" face="Arial">#Pacote MASS<br>
library(MASS)</font></p>
<p><font size="-1" face="Arial">#Crio um banco de dados simulado<br>
Features<-sort(rep(1:3,100))<br>
Resp1<-rnbinom(200,0.75,0.10) ##60 dados com distribuição
binomial negativa<br>
Resp2<-rpois(100,10) ##30 dados com distribuição de Poisson<br>
Resp<-c(Resp1,Resp2)<br>
d<-as.data.frame(cbind(Features,Resp))<br>
<br>
#Criando um loop para realizar o teste de Chi quadrado para os
três conjuntos de dados - Feature 1, 2 e 3 <br>
uniq <- unique(unlist(d$Features))<br>
res<-NULL<br>
for (i in 1:length(uniq)){<br>
data_1 <- subset(d, Features == uniq[i])<br>
k <- fitdistr(data_1$Resp,"negative binomial") #Extração
dos parâmetros média e mu da binomial negativa<br>
par <- k$estimate<br>
size <- par[1]#Parâmetro k<br>
mu <- par[2]#Média<br>
N <- length(data_1$Resp)<br>
est <-N*dnbinom(data_1$Resp,size=size,mu=mu) ## Criando
valores estimados<br>
fecdf <- ecdf(data_1$Resp) ###ecdf- Cria a distribuição
cumulativa empírica<br>
knotsX <- knots(fecdf)<br>
emp <- fecdf(c(knotsX,Inf)) # Distribuição empírica<br>
testChi<-chisq.test(table(emp),table(est),correct=TRUE)
#Teste de Qui-quadrado sobre a contagem dos dados empíricos
versus estimados<br>
result<-ifelse(testChi$p.value>0.05,"Binomial
Negative","Other distribution") ## Classifica os testes em
binomial negativa ou outra distribuição<br>
res<-rbind(res,c(uniq[i],result))<br>
colnames(res)<-c("Features","Distribution")<br>
} <br>
res</font></p>
<pre tabindex="0" class="GNKRCKGCGSB" id="rstudio_console_output" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline: none; border: none; word-break: break-all; margin: 0px; -webkit-user-select: text; white-space: pre-wrap !important; line-height: 15px; color: rgb(0, 0, 0); font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: -webkit-left; text-indent: 0px; text-transform: none; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: rgb(255, 255, 255);"><span><span class="GNKRCKGCGSB" style="outline: none; border: none; word-break: break-all; margin: 0px; -webkit-user-select: text; white-space: pre-wrap !important;"><font face="Arial">Meu output é:
</font> Features Distribution
[1,] "1" "Other distribution"
[2,] "2" "Other distribution"
[3,] "3" "Binomial Negative"
<font face="Arial">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 "</font></span></span><span><span class="GNKRCKGCGSB" style="outline: none; border: none; word-break: break-all; margin: 0px; -webkit-user-select: text; white-space: pre-wrap !important;"><font face="Arial"><span><span class="GNKRCKGCGSB" style="outline: none; border: none; word-break: break-all; margin: 0px; -webkit-user-select: text; white-space: pre-wrap !important;">Other distribution</span></span>".
Alguém teria alguma sugestão para dar?
Obrigado</font>
</span></span></pre>
<pre class="moz-signature" cols="72">--
======================================================================
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)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016
Researchgate: <a class="moz-txt-link-abbreviated" href="http://www.researchgate.net/profile/Alexandre_Santos10">www.researchgate.net/profile/Alexandre_Santos10</a>
LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635
Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/
======================================================================</pre>
</body>
</html>