<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:10pt"><div style="font-family: arial, helvetica, sans-serif; font-size: 10pt;"><span>Obrigado Vinícius,</span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 13px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><span><br></span></div><div style="background-color: transparent;"><span style="color: rgb(0, 0, 0); font-family: arial, helvetica, sans-serif; font-size: 13px; font-style: normal;"> Mas o Benilton já havia chamado a atenção para uma questão crucial: "</span><span><font color="#454545" face="Arial, Helvetica, sans-serif"><span><span style="font-size: 12px;">regressão logística modela a probabilidade de sucesso dado um conjunto de covariaveis</span></span></font></span><span style="color: rgb(0, 0, 0); font-family: arial, helvetica, sans-serif;
font-size: 13px; font-style: normal; background-color: transparent;">", então no modelo original que estava trabalhando existem erros conceituais, uma aproximação do que eu queria que era a obtenção dos coeficientes do intercepto, linear, quadrático e cúbico, ficou:</span></div><div style="background-color: transparent; color: rgb(69, 69, 69); font-size: 12px; font-family: Arial, Helvetica, sans-serif; font-style: normal;"><span style="color: rgb(0, 0, 0); font-family: arial, helvetica, sans-serif; font-size: 13px; font-style: normal; background-color: transparent;"><br></span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 10pt;"></div><div><font size="2"> pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2</font></div><div><font size="2">,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26</font></div><div><font size="2">,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de mumias</font></div><div><font
size="2">dens<-sort(rep(2^(2:6),10))#### Densidade de pulgoes oferecidos</font></div><div><font size="2">##</font></div><div><font size="2">##</font><span style="color: rgb(69, 69, 69); font-family: Arial, Helvetica, sans-serif; font-size: 12px;">Prob(Sucesso) = exp(P0+P1*N+P2*N^2+P3*N^3)/(1+exp(P0+P1*N+P2*N^2+P3*N^3))</span></div><div><font size="2">naoPred<-dens-pred</font></div><div><font size="2">#</font></div><div><font size="2">p.model2<-glm(cbind(pred, naoPred)~poly(dens, 3, raw=TRUE), family=binomial)</font></div><div><font size="2">summary(p.model2)</font></div><div><font size="2">#</font></div><div><font size="2"><br></font></div><div><font size="2">Abs</font></div><div><font size="2"><br></font></div><div><font size="2">Alexandre</font></div> <div style="font-family: arial, helvetica, sans-serif; font-size: 10pt;"> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"> <div dir="ltr"> <font size="2"
face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Vinicius Brito Rocha <viniciusbritor@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br; Alexandre Santos <alexandresantosbr@yahoo.com.br> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Terça-feira, 20 de Novembro de 2012 21:32<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] Dúvida em summary em regressão logística<br> </font> </div> <br><div id="yiv1749391770">Alexandre,<div><br></div><div>acho que vc chamou o modelo de forma errada.</div><div><br></div><div>deveria ser glm(cbind(evento,1-evento)~VA1+VA2+...VAp,data=dados,family=binomial(link = "logit"))</div><div><br>
</div><div>normalmente nos softwares o coeficiente resposta é o do modelo linear, como vc esta usando a funçao logit</div><div><br></div><div>log(P(X=evento)/P(X=nao evento))= beta0+beta1*X1+beta2*X2+...+betap*Xp</div><div>
<br></div><div><br></div><div>Abs.</div><div><br></div><div>Mas o Walmes provavelmente vai poder te responder melhor.</div><div><br></div><div>Abs</div><div><br></div><div>Vinicius Brito Rocha<br><br><div class="yiv1749391770gmail_quote">
Em 20 de novembro de 2012 16:49, Alexandre Santos <span dir="ltr"><<a rel="nofollow" ymailto="mailto:alexandresantosbr@yahoo.com.br" target="_blank" href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a>></span> escreveu:<br><blockquote class="yiv1749391770gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div style="font-size: 10pt; font-family: arial, helvetica, sans-serif;"><div align="center"><div style="text-align:left;"><span style="background-color:transparent;">Boa tarde Pessoal,</span></div><div style="text-align:left;">
<span style="background-color:transparent;"> Estou ajustando uma regressão logística e me deparei com a seguinte dúvida:</span></div><div style="text-align:left;"><span style="background-color:transparent;"> </span></div>
<div style=""><font face="arial, helvetica, sans-serif">###################################################################################################</font></div><div style=""><font face="arial, helvetica, sans-serif">#Regressão logística entre a proporção de herbívoros predados e a densidade
de herbívoros oferecidos</font></div><div style=""><font face="arial, helvetica, sans-serif">#</font></div><div style=""><font face="arial, helvetica, sans-serif">pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2</font></div>
<div style=""><font face="arial, helvetica, sans-serif">,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26</font></div><div style=""><font face="arial, helvetica, sans-serif">,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros predados</font></div>
<div style=""><font face="arial, helvetica, sans-serif">dens<-sort(rep(2^(2:6),10))#### Densidade de herbivoros oferecidos</font></div><div style=""><font face="arial, helvetica, sans-serif">##</font></div>
<div style=""><font face="arial, helvetica, sans-serif">##Regressão
pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ exp(P0+P1*N+P2*N^2+P3*N3)</font></div><div style=""><font face="arial, helvetica, sans-serif">p.model1<-glm(pred/par~dens+I(dens^2)+I(dens^3),family="binomial")</font></div>
<div style=""><font face="arial, helvetica, sans-serif">summary(p.model1)</font></div><div style=""><font face="arial, helvetica, sans-serif">#</font></div><div style="font-style: normal; font-size: 13px; background-color: transparent; font-family: arial, helvetica, sans-serif;">
<font face="arial, helvetica, sans-serif"><br></font></div><div style="font-style: normal; font-size: 13px; background-color: transparent; font-family: arial, helvetica, sans-serif;"><font face="arial, helvetica, sans-serif"><br>
</font></div><div style=""><font face="arial, helvetica, sans-serif">Minha dúvida é se o coeficiente linear que aparece no summary esta transformado em exp(x)/1+exp(x)</font></div><div style="">
<font face="arial, helvetica, sans-serif">ou trata-se do valor sem transformação?</font></div><div style=""><font face="arial, helvetica, sans-serif"><br></font></div><div style=""><font face="arial, helvetica, sans-serif">Obrigado,</font></div>
<span class="yiv1749391770HOEnZb"><font color="#888888"><div style="font-style: normal; font-size: 13px; background-color: transparent; font-family: arial, helvetica, sans-serif;"><font face="arial, helvetica, sans-serif"><br>
</font></div>
<div style="font-style: normal; font-size: 13px; background-color: transparent; font-family: arial, helvetica, sans-serif;"><font face="arial, helvetica, sans-serif"><pre>--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso - Campus Cáceres
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: <a rel="nofollow" href="">(+55) 65 8132-8112</a> (TIM) <a rel="nofollow" href="">(+55) 65 9686-6970</a> (VIVO)
<a rel="nofollow" ymailto="mailto:e-mails:alexandresantosbr@yahoo.com.br" target="_blank" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a rel="nofollow" ymailto="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
======================================================================</pre></font></div><div style="font-style: normal; font-size: 13px; background-color: transparent; font-family: arial, helvetica, sans-serif;">
<font face="arial, helvetica, sans-serif"><br></font></div></font></span></div></div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a rel="nofollow" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a rel="nofollow" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a rel="nofollow" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<i>Vinicius Brito Rocha.</i><br><i style="font-weight:bold;">Estatístico e Atuário <font size="1">(IM / UFRJ)</font></i><i style="font-weight:bold;"><br>Mestre em Pesquisa Operacional <font size="1">(COPPE / UFRJ)</font></i><br>
<br><a rel="nofollow" target="_blank" href="http://www.aplicademic.blogspot.com/">www.aplicademic.blogspot.com</a><br><a rel="nofollow" target="_blank" href="http://twitter.com/viniciusbritor">http://twitter.com/viniciusbritor</a><br><br>"Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.<br>
<br><br>
</div>
</div><br><br> </div> </div> </div></body></html>