3
$\begingroup$

I want to model the change in rate of an event in a sampled population. The populations are grouped according to some grouping variable. Each individual observation is of one organism, and the sample can either return as positive or negative (+/-). My resulting parameters are:

-number of positive samples (+)

-total number tested per grouping variable

Should this be modeled using a Poisson or binomial distribution? My sample parameters are meant to estimate the true rate value in an overall larger population.

$\endgroup$

2 Answers 2

5
$\begingroup$

As in this question, in typical count data situations, the rule of thumb is:

  • if the count has a ceiling/maximum value set by the experimental design, model the response as binomial (or some over/underdispersed variant: quasibinomial, beta-binomial, observation-level random effect ...)
  • if there is no well-defined limit (e.g. the number of trees in a 1-hectare plot can't be infinite, but we can't typically quantify the number of available "tree sites" that are available) then use a Poisson response (or some variant: quasi-Poisson, negative binomial, generalized Poisson, COM-Poisson ...)
  • if the count has a maximum value but the proportion of the maximum is always small (e.g. the number of cancer cases in a county), then a binomial and a Poisson with a log-offset term to scale the maximum value will give nearly identical results, and it's a matter of computational convenience.

If you want a literature reference I have one for you :-): Bolker, Benjamin M. Ecological Models and Data in R Princeton University Press, 2008:

p. 161:

You should only use the binomial in fitting data when there is an upper limit to the number of possible successes. When $N$ is large and $p$ is small, so that the probability of getting N successes is small, the binomial approaches the Poisson distribution, which is covered in the next section (Figure 4.17). Examples: number of surviving individuals/nests out of an initial sample; number of infested/infected animals, fruits, etc. in a sample; number of a particular class (haplotype, subspecies, etc.) in a larger population.


p. 163:

The Poisson distribution has no upper limit, although values much larger than the mean value are highly improbable. This characteristic provides a rule for choosing between the binomial and Poisson. If you expect to observe a “ceiling” on the number of counts, you should use the binomial; if you expect the number of counts to be effectively unlimited, even if it is theoretically bounded (e.g. there can’t really be an infinite number of plants in your sampling quadrat), use the Poisson.

Examples: number of seeds/seedlings falling in a gap; number of offspring produced in a season (although this might be better fit by a binomial if the number of breeding attempts is fixed); number of prey caught per unit time.


PS it turns out I am quite consistent: I wrote this answer without looking at the examples I used in 2008, and came up with almost the same wording and examples ...

$\endgroup$
2
  • $\begingroup$ This makes sense to me. Thank you. Do you know of any literature that defines these rules of thumb? $\endgroup$ Commented Jan 20, 2022 at 23:00
  • $\begingroup$ @geoscience123 If the answer makes sense to you, you might consider accepting it. If none of the answers are satisfactory to you, so that you don't want to accept any, then can you explain what aspect of your question has not yet been addressed? Looking at your profile, it appears that you hardly ever accept answers that you receive. $\endgroup$ Commented Jul 18, 2022 at 9:02
1
$\begingroup$

Poisson and Binomial distributions both model the counts of events.

The Binomial distribution is about the number of counts of events in a sequence of experiments that measure a boolean valued outcome. The events can only occur on a fixed number of cases/experiments.

The Poisson distribution is about the number of counts of events in an interval of time or space. The events can occur continuously.


See also on Wikipedia:

https://en.wikipedia.org/wiki/Poisson_distribution

In probability theory and statistics, the Poisson distribution (/ˈpwɑːsɒn/; French pronunciation: ​[pwasɔ̃]), named after French mathematician Siméon Denis Poisson, is a discrete probability distribution that expresses the probability of a given number of events occurring in a fixed interval of time or space if these events occur with a known constant mean rate and independently of the time since the last event.

https://en.wikipedia.org/wiki/Binomial_distribution

In probability theory and statistics, the binomial distribution with parameters n and p is the discrete probability distribution of the number of successes in a sequence of n independent experiments, each asking a yes–no question, and each with its own Boolean-valued outcome: success (with probability p) or failure (with probability q = 1 − p).

$\endgroup$
4
  • 1
    $\begingroup$ This doesn't really answer the question; from this reading, both are plausible, which is the reason OP asked the question in the first place. $\endgroup$
    – Noah
    Commented Jan 20, 2022 at 21:44
  • $\begingroup$ @Noah I give hints what Poisson and Binomial distributions are modelling. They both model counts. The Binomial distribution is about the number of counts in a sequence of experiments (measurements of a boolean valued outcomes). The Poisson distribution is about the number of counts in an interval of time or space. $\endgroup$ Commented Jan 20, 2022 at 23:09
  • $\begingroup$ @Noah "Each individual observation is of one organism, and the sample can either return as positive or negative (+/-). " These observations can be considered as Bernoulli trials. If the sum of the observations is given then the distribution is Binomial distributed. (Why do you believe that a Poisson distribution is plausible as well?) $\endgroup$ Commented Jan 20, 2022 at 23:15
  • $\begingroup$ Your comment completes your answer. OP could have looked up the methods on Wikipedia; it's your statistical expertise (i.e., interpreting which model is best given the study design) that makes the answer useful. $\endgroup$
    – Noah
    Commented Jan 21, 2022 at 5:34

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.