<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body 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 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="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">Warning messages:
</span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">1: </span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">In chisq.test(table(emp), table(est), correct = TRUE) :
</span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);"> </span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);"> Aproximação do qui-quadrado pode estar incorreta
</span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">2: </span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">In densfun(x, parm[1], parm[2], ...) :</span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);"> NaNs produzidos
</span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">3: </span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">In chisq.test(table(emp), table(est), correct = TRUE) :
</span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);"> </span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);"> Aproximação do qui-quadrado pode estar incorreta
</span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">4: </span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);">In chisq.test(table(emp), table(est), correct = TRUE) :
</span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);"> </span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);"> Aproximação do qui-quadrado pode estar incorreta
</span></span><span><span class="GNKRCKGCASB ace_constant" style="color: rgb(197, 6, 11);"></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>
<div class="moz-cite-prefix">Em 06/03/2019 19:41, Cesar Rabak
escreveu:<br>
</div>
<blockquote type="cite"
cite="mid:CAKrF98k=16MmZiVRBXdPbgRD38=8s3qTZgtZdVruhuc8+0boiQ@mail.gmail.com">
<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>
</body>
</html>