How can I calculate the deviance for the factors and the null and saturated models with Hurdle models?
I used the function
hurdle() from the package
pscl. The functions
deviance() are not applicable.
Cross Validated is a question and answer site for people interested in statistics, machine learning, data analysis, data mining, and data visualization. It only takes a minute to sign up.Sign up to join this community
As the hurdle and zero-inflation models are both composed of two submodels they do not have all the properties of classical GLMs (generalized linear models) like the Poisson model. That's why we decided not to provide a
anova() method for
"zeroinfl" objects. Instead there is a
logLik() method for extracting the log-likelihood (which is closely connected to the deviance in GLMs). For inference, functions like
lmtest for the likelihood ratio test can be used. Further options are discussed in
vignette("countreg", package = "pscl") (see e.g. the end of Section 2.2).