Marked Independence Series: Theory & Applications

• Marked Independence Series are defined across domains such as spatial statistics, astronomy, and combinatorics to capture full mark structures in dependency analysis.
• In spatial statistics and galaxy clustering, they quantify departures from mark independence using log-ratios and scale-indexed functional summaries like the Landy–Szalay estimator.
• In combinatorics and graph theory, they serve as generating series for marked chromatic polynomials and characterize chordal graphs via Horn hypergeometric expansions.

Across the cited literature, the expression Marked Independence Series denotes several domain-specific constructions rather than a single universal definition. In spatial statistics and astronomy, it names a scale-indexed family of diagnostics that preserves the full mark-pair structure of a marked point process while testing or quantifying departures from mark independence. In algebraic combinatorics, it denotes a multivariate generating series whose powers encode marked chromatic polynomials and related counting invariants. In graph theory, a closely related usage identifies the multivariate independence polynomial itself as a marked independence series. The common thread is that a “series” is indexed either by spatial scale or by monomial degree, and “marked” indicates that labels, types, or multiplicities are retained rather than averaged away (Takeuchi, 1 Apr 2026, Ward et al., 2024, P et al., 28 Jul 2025, Radchenko et al., 2019).

1. Product-space and dependence-resolved formulations

In the joint point-process framework for galaxy clustering, galaxies are modeled as points on the product space

$Y=\mathbb{R}^3\times\mathcal{M},$

with galaxy position $x\in\mathbb{R}^3$ and mark $m\in\mathcal{M}$. The first-order intensity $\rho^{(1)}(x,m)$ and second-order product density $\rho^{(2)}((x_1,m_1),(x_2,m_2))$ define the joint pair correlation function

$$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$

Under homogeneity and isotropy this reduces to

$$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$

where $r=|x_1-x_2|$ is the pair separation (Takeuchi, 1 Apr 2026).

The conditional mark-pair distribution at fixed separation is

$$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$

with

$$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$

The mark-independence hypothesis conditioned on $x\in\mathbb{R}^3$0 is

$x\in\mathbb{R}^3$1

equivalently

$x\in\mathbb{R}^3$2

Using the position-only pair correlation $x\in\mathbb{R}^3$3, the corresponding independence reference is

$x\in\mathbb{R}^3$4

The key diagnostic is the log-ratio

$x\in\mathbb{R}^3$5

When $x\in\mathbb{R}^3$6, the mark pair $x\in\mathbb{R}^3$7 is over-represented at separation $x\in\mathbb{R}^3$8 beyond what spatial clustering alone would produce; when $x\in\mathbb{R}^3$9, it is under-represented. The Marked Independence Series is then the collection

$m\in\mathcal{M}$0

with each $m\in\mathcal{M}$1 indexing either a continuous function over marks or a matrix over mark bins (Takeuchi, 1 Apr 2026).

This construction is explicitly designed to avoid projection. Commonly used marked summaries are treated as averages or reductions of the full joint structure. For example, weighted moments

$m\in\mathcal{M}$2

are projections of $m\in\mathcal{M}$3, while set-averaged functions and the inhomogeneous cross-$m\in\mathcal{M}$4 function discard the locality in $m\in\mathcal{M}$5 carried by $m\in\mathcal{M}$6. In this formulation, previously discussed marked effects, including assembly bias, are interpreted as projections of a non-factorizable joint structure (Takeuchi, 1 Apr 2026).

2. Scale-by-scale diagnostics in spatial point processes on spheres and convex surfaces

A second explicit use of the term arises for marked and multi-type point processes on the sphere $m\in\mathcal{M}$7 and, via a known bijection $m\in\mathcal{M}$8, on the surface of three-dimensional convex shapes. Here the geodesic distance is

$m\in\mathcal{M}$9

and the spherical cap area is

$\rho^{(1)}(x,m)$0

If data lie on a convex surface $\rho^{(1)}(x,m)$1, the image process on $\rho^{(1)}(x,m)$2 has intensity

$\rho^{(1)}(x,m)$3

where $\rho^{(1)}(x,m)$4 (Ward et al., 2024).

In the homogeneous isotropic case, for disjoint mark sets $\rho^{(1)}(x,m)$5, the paper defines the marked functional summaries

$\rho^{(1)}(x,m)$6

Under independence of $\rho^{(1)}(x,m)$7 and $\rho^{(1)}(x,m)$8,

$\rho^{(1)}(x,m)$9

The inhomogeneous analogues $\rho^{(2)}((x_1,m_1),(x_2,m_2))$0, $\rho^{(2)}((x_1,m_1),(x_2,m_2))$1, $\rho^{(2)}((x_1,m_1),(x_2,m_2))$2, and $\rho^{(2)}((x_1,m_1),(x_2,m_2))$3 are defined under intensity-reweighted moment isotropy and second-order intensity-reweighted isotropy, with the same independence baseline for $\rho^{(2)}((x_1,m_1),(x_2,m_2))$4 and the same equality $\rho^{(2)}((x_1,m_1),(x_2,m_2))$5, $\rho^{(2)}((x_1,m_1),(x_2,m_2))$6 under cross-type independence (Ward et al., 2024).

In this context, the Marked Independence Series is the coordinated family of curves

$\rho^{(2)}((x_1,m_1),(x_2,m_2))$7

interpreted jointly with global envelopes. Attraction of type $\rho^{(2)}((x_1,m_1),(x_2,m_2))$8 around type $\rho^{(2)}((x_1,m_1),(x_2,m_2))$9 at scale $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$0 corresponds to $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$1 above the null baseline and $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$2 below $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$3; repulsion corresponds to the opposite. The proposed testing mechanisms include random labeling for independent marking, parametric bootstrap of independent inhomogeneous Poisson processes for type independence, and random rotations on $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$4 when one component can be rotated while preserving the marginal structure (Ward et al., 2024).

The framework is geometric as well as inferential. Because the data may live on a convex surface rather than a plane, pair distances are computed after mapping to $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$5, and the Jacobian factor enters the intensity transformation. This permits functional summaries that account for geometry and inhomogeneity simultaneously. In the RNGC galaxy point pattern, the observed curves indicated attractive dependencies between spiral and elliptical galaxies beyond inhomogeneity: $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$6 lay above the null envelope over small to moderate $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$7, and $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$8 was below $$g((x_1,m_1),(x_2,m_2))\equiv \frac{\rho^{(2)}((x_1,m_1),(x_2,m_2))}{\rho^{(1)}(x_1,m_1)\rho^{(1)}(x_2,m_2)}.$$9 over similar scales (Ward et al., 2024).

3. Estimation, visualization, and invariance properties

The galaxy-clustering version of the Marked Independence Series is operationalized by estimating $$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$0, $$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$1, and the marginal mark law $$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$2. For the position-only component, the paper uses the Landy–Szalay estimator

$$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$3

and then sets $$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$4. For binned marks $$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$5, one estimates

$$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$6

or uses marked Landy–Szalay-style corrections, and then forms

$$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$7

For continuous marks, the paper recommends kernel smoothing or a basis expansion

$$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$8

with empirical estimator

$$g(r;m_1,m_2)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho^{(1)}(m_1)\rho^{(1)}(m_2)},$$9

to stabilize estimation (Takeuchi, 1 Apr 2026).

Visualization is integral to the construction. The cited implementations use heatmaps $r=|x_1-x_2|$0 at fixed $r=|x_1-x_2|$1, curves $r=|x_1-x_2|$2 for selected categories or quantiles, and basis-projected series based on $r=|x_1-x_2|$3. The aim is to retain scale-by-scale mark-pair locality rather than to collapse it into a single marked correlation coefficient (Takeuchi, 1 Apr 2026).

An important structural property is invariance of the joint pair correlation under mark reweighting. If

$r=|x_1-x_2|$4

then

$r=|x_1-x_2|$5

Consequently, $r=|x_1-x_2|$6 and $r=|x_1-x_2|$7 retain their interpretation regardless of importance weighting used to form summary statistics. The paper also emphasizes that invariance to monotone mark transformations is not claimed, because $r=|x_1-x_2|$8 uses $r=|x_1-x_2|$9 explicitly (Takeuchi, 1 Apr 2026).

At large scales, the binned-mark correlation matrix satisfies

$$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$0

in the large-$$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$1 limit discussed in the paper, implying an approximately rank-1 structure. This suggests that the small-scale Marked Independence Series carries information that is compressed into an effective bias description at large scales (Takeuchi, 1 Apr 2026).

4. Generating-series interpretations in hypergraphs and subspace arrangements

In algebraic combinatorics, the Marked Independence Series is a generating function attached to a hypergraph with a distinguished set of special vertices. For a hypergraph $$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$2 with special set $$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$3, a multiset $$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$4 is marked-independent if its underlying set is independent and vertices in $$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$5 appear at most once. Writing

$$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$6

the marked multivariate independence series is

$$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$7

If $$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$8, this reduces to the usual multivariate independence polynomial (P et al., 28 Jul 2025).

The corresponding marked chromatic polynomial counts marked multi-colorings. For fixed $$p(m_1,m_2\mid r)\equiv \frac{\rho^{(2)}(r;m_1,m_2)}{\rho_X^{(2)}(r)},$$9 and multiplicity vector $$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$0, a $$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$1-marked multi-coloring is a map

$$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$2

such that non-special vertices have no repeated colors, $$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$3 for all $$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$4, and for every edge $$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$5,

$$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$6

The number of such colorings is denoted $$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$7, and the paper proves that it is a polynomial in $$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$8 (P et al., 28 Jul 2025).

The central identity is

$$\rho_X^{(2)}(r)=\int\!\!\int \rho^{(2)}(r;m_1,m_2)\,dm_1\,dm_2.$$9

equivalently

$x\in\mathbb{R}^3$00

Thus the coefficients of the $x\in\mathbb{R}^3$01-th power of the Marked Independence Series coincide with the marked chromatic polynomials in $x\in\mathbb{R}^3$02 (P et al., 28 Jul 2025).

The construction extends to subspace arrangements. For an arrangement $x\in\mathbb{R}^3$03, the paper defines

$x\in\mathbb{R}^3$04

so that

$x\in\mathbb{R}^3$05

For hyperplane arrangements, the paper proves that

$x\in\mathbb{R}^3$06

has non-negative coefficients for all $x\in\mathbb{R}^3$07 (P et al., 28 Jul 2025).

A central open problem concerns hypergraphs. The paper conjectures that

$x\in\mathbb{R}^3$08

has nonnegative coefficients for all $x\in\mathbb{R}^3$09 if and only if all edges of $x\in\mathbb{R}^3$10 have even cardinality. It also proves the necessary direction in the form: if $x\in\mathbb{R}^3$11, then $x\in\mathbb{R}^3$12 must be an even hypergraph (P et al., 28 Jul 2025).

5. Graph-theoretic specialization and Horn hypergeometric expansions

For a simple graph $x\in\mathbb{R}^3$13, with $x\in\mathbb{R}^3$14, the multivariate independence polynomial is

$x\in\mathbb{R}^3$15

The paper "Independence Polynomials and Hypergeometric Series" identifies this multivariate independence polynomial itself as the marked independence series (Radchenko et al., 2019).

Its main theorem characterizes chordal graphs by the inverse expansion of this series. For a simple graph $x\in\mathbb{R}^3$16, the following are equivalent: $x\in\mathbb{R}^3$17 is chordal; the power series expansion of $x\in\mathbb{R}^3$18 at $x\in\mathbb{R}^3$19 is Horn hypergeometric; and, more generally, for every $x\in\mathbb{R}^3$20, the expansion of $x\in\mathbb{R}^3$21 is Horn hypergeometric. The proof uses a perfect elimination ordering and an upper-triangular matrix $x\in\mathbb{R}^3$22 with ones on the diagonal and $x\in\mathbb{R}^3$23 for edges $x\in\mathbb{R}^3$24 with $x\in\mathbb{R}^3$25 (Radchenko et al., 2019).

For chordal graphs, the coefficient formula is

$x\in\mathbb{R}^3$26

where

$x\in\mathbb{R}^3$27

This gives the Horn quotient recurrence

$x\in\mathbb{R}^3$28

which is rational in the multi-index $x\in\mathbb{R}^3$29 (Radchenko et al., 2019).

The converse direction reduces non-chordal graphs to induced cycles. For $x\in\mathbb{R}^3$30 with $x\in\mathbb{R}^3$31, the diagonal coefficients of $x\in\mathbb{R}^3$32 satisfy

$x\in\mathbb{R}^3$33

with asymptotic ratio

$x\in\mathbb{R}^3$34

This excludes the Horn property for cycles of length at least four and thereby for non-chordal graphs (Radchenko et al., 2019).

This graph-theoretic branch differs from the spatial-statistical uses of the term, but it preserves the same formal idea: a multivariate marked series encodes independence constraints, and structural properties of the underlying object are reflected in analytic properties of the series.

6. Related independence-testing constructions for marked or object-valued series

A plausible broader context for the term is provided by independence testing for marked, categorical, or object-valued time series, where one again studies indexed families of quantities that preserve the native data type rather than reducing to linear correlation.

For two autocorrelated time series $x\in\mathbb{R}^3$35 and $x\in\mathbb{R}^3$36, with at least one of them stationary, the shift test defines

$x\in\mathbb{R}^3$37

for an arbitrary measurable association function $x\in\mathbb{R}^3$38. If the series are independent, then the probability that the unshifted value $x\in\mathbb{R}^3$39 is among the top $x\in\mathbb{R}^3$40 values is at most

$x\in\mathbb{R}^3$41

and for large $x\in\mathbb{R}^3$42 the probability approaches

$x\in\mathbb{R}^3$43

The paper illustrates this with categorical or “marked” time series using a regularized log-odds ratio

$x\in\mathbb{R}^3$44

which peaks at $x\in\mathbb{R}^3$45 when the series share synchronized category changes (Harris, 2020).

For strictly stationary Euclidean-valued time series, a different route uses distance covariance on lag vectors

$x\in\mathbb{R}^3$46

with statistic

$x\in\mathbb{R}^3$47

Null calibration is performed by an Independent Block Bootstrap that resamples blocks from $x\in\mathbb{R}^3$48 and $x\in\mathbb{R}^3$49 separately, thereby preserving serial dependence within each series while enforcing cross-series independence in the bootstrap world (Betken et al., 2021).

For object-valued time series in a separable metric space $x\in\mathbb{R}^3$50 of strong negative type, the corresponding serial-independence problem is formulated through the auto-distance covariance

$x\in\mathbb{R}^3$51

The generalized spectral density

$x\in\mathbb{R}^3$52

and the empirical process

$x\in\mathbb{R}^3$53

lead to the Cramér–von Mises statistic

$x\in\mathbb{R}^3$54

with critical values obtained by a wild bootstrap (Jiang et al., 2023).

These constructions are not presented under the exact title “Marked Independence Series,” but they exhibit the same methodological pattern: independence is assessed through an indexed family of mark-aware or object-aware statistics, whether the index is spatial separation, graph monomial degree, time shift, or lag frequency.

7. Conceptual unification, misconceptions, and open directions

The principal misconception surrounding marked independence constructions is that classical summary statistics already preserve the relevant joint information. The spatial-point-process papers explicitly reject this view. In the galaxy-clustering formulation, marked correlation functions, weighted moments, and set-averaged summaries are projections of the underlying joint structure, whereas the Marked Independence Series retains the local dependence pattern in $x\in\mathbb{R}^3$55 at each $x\in\mathbb{R}^3$56 (Takeuchi, 1 Apr 2026). In the spherical point-process formulation, the coordinated use of $x\in\mathbb{R}^3$57, $x\in\mathbb{R}^3$58, $x\in\mathbb{R}^3$59, and $x\in\mathbb{R}^3$60 across scales serves a similar purpose: independence is not inferred from a single scalar, but from the behavior of multiple functional summaries relative to their null baselines (Ward et al., 2024).

A second misconception is that “marked independence” is synonymous with independent marking. In the cited literature, the precise null depends on context. It may mean $x\in\mathbb{R}^3$61 in a joint pair-correlation model, equality of cross-type functional summaries to their baselines, absence of hyperedges inside a marked multiset, or the vanishing of distance-covariance-based dependence at nonzero lags. The common term therefore does not fix a single hypothesis class.

The open problems are correspondingly heterogeneous. In hypergraph theory, the main unresolved question is the even-edge conjecture for the nonnegativity of $x\in\mathbb{R}^3$62 (P et al., 28 Jul 2025). In galaxy clustering, the emphasis is on estimating the full joint structure without over-aggregation and on separating observationally accessible non-factorizability from latent causal explanations (Takeuchi, 1 Apr 2026). In spherical point processes, the practical issues include intensity estimation, sensitivity to the mapping $x\in\mathbb{R}^3$63, degeneracy of rotation-based envelopes near $x\in\mathbb{R}^3$64, and variance inflation under mark imbalance (Ward et al., 2024). In time-series settings, stationarity, mixing, block-length choice, and finite-sample calibration determine whether the inferred departures from independence are trustworthy (Harris, 2020, Betken et al., 2021, Jiang et al., 2023).

Taken together, these literatures indicate that a Marked Independence Series is best understood not as a single invariant object but as a research program: preserve the mark structure, index dependence by scale or degree, and formulate independence in a way that remains visible after the indexing variable is varied rather than averaged out.