Skip to content
Prev 10414 / 20628 Next

zero Inflation in glmmadmb

I'm going to take the liberty of replying by way of r-sig-mixed-models,
as this is a generally useful question ...
On 13-07-29 04:08 AM, Bartzke, Gundula wrote:
At present glmmadmb fits only a constant-zero-inflation model,
corresponding
to a zeroinfl model of y~...|1 (zeroinfl's default is to make the
conditional mean and zero-inflation models the same). This is the
same as VGAM's default (i.e., with the 'zero' parameter excluded:
'zero=NULL' fits a model where conditional mean, dispersion, _and_
zero-inflation all use the same covariates).

  The results (see below) seem to be perfectly consistent.
The terminology is highly variable.  In this case alpha is indeed
the same as theta/k/size ...

library(pscl)
library(glmmADMB)
library(VGAM)

Owls <- transform(Owls,
                  Nest=reorder(Nest,NegPerChick),
                  logBroodSize=log(BroodSize),
                  NCalls=SiblingNegotiation)

## compare all three models, with constant zero-inflation probability
##  (default for VGAM and glmmADMB)
fit_zinb <- zeroinfl(NCalls~(FoodTreatment+ArrivalTime)*SexParent|1,
                     data=Owls, dist="negbin")
fit_zinbadmb <- glmmadmb(NCalls~(FoodTreatment+ArrivalTime)*SexParent,
                         data=Owls,
                         family="nbinom", zeroInfl=TRUE)
fit_zinegbinomial <- vglm(NCalls~(FoodTreatment+ArrivalTime)*SexParent,
                          data=Owls, zinegbinomial())

summary(fit_zinb)
summary(fit_zinbadmb)
summary(fit_zinegbinomial)

## glmmADMB fits zero-inflation and alpha on non-transformed scales
qlogis(fit_zinbadmb$pz)
log(fit_zinbadmb$alpha)

## 'alpha' (glmmADMB) = 'theta' (zeroinfl) = 'size' (VGAM)
## 'zero-inflation' = 'zero-inflation' = 'pstr0' (VGAM)
Message-ID: <51F665BF.9060905@mcmaster.ca>
In-Reply-To: <7015_1375085417_r6T8AF4m004479_34E792A3F9A41047A185428BA88DCA5497540796@NINSRV05.nina.no>