<div dir="ltr"><div>De forma geral os erros padrões obtidas por quasi verossimilhança são tidos como robustos no sentido de suportarem má especificação do modelo, tanto na estrutura de erros como propriamente na distribuição da variável resposta. Os modelos quasi não assumem uma distribuição explicitamente para os dados apenas um relacionamento entre média e variancia. Além disso, é possível mostrar que uma grande quantidade de distribuição por exemplo, a familia exponential é descrita apenas pelos primeiros dois momentos e que a função de variancia temp um comportamento do tipo potencia (power variance function), por exemplo, <br>
<br></div>V(mu) = phi*mu^power <br><div><br></div><div>Se phi e power = 1 vc tem o poisson se p = 0 vc tem a distribuição Normal power = 2 gamma e power = 3 log-normal. Todas são casos particulares da distribuição Tweedie que é descrita pelos 2 primeiros momentos. Em suma, vc precisa apenas dos dois primeiros momentos pra estimar, neste sentido a abordagem é robusta a má especificação, por consequencia os erros são tbm robustos e em geral conservadores.<br>
</div><div><br> O modelo quasi poisson suporta subdispersion (phi < 1) se vc quer saber mais sobre subdipersion veja <a href="http://www.leg.ufpr.br/doku.php/publications:papercompanions:zeviani-jas2014">http://www.leg.ufpr.br/doku.php/publications:papercompanions:zeviani-jas2014</a>.<br>
</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 8 de agosto de 2014 21:55, Manoel Galdino <span dir="ltr"><<a href="mailto:mcz.fea@gmail.com" target="_blank">mcz.fea@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Mas eu acho que o quasipoisson vai permitir é overdisperssion. No modelo de poisson, a variância é igual à média. Quasipoisson permite que a variância não seja igual à média. Mas não creio que isso signifique variância (erro padrão?) robusta (o).<div>
<br></div><div>abçs</div><div>M</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014-08-08 16:40 GMT-03:00 Leonardo Ferreira Fontenelle <span dir="ltr"><<a href="mailto:leonardof@leonardof.med.br" target="_blank">leonardof@leonardof.med.br</a>></span>:<div>
<div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>
<div><div>Obrigado, Wagner!<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 Sex 8 ago. 2014, às 04:06, Wagner Bonat escreveu:<br></div><div><div>
<blockquote type="cite"><div dir="ltr">Use family = quasipoisson ou family =quasi(link = "log", variance = "mu")<br></div>
<div dir="ltr"> </div>
<div dir="ltr"> </div>
<div><div> </div>
<div> </div>
<div><div>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></div>
<blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><u></u><br></div>
<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>
<div><span><span style="color:rgb(136,136,136)"><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>
</span></span><br></div>
</div>
<div> </div>
</div>
<div> </div>
</div>
<div><span><span style="color:rgb(136,136,136)"><div> </div>
<div> </div>
</span></span><br></div>
</div>
<div> </div>
</div>
<div> </div>
</div>
<div> </div>
</div>
<div> </div>
</div>
<div> </div>
</div>
<div> </div>
</div>
<div> </div>
</div>
<div> </div>
</blockquote></div>
<div><span><span style="color:rgb(136,136,136)"><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>
</span></span><br></div>
</div>
<div><span><span style="color:rgb(136,136,136)"><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>
</span></span><br></div>
</blockquote><div> </div>
</div>
<div> </div>
<div>_______________________________________________<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>
</blockquote></div>
<div> </div>
<div> </div>
<div> </div>
<div>-- <br></div>
<div>Wagner Hugo Bonat<br></div>
<div>
LEG - Laboratório de Estatística e Geoinformação<br></div>
<div>UFPR - Universidade Federal do Paraná<br></div>
</div>
<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>
</blockquote><div> </div>
</div></div></div>
<br>_______________________________________________<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" 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></div></div><div><div class="h5"><br><br clear="all">
<div><br></div>-- <br>
Manoel Galdino<br><a href="https://sites.google.com/site/galdinomcz/" target="_blank">https://sites.google.com/site/galdinomcz/</a><br>
</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>