20
$\begingroup$

I'm trying to get a handle on the concept of overdispersion in logistic regression. I've read that overdispersion is when observed variance of a response variable is greater than would be expected from the binomial distribution.

But if a binomial variable can only have two values (1/0), how can it have a mean and variance?

I'm fine with calculating the mean and variance of successes from x number of Bernoulli trials. But I cannot wrap my head around the concept of a mean and variance of a variable that can only have two values.

Can anyone provide an intuitive overview of:

  1. The concept of a mean and variance in a variable that can only have two values
  2. The concept of overdispersion in a variable that can only have two values
$\endgroup$
5
  • 1
    $\begingroup$ Add 20 values of $y$, where 10 are $0$ and 10 are $1$. Can you divide this by 20? Can you compute the s.d. $y$? $\endgroup$
    – Sycorax
    Commented Mar 27, 2014 at 19:53
  • $\begingroup$ Nicely put so I believe that's mean = 0.5, standard deviation = 0.11. $\endgroup$
    – luciano
    Commented Mar 27, 2014 at 20:23
  • $\begingroup$ Say my response variable had 100 successes and 5 fails. Is this likely to be overdispersed? $\endgroup$
    – luciano
    Commented Mar 27, 2014 at 20:25
  • $\begingroup$ luciano, you need more than one realization of the experiment to determine if it is overdispersed. $\endgroup$
    – Underminer
    Commented Mar 27, 2014 at 20:32
  • $\begingroup$ I recommend having a look at the comments in cran.r-project.org/web/packages/DHARMa/vignettes/… $\endgroup$ Commented Jan 27, 2021 at 15:12

2 Answers 2

15
$\begingroup$

A binomial random variable with $N$ trials and probability of success $p$ can take more than two values. The binomial random variable represents the number of successes in those $N$ trials, and can in fact take $N+1$ different values ($0,1,2,3,...,N$). So if the variance of that distribution is greater than too be expected under the binomial assumptions (perhaps there are excess zeros for instance), that is a case of overdispersion.

Overdispersion does not make sense for a Bernoulli random variable ($N = 1$)

In the context of a logistic regression curve, you can consider a "small slice", or grouping, through a narrow range of predictor value to be a realization of a binomial experiment (maybe we have 10 points in the slice with a certain number of successes and failures). Even though we do not truly have multiple trials at each predictor value and we are looking at proportions instead of raw counts, we would still expect the proportion of each of these "slices" to be close to the curve. If these "slices" have a tendency to be far away from the curve, there is too much variability in the distribution. So by grouping the observations, you create realizations of binomial random variables rather than looking at the 0/1 data individually.

The Example below is from another question on this site. Lets say the blue lines represents the expected proportion over the range of predictor variables. The blue cells indicate observed instances (in this case schools). This provides a graphical representation of how overdispersion may look. Note that there are flaws with interpreting the cells of the graph below, but it provides an idea of how overdispersion can manifest itself.

Over Dispersion Example

$\endgroup$
7
  • 2
    $\begingroup$ But I'm interested in overdispersion in the context of logistic regression. For each value of a predictor variable in logistic regression, there isn't n trials, there is only one trial. And the result of that one trial is either success or fail $\endgroup$
    – luciano
    Commented Mar 27, 2014 at 20:39
  • $\begingroup$ I just added a paragraph to address the intuition behind overdispersion in the context of linear regression. $\endgroup$
    – Underminer
    Commented Mar 27, 2014 at 21:25
  • 1
    $\begingroup$ Underminer, I'm trying to imagine what you mean by this sentence: "If these "slices" have a tendency to be far away from the curve, there is too much variability in the distribution". Here's what I think you mean: at the slice on the curve where there is say a 0.1-0.3 probability of success there are lots of successes and at the slice on the curve where there is say a 0.7-0.9 probability of success there are lots of fails. Is this what you mean and would this represent overdispersion? $\endgroup$
    – luciano
    Commented Mar 28, 2014 at 6:32
  • 1
    $\begingroup$ @luciano That is the right idea. But keep in mind there has to be a balance of "slices" that are too far above and too far below the curve in order for the fit to have occurred in the first place. So it may be more realistic to say that a slice around 0.7 has too many successes (maybe 100%) and the next slice around 0.75 has too few (50%) then 0.80 has too many (100%), etc. So there is more variance observed than would be expected. $\endgroup$
    – Underminer
    Commented Mar 28, 2014 at 14:52
  • 1
    $\begingroup$ In other words, the variance in each vertical stripe should be p(1-p), where p is the probability predicted by the logistic regression. If it is bigger, you have over-dispersion. $\endgroup$ Commented Mar 27 at 21:35
10
$\begingroup$

As already noted by others, overdispersion doesn't apply in the case of a Bernoulli (0/1) variable, since in that case, the mean necessarily determines the variance. In the context of logistic regression, this means that if your outcome is binary, you can't estimate a dispersion parameter. (N.B. This does not mean that you can ignore potential correlation between observations just because your outcome is binary!)

If, on the other hand, your outcome is a set of proportions, then you can estimate a dispersion parameter (which, although often greater than one, can also be less than one) by dividing the Pearson chi-squared statistic (or the deviance) by the residual degrees of freedom.

Remember, logistic regression with a purely binary outcome is just a special case of the more general logistic regression model in which the binomial index can exceed one (and can vary across observations). Thus, the question of whether you're fitting a logistic regression model or not is unrelated to the question of whether your data are overdispersed.

$\endgroup$

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.