PoPDivergence: A Spectrum of Divergence Methods

• PoPDivergence is a family of divergence metrics that quantify discrepancies between structured distributions across varied fields such as epidemiology, species modeling, and path analysis.
• It encompasses multiple formulations—including KL, Bregman, and Rényi divergences—highlighting its context-dependent mathematical structure and robust estimation capabilities.
• Applications span demographic studies, spatial point processes, and discrete sequence analysis, offering practical tools for comparative inference and model optimization.

PoPDivergence denotes several distinct divergence constructions across recent literature rather than a single universally fixed object. In demographic epidemiology it is a Kullback–Leibler divergence from a reference population pyramid; in species distribution modeling from presence-only data it is a robust Bregman-type divergence on Poisson point-process intensities; in path forecasting it is the divergence induced by a proper scoring rule on track space; in discrete multidimensional sequence analysis it appears as Point Divergence Gain; and in proximal optimal transport it is the authors’ name for an infimal-convolution divergence combining transport cost with an information divergence. Taken together, these usages suggest a shared role—quantifying discrepancy between structured population-level objects—while leaving the precise mathematical form application-dependent (Fonseka et al., 20 Jan 2025, Saigusa et al., 2023, Bonnier et al., 2021, Rychtáriková et al., 2017, Baptista et al., 17 May 2025).

1. Terminological scope and competing meanings

The literature does not supply a single canonical definition of PoPDivergence. Some papers define the term explicitly, whereas others provide methodological foundations for divergence-based comparison without introducing that exact name. A concise overview is therefore useful.

Usage context	Formal object	Representative source
Population pyramids	KL divergence from a reference pyramid	(Fonseka et al., 20 Jan 2025, Wijenayake et al., 17 Sep 2025)
Presence-only SDM	Bregman divergence on PPP intensities	(Saigusa et al., 2023)
Paths and tracks	Proper-scoring-rule divergence via signatures	(Bonnier et al., 2021)
Discrete multidimensional sequences	Rényi-entropy Point Divergence Gain	(Rychtáriková et al., 2017)
Probability space / OT	Proximal optimal transport divergence	(Baptista et al., 17 May 2025)

Two adjacent literatures are especially important for orientation. First, the high-dimensional hypothesis-testing paper on statistical divergences does not explicitly introduce or define a named concept “PoPDivergence”; it instead develops the general principle that model comparison can be based on divergences estimated from samples alone (Wilkinson et al., 2024). Second, the density-ratio-estimation paper similarly does not use PoPDivergence as a named construction; its central object is the Density Ratio Metric family, which interpolates between KL-divergence and integral probability metrics (Kato et al., 2022).

A separate source of confusion is acronymic rather than conceptual. POP3D, “Policy Optimization with Penalized Point Probability Distance,” is a first-order reinforcement-learning algorithm positioned as an alternative to PPO; it is not a general PoPDivergence definition and concerns a sampled-action penalty in policy optimization rather than a generic discrepancy between population distributions (Chu, 2018). This suggests that the term should always be interpreted relative to its host literature.

2. Population-pyramid PoPDivergence in epidemiology and demography

In the PoPStat framework, PoPDivergence is a scalar summary of how far a country’s population pyramid deviates from a chosen reference population pyramid. The construction treats a population pyramid as a normalized distribution over age groups, built from five-year age bins, and defines PoPDivergence by the KL form

$$D_{KL}(P\|Q)=\sum_i P_i \log\!\left(\frac{P_i}{Q_i}\right),$$

with $Q$ the reference population (Fonseka et al., 20 Jan 2025). Larger values indicate greater divergence from the selected reference, and the interpretation of that increase depends on the demographic shape of the reference itself.

The demographic workflow is explicit. Mortality or incidence is log-transformed, every country is considered in turn as a candidate reference, PoPDivergence is computed for all countries relative to that reference, and the reference is chosen by brute-force optimization to maximize the Pearson correlation between PoPDivergence and the outcome. The resulting correlation is PoPStat. Applied to 371 diseases across 204 countries, this framework is reported to outperform indicators such as median age, GDP per capita, and Human Development Index for many causes of death, with strong examples including neurological diseases, neoplasms, maternal and neonatal disorders, and neglected tropical diseases (Fonseka et al., 20 Jan 2025).

The COVID-19 adaptation uses 2019 United Nations World Population Prospects age–sex distributions and cumulative cases and deaths per million up to 5 May 2023. In that study, Malta is selected as the optimal old-skewed reference pyramid, yielding $\mathrm{PoPStat\text{-}COVID19}$ correlations of $r=-0.860$ for cases and $r=-0.821$ for deaths, both with $p<0.001$ (Wijenayake et al., 17 Sep 2025). Sensitivity tests over twenty additional reference pyramids show that old-skewed references retain strong negative correlations, while young-skewed references flip the sign but preserve significance. The paper therefore treats reference dependence as intrinsic rather than incidental.

This demographic usage is explicitly relative rather than absolute. PoPDivergence is not interpreted as “good” or “bad” in isolation; its meaning is induced by the chosen baseline and the sign of the downstream PoPStat correlation. The literature also stresses limitations: the metric is reference-dependent, the resulting correlations do not establish causation, COVID-19 outcomes are affected by under-reporting and reporting practices, and demographic structure is only one component of vulnerability (Wijenayake et al., 17 Sep 2025). A plausible implication is that PoPDivergence is best understood here as a compressed descriptor of full age-structure geometry, designed for comparative rather than mechanistic inference.

3. Robust PoPDivergence for spatial Poisson point processes

In species distribution modeling from presence-only data, PoPDivergence is introduced through a robust minimum divergence estimator for a spatial Poisson point process. The model uses a thinned PPP intensity

$$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$

with parameter $\theta=(\beta^\top,\alpha^\top)^\top$, and the usual PPP log-likelihood

$$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$

The stated motivation is that presence-only datasets often contain heterogeneous observations such as incorrect or missing coordinates, taxonomic misidentification, taxonomic shifts, records outside the species’ typical habitat, and mixed or contaminated occurrences, all of which can destabilize maximum likelihood estimation (Saigusa et al., 2023).

The divergence is constructed as a Bregman divergence on intensity functions:

$$D_\Xi(\lambda_1,\lambda_2)=\int_{\mathscr A}\left[\Xi(\lambda_1(s))-\Xi(\lambda_2(s))-\xi(\lambda_2(s))\{\lambda_1(s)-\lambda_2(s)\}\right]ds,$$

where $Q$0. The specific PoPDivergence choice defines

$Q$1

so that the estimating equations become weighted score equations with weight $Q$2. The paper uses the Pareto type II CDF

$Q$3

and in practice fixes $Q$4. Because $Q$5, observations with small fitted intensity receive small weight, which is the intended robustification.

For $Q$6, the estimating equation takes the weighted score form

$Q$7

As $Q$8, the weights converge to $Q$9, and the method reduces to the ordinary likelihood equations. The estimator is called the minimum intensity divergence estimator, denoted $\mathrm{PoPStat\text{-}COVID19}$0. The appendix establishes unbiasedness of the estimating function, consistency under regularity conditions analogous to those in Ogata (1978), Rathbun and Cressie (1994), and Assunção and Guttorp (1999), and asymptotic normality with sandwich covariance $\mathrm{PoPStat\text{-}COVID19}$1 (Saigusa et al., 2023).

Empirically, the paper reports that MIDE behaves like MLE under no contamination, but remains centered near the true target coefficients under light and heavy contamination, where MLE becomes noticeably biased. On vascular plant data from Japan, analyzed over 27 species-region combinations, MIDE AUCs are described as similar to or better than MLE AUCs across most cases. In the detailed Carpinus laxiflora example from Chugoku-Shikoku, the reported AUC improves from $\mathrm{PoPStat\text{-}COVID19}$2 for MLE to $\mathrm{PoPStat\text{-}COVID19}$3 for MIDE with selected $\mathrm{PoPStat\text{-}COVID19}$4 (Saigusa et al., 2023). In this usage, PoPDivergence is therefore a robust estimation principle specialized to point-process intensities rather than a generic information-theoretic divergence between ordinary probability vectors.

4. Geometry-aware divergences on paths, tracks, and time series

For paths and time series, PoPDivergence appears implicitly as the generalized divergence induced by a proper scoring rule on path or track space. The formal objects are probability measures on $\mathrm{PoPStat\text{-}COVID19}$5 or on $\mathrm{PoPStat\text{-}COVID19}$6, where $\mathrm{PoPStat\text{-}COVID19}$7 is a tree-like equivalence relation that identifies time-reparametrizations and certain backtracking excursions. Discrete time series are embedded as piecewise linear paths, so the framework encompasses irregular and variable-length sequences after interpolation (Bonnier et al., 2021).

The divergence is built from the standard proper-scoring-rule construction. Given a loss $\mathrm{PoPStat\text{-}COVID19}$8, the Bayes act $\mathrm{PoPStat\text{-}COVID19}$9 minimizes $r=-0.860$0, the score is $r=-0.860$1, the entropy is $r=-0.860$2, and the divergence is

$r=-0.860$3

The track-specific implementation uses the signature map

$r=-0.860$4

with $r=-0.860$5 the tensor algebra equipped with concatenation-like multiplication. For point masses, the right-loss construction yields the explicit trajectory-to-trajectory formula

$r=-0.860$6

This is the paper’s cleanest pointwise divergence formula (Bonnier et al., 2021).

Several structural properties are established. The divergence is nonnegative; it is generally asymmetric; it is invariant under time reparametrization because $r=-0.860$7 is defined on tracks; and, for symmetric losses $r=-0.860$8, left and right constructions are related by time reversal through the antipode $r=-0.860$9. A canonical experimental choice is the squared norm $r=-0.821$0. The optimization geometry is non-Euclidean: after identifying tracks with a non-linear group $r=-0.821$1, the paper uses the Pansu derivative

$r=-0.821$2

with gradient-descent update

$r=-0.821$3

This gives a divergence-respecting analogue of gradient descent on path space (Bonnier et al., 2021).

A related, looser usage appears in diffusion-based time-series generation, where the problem is framed as preserving population-level properties of a dataset rather than only individual sample realism. There the relevant divergences are the value distribution shift

$r=-0.821$4

and the functional dependency distribution shift

$r=-0.821$5

with cross-correlation as the main dependency. The proposed PaD-TS model uses

$r=-0.821$6

and an MMD-based population penalty together with same diffusion step sampling. On the reported benchmarks, PaD-TS improves FDDS by $r=-0.821$7 and VDS by $r=-0.821$8 on average relative to Diffusion-TS while maintaining comparable individual-level authenticity (Li et al., 1 Jan 2025). This suggests a broader contemporary use of “PoP” language for divergences over population-level structure in sequence data.

5. Point Divergence Gain in discrete multidimensional data

In discrete multidimensional sequence analysis, PoPDivergence appears as Point Divergence Gain, a local Rényi-entropy difference that measures the information change caused by exchanging one count of category $r=-0.821$9 for one count of category $p<0.001$0. For a discrete distribution

$p<0.001$1

and the modified distribution $p<0.001$2, the Point Divergence Gain is

$p<0.001$3

where the underlying Rényi entropy is

$p<0.001$4

An equivalent count-based form is

$p<0.001$5

with $p<0.001$6 (Rychtáriková et al., 2017).

The interpretation is explicitly local and transition-specific. $p<0.001$7 when the exchanged bins have similar frequencies; it is negative when a rare point is replaced by a frequent one; and it is positive when a frequent point is replaced by a rare one. The parameter $p<0.001$8 controls emphasis: small $p<0.001$9 accentuates rare events, while large $$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$0 emphasizes frequent events. Special cases discussed in the paper include the nearly linear form at $$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$1, the Shannon limit as $$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$2, and vanishing limits at $$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$3 and $$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$4 (Rychtáriková et al., 2017).

Two aggregate summaries extend the local quantity to whole frame pairs or multidimensional datasets. The Point Divergence Gain Entropy is

$$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$5

and the Point Divergence Gain Entropy Density is

$$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$6

where $$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$7 indicates whether a transition $$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$8 occurs at least once. $$\lambda_\theta(s)=\lambda_\beta(s)b_\alpha(s), \qquad \log \lambda_\beta(s)=\beta^\top x(s),$$9 measures total absolute information change, whereas $\theta=(\beta^\top,\alpha^\top)^\top$0 measures the density of distinct realized transitions. The paper applies these quantities to consecutive images, video-like image sequences, microscopy data, clustering and segmentation of image stacks, autofocus, and identification of stable versus changing structures (Rychtáriková et al., 2017).

This formulation differs sharply from the population-pyramid and point-process usages. Here the divergence is fundamentally local, defined by a one-point substitution in a discrete distribution, and only secondarily aggregated to a frame- or sequence-level descriptor. The commonality with other PoPDivergence usages lies in the emphasis on structurally meaningful change rather than in a shared formalism.

6. General divergence methodology, proximal OT, and adjacent frameworks

Several papers provide a methodological backdrop for PoPDivergence-like constructions even when they do not use the name. In high-dimensional hypothesis testing, statistical divergences are proposed as the basis for comparing the population distribution of observed data with competing models when likelihoods are unavailable or intractable. The paper shows that, in the large-sample limit, the usual log-likelihood-ratio statistic selects the hypothesis with smaller KL divergence to the true distribution, and then uses the variational dual form of $\theta=(\beta^\top,\alpha^\top)^\top$1-divergences to estimate divergences from samples alone by training a neural network approximation to the log-likelihood ratio. For KL, the dual objective is

$\theta=(\beta^\top,\alpha^\top)^\top$2

with a practical recipe based on train/validation splitting and validation-set lower bounds (Wilkinson et al., 2024). This is a general divergence-estimation foundation for several later application-specific “PoP” constructions.

A parallel unification appears in density-ratio estimation, where the Density Ratio Metric family

$\theta=(\beta^\top,\alpha^\top)^\top$3

interpolates between $\theta=(\beta^\top,\alpha^\top)^\top$4, $\theta=(\beta^\top,\alpha^\top)^\top$5, and $\theta=(\beta^\top,\alpha^\top)^\top$6 at $\theta=(\beta^\top,\alpha^\top)^\top$7, $\theta=(\beta^\top,\alpha^\top)^\top$8, and $\theta=(\beta^\top,\alpha^\top)^\top$9, respectively (Kato et al., 2022). This establishes a broader perspective in which divergence choice is itself a tunable modeling decision.

For Poisson laws and point patterns, a general information-theoretic framework shows that Rényi divergences between Poisson point-process laws are exactly generalized Tsallis divergences of their intensity measures:

$$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$0

The KL specialization is

$$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$1

and the Hellinger criterion characterizes absolute continuity of Poisson laws by $$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$2 together with $$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$3 (Leskelä, 2024). This result is directly relevant to point-process interpretations of PoPDivergence because it makes the divergence of Poisson laws analytically reducible to divergence of intensity measures.

Divergence-based testing frameworks provide further neighboring methodology. In logistic regression, empirical $$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$4-divergence test statistics generalize the empirical likelihood ratio test through

$$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$5

yet retain the same asymptotic $$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$6 null distribution under regularity conditions (Felipe et al., 2021). In latent class models for binary data, goodness-of-fit and nested-model tests are likewise generalized via $$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$7-divergences while preserving asymptotic chi-square limits (Felipe et al., 2014). A Bayesian alternative, co-BPM, estimates divergences directly from two sample sets by learning a coupled binary partition and computing $$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$8 on the induced piecewise-constant densities, avoiding separate density estimation (Yang et al., 2014).

A distinct explicit naming appears in proximal optimal transport. There PoPDivergence is the proximal OT divergence

$$l(\theta)=\sum_{i=1}^m \log\{\lambda_\theta(s_i)\}-\int_{\mathscr A}\lambda_\theta(s)\,ds.$$9

an infimal convolution of an OT cost with an information divergence. The paper proves

$$D_\Xi(\lambda_1,\lambda_2)=\int_{\mathscr A}\left[\Xi(\lambda_1(s))-\Xi(\lambda_2(s))-\xi(\lambda_2(s))\{\lambda_1(s)-\lambda_2(s)\}\right]ds,$$0

and establishes the limits

$$D_\Xi(\lambda_1,\lambda_2)=\int_{\mathscr A}\left[\Xi(\lambda_1(s))-\Xi(\lambda_2(s))-\xi(\lambda_2(s))\{\lambda_1(s)-\lambda_2(s)\}\right]ds,$$1

together with primal-dual and dynamic formulations (Baptista et al., 17 May 2025). This usage is mathematically distant from population-pyramid PoPDivergence, but it preserves the same central idea of comparing structured distributions through a scalar discrepancy with operational consequences for optimization.

Across these literatures, the common denominator is not a single formula but a recurring design pattern: a structured object is first encoded as a probability measure, intensity, path signature, or transition histogram; a divergence adapted to that structure is then defined; and the resulting scalar is used for testing, estimation, optimization, or comparative interpretation. This suggests that PoPDivergence is best treated as a context-bound family of divergence constructions rather than as a uniquely standardized term.