Pith. sign in

REVIEW 3 major objections 6 minor 1 cited by

Correct Estimation of Higher-Order Spectra: From Theoretical Challenges to Practical Multi-Channel Implementation in SignalSnap

T0 review · 3 major / 6 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Cumulant-based estimator removes spurious offsets from higher-order spectra

desk verdict A genuinely useful polyspectral toolbox and normalization fix, with a real but scoped i.i.d. caveat that the abstract overstates. read the letter →

arxiv 2505.01231 v2 pith:AEWQR74Q submitted 2025-05-02 physics.data-an eess.SPquant-ph

Markus Sifft , Armin Ghorbanietemad , Fabian Wagner , Daniel Hägele This is my paper · ORCID
classification physics.data-aneess.SPquant-ph
keywords higher-orderspectrapolyspectraspectralestimationunbiasedcumulantestimatorsmultichannelanalysistimeseriesk-statisticstrispectrum
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper claims that the long-standing artifacts in higher-order spectral estimation come from a solvable mismatch: existing estimators compute fourth-order spectra from moments or incomplete cumulants, leaving a spurious offset even for white Gaussian noise. It reformulates polyspectral estimation with multivariate k-statistics, which are finite-sample-unbiased estimators of cumulants, and adds window normalization factors derived from the Fourier transform of the window. The result is claimed to be the first correct software implementation of the theoretical cumulant-based trispectrum, including multi-channel cross-spectra and quasi-polyspectra for non-stationary signals. If the claim holds, higher-order spectra become a quantitative tool: experimental bispectra and trispectra can be matched to theoretical predictions and computed on datasets of hundreds of gigabytes.

What carries the argument

The load-bearing object is the multivariate k-statistics estimator: for $m$ samples, each empirical cumulant is a polynomial in sample means with $m$-dependent prefactors chosen so that its expectation equals the true cumulant at finite $m$, not only as $m$ tends to infinity. For the fourth order, the estimator combines centered fourth-, third-, and second-order products with prefactors $m+1$ and $m-1$; omitting those combinations is exactly what produces the Gaussian-noise offset. Equally important is the normalization identity connecting the finite-window estimate to the ideal spectrum: each spectral order is the Fourier-coefficient cumulant divided by a window-power integral, derived from convolving the ideal spectrum with the window transform. Together these two pieces turn raw discrete Fourier coefficients into unbiased, comparable spectral estimates.

What would settle it

Generate a stationary Gaussian process whose correlation time is longer than the analysis window and estimate its fourth-order spectrum with the proposed k-statistics estimator for increasing $m$; the unbiasedness claim predicts a statistically zero trispectrum, so a systematic nonzero or $m$-dependent result would show that it is the i.i.d. precondition, not the estimator, that removes the artifact.

Watch

Extended reading notes

Core claim

The central claim is that higher-order spectra can be estimated without bias or spurious structure by replacing natural sample-moment estimators with multivariate k-statistics. For a single channel the third- and fourth-order spectra are estimated from cumulants such as $C_3(a_k, a_l, a_{k+l}^*)$ and $C_4(a_k, a_k^*, a_l, a_l^*)$ of discrete Fourier coefficients, each normalized by $N/(T\sum_i g_i^2 g_i^*)$ or the corresponding window-power integral; the paper derives these normalizations by convolving the ideal spectrum with the window's Fourier transform. The k-statistics prefactors $m/(m-1)$, $m^2/[(m-1)(m-2)]$, and the $m$-dependent fourth-order prefactors make the estimate unbiased for every finite number $m$ of windows and consistent as $m$ grows. On white Gaussian noise, the resulting trispectrum is statistically zero, whereas the moment-based estimator used by existing tools shows a pronounced diagonal offset. The same construction is extended to multichannel cross-polyspectra of up to four channels and to quasi-polyspectra, whose dependence on $m$ flags non-stationarity.

Load-bearing premise

The entire unbiasedness guarantee rests on the assumption that Fourier coefficients from consecutive windows are independent and identically distributed, which is only approximate when the signal decorrelates within the window and is violated by unresolved sharp spectral lines.

Editorial extensions

If this is right

  • For fourth-order spectra, Gaussian noise will estimate to zero within error bars, so true non-Gaussian signatures are no longer masked by a fake offset.
  • Experimental trispectra can be compared quantitatively with theoretical cumulant predictions, for example in continuous quantum measurements where transition rates are extracted from spectral shape.
  • Multi-channel cross-polyspectra up to four channels provide a statistically correct test for inter-detector correlations, including imaginary parts that signal broken time-reversal symmetry.
  • Computing spectra with several values of $m$ gives an $m$-dependence that can reveal non-stationary, time-dependent structure in otherwise Gaussian-looking signals.
  • The GPU implementation makes bispectrum and trispectrum analysis practical on datasets exceeding hundreds of gigabytes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the i.i.d. assumption fails, say for a process with correlation time longer than the window or with an unresolved spectral line, the unbiasedness guarantee is not automatic; the paper's own caveat suggests a practical screening test: check whether the fourth-order estimate drifts with the number of windows $m$.
  • Because the artifact in existing tools is a missing lower-order correction, a simple white-Gaussian-noise test can identify any library that silently uses moment-based fourth-order estimators; a reader can run that test before trusting a trispectrum.
  • The frequency-shift trick in the appendix implies the full three-dimensional trispectrum can be assembled from the implemented two-dimensional cut by estimating spectra of artificially frequency-shifted complex signals, so the 2-D limitation is practical rather than fundamental.
  • Quasi-polyspectra could be developed into a quantitative non-stationarity statistic: under stationarity the spectrum should be $m$-independent within errors, and deviations would measure the rate of parameter drift.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper addresses long-standing estimation problems in higher-order spectral analysis. It reformulates single- and multi-channel bispectrum and trispectrum estimation using multivariate k-statistics (Eqs. 32-34), derives window-normalization factors that relate windowed discrete Fourier coefficients to Brillinger's cumulant-based polyspectra (Sec. 3 and Appendix F), and implements the resulting estimators in the open-source GPU library SignalSnap. The paper demonstrates that a moment-based estimator used in pyHOSA produces a nonzero trispectrum for white Gaussian noise, while the proposed cumulant-based estimator does not, and it presents applications to constructed signals, a two-channel telegraph-noise example, quasi-polyspectra for non-stationary signals, and a symmetry analysis of multi-channel spectra.

Significance. If the central claim holds, the paper makes an important practical contribution: a principled, finite-sample-unbiased, cumulant-based implementation of bispectra and trispectra with correct window normalization, including multi-channel cross-spectra, backed by a GPU implementation capable of handling large datasets. The derivations in Sec. 3 and Appendix F are internally consistent, the Gaussian-noise benchmark behaves as predicted, and the symmetry classification in Sec. 9 is useful. The main caveat is that the statistical guarantees are conditional on an i.i.d. assumption on window Fourier coefficients that is only approximately met for decorrelating processes; this condition is acknowledged in Sec. 4.1 but not quantified and is absent from the headline claims. The paper also provides an open-source implementation, which is a concrete strength even though the manuscript does not contain machine-checked proofs.

major comments (3)
  1. [Sec. 4.1, Eqs. (32)-(34)] The unbiasedness and consistency of the k-statistics estimators are established only for i.i.d. samples. The manuscript states in Sec. 4.1 that this requirement is “met approximately” by processes that lose memory within the window length T and is violated for sharp spectral structures, but the abstract and conclusion present the estimators as unbiased and artifact-free without this qualification. For a stationary process with correlation time comparable to T, successive Fourier-coefficient vectors are dependent; then E[c2(x,x)] differs from the marginal variance by cross-window covariance terms, and the fourth-order cancellations that make the Gaussian trispectrum vanish are no longer exact. This is a load-bearing scoping gap rather than an algebraic contradiction: the central artifact-free claim is established only for processes that decorrelate well within T. The authors should quantify the decorrelation condition, for example by giving a bound on the bias in terms of the autocorrelation at lag T or by a numerical study varying the correlation time relative to T, and should amend the abstract and conclusion to state the limitation explicitly.
  2. [Sec. 4.1, Eqs. (32)-(34)] The central estimator formulas are attributed to the authors' preprint [6], and no proof or precise statement of the underlying theorem is included in this manuscript. Because the title and abstract rest on the unbiasedness of these estimators, the paper should either reproduce the derivation, state the theorem with its conditions, or cite a peer-reviewed source; relying on an arXiv preprint for the core mathematical claim is not sufficient for a journal publication. At minimum, the authors should clarify which parts of the unbiasedness proof are new to this paper and which are taken from [6].
  3. [Sec. 3-4, Eqs. (23), (28)-(31), Appendix F] The estimators are unbiased with respect to the cumulants of windowed Fourier coefficients; the step from these cumulants to Brillinger's ideal spectrum replaces a convolution with the window kernel by the kernel area multiplied by the spectrum at the center frequency. For spectra that vary appreciably on the scale of the window mainlobe, this introduces a deterministic bias in the spectral estimator that is independent of m. The paper should state explicitly that the “unbiased and consistent” claim refers to the cumulant part of the estimator, and that the window-convolution approximation is controlled only when the true spectrum is sufficiently smooth on the frequency-resolution scale. This distinction is important for the claim that the estimates “precisely align with Brillinger’s theoretical definitions.”
minor comments (6)
  1. [Throughout] There are numerous typographical errors, including “o ffer” in the abstract, “e fficiently” in the introduction, “a such an offset” in Sec. 1, and “negtive maginary part” in the caption of Fig. 7. The manuscript should be carefully proofread.
  2. [References] Reference [40] has an incomplete URL (“https://.com/arrayfire/arrayfire”); the correct repository path should be provided.
  3. [Eq. (37)] Equation (37) contains a misplaced closing parenthesis and a stray comma in the integral; the equation should be typeset correctly.
  4. [Figure 4 caption] The caption for Fig. 4 says “ax1(t)+ax3(t)” where the text defines y3(t) = x1(t)x2(t)+ax1(t)+ax2(t); the caption should be corrected.
  5. [Sec. 4.3] The claim that “no existing software library correctly implements Brillinger’s cumulant-based trispectrum” is stronger than the evidence presented, since the review covers HOSA and pyHOSA only. The sentence should be qualified as referring to the libraries examined.
  6. [Sec. 4.2, Eq. (36)] The standard error formula in Eq. (36) would benefit from an explicit statement that x denotes the sample mean over the Np parts and that Np is the number of independent parts; currently this is only implicit in the surrounding text.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the estimators derive from Brillinger's cumulant definition and are benchmarked externally; the one self-citation to the k-statistics formulas is parameter-free and not load-bearing.

full rationale

I walked the derivation chain from Brillinger's definition, Eq. (24), through the windowing and convolution relations, Eqs. (21) and (F.2)-(F.7), to the normalized estimators in Eqs. (28)-(31) and the multi-channel forms in Eqs. (47)-(50). The estimators are not fitted to the benchmark outputs: the Gaussian-noise zero trispectrum in Fig. 3 and the constructed-signal spectra in Figs. 4-6 are consequences of the cumulant definitions and are verified numerically, not used as fit targets. The only self-citation is to the authors' prior work [6] for the explicit multivariate k-statistics in Eqs. (32)-(34) and their unbiasedness. That cited result is parameter-free, its stated i.i.d. assumption is disclosed in Sec. 4.1, and it does not already contain the spectral application or the claimed artifact-free behavior; the formulas themselves appear in the present paper and are externally checkable. The acknowledged limitation that the i.i.d. condition is only approximate for processes with memory comparable to the window length is a scoping caveat, not a circular step. No equation reduces to its input by construction, and no fitted parameter is renamed as a prediction. The derivation is therefore self-contained with respect to the circularity concerns considered here.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The ledger has no fitted free parameters: m, the number of windows, and the window width are implementation defaults, not coefficients adjusted to make the demonstration succeed. The load-bearing inputs are the multivariate k-statistics estimators from the authors' earlier preprint and the domain assumption that windowed Fourier coefficients are i.i.d.; both are identified and discussed in the text.

assumptions (4)
  • standard math The multivariate k-statistics estimators c3 and c4 in Eqs. (33)-(34) and (H.1) are unbiased and consistent for the corresponding cumulants.
    Invoked in Sec. 4.1; formulas are quoted from the authors' earlier preprint [6] and not re-derived in this paper. They are standard in form, but the paper itself does not supply the proof.
  • domain assumption Fourier coefficients a_k from consecutive time windows are independent and identically distributed.
    Stated in Sec. 4.1: this requires the process to lose memory within the window length T and to be stationary; sharp spectral structures are explicitly acknowledged to violate the assumption.
  • domain assumption The discrete windowed Fourier coefficient a_k is well approximated by the convolution of the continuous Fourier transform of z with the window.
    Used in Eq. (20) and Appendix F; requires smooth z(t) and negligible aliasing between the T/N grid and the continuous integral.
  • domain assumption Signals are stationary with finite memory for standard polyspectra; quasi-polyspectra relax this heuristically.
    Stationarity is the basis of the delta-function in Eq. (24); Sec. 8 deliberately applies the estimators to non-stationary signals and labels the result quasi-polyspectra.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Correct Estimation of Higher-Order Spectra: From Theoretical Challenges to Practical Multi-Channel Implementation in SignalSnap." pith.science (2026). https://pith.science/paper/AEWQR74Q

@misc{pith2026250501231,
  author       = {Pith},
  title        = {Pith review of: Correct Estimation of Higher-Order Spectra: From Theoretical Challenges to Practical Multi-Channel Implementation in SignalSnap},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AEWQR74Q}},
  note         = {Machine review of arXiv:2505.01231}
}
read the original abstract

Higher-order spectra (Brillinger's polyspectra) offer powerful methods for solving critical problems in signal processing and data analysis. Despite their significant potential, their practical use has remained limited due to unresolved mathematical issues in spectral estimation, including the absence of unbiased and consistent estimators and the high computational cost associated with evaluating multidimensional spectra. Consequently, existing tools frequently produce artifacts, no existing software library correctly implements Brillinger's cumulant-based trispectrum, or fail to scale effectively to real-world data volumes, leaving crucial applications like multi-detector spectral analysis largely unexplored. In this paper, we revisit higher-order spectra from a modern perspective, addressing the root causes of their historical underuse. We reformulate higher-order spectral estimation using recently derived multivariate k-statistics, yielding unbiased and consistent estimators that eliminate spurious artifacts and precisely align with Brillinger's theoretical definitions. Our methodology covers single- and multi-channel spectral analysis up to the bispectrum (third order) and trispectrum (fourth order), enabling robust investigations of inter-frequency coupling, non-Gaussian behavior, and time-reversal symmetry breaking. Additionally, we introduce quasi-polyspectra to uncover non-stationary, time-dependent higher-order features. We implement these new estimators in SignalSnap, an open-source GPU-accelerated library capable of efficiently analyzing datasets exceeding hundreds of gigabytes within minutes. In applications such as continuous quantum measurements, SignalSnap's rigorous estimators enable precise quantitative matching between experimental data and theoretical models.

Figures

Figures reproduced from arXiv: 2505.01231 by the authors.

Figure 1
Figure 1. Comparison of RC-filtered white noise (upper row) and random telegraph noise (lower row). (a) Samples of the time-dependent signals are clearly distinct. [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. SignalSnap applies the Approximate Confined Gaussian Window [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Comparison of trispectra estimated for white Gaussian noise (length [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: Comparison of signals y1(t) to y4(t) and their spectra S (2) , S (3) , S (4). The signals y j(t) are constructed from Gaussian signals x1(t) and x2(t) via linear and nonlinear operations. Row 1: The signal of the linear superposition y1(t) = x1(t) + x2(t) shows Lorentz…
Figure 5
Figure 5. Figure 5: Polyspectra of a non-Gaussian signal z(t) [see Eq. (45)]. Spectra from unbiased and biased estimators are compared. The unbiased spectra S (2) z and S (3) z are scaled down by a constant prefactor m/(m − 1) and m 2 /[(m − 1)(m − 2)]. The biased spectrum S (4) z exhibit…
Figure 6
Figure 6. Figure 6: (a) Samples of signals u(t) and x(t), where u(t) is random telegraph noise. The signal x(t) represents the position of a stochastically driven harmonic oscillators whose frequency switches between two frequencies as u(t) changes values. (b) The power-spectra S (2) uu a…
Figure 7
Figure 7. Figure 7: (a) Position x(t) and velocity v(t) of the process discussed in the text. (b) Power-spectra of x(t) (dashed line) and v(t) (solid line). (c) The cross-correlation spectrum S (2) vx has a zero real part and a negtive maginary part with peaks at 1 kHz and 2 kHz showing t…
Figure 8
Figure 8. Figure 8: Two-channel bispectra S (3) and trispectra S (4) of the signals u(t) and x(t) (see [PITH_FULL_IMAGE:figures/full_fig_p014_8.png]
Figure 9
Figure 9. Figure 9: (a) Quasi-polyspectra S (4) of a non-stationary signal x(t) that corresponds to a stochastically driven harmonic oscillator with a frequency that slowly drifts from 0.3 kHz to 0.7 kHz (see text). The spectra are an average of short-time estimates that are calculated ea…
Figure 10
Figure 10. Figure 10: Symmetries of polyspectra S (3) and S (4) regarding their dependency on ω1 and ω2. The red and black ”F” labels sectors where the spectral values are the complex conjugate of each other. Spectra (b) and (d) are equivalent and can be mapped onto each other. S (4) xyxy(…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Bandable Cumulant Tensors: Optimal Estimation and Applications in Non-Gaussian Data Modeling

    stat.ME 2026-08 conditional novelty 7.0 of 10

    A tapered sample cumulant estimator is minimax rate-optimal for bandable higher-order cumulant tensors under spectral norm when the bandwidth is near the oracle value and the sample size is sufficiently large.

Reference graph

Works this paper leans on

42 extracted references · 31 canonical work pages · cited by 1 Pith paper

  1. [6]

    Schefczik, D

    F. Schefczik, D. H¨agele, Ready-to-use unbiased estimators for multivariate cumulants including one that outperforms x3 (2019). arXiv:1904.12154. URL https://arxiv.org/abs/1904.12154

    arXiv 2019

  2. [1]

    Birkelund, A

    Y . Birkelund, A. Hanssen, E. J. Powers, Multitaper estimators of polyspec- tra, Signal Process. 83 (2003) 545

    work page 2003

  3. [2]

    D. R. Brillinger, An Introduction to Polyspectra, The Annals of Mathemat- ical Statistics 36 (5) (1965) 1351–1374. doi:10.1214/aoms/1177699896

    arXiv 1965

  4. [3]

    Trapp, P

    A. Trapp, P. Wolfsteiner, Estimating higher-order spectra via filtering-averaging, Mech. Syst. Signal Process 150 (2021) 107256. doi:10.1016/j.ymssp.2020.107256

    arXiv 2021

  5. [4]

    J. F. Kenney, E. S. Keeping, Mathematics of Statistics: Part two, 2nd Edition, Van Nostrand, Princeton, New Jersey, 1951

    work page 1951

  6. [5]

    Starosielec, R

    S. Starosielec, R. Fainblat, J. Rudolph, D. H¨agele, Two-dimensional higher order noise spectroscopy up to radio frequencies, Rev. Scientific Instrum. 81 (2010) 125101

    work page 2010

  7. [7]

    Harris, On the use of windows for harmonic analysis with the dis- crete fourier transform, Proceedings of the IEEE 66 (1) (1978) 51–83

    F. Harris, On the use of windows for harmonic analysis with the dis- crete fourier transform, Proceedings of the IEEE 66 (1) (1978) 51–83. doi:10.1109/PROC.1978.10837

    arXiv 1978

  8. [8]

    Brillinger, Time series

    D. Brillinger, Time series. data analysis and theory. repr, Time Series: Data Analysis and Theory V ol. 36 (01 2001). doi:10.1137/1.9780898719246

    work page doi:10.1137/1.9780898719246 2001

Show all 42 references
  1. [9]

    Papoulis, S

    A. Papoulis, S. U. Pillai, Probability, Random Variables, and Stochastic Processes, 4th Edition, McGraw-Hill, 2002

    2002

  2. [10]

    Oppenheim, R

    A. Oppenheim, R. Schafer, J. Buck, Discrete-Time Signal Processing, V ol. 23, 1999

    1999

  3. [11]

    Black, M

    F. Black, M. Scholes, The pricing of options and corporate liabilities, Journal of Political Economy 81 (3) (1973) 637–654. URL http://www.jstor.org/stable/1831029

    1973

  4. [12]

    J ¨ager, A

    M. J ¨ager, A. Kiel, D.-P. Herten, F. A. Hamprecht, Anal- ysis of single-molecule fluorescence spectroscopic data with a markov-modulated poisson process, ChemPhysChem 10 (14) (2009) 2486–2495. arXiv:https: //chemistry- europe.onlinelibrary.wiley.com/doi/pdf/10.1002/cphc.20090...

    2009 doi

  5. [13]

    A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, R. J. Schoelkopf, Introduction to quantum noise, measurement, and amplification, Rev. Mod. Phys. 82 (2010) 1155. doi:10.1103/RevModPhys.82.1155

    2010 doi

  6. [14]

    H¨agele, F

    D. H¨agele, F. Schefczik, Higher-order moments, cumulants, and spectra of continuous quantum noise measurements, Phys. Rev. B 98 (2018) 205143

    2018

  7. [15]

    Sifft, A

    M. Sifft, A. Kurzmann, J. Kerski, R. Schott, A. Ludwig, A. D. Wieck, A. Lorke, M. Geller, D. H¨agele, Quantum polyspectra for modeling and evaluating quantum transport measurements: A unifying approach to the strong and weak measurement regime, Phys. Rev. Res. 3 (2021) 033123....

    2021 doi

  8. [16]

    Sifft, D

    M. Sifft, D. H¨agele, Random-time quantum measurements, Phys. Rev. A 107 (2023) 052203. doi:10.1103/PhysRevA.107.052203. URL https://link.aps.org/doi/10.1103/PhysRevA.107.052203

    2023 doi

  9. [17]

    Sifft, A

    M. Sifft, A. Kurzmann, J. Kerski, R. Schott, A. Ludwig, A. D. Wieck, A. Lorke, M. Geller, D. H ¨agele, Quantum polyspectra approach to the dynamics of blinking quantum emitters at low photon rates without binning: Making every photon count, Phys. Rev. A 109 (2024) 062210. doi:...

    2024 doi

  10. [18]

    Dutta, P

    P. Dutta, P. M. Horn, Low-frequency fluctuations in solids: 1 f noise, Rev. Mod. Phys. 53 (1981) 497–516. doi:10.1103/RevModPhys.53.497. URL https://link.aps.org/doi/10.1103/RevModPhys.53.497

    1981 doi

  11. [19]

    Blanc-Lapierre, R

    A. Blanc-Lapierre, R. Fortet, Th´eorie des fonctions al´eatoires. Applications a divers ph´enom`enes de fluctuation, 1st Edition, Masson, Paris, 1953

    1953

  12. [20]

    J. M. Mendel, Tutorial on Higher-Order Statistics (Spectra) in Signal Pro- cessing and System Theory: Theoretical Results and Some Applications, Proc. IEEE 19 (1991) 278

    1991

  13. [21]

    C. L. Nikias, J. M. Mendel, Signal processing with higher-order spectra, IEEE Signal Proc. Mag. 10 (1993) 10

    1993

  14. [22]

    Birkelund, A

    Y . Birkelund, A. Hanssen, Improved bispectrum based tests for Gaussianity and linearity, Signal Process. 89 (2009) 2537. doi:10.1016/j.sigpro.2009.04.013

    2009 doi

  15. [23]

    Hasselmann, W

    K. Hasselmann, W. Munk, G. J. MacDonald, Bispectra of ocean waves, in: Symposium on time series analysis, New York, 1963, pp. 125–139

    1963

  16. [24]

    Si fft, J

    M. Si fft, J. C. Bayer, D. H ¨agele, R. J. Haug, Revealing hidden states in quantum dot array dynamics: Quantum polyspectra versus waiting time analysis, Phys. Rev. B 111 (2025) 205412. doi:10.1103/PhysRevB.111.205412. URL https://link.aps.org/doi/10.1103/PhysRevB.111.205412

    2025 doi

  17. [25]

    Swami, HOSA - Higher Order Spectral Analysis Toolbox, MATLAB Central File Exchange

    A. Swami, HOSA - Higher Order Spectral Analysis Toolbox, MATLAB Central File Exchange. Available at https://www.mathworks.com/matlabcentral/fileexchange/ 3013-hosa-higher-order-spectral-analysis-toolbox, re- trieved January 10, 2025 (2025)

    2025

  18. [26]

    Gardiner, Stochastic Methods, 4th Edition, Springer, Berlin Heidelberg, 2009

    C. Gardiner, Stochastic Methods, 4th Edition, Springer, Berlin Heidelberg, 2009

    2009

  19. [27]

    Starosielec, D

    S. Starosielec, D. H ¨agele, Discrete-time windows with minimal RMS bandwidth for given RMS temporal width, Signal Process. 102 (2014) 240

    2014

  20. [28]

    F. Schefczik, Theorie kontinuierlicher Quantenmessungen h¨oherer Ord- nung und Realisierung eines GPU-GHz-Korrelationsspektrometers, doc- toral thesis, Ruhr-Universit ¨at Bochum, Universit ¨atsbibliothek (2020). doi:10.13154/294-7525

    2020 doi

  21. [29]

    M. A. Wolinsky, Invitation to the bispectrum, Rep. APL-TR-88-7 (1988)

    1988

  22. [30]

    Huber, B

    P. Huber, B. Kleiner, T. Gasser, G. Dumermuth, Statistical meth- ods for investigating phase relations in stationary stochastic pro- cesses, IEEE Trans. Audio Electroacoust. 19 (1) (1971) 78–86. doi:10.1109/TAU.1971.1162163

    1971

  23. [31]

    R. A. Fisher, Moments and product moments of sampling distributions, Proceedings of the London Mathematical Society s2-30 (1) (1928) 199– 238

    1928

  24. [32]

    M. G. Kendall, Advanced theory of statistics. V ol. 1, 1st Edition, Charles Griffin and Co., London, 1943

    1943

  25. [33]

    M. B. Cook, Bi-variate k-statistics and cumulants of their joint sampling distribution, Biometrika 38 (1951) 179

    1951

  26. [34]

    Gardner, C

    W. Gardner, C. Spooner, The cumulant theory of cyclostationary time- series. i. foundation, IEEE Transactions on Signal Processing 42 (12) (1994) 3387. doi:10.1109/78.340775

    1994 doi

  27. [35]

    F. J. Anscombe, Examination of residuals, Proc. Fourth Berkeley Symp. on Math. Statist. and Prob. 1 (1961) 1

    1961

  28. [36]

    L. H. Chan, K. Chen, C. Li, C. W. Wong, C. Y . Yau, On higher-order moment and cumulant estimation, J. Stat. Comput. Simul. 90 (2020) 747. doi:10.1080/00949655.2019.1700987

    2020

  29. [37]

    Chatterjee, O

    R. Chatterjee, O. Grover, Higher order spectral analysis toolkit, https://github.com/synergetics/spectrum, accessed: July 15, 2025 (2025)

    2025

  30. [38]

    Rosenblatt, J

    M. Rosenblatt, J. W. Van Ness, Estimation of the bispectrum, Ann. Math. Statist. 36 (1965) 1120

    1965

  31. [39]

    L. A. Pflug, J. W. Ioup, R. L. Field, Properties of higher-order correlations and spectra for bandlimited, deterministic transients, J. Acoust. Soc. Am. 91 (1992) 975. doi:10.1121/1.402503

    1992 doi

  32. [40]

    Yalamanchili, U

    P. Yalamanchili, U. Arshad, Z. Mohammed, P. Garigipati, P. Entschev, 19 B. Kloppenborg, J. Malcolm, J. Melonakos, ArrayFire - A high perfor- mance software library for parallel computing with an easy-to-use API (2015). URL https://.com/arrayfire/arrayfire

    2015

  33. [41]

    Frigo, S

    M. Frigo, S. Johnson, The design and implementation of FFTW3, Proceed- ings of the IEEE 93 (2) (2005) 216. doi:10.1109/JPROC.2004.840301

    2005

  34. [42]

    V . V . Efimov, M. V . Shokurov, Bispectral analysis of geopotential-height anomalies for synoptic-intraseasonal time-scale, Q. J. R. Meteorol. Soc. 127 (2001) 1707. 20

    2001

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.