<div dir="auto"><i style="font-family:verdana,geneva,helvetica,arial,sans-serif;font-size:13.44px;background-color:rgb(255,255,255)"><span style="font-family:'arial'">> result<-ifelse(testChi$p.</span></i><span style="font-size:13.44px;background-color:rgb(255,255,255);font-family:arial"><i>value>0.05,"Binomial Negative","Other distribution")</i> ## Classifica os testes em binomial negativa ou outra distribuição</span><span style="font-family:verdana,geneva,helvetica,arial,sans-serif;font-size:13.44px;background-color:rgb(255,255,255)"> </span><div dir="auto"><span style="font-family:verdana,geneva,helvetica,arial,sans-serif;font-size:13.44px;background-color:rgb(255,255,255)"><br></span></div><div dir="auto"><span style="font-family:verdana,geneva,helvetica,arial,sans-serif;font-size:13.44px;background-color:rgb(255,255,255)">Se a distribuição estimada table(est) for igual a empírica table (emp) pelo teste de independência de Qui-quadrado, ou seja, </span><i style="font-family:verdana,geneva,helvetica,arial,sans-serif;font-size:13.44px;background-color:rgb(255,255,255)"><span style="font-family:'arial'">testChi$p.</span></i><span style="font-size:13.44px;background-color:rgb(255,255,255);font-family:arial"><i>value>0.05, então os dados tem distribuição binomial negativa, caso contrário, </i></span><i style="font-family:verdana,geneva,helvetica,arial,sans-serif;font-size:13.44px;background-color:rgb(255,255,255)"><span style="font-family:'arial'">testChi$p.</span></i><span style="font-size:13.44px;background-color:rgb(255,255,255);font-family:arial"><i>value<0.05, então os dados tem outra distribuição de probabilidade.</i></span></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Em qua, 6 de mar de 2019 22:24, ASANTOS <<a href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>> escreveu:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF">
<p>Boa noite Cesar,</p>
<p> Obrigado pela ajuda. Para mim, quando eu rodo aparece também
alguns alertas:</p>
<pre class="m_-5677666329908503353GNKRCKGCGSB" id="m_-5677666329908503353rstudio_console_output" style="font-family:'Lucida Console';font-size:10pt!important;outline:none;border:none;word-break:break-all;margin:0px;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;text-align:-webkit-left;text-indent:0px;text-transform:none;word-spacing:0px;background-color:rgb(255,255,255)"><span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">Warning messages:
</span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">1: </span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">In chisq.test(table(emp), table(est), correct = TRUE) :
</span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)"> </span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)"> Aproximação do qui-quadrado pode estar incorreta
</span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">2: </span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">In densfun(x, parm[1], parm[2], ...) :</span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)"> NaNs produzidos
</span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">3: </span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">In chisq.test(table(emp), table(est), correct = TRUE) :
</span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)"> </span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)"> Aproximação do qui-quadrado pode estar incorreta
</span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">4: </span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)">In chisq.test(table(emp), table(est), correct = TRUE) :
</span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)"> </span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)"> Aproximação do qui-quadrado pode estar incorreta
</span></span><span><span class="m_-5677666329908503353GNKRCKGCASB m_-5677666329908503353ace_constant" style="color:rgb(197,6,11)"></span></span></pre>
<pre class="m_-5677666329908503353moz-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="m_-5677666329908503353moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br" target="_blank" rel="noreferrer">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="m_-5677666329908503353moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank" rel="noreferrer">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="m_-5677666329908503353moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680" target="_blank" rel="noreferrer">http://lattes.cnpq.br/1360403201088680</a>
OrcID: <a href="http://orcid.org/0000-0001-8232-6722" target="_blank" rel="noreferrer">orcid.org/0000-0001-8232-6722</a> - ResearcherID: A-5790-2016
Researchgate: <a class="m_-5677666329908503353moz-txt-link-abbreviated" href="http://www.researchgate.net/profile/Alexandre_Santos10" target="_blank" rel="noreferrer">www.researchgate.net/profile/Alexandre_Santos10</a>
LinkedIn: <a href="http://br.linkedin.com/in/alexandre-dos-santos-87961635" target="_blank" rel="noreferrer">br.linkedin.com/in/alexandre-dos-santos-87961635</a>
Mendeley:<a href="http://www.mendeley.com/profiles/alexandre-dos-santos6/" target="_blank" rel="noreferrer">www.mendeley.com/profiles/alexandre-dos-santos6/</a>
======================================================================</pre>
<div class="m_-5677666329908503353moz-cite-prefix">Em 06/03/2019 19:41, Cesar Rabak
escreveu:<br>
</div>
<blockquote type="cite">
<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>
</blockquote>
</div>
</blockquote></div>