[R] Zinb for Non-interger data
Alain Zuur
highstat at highstat.com
Mon Jul 20 13:47:42 CEST 2009
JPS2009 wrote:
>
> Sorry bit of a Newbie question, and I promise I have searched the forum
> already, but I'm getting a bit desperate!
>
> I have over-dispersed, zero inflated data, with variance greater than the
> mean, suggesting Zero-Inflated Negative Binomial - which I attempted in R
> with the pscl package suggested on
> http://www.ats.ucla.edu/stat/R/dae/zinbreg.htm
>
> However my data is non-integer with some pesky decimals (i.e. 33.12) and
> zinb / pscl doesn't like that - not surprising as zinb is for count data,
> normally whole integers etc.
>
> Does anyone know of a different zinb package that will allow non-integers
> or and equivalent test/ model to zinb for non-integer data? Or should I
> try something else like a quasi-Poisson GLM?
>
>
> Apologies for the Newbie question! Any help much appreciated!
> Thanks!
>
Is it really non-integer...or is it a density (in which case you could use
NB + offset)?
The quasi-Poisson will not help you with the zero inflation.
I'm afraid you will have to do some hard programming by combining the 0-1
binomial part with a continuous distribution on the second part of the
data......and I guess the easiest is to do this in MCMC. Perhaps the Gamma
distribution can be used? You would have to adjust all likelihood equations
as Gamma doesn't allow for zeros. But perhaps another continuous
distribution is more appropriate...depends on your data.
Alain Zuur
--
View this message in context: http://www.nabble.com/Zinb-for-Non-interger-data-tp24550044p24568326.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list