[R-br] Regresssão poisson
Manoel Galdino
mcz.fea em gmail.com
Sexta Agosto 8 16:55:30 BRT 2014
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 em 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 em 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 em 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 em 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 em 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 em 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 em 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 em 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 em 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 em 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/
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140808/b2a6d2f8/attachment.html>
Mais detalhes sobre a lista de discussão R-br