<div dir="ltr">Use family = quasipoisson ou family =quasi(link = "log", variance = "mu")<br><br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 8 de agosto de 2014 03:51, Leonardo Ferreira Fontenelle <span dir="ltr"><<a href="mailto:leonardof@leonardof.med.br" target="_blank">leonardof@leonardof.med.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>




<div><div>Como se consegue no R uma regressão Poisson com variância robusta? Quasipoisson?<br></div>
<div> </div>
<div>Procurando a internet encontrei sugestões de como encontrar a variância robusta em si; seria necessário uma trabalheira para recalcular os intervalos de confiança e/ou valores p.<br></div>
<div> </div>
<div><div><a href="http://lattes.cnpq.br/9234772336296638" title="Currículo Lattes" target="_blank">Leonardo Ferreira Fontenelle</a><br></div>
</div>
<div> </div>
<div> </div>
<div>Em Qua 6 ago. 2014, às 17:03, Manoel Galdino escreveu:<br></div>
<blockquote type="cite"><div dir="ltr"><div>Mas o fato da regressão ser de poisson não tem relação com erro padrão ser robusto ou não. O que a regressão de poisson faz é mudar a forma funcional do modelo. Mas o erro padrão, se robusto ou não, não tem a ver com isso.<br>
</div>
<div>Basicamente, temos:<br></div>
<div> E[y] = f^-1(u), u = BX<br></div>
<div>na linear: f^-1(u) = u (identidade)<br></div>
<div>na logistica: f^-1(u) = invlogit(u)<br></div>
<div>na poisson: f^-1(u) =  exp(u)<br></div>
<div> </div>
<div>E assim por diante. <a href="http://en.wikipedia.org/wiki/Generalized_linear_model" target="_blank">A wikipedia explica isso em detalhes</a>.<br></div>
<div> </div>
<div>abçs<br></div>
<div>M<br></div>
</div>
<div><div> </div>
<div> </div>
<div><div>2014-08-06 16:21 GMT-03:00 geovane barbosa <span dir="ltr"><<a href="mailto:geovanecb@yahoo.com.br" target="_blank">geovanecb@yahoo.com.br</a>></span>:<br></div>
<blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:10pt">
<div><span>Foi uma sugestão do revisor para termos erros padrões robustos. Eu também fique com a mesma dúvida. Realizei a regressão logsitica e deu tudo certo. Só que na regressão de poisson esta dando esse erro. vou verificar .</span><br>
</div>
<div> </div>
<div> </div>
<div> </div>
<div> </div>
<div> </div>
<div style="display:block"><div style="font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif;font-size:10pt"><div style="font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif;font-size:12pt">
<div dir="ltr"><span style="font-family:Arial"> Em Quarta-feira, 6 de Agosto de 2014 16:18, Manoel Galdino <<a href="mailto:mcz.fea@gmail.com" target="_blank">mcz.fea@gmail.com</a>> escreveu:</span><br></div>
<div><div><div> </div>
<div> </div>
<div><div><div><div dir="ltr"><div>Mas não faz sentido aplicar regressão de poisson para uma vd binária, faz?<br></div>
<div><div> </div>
<div>De todo modo, pelo, erro, AF tem valor não permitido. Dê um summary em AF. <br></div>
<div>summary(dados$AF) ou unique(dados$AF) e você deve encontrar o erro.<br></div>
<div> </div>
<div>abç<br></div>
<div>M<br></div>
</div>
</div>
<div><div> </div>
<div> </div>
<div><div>2014-08-06 16:11 GMT-03:00 geovane barbosa <span dir="ltr"><<a shape="rect" href="mailto:geovanecb@yahoo.com.br" target="_blank">geovanecb@yahoo.com.br</a>></span>:<br></div>
<blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div style="color:rgb(0,0,0);font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif;font-size:10pt;background-color:rgb(255,255,255)">
<div> </div>
<div> </div>
<div> </div>
<div>Olá pessoal boa tarde. Estou rodando um modelo de regressão de poisson, porém esta dando o seguinte erro.<br></div>
<div> </div>
<div><div>model2<- glm(AF ~ Sexo + Idade + Cor + Esc.mae + Escolinha + Socio.eco + EP,  data=dados, family="poisson",na.action="na.omit")<br></div>
<div> </div>
</div>
<div> </div>
<div>AF é binária.<br></div>
<div> </div>
<div><div>Erro em if (any(y < 0)) stop("negative values not allowed
 for the 'Poisson' family") : <br></div>
<div>  valor ausente onde TRUE/FALSE necessário<br></div>
<div>Além disso: Mensagens de aviso perdidas:<br></div>
<div>In Ops.factor(y, 0) : < not meaningful for factors<br></div>
<div>> <br></div>
<div> </div>
<div>grato<br></div>
</div>
</div>
</div>
</div>
<div> </div>
<div>_______________________________________________<br></div>
<div>
R-br mailing list<br></div>
<div> <a shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br></div>
<div> <a shape="rect" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br></div>
<div>
Leia o guia de postagem (<a shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></div>
</blockquote></div>
<div> </div>
<div> </div>
<div> </div><span class="HOEnZb"><font color="#888888">
<div>-- <br></div>
<div>

Manoel Galdino<br></div>
<div><a shape="rect" href="https://sites.google.com/site/galdinomcz/" target="_blank">https://sites.google.com/site/galdinomcz/</a><br></div>
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
<div> </div>
<div> </div>
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="HOEnZb"><font color="#888888">
<div> </div>
<div> </div>
<div> </div>
<div>-- <br></div>
<div>Manoel Galdino<br></div>
<div><a href="https://sites.google.com/site/galdinomcz/" target="_blank">https://sites.google.com/site/galdinomcz/</a><br></div>
</font></span></div><span class="HOEnZb"><font color="#888888">
<div><u>_______________________________________________</u><br></div>
<div>R-br mailing list<br></div>
<div><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br></div>
<div><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br></div>
<div>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></div>
</font></span></blockquote><div> </div>
</div>

<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">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"><br>-- <br>Wagner Hugo Bonat<br>
LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná
</div>