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çsM2014-08-08 16:40 GMT-03:00 Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br>:
Obrigado, Wagner!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.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 = BXna 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.abçsM2014-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çM2014-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árioAlém disso: Mensagens de aviso perdidas:In Ops.factor(y, 0) : < not meaningful for factors>grato_______________________________________________R-br mailing listLeia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.--Manoel Galdino--Manoel Galdino_______________________________________________R-br mailing listLeia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível._______________________________________________R-br mailing listLeia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.--Wagner Hugo BonatLEG - Laboratório de Estatística e GeoinformaçãoUFPR - Universidade Federal do Paraná_______________________________________________R-br mailing listLeia 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.
_______________________________________________
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.