Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper introduces two novel kernel independence tests, $m\mathrm{HSIC}$ and $md\mathrm{HSIC}$, which address the computational bottleneck of permutation calibration required by traditional Hilbert-Schmidt Independence Criterion (HSIC) based tests. By adapting a martingale MMD construction, the proposed tests yield studentized statistics with standard normal null distributions, eliminating the need for computationally expensive permutations. Empirical results demonstrate that $m\mathrm{HSIC}$ and $md\mathrm{HSIC}$ achieve comparable type-I error rates and test power to permutation-calibrated baselines while offering a 25-60x speedup.
Kernel independence testing just got 60x faster: martingale methods sidestep the permutation bottleneck of HSIC, enabling efficient independence checks even with high-dimensional data.
The Hilbert-Schmidt Independence Criterion (HSIC) and its joint-independence extension $d\mathrm{HSIC}$ are degenerate $V$-statistics whose data-dependent weighted-$χ^2$ null limits force a permutation calibration that multiplies the per-test cost by the number of permutations, in practice two orders of magnitude. Adapting the recent martingale MMD construction for two-sample testing to the (joint) independence problem, we introduce two studentised statistics whose null distributions are standard normal regardless of the data law, so that a single normal-quantile lookup replaces the permutation step entirely. The first, $m\mathrm{HSIC}$, is a self-normalised lower-triangular sum of the Hadamard product of two empirically centred Gram matrices. Under independence and bounded-fourth-moment kernels it converges to a standard normal. It is consistent against every fixed alternative, and runs at quadratic cost in the sample size without any sample split, matching the biased HSIC $V$-statistic. Our second statistic, $md\mathrm{HSIC}$, achieves finite-sample consistency with a single half-sample split: the centring is estimated on one half and the lower-triangular self-normalised martingale is run on the other, shrinking the conditional-mean residual to a quantity that is exponentially small in $d$, so the statistic is asymptotically standard normal at every fixed number of jointly tested variables, with a per-test cost that grows only linearly in $d$. On synthetic data with per-variable input dimension from $1$ to $500$ and between $2$ and $10$ jointly tested variables, both statistics match the empirical type-I error rate and test power of permutation-calibrated baselines while running $25$ to $60\times$ faster.