everyone. I am fitting a glmm model using the R library glmmTMB
for predicting a count response variable with excess-zeros and overdispersion (nbinom2
> Poisson
).
Additionally, I am insterested in exploring the zero-inflation component of my sample data, therefore a zero-inflation
vs hurdle
question arises. Even though both approaches are two-component models with finite mixture, two base assumptions between them arises
hurdle
models are those in which the zero vs non-zero component of the count are mechanistically splitted (i.e. you assume mechanism producing the zero part are different from the non-zero part), which mean assuming all your sample zeros are true or structural zeros. Then, a logit model for the zero vs no-zero plus a count model for counts > 0 are presented together.zero-inflated
models, on the other hand, while also splitting the zero vs non-zero components, allow the count model to produce zero observations (i.e., true zeros), while the zero-inflation component becomes a point-mass extra-zero generator (i.e., false or excess zeros). Those zero counts would represent false abscenses e.g. a sampling error or a mistake in the measurement (maybe you were not paying enough attention...)
When implementing those in glmmTMB
, the argument ziformula =
enables the zero-inflation model, while family = truncated_nbinom2()
defines the hurdle model
However, I feel puzzled since attempting to fit a hurdle truncated_nbinom2
model yield the following error message:
'y' contains zeros (or values close to zero). Zeros are compatible with a truncated distribution only when zero-inflation is added
Then my model can be defined as follows:
m0 <- glmmTMB(formula + (1 | id),
offset = log(offset),
data = df,
family = truncated_nbinom2(),
ziformula= ~ .)
which is producing a nice model that meets my requirements. In fact, it is the only model that i have been able to affectively fit.
I have seen these kind of models in other examples, and it is obvious glmmTMB
has specific implementations of zero-inflation with truncated families, but I have never read in the literature about hurdle + zi models as they are always presented as mutually exclussive, focusing on the criteria for selecting and comparing between them
What am I missing here??
- Is this approach statistically valid?
- How is the result interpreted? I guess the zero-inflation component is the same as any other zi model. But my conditional model, as a truncated family, should no be able to account for true zero counts and therefore the count model might be biased
I'd love if someone could shed some light into this issue
P.S: i also found some questions regarding prediction with these kind of models, truncated with zero-inflation here