Higher-Order Statistics

Skewness, kurtosis, tail asymmetry, non-Gaussianity
dynamicaldim mixed (3rd/4th order)6 metrics

What It Measures

The statistical fingerprint beyond mean and variance — skewness, kurtosis, and nonlinear frequency coupling.

Projects pairs of consecutive values onto 20 directions in 2D and computes the skewness and kurtosis along each projection. Also computes permutation entropy (how many ordinal patterns appear) and bicoherence (whether triplets of frequencies are phase-locked). Together these detect non-Gaussianity from four independent angles: directional asymmetry, tail heaviness, ordinal structure, and nonlinear coupling.

Metrics

bicoherence_max

Are any frequency triplets phase-locked? Measures whether the phases of frequencies f1, f2, and f1+f2 are correlated across segments, normalized by sqrt(P1*P2*P3) — a non-standard normalization that is unbounded above 1.0 (standard bicoherence uses a different denominator). Devil's Staircase scores 7.92: its self-similar staircase structure creates strong quadratic phase coupling between harmonics. Quantum Walk (5.54) is next — interference creates phase-locked frequency relationships. Linear systems and noise score near zero because their phases are independent.

kurt_max

Maximum projection kurtosis across 20 directions. Heavy tails in any direction push this up. Rainfall dominates at 1373: most hours are dry (near zero), but rare downpours create extreme outliers. Forest Fire (886) has similar bursty dynamics. Gaussian noise scores near 0 (the excess kurtosis baseline — Fisher's definition, not Pearson's). Constants and periodic orbits score 0.0 (degenerate single-point projections).

perm_entropy

What fraction of ordinal patterns are observed, weighted by their frequencies? Pi Digits and XorShift32 score 0.999 (all 120 possible 5-element rank orderings appear, roughly equally). Forest Fire scores 0.008: its avalanche dynamics visit almost no ordinal patterns. This is encoding-invariant — it ignores absolute values entirely.

perm_forbidden

What fraction of ordinal patterns never appear? Logistic Period-2 scores 0.983 (only 2 of 120 patterns observed — alternating between two values constrains the rank structure almost completely). Noise, prime gaps, and pink noise score 0.0 (all patterns appear). This is a zero-training determinism detector: forbidden patterns are the fingerprint of low-dimensional dynamics.

skew_mean

Average absolute skewness across projections. Rainfall (20.8) and Forest Fire (8.32) have the most asymmetric distributions: rare large events create long right tails in nearly every direction. OpenBSD ELF (6.28) is skewed because binary executables have non-uniform byte distributions. Symmetric distributions (Gaussian noise, sine waves) score near zero.

c3_energy

Energy of the third-order cumulant (bispectrum). Measures nonlinear coupling strength beyond what bicoherence_max captures — the total energy rather than just the peak. Not yet in the atlas — recently added.

Atlas Rankings

bicoherence_max
SourceDomainValue
De Bruijn Sequencenumber_theory1.0000
Logistic Edge-of-Chaoschaos1.0000
L-System (Dragon Curve)exotic0.9999
···
Logistic r=3.5 (Period-4)chaos0.0000
Logistic r=3.2 (Period-2)chaos0.0000
Thue-Morseexotic0.0163
c3_energy
SourceDomainValue
Lotka-Volterrabio3.7755
Forest Fireexotic3.3996
Hodgkin-Huxleybio2.7348
···
Logistic r=3.2 (Period-2)chaos0.0000
Thue-Morseexotic0.0001
Noisy Period-2chaos0.0016
kurt_max
SourceDomainValue
Rainfall (ORD Hourly)climate1400.3451
Forest Fireexotic885.8015
Accel Sitmotion647.7791
···
Logistic r=3.2 (Period-2)chaos0.0000
Neural Net (Dense)binary0.0787
ARMA(2,1)noise0.0827
perm_entropy
SourceDomainValue
Pi Digitsnumber_theory0.9994
glibc LCGbinary0.9992
XorShift32binary0.9992
···
Forest Fireexotic0.0084
Square Wavewaveform0.0120
Sawtooth Wavewaveform0.0128
perm_forbidden
SourceDomainValue
Logistic r=3.2 (Period-2)chaos0.9833
Logistic r=3.83 (Period-3 Window)chaos0.9750
Logistic r=3.5 (Period-4)chaos0.9667
···
Sunspot Numberastro0.0000
Sandpileexotic0.0000
Coupled Map Latticechaos0.0000
skew_mean
SourceDomainValue
Rainfall (ORD Hourly)climate20.6042
Forest Fireexotic8.3203
OpenBSD ELF x86-64binary6.2805
···
Logistic r=3.5 (Period-4)chaos0.0000
Logistic r=3.2 (Period-2)chaos0.0000
Thue-Morseexotic0.0001

When It Lights Up

Higher-Order Statistics separates signals that look identical under mean-and-variance analysis. Two sources with the same entropy and autocorrelation can have completely different kurtosis profiles if one has heavy tails and the other doesn't. Bicoherence is the only metric in the framework that detects nonlinear frequency coupling — linear filters destroy it, so high bicoherence is a direct signature of nonlinear dynamics. In the atlas, the combination of kurt_max and perm_forbidden separates bursty processes (high kurtosis, low forbidden) from low-dimensional chaos (moderate kurtosis, high forbidden).

Open in Atlas
← Spectral AnalysisRecurrence Quantification →