[R-sig-ME] Is glmmADMB the only package in R which can handle, a zero-inflated Poisson model for repeated measures? (Highland Statistics Ltd)
Mueller,John Martin
john.mueller at louisville.edu
Wed May 30 17:03:51 CEST 2012
Alain, thanks for the reply and links to your book. I am going to try the glmmADMB option first. Then if that doesn't work, I will look at the Markov Chain Monte Carlo / BUGS option. - JM
----------------------------------------------------------------------
Message: 1
Date: Sun, 27 May 2012 11:36:54 +0100
From: Highland Statistics Ltd <highstat at highstat.com>
To: "r-sig-mixed-models at r-project.org"
<r-sig-mixed-models at r-project.org>
Subject: Re: [R-sig-ME] Is glmmADMB the only package in R which can
handle, a zero-inflated Poisson model for repeated measures?
Message-ID: <4FC203C6.6030303 at highstat.com>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
------------------------------
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
More information about the R-sig-mixed-models
mailing list