[R-sig-ME] Is glmmADMB the only package in R which can handle, a zero-inflated Poisson model for repeated measures?

[Whit Armstrong]

If you care to upload an example BUGS model, I'll re-write the example
for rcppbugs and post it back.

Here is one example using the cbpp dataset:
https://github.com/armstrtw/rcppbugs.examples/blob/master/herd.r

-Whit


On Sun, May 27, 2012 at 6:36 AM, Highland Statistics Ltd
<highstat at highstat.com> wrote:
>
>
>
>
> ------------------------------
>
> Message: 2
> Date: Sat, 26 May 2012 17:59:28 +0000
> From: "Mueller,John Martin"<john.mueller at louisville.edu>
> To: "r-sig-mixed-models at r-project.org"
>        <r-sig-mixed-models at r-project.org>
> Subject: [R-sig-ME] Is glmmADMB the only package in R which can handle
>        a zero-inflated Poisson model for repeated measures?
> Message-ID:
>        <109E4F6374337D4C9153497DE5F3EEC51CA2A3E7 at EXMBX07.ad.louisville.edu>
> Content-Type: text/plain
>
> I am looking into which package to use in R to handle a zero-inflated
> Poisson model for repeated measures (growth curve analysis through multiple
> level modeling).
>
>> From what I can tell, the lme4 package can handle the Poisson model and
>> repeated measures but can't handle the zero-inflated Poisson model.  The
>> pscl package can handle the zero-inflate Poisson model but can't handle
>> repeated measures.  And that leaves the glmmADMB package as being the only
>> package in R that can handle a zero-inflated Poisson model for repeated
>> measure.
>
>
> Is this correct?  Or is there a package in R other than the glmmADMB package
> that I should be considering?
>
> And if the glmmADMB package is the only package for this purpose, can the
> glmmADMB handle a 3-level multi-level model (hierarchical linear model)?
>
> The model that I am working on has TIME at level 1, FIRMS at level 2, and
> LOCATIONS (metropolitan statistical areas) at level 3.  The dependent
> variable at level 1 is a count variable (number of credit sources used by a
> firm) which has an over dispersion of zeros at various points in time.
>  Thus, the reason why I am having to use a zero-inflated Poisson model (or I
> could use a hurdle model).
>
> Thanks for your help.
>
> - JM
>
>        [[alternative HTML version deleted]]
>
>
>
> ------------------------------
>
>
>
> John,
> You can easily do this with R2WinBUGS and then in WinBUGS. See Chapter 4 of
> our 2012 book
> in which we analyse 2-way nested data; R code is include. Extension to your
> 3-way nested data would be
> an extra loop, plus an extra nesting level of the RE. It is also explained
> how you can calculate an
> intraclass correlation for each level.
>
> The choice between ZIP and hurdle depends on your underlying question.
>
> Alain
>
>
> --
>
> Dr. Alain F. Zuur
> First author of:
>
> 1. Analysing Ecological Data (2007).
> Zuur, AF, Ieno, EN and Smith, GM. Springer. 680 p.
> URL: www.springer.com/0-387-45967-7
>
>
> 2. Mixed effects models and extensions in ecology with R. (2009).
> Zuur, AF, Ieno, EN, Walker, N, Saveliev, AA, and Smith, GM. Springer.
> http://www.springer.com/life+sci/ecology/book/978-0-387-87457-9
>
>
> 3. A Beginner's Guide to R (2009).
> Zuur, AF, Ieno, EN, Meesters, EHWG. Springer
> http://www.springer.com/statistics/computational/book/978-0-387-93836-3
>
>
> 4. Zero Inflated Models and Generalized Linear Mixed Models with R. (2012)
> Zuur, Saveliev, Ieno.
> http://www.highstat.com/book4.htm
>
> Other books: http://www.highstat.com/books.htm
>
>
> Statistical consultancy, courses, data analysis and software
> Highland Statistics Ltd.
> 6 Laverock road
> UK - AB41 6FN Newburgh
> Tel: 0044 1358 788177
> Email: highstat at highstat.com
> URL: www.highstat.com
> URL: www.brodgar.com
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models