My dataset (with the dependent variable being a count variable) has two features that are different from a standard Poisson regression. I haven't been able to find a more advanced model that accounts for both irregularities, but have found alternatives that can deal with one or the other irregularity. I'll outline the situation below, but my questions are:
Given the features below, is there some hybrid model I could use that simultaneously accounts for both issues?
If not, is my strategy below reasonable in the circumstances?
I have an experiment where people are given a cue word and have to respond with another word that they think other people will pick too. I.e. they have to try coordinate. The dependent count variable is simply the number of people picking each word. The minimum value of this variable, then, is 1. However, while some people were able to coordinate well, many more were not. This means there are a lot of data points at y=1.
A zero-truncated poisson (e.g. in R,
vglm(y~x,family=pospoisson)) explicitly models the fact that the minimum count is 1 rather than 0 but doesn't address over-inflation of this minimum count. On the other hand, a hurdle model (e.g. in R
hurdle(y~x,dist=poisson)) explicitly models inflated minimum counts, but it expects a minimum count of 0 rather than 1, so to use it for my data, I'd have to use
hurdle(y-1~x), so it doesn't address the fact that I have a positive poisson distribution. This doesn't seem conceptually bothersome to me, since a count of 1 (only one person picked the word) means there has been a failure to coordinate.
I don't know of a model that can address both issues simultaneously. I have asked another question here about whether shifting y-1 so that I can use the hurdle model is in principle problematic, or similar enough for this to be informative.
My strategy so far has been to do all three (regular poisson, zero-truncated poisson, hurdle model with shifted values), get boot strapped confidence intervals for the model parameters to show that they're not significantly different, and conclude that my predictor consistently has a significant effect that is similar across these models, and that over-inflation and having a positive poisson distribution thus do not undermine the results of the basic model. Is that a reasonable strategy, or can you recommend something better?
The dependent count variable ranges from 1 to 18 (20 people were given each cue word, but there were no cases where all 20 gave the same response), with most of the weight at 1. There is one predictor, which ranges continuously from 0 to 1. For extremely low values of that predictor, I expect no coordination (y=1 meaning only one person generated the response). For low values of the predictor, I expect small amounts of coordination, and more people should coordinate over a response as the value of the predictor increases.