Helminth Burden Distribution Simulator: Negative Binomial Worm Aggregation

simulator intermediate ~10 min
Loading simulation...
15.2% heavily infected (>=100 worms), harboring 78% of the total parasite population

With mean burden 30, k=0.3, and heavy infection threshold of 100 worms, approximately 15.2% of the population is heavily infected, yet this group carries about 78% of all worms in the community. The variance (3030) far exceeds the mean (30), confirming extreme overdispersion — the hallmark of macroparasite distributions.

Formula

P(X=j) = C(j+k-1, j) * (mu/(mu+k))^j * (k/(mu+k))^k
Variance = mean + mean^2 / k
Variance-to-mean ratio = 1 + mean/k

The Overdispersion Phenomenon

One of the most robust findings in parasitology is that macroparasites (helminths, ectoparasites) are never randomly distributed among hosts. Instead, they follow a characteristically overdispersed pattern: the variance of parasite counts far exceeds the mean. In practical terms, most individuals in an infected community harbor few or no worms, while a small minority carry extraordinarily heavy burdens — often hundreds or thousands of worms.

The Negative Binomial Model

The negative binomial distribution with parameters mean (mu) and aggregation (k) provides the standard mathematical description of this overdispersion. The aggregation parameter k controls the degree of clumping: k < 1 indicates extreme overdispersion typical of most helminth-host systems, k = 1 gives a geometric distribution, and k approaching infinity recovers the Poisson (random) distribution. The variance equals mu + mu^2/k, so for Ascaris with mu=30 and k=0.3, the variance is 3030 — a hundred-fold the mean.

Epidemiological Consequences

Overdispersion has profound implications for disease control. Since morbidity is primarily determined by worm burden (not merely infection status), the small proportion of heavily infected individuals suffers most disease and also contributes disproportionately to environmental contamination and ongoing transmission. This principle underlies the WHO strategy of targeting school-age children with mass drug administration in endemic areas.

Visualizing the Distribution

Adjust the mean burden and k parameter to explore different helminth species: Ascaris lumbricoides (k ~ 0.3-0.5), hookworm (k ~ 0.3-0.4), and Trichuris trichiura (k ~ 0.2-0.4). Move the pathology threshold to see what fraction of the population exceeds clinically significant burdens. The histogram visualization clearly shows the long right tail characteristic of overdispersed distributions.

FAQ

Why are helminths overdispersed in host populations?

Overdispersion arises from heterogeneity in exposure (some individuals contact more infective stages), susceptibility (genetic and immune variation), and predisposition (previously infected individuals tend to be reinfected). This creates the characteristic pattern where most people have few worms while a minority are heavily infected.

What is the negative binomial distribution?

The negative binomial is a probability distribution with two parameters: the mean and aggregation parameter k. When k is small (<1), the distribution is highly overdispersed (right-skewed). As k approaches infinity, it converges to a Poisson distribution. It is the standard statistical model for macroparasite burdens.

Why does aggregation matter for control programs?

When parasites are highly aggregated, treating the entire population is inefficient because most people have few worms. Targeted treatment of heavily infected individuals (who harbor most parasites) removes the majority of the transmission reservoir at a fraction of the cost. This is the rationale for school-based deworming programs.

What is the 20/80 rule in parasitology?

Empirically, roughly 20% of hosts harbor approximately 80% of parasites. This Pareto-like distribution is a direct consequence of overdispersion (low k values). It implies that control programs achieving 80% worm reduction may only need to treat the most heavily infected 20% of the population.

Sources

Embed

<iframe src="https://homo-deus.com/lab/parasitology/helminth-burden/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub