[R] R wont accept my zero count values in the GLM with quasi_poisson dsitribution

Göran Broström goran.brostrom at umu.se
Tue Jul 28 09:05:38 CEST 2015



On 28/07/15 08:33, Charlotte wrote:
> Hello
>
> I have count values for abundance which follow a pattern of
> over-dispersal with many zero values.  I have read a number of
> documents which suggest that I don't use data transforming methods
> but rather than I run the GLM with the quasi poisson distribution.
> So I have written my script and R is telling me that Y should be more
> than 0.

No,  R  is telling you that you must have 0 <= y <= 1 (see below).
For count data you should not use the binomial family, but rather 
'poisson', or 'quasipoisson'.

?family

Göran

> Everything I read tells me to do it this way but I can't get R to
> agree. Did I need to add something else to my script to get it to
> work and keep my data untransformed? The script I wrote is as
> follows:
>
>> fit <- glm(abundance~Gender,data=teminfest,family=binomial())
>
> then I get this error Error in eval(expr, envir, enclos) : y values
> must be 0 <= y <= 1
>
> I don't use R a lot so I am having trouble figuring out what to do
> next.
>
> I would appreciate some help
>
> Many Thanks Charlotte
>
>
>
>
>
>
> -- View this message in context:
> http://r.789695.n4.nabble.com/R-wont-accept-my-zero-count-values-in-the-GLM-with-quasi-poisson-dsitribution-tp4710462.html
>
>
Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________ R-help at r-project.org
> mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the
> posting guide http://www.R-project.org/posting-guide.html and provide
> commented, minimal, self-contained, reproducible code.
>



More information about the R-help mailing list