← Back to arXiv
arXivProbabilityarXiv:2607.21931

Lorenz and convex ordering of parasite burden distributions with density-dependent deaths

The paper studies how parasites accumulate in a host when the death rate of parasites depends on how many are already present. This is a realistic biological assumption: crowding effects, immune responses, or resource competition often mean that each additional parasite faces a higher (or lower) chance of dying depending on the total burden already carried. The mathematical object of interest is the probability distribution describing how many parasites a randomly chosen host ends up carrying at equilibrium under these competing forces of acquisition and density-dependent death.

The central result concerns a way of comparing two such distributions called "convex ordering." One distribution is said to be less than another in the convex order if it is less spread out in a precise probabilistic sense that goes beyond just comparing variances. The paper proves that if you take two different density-dependent death regimes and the ratio of their death rates is non-increasing as parasite counts grow, then the resulting equilibrium distributions can be ranked in the convex order. Intuitively, a death rate that falls off more slowly with crowding keeps parasite numbers more tightly clustered, while one that eases off quickly allows more extreme variation in burden across hosts.

A notable application of this general result is a cleaner understanding of the Conway-Maxwell-Poisson distribution, a flexible statistical model that can describe both under- and over-dispersed count data relative to the standard Poisson. The paper shows that as the dispersion parameter of this distribution changes, the family of distributions with the same mean is convex-ordered by that parameter. This gives a rigorous mathematical foundation for the intuition that the dispersion parameter genuinely controls spread in this strong ordering sense, which has practical value for anyone using this distribution to model parasite loads or other biological count data.

Read original →