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). abçs M 2014-08-08 16:40 GMT-03:00 Leonardo Ferreira Fontenelle < leonardof@leonardof.med.br>:
Obrigado, Wagner!
Leonardo Ferreira Fontenelle <http://lattes.cnpq.br/9234772336296638>
Em Sex 8 ago. 2014, às 04:06, Wagner Bonat escreveu:
Use family = quasipoisson ou family =quasi(link = "log", variance = "mu")
Em 8 de agosto de 2014 03:51, Leonardo Ferreira Fontenelle < leonardof@leonardof.med.br> escreveu:
Como se consegue no R uma regressão Poisson com variância robusta? Quasipoisson?
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.
Leonardo Ferreira Fontenelle <http://lattes.cnpq.br/9234772336296638>
Em Qua 6 ago. 2014, às 17:03, Manoel Galdino escreveu:
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. Basicamente, temos: E[y] = f^-1(u), u = BX na linear: f^-1(u) = u (identidade) na logistica: f^-1(u) = invlogit(u) na poisson: f^-1(u) = exp(u)
E assim por diante. A wikipedia explica isso em detalhes <http://en.wikipedia.org/wiki/Generalized_linear_model>.
abçs M
2014-08-06 16:21 GMT-03:00 geovane barbosa <geovanecb@yahoo.com.br>:
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 .
Em Quarta-feira, 6 de Agosto de 2014 16:18, Manoel Galdino < mcz.fea@gmail.com> escreveu:
Mas não faz sentido aplicar regressão de poisson para uma vd binária, faz?
De todo modo, pelo, erro, AF tem valor não permitido. Dê um summary em AF. summary(dados$AF) ou unique(dados$AF) e você deve encontrar o erro.
abç M
2014-08-06 16:11 GMT-03:00 geovane barbosa <geovanecb@yahoo.com.br>:
Olá pessoal boa tarde. Estou rodando um modelo de regressão de poisson, porém esta dando o seguinte erro.
model2<- glm(AF ~ Sexo + Idade + Cor + Esc.mae + Escolinha + Socio.eco + EP, data=dados, family="poisson",na.action="na.omit")
AF é binária.
Erro em if (any(y < 0)) stop("negative values not allowed for the 'Poisson' family") : valor ausente onde TRUE/FALSE necessário Além disso: Mensagens de aviso perdidas: In Ops.factor(y, 0) : < not meaningful for factors
grato
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Manoel Galdino https://sites.google.com/site/galdinomcz/
-- Manoel Galdino https://sites.google.com/site/galdinomcz/
*_______________________________________________* R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Wagner Hugo Bonat LEG - Laboratório de Estatística e Geoinformação UFPR - Universidade Federal do Paraná *_______________________________________________* R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Manoel Galdino https://sites.google.com/site/galdinomcz/