[R-br] glm.fit: fitted probabilities numerically 0 or 1 occurred
Andre Oliveira
andreolsouza em yahoo.com.br
Quinta Outubro 13 05:50:20 BRT 2011
Uma possível solução para dados deste tipo.
require(brglm)
Carregando pacotes exigidos: brglm
Carregando pacotes exigidos: profileModel
m1<-brglm(y~ x1+x2, family=binomial)
summary(m1)
André Oliveira Souza
________________________________
De: Benilton Carvalho <beniltoncarvalho em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; Marcelo Luiz de Laia <marcelolaia em gmail.com>
Enviadas: Quarta-feira, 12 de Outubro de 2011 20:10
Assunto: Re: [R-br] glm.fit: fitted probabilities numerically 0 or 1 occurred
Imagino que faltou clareza na minha explicacao, minhas desculpas. Com
regressao logistica, vc esta' modelando
P(Y = 1 | X) e, consequentemente, P(Y = 0 | X),
na qual X e' o seu conjunto de covariaveis. Assim, o que eu tive a
intencao de dizer e' que existe uma combinacao de X para a qual vc
consegue prever perfeitamente a resposta (coisa como nenhuma semente
germinou no dia 2 e 3).
O que fazer? Encontrar exatamente qual eh esse subespaco para o qual a
predicao perfeita acontece e perguntar-se se isso faz mesmo sentido.
Na pratica, eu acho pouco provavel que um experimento (bem) delineado
gere resultados assim (note que eventos de probabilidade muito baixa
podem acontecer, caso contrario ninguem ganharia na mega-sena)... e,
dito isso, "minha solucao" seria gerar mais dados (estou certo de que
outros tem outras sugestoes, como considerar modelos alternativos, que
podem ser de fato a melhor solucao, dependendo do seu conhecimento
previo a respeito dos dados).
Tudo isso soa como algo para o qual Walmes possua solucoes mais elegantes.
b
_______________________________________________
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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111013/d18a2169/attachment.html>
Mais detalhes sobre a lista de discussão R-br