<div dir="ltr">Alexandre:<div><br></div><div>Por favor explique-no o seu raciocínio para fazer este teste aqui:<br><br><i><span style="font-family:Arial">> result<-ifelse(testChi$p.</span></i><span style="font-family:Arial"><i>value>0.05,"Binomial Negative","Other distribution")</i> ## Classifica os testes em binomial negativa ou outra distribuição</span> </div><div><br></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Mar 5, 2019 at 5:07 PM ASANTOS por (R-br) <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div 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 class="gmail-m_-2788180973372379985GNKRCKGCGSB" id="gmail-m_-2788180973372379985rstudio_console_output" style="font-family:"Lucida Console";outline:none;border:none;word-break:break-all;margin:0px;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);font-size:10pt;white-space:pre-wrap"><span><span class="gmail-m_-2788180973372379985GNKRCKGCGSB" style="outline:none;border:none;word-break:break-all;margin:0px;white-space:pre-wrap"><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="gmail-m_-2788180973372379985GNKRCKGCGSB" style="outline:none;border:none;word-break:break-all;margin:0px;white-space:pre-wrap"><font face="Arial"><span><span class="gmail-m_-2788180973372379985GNKRCKGCGSB" style="outline:none;border:none;word-break:break-all;margin:0px;white-space:pre-wrap">Other distribution</span></span>".
Alguém teria alguma sugestão para dar?
Obrigado</font>
</span></span></pre>
<pre class="gmail-m_-2788180973372379985moz-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="gmail-m_-2788180973372379985moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="gmail-m_-2788180973372379985moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="gmail-m_-2788180973372379985moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/1360403201088680</a>
OrcID: <a href="http://orcid.org/0000-0001-8232-6722" target="_blank">orcid.org/0000-0001-8232-6722</a> - ResearcherID: A-5790-2016
Researchgate: <a class="gmail-m_-2788180973372379985moz-txt-link-abbreviated" href="http://www.researchgate.net/profile/Alexandre_Santos10" target="_blank">www.researchgate.net/profile/Alexandre_Santos10</a>
LinkedIn: <a href="http://br.linkedin.com/in/alexandre-dos-santos-87961635" target="_blank">br.linkedin.com/in/alexandre-dos-santos-87961635</a>
Mendeley:<a href="http://www.mendeley.com/profiles/alexandre-dos-santos6/" target="_blank">www.mendeley.com/profiles/alexandre-dos-santos6/</a>
======================================================================</pre>
</div>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</blockquote></div>