Maybe the dumbiest question ever posted on CV:
I want to analyze the relationship between proportion data and some covariates in a generalized linear (mixed, but I don't think it matters) model.
The natural distribution to use would be binomial but most of the proportions are close to 0 and their distribution looks like a (slightly overdispersed) Poisson.
My proportions are by definition not integer values, and so cannot be directly modeled using some implementations of gl(m)m with Poisson errors. For instance I'm using MCMCglmm in R which (rightfully) only allows integer values (and Beta-binomial is not implemented).
Could I simply multiply my proportion by let's say 1000 (to keep some accuracy), round them (to make them integer) and do the glmm with Poisson errors on them. Then I could rescale the parameter estimates (and in glmm the relative importance of variance components will stay the same)?
I assume (maybe wrongly) that there would be not problem with this, but as I've been unable to find it discussed somewhere, I'm starting to doubt, and would be pleased if others thinking about this later can find an answer (whatever it is) here.
Many thanks for the support.