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@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.
______________________________________________
R-help@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.
