Tutorials › Generalized Linear Mixed Models (GLMM) › Structural Zeros and Sampling Zeros

Overdispersion and Zero Inflation · Tutorial 540 of 698

Structural Zeros and Sampling Zeros

Learn to define structural and sampling zeros, represent their roles in a count mixture, and judge what evidence can—and cannot—distinguish them.

Advanced 11 min read

What You'll Learn

  • Define structural zeros relative to a specified count outcome, population, and observation period.
  • Explain how sampling zeros arise when an eligible unit has a positive event rate.
  • Use the law of total probability to describe zero probability in a latent-state count mixture.
  • Distinguish zero provenance using ecological and healthcare data-generating examples.
  • Recognize why extra observed zeros do not, by themselves, establish a structural-zero process.
  • Compare structural-zero explanations with ordinary Poisson variation, overdispersion, and hurdle interpretations.

Why the Meaning of a Zero Matters

A recorded zero is an observation, not an explanation. It may mean that the event-generating process could not operate for that unit during the observation period, or that the process could operate but happened not to produce an event. Those explanations have different implications for the population being described and for a count model.

For example, a species count of zero at a plot could reflect habitat in which the species cannot occur, or a suitable plot where the species was present but missed during sampling. A zero count of medication-related adverse events could occur in someone who was not exposed to the medication, or in an exposed patient who experienced no event. In both settings, the observed value is identical while its data-generating interpretation differs.

Definition: A structural zero is a zero that is guaranteed by the unit's state or circumstances: under the defined outcome, exposure, and observation period, the count-generating process cannot produce a positive count for that unit. A sampling zero is a zero produced by chance when the unit is eligible for the process and a positive count is possible, but no events occur in the observed period.

These definitions are relative to a precise outcome and time window. A plot may be structurally unable to produce a count of one particular species under the habitat definition being used, but not structurally unable to produce counts of other species. A patient not exposed to a drug may be structurally unable to have an adverse event caused by that drug during the specified period, but can still have other adverse events. “Structural” does not mean permanent or intrinsic; a unit's state can change across periods or under a different definition of the outcome.

A Latent-State Description of Zeros

In a count mixture, units may belong to different data-generating states. Let \(S_i=0\) indicate that unit \(i\) is in a state that cannot generate the specified event, and \(S_i=1\) indicate that it is eligible for the event-generating process. The state may be known from design information, or it may be unobserved. A structural zero occurs when \(S_i=0\); when \(S_i=1\), a zero can still occur by sampling.

To make the distinction precise, suppose that for eligible units the count \(Y_i\), conditional on the unit's rate \(\lambda_i>0\), follows a Poisson distribution. A Poisson count has probability \(e^{-\lambda_i}\) of being zero. If ineligible units have \(Y_i=0\) with certainty, the law of total probability gives:

$$ P(Y_i=0) = P(S_i=0) + P(S_i=1)P(Y_i=0\mid S_i=1) = P(S_i=0)+P(S_i=1)e^{-\lambda_i}. $$

The first term is the structural-zero contribution; the second is the sampling-zero contribution. This equation is a probability accounting identity under the stated assumptions, not a way to label an individual observed zero. It also shows why the phrase “extra zeros” needs care: even when every unit is eligible, some zeros are expected whenever counts are stochastic. If the eligible-unit count distribution is not Poisson, replace \(e^{-\lambda_i}\) by that distribution's zero probability; the distinction between states still applies.

Assumptions behind the Poisson illustration: The structural state guarantees a zero for the specified count and period. For eligible units, the rate is positive and the count is Poisson conditional on that rate and relevant covariates. If units have different rates or random effects, the zero probability must also account for that heterogeneity; a single common \(\lambda_i\) is not automatically adequate.

Worked Examples

Worked Example: Zeros in an Ecological Survey

A hypothetical survey records the number of a plant species in 12 plots. Independent habitat criteria classify 5 plots as unsuitable for the species during the survey, while 7 plots are suitable. Suppose, for illustration, that counts in suitable plots follow a Poisson distribution with mean \(\lambda=1.5\) plants per plot. Explain the two sources of zeros and calculate the expected number of sampling zeros among the 7 suitable plots.

For each unsuitable plot, the model's state definition makes the count zero: these are structural zeros for this species and survey. In a suitable plot, zero is possible but not guaranteed. Its probability is:

$$ P(Y=0\mid S=1)=e^{-1.5}=0.2231. $$

Thus the expected number of zeros among the 7 suitable plots is:

$$ 7(0.2231)=1.562 \text{ plots}. $$

This expectation does not mean that 1.562 plots must be zero; the actual number is an integer and varies from survey to survey. It does mean that sampling zeros are quite plausible even in plots known to be suitable. The habitat information supports the structural interpretation for the 5 unsuitable plots; the zero counts alone would not establish that interpretation for the other plots.

Worked Example: Medication-Related Events in Healthcare

In a hypothetical 30-day record, 100 patients are considered: 80 received a specified medication throughout the period and 20 did not receive it. The outcome is the count of adverse events caused by that medication. Suppose counts among exposed patients follow a Poisson distribution with mean \(\lambda=0.08\) per patient over 30 days. Under the stated outcome definition, patients with no exposure cannot have an event caused by this medication during the period. Calculate the expected total number of zero counts.

The 20 unexposed patients have structural zeros for this specific outcome and time window. For an exposed patient, the chance of a sampling zero is:

$$ P(Y=0\mid \text{exposed})=e^{-0.08}=0.9231. $$

The expected number of zero counts among the 80 exposed patients is \(80(0.9231)=73.85\). Adding the 20 guaranteed zeros gives:

$$ 20+80(0.9231)=93.85 \text{ zero counts, in expectation}. $$

This calculation is conditional on the hypothetical rate model and the stated exposure information; it is not a claim that the observed data must contain that number of zeros. Also, the structural-zero label applies only to medication-caused events. If the outcome were counts of all adverse events, no medication exposure would not guarantee a zero.

Worked Example: Distinguishing Low Exposure From Structural Ineligibility

A hypothetical ecological team counts insects during fixed-length visits. Both of two sites are suitable for the insect, but the count rate is 0.4 insects per hour at each site. One visit lasts 2 hours and the other lasts 10 hours. Assuming a Poisson count conditional on rate and visit duration, calculate the probability of observing zero at each site and interpret the results.

For a rate \(r\) per hour and exposure duration \(t\) hours, the Poisson mean for a visit is \(\lambda=rt\). At the 2-hour site, \(\lambda=(0.4)(2)=0.8\), so:

$$ P(Y=0)=e^{-0.8}=0.4493. $$

At the 10-hour site, \(\lambda=(0.4)(10)=4\), so:

$$ P(Y=0)=e^{-4}=0.01832. $$

Both zeros, if observed, would be sampling zeros under the assumptions: both sites are suitable and have positive rates. The shorter visit has a much higher chance of yielding a zero because it provides less observation time. A zero at the shorter visit is not evidence by itself that the site is structurally unsuitable.

What the Zero Count Can and Cannot Establish

The derivation above follows directly by conditioning on the latent state: the events \(S_i=0\) and \(S_i=1\) partition the possible states, so their contributions to \(P(Y_i=0)\) add. But a data table usually records \(Y_i\), not \(S_i\). When \(Y_i=0\), the observed value alone cannot reveal which state generated it. A structural explanation needs substantive information about eligibility, exposure, or the process—not just a high proportion of zeros.

The same caution applies at the group level. If more zeros are observed than predicted by a fitted Poisson model, that discrepancy can arise from structural zeros, but also from omitted covariates, an incorrect exposure measure, unmodeled differences in rates, or other extra-Poisson variation. The earlier material on overdispersion establishes that count variation can exceed the Poisson pattern for multiple reasons. The zero count is therefore a diagnostic clue, not a unique diagnosis.

A useful check is to ask whether the proposed structural state is genuinely unable to produce a positive value under the exact outcome definition. If it merely has a low chance of a positive value, its zeros are sampling zeros. If “unsuitable” is defined only after seeing that a plot had zero plants, the explanation risks being circular: it classifies the observed outcome rather than identifying an independent state that could have generated it.

1
Define the count precisely.
Specify the event, unit, exposure, and observation period. A zero is structural only relative to that definition.
2
Identify eligibility information.
Use design records, scientific constraints, or measured covariates to identify states in which positive counts are impossible. Do not infer impossibility merely from a zero.
3
Model eligible-unit variation.
Check the rate, exposure, and count distribution for eligible units. A positive rate does not mean a positive count is guaranteed.
4
Compare predicted and observed patterns.
Assess zeros alongside the positive counts, covariates, exposure, and group structure. A mismatch can flag model inadequacy without identifying its cause.

Alternatives, Diagnostics, and Limitations

The main alternatives differ in what they say about the source of zeros. A single Poisson model treats all units as governed by the same eligible count process, conditional on its predictors and any included random effects. A negative binomial model allows extra count variation relative to Poisson; it can consequently produce more zeros without positing a separate class that must be zero. A hurdle model separates the probability of zero from the distribution of positive counts, but that separation does not by itself establish that zeros are caused by ineligibility. A structural-state explanation instead asserts that some units cannot produce a positive count for substantive reasons. These descriptions are not interchangeable.

Approach or interpretationWhat it says about zero countsKey limitation
Poisson count processZeros occur by chance under the modeled rate and exposure.Can understate variation if rates or the count distribution are misspecified.
Negative binomial count processExtra variation can increase the frequency of zeros as well as large counts.Extra zeros do not prove a distinct structurally ineligible group.
Structural-state interpretationSome units are unable to generate a positive count for the defined outcome.Requires defensible information or assumptions about that state; zero values alone do not identify it.
Hurdle interpretationModels zero versus positive occurrence separately from positive counts.A separate zero process is not, by itself, evidence that zeros are structurally impossible in one class.

Diagnostics should therefore examine the whole count distribution, not only the zero proportion. Check whether exposure is recorded correctly, whether eligible units have comparable rates, whether zeros cluster by site or patient group, and whether the model reproduces both the zeros and the positive counts. For a GLMM, conditional assumptions about counts given random effects do not automatically guarantee that the marginal pattern across units is well represented. Persistent discrepancies call for investigating the data-generating assumptions rather than automatically declaring structural zeros.

There are also limits to what can be learned from observational counts. If eligibility is unobserved and the eligible units can themselves have very small rates, different latent explanations may fit the same observed zero frequency. Structural and sampling zero probabilities can be difficult to separate without measured state indicators, repeated observations, design information, or defensible restrictions. Even when a mixture is appropriate, an estimated latent class is a model-based description; it is not automatically a verified biological or clinical category.

Common Mistake

Common mistake: “This group has many zero counts, so its zeros are structural.” That reasoning confuses an observed frequency with a data-generating impossibility. Many sampling zeros can arise from low rates, short exposure, heterogeneity, or a misspecified count distribution. The correction is to define a state in which positive counts truly cannot occur, justify that state independently of the observed zeros, and assess sampling-zero expectations among eligible units.
Key takeaway: A structural zero is guaranteed by a clearly defined ineligible state; a sampling zero is a chance outcome from an eligible count process. An observed zero does not disclose which explanation applies.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A species is capable of living in a plot, but a short survey records none. Is this necessarily a structural zero? Explain.
  2. For eligible units with a Poisson mean of \(0.5\), calculate the sampling-zero probability.
  3. Patients without exposure to a medication have zero counts of adverse events caused by that medication. When is it appropriate to call these structural zeros?
  4. An analysis finds more zeros than a Poisson model predicts. Give two explanations other than structural zeros.
  5. Why is a hurdle model's separate zero process not sufficient evidence for structural ineligibility?

Solutions. (1) No. The plot is eligible, and the species could have been missed or absent by chance; this is a sampling zero under the stated circumstances. (2) \(P(Y=0)=e^{-0.5}=0.6065\), or about 60.65%. A positive mean does not rule out zero. (3) Only when the outcome is specifically medication-caused events and the relevant period contains no medication exposure, so such an event is impossible by definition. For all adverse events, no exposure would not guarantee zero. (4) For example, omitted differences in rates or extra-Poisson variation such as overdispersion; inaccurate exposure measurement is another possibility. (5) The hurdle model describes a statistical separation of zero and positive outcomes; that construction does not establish that any units are scientifically unable to generate positive counts.

Clinical Biostats