[R-es] Warning en modelo ZINB

[Carlos J. Gil Bellosta]

¿Podría ser que para algún nivel de alguna variable independiente
categórica solo hubiese ceros? En ese caso, casi seguro, aparecería ese
tipo de warning.

El lun., 9 abr. 2018 a las 19:00, <miriam.alzate en unavarra.es> escribió:

> Muchas gracias por la respuesta. He mirado y los coeficientes no son altos
> pero sí tengo una gran cantidad de ceros en la variable dependiente (más
> del 90%). Sin embargo, al incluir otro tipo de variables independientes no
> me da ese aviso, dejando la misma variable dependiente.
>
> ¿Cómo podría utilizar stan/rstan de forma sencilla para diagnosticar el
> modelo?
>
> Muchas gracias
>
> El Lun, 9 de Abril de 2018, 18:48, Carlos J. Gil Bellosta escribió:
> > Hola, ¿qué tal?
> >
> > El "warning" que comentas aparece en glm.fit precisamente cuando un
> > coeficiente diverge.
> >
> > El aviso puede ser malo o irrelevante, depende. Puede que haya sido
> > emitido
> > en algún paso intermedio del ajuste (por lo que no habría mayor
> problema).
> > O que afecte al ajuste entero. Uno de los síntomas de que el ajuste es
> > malo
> > es que algún coeficiente de tu modelo es excesivamente grande. O que el
> > peso correspondiente a tus ceros esté muy próximo a 0 o 1 (que es lo
> mismo
> > en el contexto de un modelo inflado). Mira a ver qué pinta tienen los
> > coeficientes y si tienen cierto sentido a la vista de tus datos. Vigila
> > los
> > extremadamente altos.
> >
> > Si no estás segura de tus coeficientes, puedes implementar el modelo
> > generativamente vía stan/rstan. Es sencillo y tendrías buenos
> diagnósticos
> > de lo que esta ocurriendo.
> >
> > Un saludo,
> >
> > Carlos J. Gil Bellosta
> > http://www.datanalytics.com
> >
> >
> >
> >
> > El lun., 9 abr. 2018 a las 18:34, <miriam.alzate en unavarra.es> escribió:
> >
> >> Buenas tardes,
> >>
> >> Estoy estimando un modelo binomial negativo de ceros inflados (ZINB)
> >> utilizando el comando zeroinfl() del paquete pscl. Al ejecutarlo me da
> >> el
> >> siguiente aviso:
> >> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
> >>
> >> ¿Sabéis que significa y si puedo usar el modelo aún con ese aviso? ¿Los
> >> coeficientes son fiables?
> >>
> >> Muchas gracias,
> >>
> >> Miriam
> >>
> >> _______________________________________________
> >> R-help-es mailing list
> >> R-help-es en r-project.org
> >> https://stat.ethz.ch/mailman/listinfo/r-help-es
> >>
> >
>
>
>

	[[alternative HTML version deleted]]