---
title: Intersection Supermartingale Products
url: https://www.emergentmind.com/topics/intersection-supermartingale-products
type: topic
---

# Intersection Supermartingale Products

Intersection supermartingale products are a rigorous probabilistic apparatus combining single-stratum nonnegative supermartingale (NNSM) processes to enable composite union-intersection hypothesis tests, particularly relevant in stratified risk-limiting audits (RLAs). This construction controls Type I error guarantees nonparametrically, applies under any predictable sampling strategy, and supports adaptive allocation of sampling resources across strata. Intersection supermartingale products yield anytime-valid $P$-values and can substantially increase efficiency in practical RLA settings, offering improvements over earlier approaches such as SUITE and SHANGRLA by exploiting the generality and flexibility of NNSMs [2207.03379].

## 1. Single-Stratum Nonnegative Supermartingales

The fundamental component is the single-stratum NNSM, defined for a stratum $s$ of size $N_s$ with bounded nonnegative assorter values $X_{1s},\dots,X_{N_s s}\in[0,u_s]$. To test the one-sided null hypothesis $H_{0s}: \mu_s \leq \theta_s$ for some $\theta_s\in[0,u_s]$, samples are drawn sequentially. At each sample, a predictable, typically upward-biased estimator $\tau_{i-1,s}$ of the mean $\mu_s$ is used to define a nonnegative factor $Y_{is}(\theta_s)$ satisfying $\mathbb{E}[Y_{is}(\theta_s)\mid \text{history}]\leq 1$ under $H_{0s}$.

The single-stratum NNSM process is then
$$
M_{n_s}^{(s)}(\theta_s) = \prod_{i=1}^{n_s} Y_{is}(\theta_s),\quad M_0^{(s)} = 1.
$$
Ville’s inequality ensures $\Pr(\exists\,n_s: M_{n_s}^{(s)}\geq 1/\alpha)\leq \alpha$ under $H_{0s}$ for any $n_s$, granting finite-sample error control and anytime validity.

## 2. Intersection Supermartingale Construction

For testing the global null hypothesis $H_0^{\text{int}}: \mu_s \leq \theta_s$ for all strata $s=1,...,S$, intersection supermartingale products combine the single-stratum NNSMs into a global process. At each sampling epoch $t$, a (predictable) stratum selector $\kappa_t\in\{1,\dots,S\}$ identifies the stratum to be sampled, measurable with respect to the filtration $\mathcal{F}_{t-1}$ generated by the previous samples and stratum selections. Let $\nu^s_t$ be the total draws from stratum $s$ up to time $t$.

The update at time $t$ is:
$$
Z_t = Y_{\nu^s_t, s}(\theta_s;\, X_t)
$$
where $X_t$ is the new sample from stratum $s = \kappa_t$. The global intersection supermartingale is
$$
M_t = \prod_{i=1}^t Z_i,\quad M_0=1.
$$
By construction, $\mathbb{E}[M_t\mid\mathcal{F}_{t-1}]\leq M_{t-1}$ under $H_0^{\text{int}}$, and Ville's inequality immediately applies:
$$
\Pr\left( \exists\, t : M_t \geq 1/\alpha \right) \leq \alpha.
$$
Therefore, the first crossing time $t$ such that $M_t \geq 1/\alpha$ achieves strong control of the familywise Type I error for the intersection hypothesis.

## 3. Adaptive Stratum Selection and Efficiency

The validity of the global process for any predictable sampling rule enables adaptive allocation. Two canonical strategies are:

- **Proportional Allocation**: Sampling strata in proportion to their sizes ($N_s$), a nonadaptive but valid rule.
- **Lower-Sided Screening (Hard Stop)**: In parallel, a one-sided test is run for each stratum for the alternative $H'_{0s}:\mu_s\ge\theta_s$, using reversed bets. When this lower-sided supermartingale exceeds $1/\alpha$, sampling for that stratum ceases; subsequent $\kappa_t$ selections are restricted to remaining strata.

The lower-sided screening approach can reduce expected sample size by up to approximately $70\%$ in illustrative two-stratum cases. Once a stratum ceases to provide evidence for $H_0^{\text{int}}$, it is efficiently excluded from further sampling [2207.03379].

## 4. Empirical Results and Case Studies

Application to both real-world and simulated audits demonstrates marked efficiency gains. In the 2018 Whitmer–Schuette race pilot audit in Kalamazoo, MI (two strata: CVR $N_1=5\,294$, no-CVR $N_2=22\,732$), the SUITE approach yielded a $P$-value of $0.037$. Re-analysis with intersection supermartingales reduced this to approximately $0.003$, representing a more than order-of-magnitude reduction in measured risk, even after conditioning on observed samples.

In a large-scale simulation across 58 California counties (total $N\approx 1.75\times10^7$), employing countywise Empirical-Bernstein supermartingales combined by the intersection product, sampling proportional to county turnout and stopping at $M_t\geq 1/\alpha$ with $\alpha=5\%$, $91\%$ of $300$ replicates completed after $t=70,580$ draws. The global $P$-value computation, requiring a small LP (since the EB log-$P$ is linear in county means), took approximately $3.5$ seconds per instance on commodity hardware [2207.03379].

## 5. Practical Algorithm and Implementation

The intersection supermartingale audit can be expressed concisely as:

```python
Input: strata sizes N[1..S], upper-bounds u[1..S],
       null-means θ[1..S], risk α,
       bet-rule τ_s(i−1), selector κ_t(⋅).

for each s=1..S:
    M_s = 1.0      # stratum-wise supermartingale
    n_s = 0        # draws so far
Global M = 1.0     # intersection supermartingale

for t = 1,2,3,…:
    # Select stratum
    s = κ_t({M_s, n_s: s=1..S})

    # Draw new ballot, record assorter X
    X = sample_stratum(s)
    n_s += 1

    # Compute next factor
    θs = θ[s]; us = u[s]
    τ = τ_s(n_s - 1)
    Y = compute_factor(X, θs, us, τ)
    M_s *= Y
    M *= Y

    # Check stopping
    if M >= 1/α:
        reject intersection-null; break
# If loop ends without crossing, fail to reject
```

The function `compute_factor` is required to ensure $\mathbb{E}[Y\mid\text{past}]\leq 1$ under $\mu_s\leq\theta_s$, guaranteeing the supermartingale property. For numerical safety, truncation of $Y$-factors is allowed provided it only decreases the value, preserving the error control.

A valid omnibus $P$-value at time $t$ is $P_t = 1 \wedge \frac{1}{M_t}$, forming an anytime-valid $P$-process.

## 6. Theoretical Guarantees and Connections

Intersection supermartingale products deliver Type I error control for global intersection nulls, generalizing classical sequential tests. The supermartingale property is robust to arbitrary predictable sampling rules, encompassing both non-adaptive and adaptive stratum selection. This universal validity supports optimal resource allocation in practical RLAs and ensures computational tractability even in high-dimensional stratification (e.g., many counties).

The intersection product structure provides a rigorous foundation for union–intersection testing in stratified designs and generalizes to any family of one-sided hypotheses with stratumwise NNSMs. The compatibility with Ville’s inequality and anytime $P$-value construction supports efficient, sequential evidence accumulation. A plausible implication is that this approach could transfer to other multi-group sequential analysis tasks in statistics or machine learning involving stratified or dependent data [2207.03379].

## 7. Related Methodologies and Generalizations

Intersection supermartingale products substantially synergize techniques from SHANGRLA, SUITE, and ALPHA. They inherit the computational tractability and modularity of SHANGRLA’s reduction to mean testing, the union-intersection logic of SUITE, and the nonparametric flexibility of ALPHA’s NNSM-based construction. Key references include Stark (2022) on intersection supermartingales for union–intersection tests, Waudby‐Smith et al. (2021) and Stark (2020) on betting martingales and the SHANGRLA framework.

The underlying principles are not reliant on specifics of ballot sampling and are extensible to other contexts—anywhere nonnegative supermartingale tests for intersection nulls, stratified adaptivity, or anytime error control are vital. This suggests broad applicability wherever composite hypotheses are tested under sequential, adaptive, or stratified sample collection regimes. 

**References:**
- "Sweeter than SUITE: Supermartingale Stratified Union-Intersection Tests of Elections" [2207.03379]
- Stark (2022), “Intersection supermartingales for union–intersection tests,” Stat Sci.
- Waudby-Smith et al. (2021), “Betting martingales for RLAs.”
- Stark (2020), “SHANGRLA: A framework for RLAs.”

Source: https://www.emergentmind.com/topics/intersection-supermartingale-products