pith. sign in

arxiv: 2604.16219 · v1 · submitted 2026-04-17 · 🧮 math.ST · stat.ME· stat.TH

Simultaneous Inference for Covariance and Precision Matrices of Long-Range Dependent Time Series

Pith reviewed 2026-05-10 07:02 UTC · model grok-4.3

classification 🧮 math.ST stat.MEstat.TH
keywords long-range dependencecovariance matrixprecision matrixBerry-Esseen boundGaussian approximationblock bootstraphigh-dimensional time seriessimultaneous inference
0
0 comments X

The pith

A Berry-Esseen bound controls the maximum error in sample covariance for long-range dependent time series in ultra-high dimensions.

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

The paper establishes a finite-sample Gaussian approximation for the infinity-norm error of the sample covariance matrix under long-range temporal dependence. It derives a Berry-Esseen type bound on the Kolmogorov distance between this error distribution and its Gaussian limit by combining martingale approximations, m-dependent approximations, and a triadic block construction. A block bootstrap procedure is shown to remain valid under the same dependence structure. These results hold without any sparsity or other structural assumptions on the matrices, and they permit the dimension to grow sub-exponentially with sample size for covariance inference while allowing an analogous treatment of precision matrices in low-dimensional regimes.

Core claim

We propose a Berry-Esseen type Gaussian approximation result that gives a finite-sample bound for the Kolmogorov distance between the infinity norms of the estimation error of sample covariance matrix and the corresponding Gaussian approximation. The method utilizes martingale and m-dependent approximation and relies on constructing triadic blocks. We also establish a bootstrapping result with block sampling method, which preserves validity despite strong temporal dependence. Our results on covariance allow ultra-high-dimensional settings where the dimension of time series can grow sub-exponentially with sample size. Similar results can be built for precision matrix under low-dimensional tr

What carries the argument

Berry-Esseen Gaussian approximation for the infinity norm of covariance estimation error, obtained via martingale and m-dependent approximations together with triadic block construction.

If this is right

  • Covariance matrix inference remains valid in ultra-high dimensions where the number of variables grows sub-exponentially with sample size.
  • Block bootstrap inference continues to be consistent even when temporal dependence is strong.
  • Precision matrix inference is feasible under low-dimensional regimes with the same lack of structural assumptions.
  • Simultaneous inference on all entries is possible without requiring sparsity or other matrix structure.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same approximation technique could be tested on other matrix-valued functionals of long-memory processes, such as correlation matrices or autocovariance operators at multiple lags.
  • Applications in fields with natural long-range dependence, including financial volatility modeling and climate time series, would gain finite-sample error bounds that do not rely on asymptotic regimes.
  • Extension to settings where both dimension and dependence strength vary jointly might be feasible by adjusting the block construction parameters.

Load-bearing premise

The given long-range dependent process admits martingale and m-dependent approximations that justify the triadic block construction without further explicit decay rates or moment conditions.

What would settle it

A simulation study on fractional Gaussian noise with Hurst parameter above 0.5 in which the empirical Kolmogorov distance between the infinity-norm error distribution and the Gaussian approximation exceeds the paper's derived finite-sample bound for moderate sample sizes.

Figures

Figures reproduced from arXiv: 2604.16219 by Mladen Kolar, Percy S. Zhai, Wei Biao Wu.

Figure 1
Figure 1. Figure 1: QQ-plots for covariance matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p019_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: QQ-plots for covariance matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p019_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: QQ-plots for covariance matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p020_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: QQ-plots for precision matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p021_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: QQ-plots for precision matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p021_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: QQ-plots for precision matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p022_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: QQ-plots for covariance matrix in the high-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p023_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: QQ-plots for covariance matrix in the high-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p023_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: QQ-plots for covariance matrix in the high-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p024_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Histograms of Hurst exponents of all time series within each group. The red dashed line is [PITH_FULL_IMAGE:figures/full_fig_p024_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: ACF plots of select time series. The blue dashed lines are [PITH_FULL_IMAGE:figures/full_fig_p025_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Conditional independence graph (CIG) based on the precision matrix of brain fMRI signals over [PITH_FULL_IMAGE:figures/full_fig_p025_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Conditional independence graph (CIG) based on the precision matrix of brain fMRI signals over [PITH_FULL_IMAGE:figures/full_fig_p026_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Empirical CDF in the low-dimensional regime, which contains the covariance error [PITH_FULL_IMAGE:figures/full_fig_p055_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Empirical CDF in the low-dimensional regime, which contains the precision error [PITH_FULL_IMAGE:figures/full_fig_p055_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Empirical CDF in the high-dimensional regime, containing the covariance error [PITH_FULL_IMAGE:figures/full_fig_p056_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: QQ-plots for covariance matrix in the low-dimensional regime, [PITH_FULL_IMAGE:figures/full_fig_p059_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: QQ-plots for covariance matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p059_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: QQ-plots for covariance matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p060_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Empirical CDF in the low-dimensional regime under the banded setup, which contains the covariance error [PITH_FULL_IMAGE:figures/full_fig_p060_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: QQ-plots for precision matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p061_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: QQ-plots for precision matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p061_22.png] view at source ↗
Figure 23
Figure 23. Figure 23: QQ-plots for precision matrix in the low-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p062_23.png] view at source ↗
Figure 24
Figure 24. Figure 24: Empirical CDF in the low-dimensional regime under the banded setup, which contains the precision error [PITH_FULL_IMAGE:figures/full_fig_p062_24.png] view at source ↗
Figure 25
Figure 25. Figure 25: QQ-plots for covariance matrix in the high-dimensional regime, [PITH_FULL_IMAGE:figures/full_fig_p065_25.png] view at source ↗
Figure 26
Figure 26. Figure 26: QQ-plots for covariance matrix in the high-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p065_26.png] view at source ↗
Figure 27
Figure 27. Figure 27: QQ-plots for covariance matrix in the high-dimensional regime, ˆ [PITH_FULL_IMAGE:figures/full_fig_p066_27.png] view at source ↗
Figure 28
Figure 28. Figure 28: Empirical CDF in the high-dimensional regime under the banded setup, which contains the covariance error [PITH_FULL_IMAGE:figures/full_fig_p066_28.png] view at source ↗
read the original abstract

For time series with long-range temporal dependence, inference for covariance and precision matrices is non-trivial. We propose a Berry-Esseen type Gaussian approximation result that gives a finite-sample bound for the Kolmogorov distance between the infinity norms of the estimation error of sample covariance matrix and the corresponding Gaussian approximation. The method utilizes martingale and m-dependent approximation and relies on constructing triadic blocks. We also establish a bootstrapping result with block sampling method, which preserves validity despite strong temporal dependence. Our results on covariance allow ultra-high-dimensional settings where the dimension of time series can grow sub-exponentially with sample size. Similar results can be built for precision matrix under low-dimensional settings. No assumption is required on the structure of covariance and precision matrices.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a Berry-Esseen-type Gaussian approximation result for the Kolmogorov distance between the infinity-norm error of the sample covariance estimator and its Gaussian limit for long-range dependent time series. The approach relies on martingale and m-dependent approximations constructed via triadic blocks, yielding finite-sample bounds that permit ultra-high-dimensional regimes in which the dimension p grows sub-exponentially with sample size n. A block-sampling bootstrap is shown to be consistent under the same dependence, and analogous (lower-dimensional) results are indicated for precision-matrix estimation. No structural assumptions are imposed on the covariance or precision matrices themselves.

Significance. If the dependence and moment conditions are stated explicitly and the approximation errors are controlled uniformly over the p^2 entries, the results would meaningfully extend simultaneous inference tools to LRD settings that arise in finance, neuroscience, and climate data, where sub-exponential dimension growth is practically relevant. The triadic-block construction and martingale approximation are standard devices, but their successful application here to obtain dimension-free rates would be a technical contribution.

major comments (2)
  1. [Abstract and §1] The abstract and introduction assert results for general long-range dependent processes without stating explicit decay rates (e.g., |Cov(X_t,X_{t+k})| = O(k^{-α}) for α > 0) or moment conditions (e.g., E|X_t|^{2+δ}). These rates are load-bearing for the claim that the m-dependent approximation error is o((log p)^{-1/2}) uniformly when p grows sub-exponentially; without them the Kolmogorov bound in the main theorem cannot be verified for the full range of LRD processes.
  2. [§3 (Martingale approximation)] The triadic-block construction is invoked to control dependence within and across blocks, yet no quantitative bound is given on the block-size choice m = m(n) that simultaneously keeps the number of blocks ~ n/m and makes the approximation error negligible relative to the Gaussian tail when p = exp(o(n^β)). This choice is central to the ultra-high-dimensional regime.
minor comments (2)
  1. [§3] Notation for the triadic blocks and the associated filtration should be introduced with a diagram or explicit indexing to improve readability.
  2. [§5] The bootstrap section would benefit from a clear statement of the block length choice and its relation to the dependence parameter α.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and will make the suggested revisions to improve clarity and rigor.

read point-by-point responses
  1. Referee: [Abstract and §1] The abstract and introduction assert results for general long-range dependent processes without stating explicit decay rates (e.g., |Cov(X_t,X_{t+k})| = O(k^{-α}) for α > 0) or moment conditions (e.g., E|X_t|^{2+δ}). These rates are load-bearing for the claim that the m-dependent approximation error is o((log p)^{-1/2}) uniformly when p grows sub-exponentially; without them the Kolmogorov bound in the main theorem cannot be verified for the full range of LRD processes.

    Authors: We agree that the abstract and introduction should explicitly state the decay rates and moment conditions. The theorems in the manuscript are proved under the assumptions that the covariance satisfies |Cov(X_t, X_{t+k})| = O(k^{-α}) for some α > 0 and that E|X_t|^{2+δ} < ∞ for δ > 0, which ensure the m-dependent approximation error is o((log p)^{-1/2}) uniformly. These conditions were implicit in the technical sections but not highlighted upfront. In the revision we will update the abstract and Section 1 to state them explicitly, thereby clarifying the precise scope of the sub-exponential dimension growth. revision: yes

  2. Referee: [§3 (Martingale approximation)] The triadic-block construction is invoked to control dependence within and across blocks, yet no quantitative bound is given on the block-size choice m = m(n) that simultaneously keeps the number of blocks ~ n/m and makes the approximation error negligible relative to the Gaussian tail when p = exp(o(n^β)). This choice is central to the ultra-high-dimensional regime.

    Authors: We acknowledge that an explicit quantitative range for the block size m(n) is needed to justify the ultra-high-dimensional regime. The triadic-block construction is used to obtain a martingale approximation whose error is controlled relative to the Gaussian tail, but the manuscript does not spell out the admissible exponents for m = n^γ in terms of the dependence parameter α and the dimension-growth exponent β. In the revision we will add a supporting lemma or remark in Section 3 that specifies the allowable range for γ (for example, γ small enough to keep the approximation error o((log p)^{-1/2}) yet large enough to preserve the block count ~ n/m) under the stated moment and decay conditions when p = exp(o(n^β)). revision: yes

Circularity Check

0 steps flagged

No circularity: derivation uses standard martingale/m-dependent approximations and triadic blocking without self-referential definitions or fitted inputs renamed as predictions

full rationale

The abstract and provided text describe a Berry-Esseen Gaussian approximation for the infinity-norm error of the sample covariance under long-range dependence, obtained via martingale and m-dependent approximations together with triadic block construction. No equation or step is shown to reduce by construction to its own inputs (e.g., no parameter fitted on a subset and then called a prediction, no self-definition of the target quantity, no uniqueness theorem imported solely from the authors' prior work). The central claims rest on established blocking techniques for dependent processes rather than on any load-bearing self-citation chain or ansatz smuggled via citation. The absence of explicit decay rates or moment conditions in the abstract is a potential gap in stated assumptions, but does not constitute circularity in the derivation itself. The result is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the applicability of martingale and m-dependent approximations to long-range dependent processes and on the construction of triadic blocks; these are domain-specific technical assumptions rather than free parameters or new entities.

axioms (2)
  • domain assumption Martingale and m-dependent approximations can be applied to the long-range dependent time series under study
    Invoked to obtain the Gaussian approximation and Kolmogorov-distance bound.
  • ad hoc to paper Triadic blocks can be constructed so that the dependence within and across blocks is controllable
    Central technical device used to handle strong temporal dependence.

pith-pipeline@v0.9.0 · 5428 in / 1459 out tokens · 57117 ms · 2026-05-10T07:02:07.540281+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

89 extracted references · 89 canonical work pages

  1. [1]

    Gaussian approximations and multiplier bootstrap for maxima of sums of high-dimensional random vectors,

    V . Chernozhukov, D. Chetverikov, and K. Kato, “Gaussian approximations and multiplier bootstrap for maxima of sums of high-dimensional random vectors,”The Annals of Statistics, vol. 41, no. 6, pp. 2786–2819, 2013

  2. [2]

    Comparison and anti-concentration bounds for maxima of gaussian random vectors,

    ——, “Comparison and anti-concentration bounds for maxima of gaussian random vectors,”Probability Theory and Related Fields, vol. 162, no. 1, pp. 47–70, 2015

  3. [3]

    Central limit theorems and bootstrap in high dimensions,

    ——, “Central limit theorems and bootstrap in high dimensions,”The Annals of Probability, vol. 45, no. 4, pp. 2309–2352, 2017

  4. [4]

    Gaussian approximation for high dimensional time series,

    D. Zhang and W. B. Wu, “Gaussian approximation for high dimensional time series,”The Annals of Statistics, vol. 45, no. 5, pp. 1895–1919, 2017

  5. [5]

    The asymptotic distributions of the largest entries of sample correlation matrices,

    T. Jiang, “The asymptotic distributions of the largest entries of sample correlation matrices,”The Annals of Applied Probability, vol. 14, no. 2, pp. 865–880, 2004

  6. [6]

    Some strong limit theorems for the largest entries of sample correlation matrices,

    D. Li and A. Rosalsky, “Some strong limit theorems for the largest entries of sample correlation matrices,”The Annals of Applied Probability, vol. 16, no. 1, pp. 423–447, 2006

  7. [7]

    Asymptotic distribution of the largest off-diagonal entry of correlation matrices,

    W. Zhou, “Asymptotic distribution of the largest off-diagonal entry of correlation matrices,”Transactions of the American Mathematical Society, vol. 359, no. 11, pp. 5345–5363, 2007

  8. [8]

    The asymptotic distribution and berry–esseen bound of a new test for independence in high dimension with an application to stochastic optimization,

    W.-D. Liu, Z. Lin, and Q.-M. Shao, “The asymptotic distribution and berry–esseen bound of a new test for independence in high dimension with an application to stochastic optimization,”The Annals of Applied Probability, vol. 18, no. 6, pp. 2337–2366, 2008

  9. [9]

    Necessary and sufficient conditions for the asymptotic distribution of the largest entry of a sample correlation matrix,

    D. Li, W.-D. Liu, and A. Rosalsky, “Necessary and sufficient conditions for the asymptotic distribution of the largest entry of a sample correlation matrix,”Probability Theory and Related Fields, vol. 148, no. 1-2, pp. 5–35, 2010

  10. [10]

    Limiting laws of coherence of random matrices with applications to testing covariance structure and construction of compressed sensing matrices,

    T. T. Cai and T. Jiang, “Limiting laws of coherence of random matrices with applications to testing covariance structure and construction of compressed sensing matrices,”The Annals of Statistics, vol. 39, no. 3, pp. 1496–1525, 2011

  11. [11]

    Asymptotic theory for maximum deviations of sample covariance matrix estimates,

    H. Xiao and W. B. Wu, “Asymptotic theory for maximum deviations of sample covariance matrix estimates,”Stochastic Processes and their Applications, vol. 123, no. 7, pp. 2899–2920, 2013

  12. [12]

    Largest entries of sample correlation matrices from equi-correlated normal populations,

    J. Fan and T. Jiang, “Largest entries of sample correlation matrices from equi-correlated normal populations,”The Annals of Probability, vol. 47, no. 5, pp. 3321–3374, 2019

  13. [13]

    Gaussian and bootstrap approximations for high-dimensionalu-statistics and their applications,

    X. Chen, “Gaussian and bootstrap approximations for high-dimensionalu-statistics and their applications,”The Annals of Statistics, vol. 46, no. 2, pp. 642–678, 2018

  14. [14]

    Asymptotic inference of autocovariances of stationary processes,

    H. Xiao and W. B. Wu, “Asymptotic inference of autocovariances of stationary processes,”arXiv preprint arXiv:1105.3423, 2011

  15. [15]

    Estimation of large covariance and precision matrices from temporally dependent observations,

    H. Shu and B. Nan, “Estimation of large covariance and precision matrices from temporally dependent observations,”The Annals of Statistics, vol. 47, no. 3, pp. 1321 – 1350, 2019

  16. [16]

    Local Whittle estimation of high-dimensional long-run variance and precision matrices,

    C. Baek, M.-C. D ¨uker, and V . Pipiras, “Local Whittle estimation of high-dimensional long-run variance and precision matrices,”The Annals of Statistics, vol. 51, no. 6, pp. 2386–2414, 2023

  17. [17]

    On asymptotically optimal confidence regions and tests for high-dimensional models,

    S. van de Geer, P. B ¨uhlmann, Y . Ritov, and R. Dezeure, “On asymptotically optimal confidence regions and tests for high-dimensional models,”The Annals of Statistics, vol. 42, no. 3, pp. 1166 – 1202, 2014

  18. [18]

    Confidence intervals and hypothesis testing for high-dimensional regression,

    A. Javanmard and A. Montanari, “Confidence intervals and hypothesis testing for high-dimensional regression,”Journal of Machine Learning Research, vol. 15, no. 82, pp. 2869–2909, 2014

  19. [19]

    Statistical inference of high-dimensional vector autoregressive time series with non-i.i.d. innovations,

    Y . Zhang, “Statistical inference of high-dimensional vector autoregressive time series with non-i.i.d. innovations,”arXiv preprint arXiv:2310.07364, 2023

  20. [20]

    On statistical inference for high-dimensional binary time series,

    D. Dai and Y . Zhang, “On statistical inference for high-dimensional binary time series,”arXiv preprint arXiv:2512.00338, 2025

  21. [21]

    Bootstrapping High Dimensional Time Series

    X. Zhang and G. Cheng, “Bootstrapping high dimensional time series,”arXiv preprint arXiv:1406.1037, 2014

  22. [22]

    Block sampling under strong dependence,

    T. Zhang, H.-C. Ho, M. Wendler, and W. B. Wu, “Block sampling under strong dependence,”Stochastic Processes and their Applications, vol. 123, no. 6, pp. 2323–2339, 2013

  23. [23]

    Improved central limit theorem and bootstrap approximations in high dimensions,

    V . Chernozhuokov, D. Chetverikov, K. Kato, and Y . Koike, “Improved central limit theorem and bootstrap approximations in high dimensions,”The Annals of Statistics, vol. 50, no. 5, pp. 2562 – 2586, 2022

  24. [24]

    High-dimensional central limit theorems by stein’s method,

    X. Fang and Y . Koike, “High-dimensional central limit theorems by stein’s method,”The Annals of Applied Probability, vol. 31, no. 4, pp. 1660–1686, 2021

  25. [25]

    Nearly optimal central limit theorem and bootstrap approximations in high dimensions,

    V . Chernozhukov, D. Chetverikov, and Y . Koike, “Nearly optimal central limit theorem and bootstrap approximations in high dimensions,”The Annals of Applied Probability, vol. 33, no. 3, pp. 2374 – 2425, 2023

  26. [26]

    Approximating high-dimensional infinite-orderu-statistics: statistical and computational guarantees,

    Y . Song, X. Chen, and K. Kato, “Approximating high-dimensional infinite-orderu-statistics: statistical and computational guarantees,”Electronic Journal of Statistics, vol. 13, no. 2, pp. 4794–4848, 2019

  27. [27]

    Randomized incompleteu-statistics in high dimensions,

    X. Chen and K. Kato, “Randomized incompleteu-statistics in high dimensions,”The Annals of Statistics, vol. 47, no. 6, pp. 3127–3156, 2019

  28. [28]

    Jackknife multiplier bootstrap: finite sample approximations to theu-process supremum with applications,

    ——, “Jackknife multiplier bootstrap: finite sample approximations to theu-process supremum with applications,”Probability Theory and Related Fields, vol. 176, no. 3, pp. 1097–1163, 2020

  29. [29]

    Gaussian approximation of maxima of wiener functionals and its application to high-frequency data,

    Y . Koike, “Gaussian approximation of maxima of wiener functionals and its application to high-frequency data,”The Annals of Statistics, vol. 47, no. 3, pp. 1663–1687, 2019

  30. [30]

    Inference on causal and structural parameters using many moment inequalities,

    V . Chernozhukov, D. Chetverikov, and K. Kato, “Inference on causal and structural parameters using many moment inequalities,”The Review of Economic Studies, vol. 86, no. 5, pp. 1867–1900, 2019

  31. [31]

    Lasso-driven inference in time and space,

    V . Chernozhukov, W. K. H ¨ardle, C. Huang, and W. Wang, “Lasso-driven inference in time and space,”The Annals of Statistics, vol. 49, no. 3, pp. 1702–1735, 2021

  32. [32]

    Koml ´os–Major–Tusn´ady approximation under dependence,

    I. Berkes, W. Liu, and W. B. Wu, “Koml ´os–Major–Tusn´ady approximation under dependence,”The Annals of Probability, vol. 42, no. 2, pp. 794–817, 2014

  33. [33]

    High-dimensional graphs and variable selection with the Lasso,

    N. Meinshausen and P. B ¨uhlmann, “High-dimensional graphs and variable selection with the Lasso,”The Annals of Statistics, vol. 34, no. 3, pp. 1436 – 1462, 2006

  34. [34]

    A constrainedℓ 1-minimization approach to sparse precision matrix estimation,

    T. Cai, W. Liu, and X. Luo, “A constrainedℓ 1-minimization approach to sparse precision matrix estimation,”Journal of the American Statistical Association, vol. 106, no. 494, pp. 594–607, 2011

  35. [35]

    Test for high-dimensional covariance matrices,

    Y . Han and W. B. Wu, “Test for high-dimensional covariance matrices,”The Annals of Statistics, vol. 48, no. 6, pp. 3565 – 3588, 2020. ACCEPTED FOR PUBLICATION IN IEEE TRANSACTIONS ON INFORMATION THEORY 28

  36. [36]

    High-dimensional econometrics and regularized GMM,

    A. Belloni, V . Chernozhukov, D. Chetverikov, C. Hansen, and K. Kato, “High-dimensional econometrics and regularized GMM,” Centre for Microdata Methods and Practice (cemmap), Tech. Rep. CWP35/18, 2018, contributed chapter for the Handbook of Econometrics

  37. [37]

    M. L. Stein,Interpolation of spatial data: some theory for kriging, ser. Springer Series in Statistics. New York: Springer, 1999

  38. [38]

    Simulating sample paths of linear fractional stable motion,

    W. B. Wu, G. Michailidis, and D. Zhang, “Simulating sample paths of linear fractional stable motion,”IEEE Transactions on Information Theory, vol. 50, no. 6, pp. 1086–1096, 2004

  39. [39]

    Resting brains never rest: computational insights into potential cognitive architectures,

    G. Deco, V . K. Jirsa, and A. R. McIntosh, “Resting brains never rest: computational insights into potential cognitive architectures,”Trends in Neurosciences, vol. 36, no. 5, pp. 268–274, 2013

  40. [40]

    Large-scale cortical correlation structure of spontaneous oscillatory activity,

    J. F. Hipp, D. J. Hawellek, M. Corbetta, M. Siegel, and A. K. Engel, “Large-scale cortical correlation structure of spontaneous oscillatory activity,” Nature Neuroscience, vol. 15, no. 6, pp. 884–890, 2012

  41. [41]

    Dynamic functional connectivity: promise, issues, and interpretations,

    R. M. Hutchison, T. Womelsdorf, E. A. Allen, P. A. Bandettini, V . D. Calhoun, M. Corbetta, S. Della Penna, J. H. Duyn, G. H. Glover, J. Gonzalez-Castillo et al., “Dynamic functional connectivity: promise, issues, and interpretations,”NeuroImage, vol. 80, pp. 360–378, 2013

  42. [42]

    Adding dynamics to the Human Connectome Project with MEG,

    L. J. Larson-Prior, R. Oostenveld, S. Della Penna, G. Michalareas, F. Prior, A. Babajani-Feremi, J.-M. Schoffelen, L. Marzetti, F. de Pasquale, F. Di Pompeo et al., “Adding dynamics to the Human Connectome Project with MEG,”NeuroImage, vol. 80, pp. 190–201, 2013

  43. [43]

    High-dimensional functional graphical model structure learning via neighborhood selection approach,

    B. Zhao, P. S. Zhai, Y . S. Wang, and M. Kolar, “High-dimensional functional graphical model structure learning via neighborhood selection approach,” Electronic Journal of Statistics, vol. 18, no. 1, pp. 1042–1129, 2024

  44. [44]

    The variation of some other speculative prices,

    B. Mandelbrot, “The variation of some other speculative prices,”The Journal of Business, vol. 40, no. 4, pp. 393–413, 1967

  45. [45]

    A long memory property of stock market returns and a new model,

    Z. Ding, C. W. Granger, and R. F. Engle, “A long memory property of stock market returns and a new model,”Journal of Empirical Finance, vol. 1, no. 1, pp. 83–106, 1993

  46. [46]

    Beran,Statistics for long-memory processes, ser

    J. Beran,Statistics for long-memory processes, ser. Monographs on Statistics and Applied Probability. New York: Chapman & Hall, 1994, vol. 61

  47. [47]

    M. S. Bartlett,An introduction to stochastic processes: with special reference to methods and applications, 3rd ed. Cambridge and New York: Cambridge University Press, 1978

  48. [48]

    Concentration inequalities for polynomials inα-sub-exponential random variables,

    F. G ¨otze, H. Sambale, and A. Sinulis, “Concentration inequalities for polynomials inα-sub-exponential random variables,”Electronic Journal of Probability, vol. 26, no. 48, pp. 1–22, 2021

  49. [49]

    Block length selection in the bootstrap for time series,

    P. B ¨uhlmann and H. R. K ¨unsch, “Block length selection in the bootstrap for time series,”Computational Statistics & Data Analysis, vol. 31, no. 3, pp. 295–310, 1999

  50. [50]

    Automatic block-length selection for the dependent bootstrap,

    D. N. Politis and H. White, “Automatic block-length selection for the dependent bootstrap,”Econometric Reviews, vol. 23, no. 1, pp. 53–70, 2004

  51. [51]

    Strong invariance principles for dependent random variables,

    W. B. Wu, “Strong invariance principles for dependent random variables,”The Annals of Probability, vol. 35, no. 6, pp. 2294 – 2320, 2007

  52. [52]

    On the sampling window method for long-range dependent data,

    P. Hall, B.-Y . Jing, and S. N. Lahiri, “On the sampling window method for long-range dependent data,”Statistica Sinica, vol. 8, no. 4, pp. 1189–1204, 1998

  53. [53]

    The tight constant in the Dvoretzky–Kiefer–Wolfowitz inequality,

    P. Massart, “The tight constant in the Dvoretzky–Kiefer–Wolfowitz inequality,”The Annals of Probability, vol. 18, no. 3, pp. 1269–1283, 1990

  54. [54]

    Concentration inequalities and moment bounds for sample covariance operators,

    V . Koltchinskii and K. Lounici, “Concentration inequalities and moment bounds for sample covariance operators,”Bernoulli, vol. 23, no. 1, pp. 110–133, 2017

  55. [55]

    Introduction to the non-asymptotic analysis of random matrices,

    R. Vershynin, “Introduction to the non-asymptotic analysis of random matrices,” inCompressed Sensing: Theory and Applications, Y . C. Eldar and G. Kutyniok, Eds. Cambridge: Cambridge University Press, 2012, pp. 210–268

  56. [56]

    Simulation of stationary gaussian processes in[0,1] d,

    A. T. A. Wood and G. Chan, “Simulation of stationary gaussian processes in[0,1] d,”Journal of Computational and Graphical Statistics, vol. 3, no. 4, pp. 409–432, 1994

  57. [57]

    The Autism Brain Imaging Data Exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism,

    A. Di Martino, C.-G. Yan, Q. Li, E. Denio, F. X. Castellanos, K. Alaerts, J. S. Anderson, M. Assaf, S. Y . Bookheimer, M. Daprettoet al., “The Autism Brain Imaging Data Exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism,”Molecular Psychiatry, vol. 19, no. 6, pp. 659–667, 2014

  58. [58]

    Neuroimaging in autism spectrum disorder: brain structure and function across the lifespan,

    C. Ecker, S. Y . Bookheimer, and D. G. Murphy, “Neuroimaging in autism spectrum disorder: brain structure and function across the lifespan,”The Lancet Neurology, vol. 14, no. 11, pp. 1121–1134, 2015

  59. [59]

    Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain,

    N. Tzourio-Mazoyer, B. Landeau, D. Papathanassiou, F. Crivello, O. Etard, N. Delcroix, B. Mazoyer, and M. Joliot, “Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain,”NeuroImage, vol. 15, no. 1, pp. 273–289, 2002

  60. [60]

    The precuneus/posterior cingulate cortex plays a pivotal role in the default mode network: Evidence from a partial correlation network analysis,

    P. Fransson and G. Marrelec, “The precuneus/posterior cingulate cortex plays a pivotal role in the default mode network: Evidence from a partial correlation network analysis,”NeuroImage, vol. 42, no. 3, pp. 1178–1184, 2008

  61. [61]

    How ongoing neuronal oscillations account for evoked fMRI variability,

    R. Becker, M. Reinacher, F. Freyer, A. Villringer, and P. Ritter, “How ongoing neuronal oscillations account for evoked fMRI variability,”The Journal of Neuroscience, vol. 31, no. 30, pp. 11 016–11 027, 2011

  62. [62]

    Event-related fMRI: characterizing differential responses,

    K. J. Friston, P. Fletcher, O. Josephs, A. Holmes, M. Rugg, and R. Turner, “Event-related fMRI: characterizing differential responses,”NeuroImage, vol. 7, no. 1, pp. 30–40, 1998

  63. [63]

    Growing together and growing apart: regional and sex differences in the lifespan developmental trajectories of functional homotopy,

    X.-N. Zuo, C. Kelly, A. Di Martino, M. Mennes, D. S. Margulies, S. Bangaru, R. Grzadzinski, A. C. Evans, Y .-F. Zang, F. X. Castellanoset al., “Growing together and growing apart: regional and sex differences in the lifespan developmental trajectories of functional homotopy,”The Journal of Neuroscience, vol. 30, no. 45, pp. 15 034–15 043, 2010

  64. [64]

    The organization of local and distant functional connectivity in the human brain,

    J. Sepulcre, H. Liu, T. Talukdar, I. Martincorena, B. T. Yeo, and R. L. Buckner, “The organization of local and distant functional connectivity in the human brain,”PLOS Computational Biology, vol. 6, no. 6, p. e1000808, 2010

  65. [65]

    Atypical cross talk between mentalizing and mirror neuron networks in autism spectrum disorder,

    I. Fishman, C. L. Keown, A. J. Lincoln, J. A. Pineda, and R.-A. M ¨uller, “Atypical cross talk between mentalizing and mirror neuron networks in autism spectrum disorder,”JAMA Psychiatry, vol. 71, no. 7, pp. 751–760, 2014

  66. [66]

    Executive dysfunction in autism,

    E. L. Hill, “Executive dysfunction in autism,”Trends in Cognitive Sciences, vol. 8, no. 1, pp. 26–32, 2004

  67. [67]

    Abnormal functional connectivity of default mode sub-networks in autism spectrum disorder patients,

    M. Assaf, K. Jagannathan, V . D. Calhoun, L. Miller, M. C. Stevens, R. Sahl, J. G. O’Boyle, R. T. Schultz, and G. D. Pearlson, “Abnormal functional connectivity of default mode sub-networks in autism spectrum disorder patients,”NeuroImage, vol. 53, no. 1, pp. 247–256, 2010

  68. [68]

    Disruption of functional organization within the primary motor cortex in children with autism,

    M. B. Nebel, S. E. Joel, J. Muschelli, A. D. Barber, B. S. Caffo, J. J. Pekar, and S. H. Mostofsky, “Disruption of functional organization within the primary motor cortex in children with autism,”Human Brain Mapping, vol. 35, no. 2, pp. 567–580, 2014

  69. [69]

    Deviant functional magnetic resonance imaging patterns of brain activity to speech in 2–3-year-old children with autism spectrum disorder,

    E. Redcay and E. Courchesne, “Deviant functional magnetic resonance imaging patterns of brain activity to speech in 2–3-year-old children with autism spectrum disorder,”Biological Psychiatry, vol. 64, no. 7, pp. 589–598, 2008

  70. [70]

    Altered network topologies and hub organization in adults with autism: a resting-state fMRI study,

    T. Itahashi, T. Yamada, H. Watanabe, M. Nakamura, D. Jimbo, S. Shioda, K. Toriizuka, N. Kato, and R. Hashimoto, “Altered network topologies and hub organization in adults with autism: a resting-state fMRI study,”PLOS ONE, vol. 9, no. 4, p. e94115, 2014

  71. [71]

    A compensatory role for declarative memory in neurodevelopmental disorders,

    M. T. Ullman and M. Y . Pullman, “A compensatory role for declarative memory in neurodevelopmental disorders,”Neuroscience & Biobehavioral Reviews, vol. 51, pp. 205–222, 2015

  72. [72]

    Atypically diffuse functional connectivity between caudate nuclei and cerebral cortex in autism,

    K. C. Turner, L. Frost, D. Linsenbardt, J. R. McIlroy, and R.-A. M ¨uller, “Atypically diffuse functional connectivity between caudate nuclei and cerebral cortex in autism,”Behavioral and Brain Functions, vol. 2, p. 34, 2006

  73. [73]

    Cortical activation and synchronization during sentence comprehension in high-functioning autism: evidence of underconnectivity,

    M. A. Just, V . L. Cherkassky, T. A. Keller, and N. J. Minshew, “Cortical activation and synchronization during sentence comprehension in high-functioning autism: evidence of underconnectivity,”Brain, vol. 127, no. 8, pp. 1811–1821, 2004

  74. [74]

    Functional graphical models,

    X. Qiao, S. Guo, and G. M. James, “Functional graphical models,”Journal of the American Statistical Association, vol. 114, no. 525, pp. 211–222, 2019

  75. [75]

    A bound on tail probabilities for quadratic forms in independent random variables,

    D. L. Hanson and F. T. Wright, “A bound on tail probabilities for quadratic forms in independent random variables,”The Annals of Mathematical Statistics, vol. 42, no. 3, pp. 1079–1083, 1971

  76. [76]

    Hanson–Wright inequality and sub-Gaussian concentration,

    M. Rudelson and R. Vershynin, “Hanson–Wright inequality and sub-Gaussian concentration,”Electronic Communications in Probability, vol. 18, no. 82, pp. 1–9, 2013. ACCEPTED FOR PUBLICATION IN IEEE TRANSACTIONS ON INFORMATION THEORY 29

  77. [77]

    Quelques in ´egalit´es pour les op ´erations lin ´eaires,

    J. Marcinkiewicz and A. Zygmund, “Quelques in ´egalit´es pour les op ´erations lin ´eaires,”Fundamenta Mathematicae, vol. 32, no. 1, pp. 115–121, 1939

  78. [78]

    On direct product matrices,

    W. E. Roth, “On direct product matrices,”Bulletin of the American Mathematical Society, vol. 40, no. 6, pp. 461–468, 1934. Percy S. Zhaiis currently a Ph.D. student in econometrics and statistic at the University of Chicago Booth School of Business. His research interests include time series analysis, graphical models, uncertainty quantification, statisti...

  79. [79]

    Since ACCEPTED FOR PUBLICATION IN IEEE TRANSACTIONS ON INFORMATION THEORY 33 |G(s,t)|op ≤ |G (s,t)|F, we focus on the latter

    We need to bound the operator and the Frobenius norms ofG (s,t). Since ACCEPTED FOR PUBLICATION IN IEEE TRANSACTIONS ON INFORMATION THEORY 33 |G(s,t)|op ≤ |G (s,t)|F, we focus on the latter. It suffices to compute P j≤n P k≤n |G(s,t) jk |2 F. By Condition 1 and Cauchy- Schwarz inequality, up to a constant related toC 0, nX j=1 nX k=1 |G(s,t) jk |2 F ≲ nX ...

  80. [80]

    We also have (Fj,1)(s,t) = (3j+2)mX i=(3j+1)m+1 " ∞X l=0 (Al)s·ϵi(Al)t·ϵi −Σ st + m−1X k=1 ∞X l=0 (Al+k)s·ϵi−k(Al)t·ϵi + m−1X k=1 ∞X l=0 (Al)s·ϵi(Al+k)t·ϵi−k # , (A13) which is mean-zero, and unrelated toϵ ∗

Showing first 80 references.