[R-sig-ME] Mixed model with zero truncated Poisson distributed data
Ben Bolker
bbolker at gmail.com
Wed Jan 11 19:46:46 CET 2012
Helen Ward <h.l.ward at ...> writes:
> I would like to describe the relationship between age and male
> reproductive success in a population of greater horseshoe bats.
>
> My data consists of three columns: MaleID, Age, NumberofPups (at that
> age). Many of the males appear multiple times in the data set, so I
> believe I need to derive a mixed model with MaleID as a random variable.
>
> The data is Poisson distributed, but zero-truncated. So far I have only
> succeeded in making a mixed model with a poisson distribution (using
> glmmPQL in the MASS package), and a zero truncated poisson model (using
> vglm in the VGAM package), but not a mixed model capable of handling
> zero truncated Poisson data.
>
> It has been suggested that I could just minus 1 from each value in the
> NumberofPups column to make a more usual Poisson distribution, so I can
> ignore the zero truncated bit. I have tried this and it changes the
> results of the model, but is this an acceptable transformation?
>
> If not, can anyone advise me on a mixed model that can handle zero
> truncated Poisson data please?
>
Thanks for letting us know about cross-posting.
You should be able to do this in either the MCMCglmm package or
(recent versions of) the glmmADMB package. In MCMCglmm, use
family="ztpoisson"; in glmmADMB, use family="truncpoiss" ...
Ben Bolker
More information about the R-sig-mixed-models
mailing list