From sparse to dense functional time series: phase transitions of detecting structural breaks and beyond

Leheng Cai; Qirui Hu

REVIEW 3 major objections 6 minor 1 cited by

From sparse to dense functional time series: phase transitions of detecting structural breaks and beyond

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

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

Pith's one-line read A B-spline-smoothed CUMSUM process yields functional break-point tests whose null distributions hold uniformly across sparse, semi-dense, and dense sampling, alongside a phase-transition boundary and post-break inference.

desk verdict Unified break-point test for functional time series from sparse to dense; real contribution, but the automatic knot choice runs ahead of the theory and the fixed-design applications are outside the theorem. read the letter →

arxiv 2412.20858 v1 pith:BJHEEFHM submitted 2024-12-30 stat.ME math.STstat.TH

Leheng Cai , Qirui Hu This is my paper · ORCID

classification stat.MEmath.STstat.TH MSC 62G1062G2062G0862M1062R10

keywords functionaltimeseriesstructuralbreakschangepointdetectionB-splinesmoothingCUMSUMprocesssparsetodensedataphasetransitionGaussianapproximation

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

This paper tries to establish that the same change-point tests for the mean of a functional time series work no matter how thinly or densely each curve is sampled. Previous tests either require fully observed trajectories, regular grids, or sampling frequencies that grow at a prescribed rate relative to the number of curves. The authors construct L∞ and L2 statistics from a B-spline-smoothed CUMSUM process and prove, under a common set of assumptions, that both have asymptotically correct size across all sampling regimes and a Gaussian-field null limit. They also derive the boundary where the effective limit switches from a noise-dominated sparse regime to a dependence-dominated dense regime, give convergence rates for the estimated break date, and build simultaneous confidence bands for the jump function. If true, one procedure could replace a patchwork of regime-specific tests.

What carries the argument

The load-bearing object is the standardized smoothed CUMSUM process $Q_n(t,x)=([nt]\hat m(t,x)/\sqrt n - \sqrt n\,t\,\hat m(x))/\sqrt{B^{\top}(x)\Sigma B(x)}$, where $\hat m(t,\cdot)$ and $\hat m(\cdot)$ are B-spline least-squares estimates of the partial and global mean with per-curve weight $N_i^{-1}$. Under the null this process is approximated by a zero-mean Gaussian field $Z_n(t,x)$ with covariance $(\min\{t,t'\}-tt')B^{\top}(x)\Sigma B(x')$, and the normalizer $\sqrt{B^{\top}\Sigma B}$ keeps the statistic from diverging when few points per curve make the noise term dominate. The proof machinery combines sequential Gaussian approximation for dependent vector-valued series, anti-concentration bounds for maxima of Gaussian vectors, and spline basis properties; the decomposition $\Sigma=\Sigma_1+\Sigma_2$ is what turns the sparse-to-dense phase transition into a comparison of orders.

What would settle it

Simulate the null model (1.1) in a sparse regime where (A6) is violated—for example $E(N^{-1})J_n \gg 1$ while $J_n^{-q^*-1/2} n^{1/2}E^{-1/2}(N^{-1})a_n$ diverges—and compare the empirical distribution of $T_n$ with the quantile of $\sup_{\epsilon\le t\le1-\epsilon}\sup_{x\in[0,1]}|Z_n(t,x)|/\sqrt{B^{\top}(x)\Sigma B(x)}$; a persistent size distortion would show these rate conditions are doing real work.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that, under assumptions (A1)-(A6), $\sup_z |P(T_n \le z) - P(\sup_{\epsilon\le t\le1-\epsilon}\sup_{x\in[0,1]} |Z_n(t,x)|/\sqrt{B^{\top}(x)\Sigma B(x)} \le z)| \to 0$ with no restriction on how $n$ and the per-curve counts $N_i$ are related; the same holds for $S_n$ with the integral of the squared field under the slightly lighter (A6'). Corollary 3.1 then splits the regimes by comparing $E(N^{-1})$ with $n^{1/(2q^*)}$: sparse sampling is dominated by the noise term $\Sigma_2$, dense by the dependence term $\Sigma_1$, and the overlap is a mixed semi-dense regime. Under local alternatives the minimal detectable jump is of order $\max\{\sqrt{\log n/n}, \sqrt{E(N^{-1})J_n\log n/n}\}$ in sup-norm and its $L^2$ analogue, so the tests are consistent whenever the jump exceeds the regime-dependent threshold. The same construction yields break-point estimators and an asymptotically valid simultaneous confidence band for the jump magnitude $\Delta(\cdot)$.

Load-bearing premise

The theory holds only if the spline order, moment conditions, and per-curve sampling intensities satisfy a list of coupled rate restrictions, and the data-driven BIC knot selection used in the applications is not shown to meet them.

Editorial extensions

If this is right

The same two test statistics can be applied without first classifying the sampling scheme; the normalizer $\sqrt{B^{\top}(x)\Sigma B(x)}$ automatically tunes the limit from noise-dominated to dependence-dominated.
The phase-transition boundary separates sparse ($E^{-1}(N^{-1}) \ll n^{1/(2q^*)}$) from dense ($E^{-1}(N^{-1}) \gg n^{1/(2q^*)}$) sampling, up to logarithmic factors, with a mixed semi-dense regime in between.
Under local alternatives, the tests have nontrivial power once the jump exceeds $n^{-1/2}$ in dense or semi-dense settings and $\sqrt{E(N^{-1})J_n/n}$ in sparse settings, and power tends to 1 when the jump is larger.
The $L^2$-based break-point estimator is within $O_p(\max\{1, J_n E(N^{-1})\})$ of the true break index, which is optimal in the dense case; the $L^\infty$ estimator has a slower proven rate but better finite-sample accuracy for sharp jumps.
The two-step simultaneous confidence band for the jump magnitude has asymptotically correct coverage across all sampling regimes, so post-break inference does not require a separate dense-data method.

Reading between the lines

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

Editorial inference — The BIC-based knot choice used in the simulations and applications is not proved to satisfy (A6), so the guarantees proven here may not cover the exact procedure as implemented.
Editorial inference — The threshold $E^{-1}(N^{-1}) \asymp n^{1/(2q^*)}$ coincides with the known estimation phase transition for mean functions, suggesting a general rule: change-point detection pays the nonparametric smoothing price when curves are sparse and recovers the parametric $n^{-1/2}$ rate once curves are dense.
Editorial inference — The same normalized smoothed CUMSUM construction could be adapted to test breaks in covariance operators or locally stationary functional series; the authors list both as future directions rather than established results.
Editorial inference — Taken together, the simulations imply a practical rule the paper does not formally state: use the L∞ test for spiky jump functions and the L2 test for flat ones, and interpret a rejection by either as evidence of a break.

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.

Desk Editor's Note

read the letter

Here's the short take: this is a genuinely useful unification of break-point testing for functional time series under sparse-to-dense sampling. One B-spline smoothed CUMSUM framework, both L∞ and L2 tests, a phase-transition corollary, break-point estimators, and a simultaneous band for the jump. If the supplement is right, it goes beyond Aue et al. (2018) and Bai et al. (2024) to cover irregular, unequally sampled designs. That is a real contribution.

What the paper does well: the theoretical architecture is sensible. It uses Mies and Steland's sequential Gaussian approximation and Chernozhukov-type anti-concentration, and the derived limits reduce to known results in fully observed and semi-dense special cases—that is a good sign. The simulations are extensive: four sampling schemes, three jump shapes, three error distributions, n=200 and 400. The empirical size is near nominal, power patterns make sense, and the L∞/L2 trade-off for spike-like versus smooth jumps is a nice practical insight.

Soft spots, in order of importance. First, proofs are in the supplemental material, which is not part of the posted arXiv v1. So the central derivations cannot be checked from what's here. Second, the stress-test note is on target: the unified guarantee in Theorem 3.1 is conditional on (A6), a coupled knot-moment-sampling condition, and the BIC knot selector in Section 5.1 is not shown to satisfy it. The theorem as written covers an oracle J_n in an admissible window; the tables report an automatic BIC choice. That's a genuine theory-practice gap. It's probably patchable—for the smooth mean functions and normal errors in the simulations, one can take the moment order r large enough—but as written, the size guarantees in Tables 1–2 don't directly follow from the theorem. Third, the abstract says 'arbitrary sampling schemes,' but the main assumptions require random i.i.d. design (X_{ij} with a density, N_i i.i.d.). Both applications use fixed designs (daily electricity demand, day-of-year temperature), which fall outside the theorem. The simplified fixed-design estimator in Section 4.1 has no supporting theorem. This is an overreach, not a fatal flaw.

Bottom line: solid, significant methodology, honestly presented, with a couple of gaps between the theory and the implementation. I'd send it to a serious referee and expect they'll want to see the supplement and a fix or explicit discussion of the BIC condition. I'd also bring it to a reading group.

Referee Report

3 major / 6 minor

Summary. The paper develops a unified methodology for testing and estimating a structural break in the mean function of functional time series whose trajectories are observed at irregular, possibly sparse locations. The authors construct a smoothed CUMSUM process via B-spline estimation of partial and global mean functions, then define L-infinity and L-2 test statistics. The main theoretical results are a Gaussian approximation for the null distribution (Theorem 3.1), local alternative power analysis (Theorem 3.2), convergence rates for break-point estimators (Theorem 3.3), and a simultaneous confidence band for the jump magnitude (Theorem 3.4). Corollary 3.1 gives phase-transition boundaries in terms of the sampling intensity E(N^{-1}). The paper also contains an extensive simulation study across four sampling schemes and three jump shapes, plus applications to German electricity price data and Sydney temperature data.

Significance. If the theoretical results are correct, the paper would provide a genuinely unified treatment of structural-break testing from sparse to dense functional data, and the L-infinity statistic is a useful complement to the existing L-2-based procedures. The phase-transition analysis in Corollary 3.1 is a substantive contribution, and the simulation study is unusually thorough, covering several sampling regimes, error distributions, and jump shapes. The main novelty—Gaussian approximation of a smoothed CUMSUM process without a restriction on the relationship between n and Ni—is attractive and would be of broad interest to the functional data and change-point communities. However, the posted arXiv version does not include the supplementary proofs, and the gap between the theoretical Assumption (A6) and the implemented BIC knot-selection procedure means that the paper as written does not yet establish the advertised guarantees for the practical algorithm.

major comments (3)

[Section 5.1, Assumption (A6), Theorem 3.1] The null-distribution result is stated for a user-specified number of knots J_n satisfying the coupled rate conditions in (A6), but the implementation selects J_n by minimizing BIC over the range [min{0.5(n\bar N)^{1/9}, 0.5n^{1/8}}, max{(n\bar N)^{1/7}, n^{1/6}}]. No lemma or argument shows that the BIC-selected J_n satisfies either branch of (A6), nor that the moment-bias product in (A6) tends to zero. In sparse Setting (1), E(N^{-1}) is bounded away from zero and, for n=200, the BIC range yields J_n of order 1--2; whether that lies in the admissible region depends on the effective smoothness q* and the available moment order r. If r is small (for instance r=3), the displayed moment-bias expression can grow, so the Gaussian approximation in Theorem 3.1 can fail exactly in the sparse regime the paper claims to unify. Consequently, the empirical sizes and powers in Tables 1--3 are reported for a procedure whose asymptotic guarantees are not established; the paper proves the result for an oracle statistic with J_n in (A6), not for the automatic BIC procedure actually used.
[Theorem 3.1 and the paragraph following it] The text says that 'Since no restrictions are imposed on the relationship between n and Ni in Assumptions (A1)-(A6), Theorem 3.1 contains scenarios with arbitrary sampling schemes from sparse to dense.' This is inaccurate: Assumption (A6) restricts J_n, n, E(N^{-1}), and r jointly, for example by requiring J_n^{-q*-1/2} n^{1/2} E^{-1/2}(N^{-1}) log^{1/2} n = O(1) when E(N^{-1}) J_n \gg 1. Thus 'arbitrary sampling schemes' is too strong; the theorem holds only for sampling intensities and knot numbers lying in the coupled window described by (A6). The abstract's claim of being 'adaptable to arbitrary sampling schemes' should be weakened, or the authors should show that (A6) is automatically satisfied under the stated model assumptions, which is not done.
[Supplement (end of Section 1 and beginning of Section 3)] The technical proofs of Theorems 3.1--3.4 are deferred to a supplementary file that is not included in the posted arXiv v1. Because the paper's central claims rest entirely on these Gaussian approximation results—and the reader can neither check the derivations, the precise conditions under which the phase-transition boundaries in Corollary 3.1 hold, nor the treatment of the estimated covariance matrix \hat\Sigma—the posted version is not self-contained. For a revised submission, the supplement should be made available, or the main text should summarize the key intermediate lemmas and indicate where each part is proven.

minor comments (6)

[Abstract and Section 1] The abstract contains the duplicated phrase 'test statistics statistics'; also in Section 1, 'the discrete girds' should read 'the discrete grids.'
[Theorem 3.2] In the local alternative display, the expression 'max{n√ n−1, p n−1E(N −1)Jn}' is confusing because √ n^{-1} equals n^{-1/2}; writing n^{-1/2} explicitly would improve readability.
[Section 3.2] The final sentence says 'bkn,∞ becomes more accurate than bkn,∞ with finite sample sizes,' which should presumably read 'bkn,∞ becomes more accurate than bkn,2.'
[Table 3 caption] The caption refers to 'SCRs' but the text and method name use 'SCB'; please correct the abbreviation.
[References] Several references are incomplete or informal, for instance 'Bai, L., Hu, Q., and Wu, W. (2024+)' and 'Dette, H. and Wu, W. (2024+)' are listed as 'Manuscript' without an arXiv or DOI identifier; please update these entries.
[Table 2 header] The column header 'Linfinity' should be typeset as 'L∞' for consistency with the rest of the paper.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: central Gaussian approximation and phase transition are rate-based, not fitted; self-citations are non-load-bearing.

full rationale

The derivation chain is self-contained rather than circular. Theorem 3.1 is a Gaussian approximation result built from external ingredients (Mies and Steland 2023; Chernozhukov et al. 2015; spline approximation theory from De Boor 1978 and Schumaker 2007), and the covariance of the limiting Gaussian field is computed from the smoothed CUMSUM process, not fitted. The phase-transition boundaries in Corollary 3.1 follow from comparing the orders of B^T Sigma_1 B and B^T Sigma_2 B as functions of J_n and E(N^{-1}); this is a rate calculation, not a fitted prediction. Self-citations to Cai and Hu (2024a, 2024b) supply only a parametrization idea and an oracle-efficiency reference; no theorem in the paper is justified solely by those citations, and the main proof ingredients are external. The empirical covariance estimator bSigma and the break-point plug-in are standard nuisance-parameter estimation, not circular. The implementation's BIC selection of J_n in Section 5.1 is not shown to satisfy Assumption (A6), so the practical procedure may operate outside the conditions of Theorem 3.1; this is an assumption-implementation gap and a correctness risk, not circularity, because (A6) is an explicit stated assumption rather than an output derived from the BIC selector. The paper also openly notes sub-optimality of the L_infty break-point rate as a proof artifact. Overall, no step reduces by construction to its own input.

Assumptions & free parameters 5 free parameters · 8 assumptions · 0 invented entities

The central claim rests on the structural model A1-A6 and on cited Gaussian approximation and spline theorems. No constant is fitted to make the asymptotic distributions match data; the free parameters are tuning choices for implementation. The main unverified ingredient is Assumption A6, whose rate conditions the data-driven implementation does not certify, and the treatment of fixed designs in the application is not covered by the main theorem.

free parameters (5)

Number of interior knots J_n = Data-driven via BIC; search range [min{0.5(n Nbar)^{1/9},0.5 n^{1/8}}, max{(n Nbar)^{1/7}, n^{1/6}}]
Smoothing parameter; the theory constrains its asymptotic order in Assumption (A6), and the implementation selects it by BIC in Section 5.1.
Long-run covariance bandwidth L = L = [n^{1/5}] chosen as default from range [n^{1/5}, n^{1/5} log log n]
Tuning parameter in Section 4.1 for spectral estimation of Sigma; the paper states results were similar across the range.
Trimming parameter epsilon = Any fixed value in (0,1/2); concrete value not given in simulations
Excludes change points near the boundary so partial mean estimators are well defined, as described in Section 2.
Spline order p = p=4 (cubic) in numerical experiments
Must satisfy p >= q* per Assumption (A6); cubic is the default implementation choice.
Explained-variance threshold for critical-value approximation = 0.99
Used in Section 4.1 to truncate the Karhunen-Loeve expansion when simulating Gaussian critical values.

assumptions (8)

domain assumption Model (1.1): Y_ij = m_i(X_ij) + xi_i(X_ij) + sigma(X_ij) epsilon_ij with mutually independent {xi_i}, {X_ij}, {epsilon_ij} and i.i.d. N_i.
Structural model defined in Section 1; every theorem operates inside this model.
domain assumption A1: mean functions m_i lie in H^{q,nu}[0,1] with q* = q + nu.
Smoothness assumption used for B-spline approximation error.
domain assumption A2: design density f is bounded above and away from zero.
Needed for identifiability of B-spline weighted least squares; contradicts a literal reading of arbitrary sampling schemes.
domain assumption A3-A4: sigma^2 is bounded, moments of epsilon and eta of order r >= 3 exist, and long-run covariance G is Holder with inf_x G(x,x) > 0.
Moment and regularity conditions used in the Gaussian approximation.
domain assumption A5: physical dependence measure delta_H(i,r2) = O(i^{-beta}) for beta >= 3.
Weak dependence condition on the functional time series, following Zhou and Dette (2023).
ad hoc to paper A6: coupled rate conditions on J_n, n, E(N^{-1}), and moment index r; A6' sets a_n = 1 for the L2 statistic.
Tailored high-level conditions that make the proofs go through; they are not verified in the data applications.
standard math Cited external tools: sequential Gaussian approximation of Mies and Steland (2023), comparison and anti-concentration bounds of Chernozhukov et al. (2015), and B-spline approximation properties of De Boor (1978) and Schumaker (2007).
The paper interweaves these theorems to obtain the approximating Gaussian fields, as stated in Section 1.
domain assumption B1: the break-point estimator used for residuals and jump SCB satisfies the stated O_p rate in Section 3.3.
Theorem 3.4 and the residual-based covariance estimator in Section 4.1 require a consistent break-point estimator satisfying B1.

how reviews work

0 comments

Cite this review

Pith. "Pith review of From sparse to dense functional time series: phase transitions of detecting structural breaks and beyond." pith.science (2026). https://pith.science/paper/BJHEEFHM

@misc{pith2026241220858,
  author       = {Pith},
  title        = {Pith review of: From sparse to dense functional time series: phase transitions of detecting structural breaks and beyond},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BJHEEFHM}},
  note         = {Machine review of arXiv:2412.20858}
}

abstract

We develop a novel methodology for detecting abrupt break points in mean functions of functional time series, adaptable to arbitrary sampling schemes. By employing B-spline smoothing, we introduce $\mathcal L_{\infty}$ and $\mathcal L_2$ test statistics statistics based on a smoothed cumulative summation (CUMSUM) process, and derive the corresponding asymptotic distributions under the null and local alternative hypothesis, as well as the phase transition boundary from sparse to dense. We further establish the convergence rate of the proposed break point estimators and conduct statistical inference on the jump magnitude based on the estimated break point, also applicable across sparsely, semi-densely, and densely, observed random functions. Extensive numerical experiments validate the effectiveness of the proposed procedures. To illustrate the practical relevance, we apply the developed methods to analyze electricity price data and temperature data.

Figures

Figures reproduced from arXiv: 2412.20858 by the authors.

**Figure 2.** Examples of empirical performance of estimator [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗

**Figure 3.** Left panel: The raw data before and after March 15, 2011 (before: triangle [PITH_FULL_IMAGE:figures/full_fig_p028_3.png] view at source ↗

Figures from the paper (1 more)

**Figure 4.** Figure 4: Top-left panel: of the maximum temperature before and after 1988 in Sydney, [PITH_FULL_IMAGE:figures/full_fig_p029_4.png]

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

Selfnormalization for relevant inference with supremum-type statistics
math.ST 2026-07 conditional novelty 7.0 of 10

A bias-corrected, selfnormalized statistic based on a log-sum-exp smoothing of the supremum norm yields an asymptotically pivotal test for relevant changes in functional time series.

Reference graph

Works this paper leans on

43 extracted references · 39 canonical work pages · cited by 1 Pith paper

[1]

Aue, A., Gabrys, R., Horv \'a th, L., and Kokoszka, P. (2009). Estimation of a change-point in the mean function of functional data. Journal of Multivariate Analysis , 100(10):2254--2269

work page 2009
[2]

Aue, A., Rice, G., and S \"o nmez, O. (2018). Detecting and dating structural breaks in functional data without dimension reduction. Journal of the Royal Statistical Society Series B: Statistical Methodology , 80(3):509--529

work page 2018
[3]

Baek, C., Kokoszka, P., and Meng, X. (2024). Test of change point versus long-range dependence in functional time series. Journal of Time Series Analysis , 45(4):497--512

work page 2024
[4]

Bai, L., Hu, Q., and Wu, W. (2024+). Inference for structural changes in nonstationary functional time series with partial measurement error. Manuscript

work page 2024
[5]

Bastian, P., Basu, R., and Dette, H. (2023). Multiple change point detection in functional data with applications to biomechanical fatigue data. arXiv preprint arXiv:2312.11108

work page Pith review arXiv 2023
[6]

Berger, M., Hermann, P., and Holzmann, H. (2023). From dense to sparse design: Optimal rates under the supremum norm for estimating the mean function in functional data analysis. arXiv preprint arXiv:2306.04550

arXiv 2023
[7]

Berkes, I., Gabrys, R., Horv \'a th, L., and Kokoszka, P. (2009). Detecting changes in the mean of functional observations. Journal of the Royal Statistical Society Series B: Statistical Methodology , 71(5):927--946

work page 2009
[8]

Bosq, D. (2000). Linear processes in function spaces: theory and applications , volume 149. Springer Science & Business Media

work page 2000

Show all 43 references

[9]

and Hu, Q

Cai, L. and Hu, Q. (2024a). From sparse to dense functional data: Phase transitions from a simultaneous inference perspective. arXiv preprint arXiv:2401.17646

2024 arXiv
[10]

and Hu, Q

Cai, L. and Hu, Q. (2024+b). Simultaneous inference for distribution of functional principal scores. Statist. Sinica . DOI: 10.5705/ss.202023.0246

2024
[11]

Cai, T. T. and Yuan, M. (2011). Optimal estimation of the mean function based on discretely sampled functional data: Phase transition . The Annals of Statistics , 39(5):2330--2355

2011
[12]

Cao, G., Yang, L., and Todem, D. (2012). Simultaneous inference for the mean function based on dense functional data. Journal of nonparametric statistics , 24(2):359--377

2012
[13]

Chernozhukov, V., Chetverikov, D., and Kato, K. (2015). Comparison and anti-concentration bounds for maxima of gaussian random vectors. Probability Theory and Related Fields , 162:47--70

2015
[14]

and Zhou, Z

Cui, Y. and Zhou, Z. (2023). Simultaneous inference for time series functional linear regression

2023
[15]

De Boor, C. (1978). A practical guide to splines , volume 27. springer-verlag New York

1978
[16]

and Wu, W

Dette, H. and Wu, W. (2024+). Confidence surfaces for the mean of locally stationary functional time series. Statistica sinica . DOI: 10.5705/ss.202023.0150

2024
[17]

Guo, S., Li, D., Qiao, X., and Wang, Y. (2023). From sparse to dense functional data in high dimensions: Revisiting phase transitions from a non-asymptotic perspective. arXiv preprint arXiv:2306.00476

2023 arXiv
[18]

Harris, T., Li, B., and Tucker, J. D. (2022). Scalable multiple changepoint detection for functional data sequences. Environmetrics , 33(2):e2710

2022
[19]

and Kokoszka, P

H \"o rmann, S. and Kokoszka, P. (2010). Weakly dependent functional data . The Annals of Statistics , 38(3):1845--1884

2010
[20]

and Kokoszka, P

Horv \'a th, L. and Kokoszka, P. (2012). Inference for functional data with applications , volume 200. Springer Science & Business Media

2012
[21]

Horv \'a th, L., Kokoszka, P., and Reeder, R. (2013). Estimation of the mean of functional time series and a two-sample problem. Journal of the Royal Statistical Society Series B: Statistical Methodology , 75(1):103--122

2013
[22]

and Eubank, R

Hsing, T. and Eubank, R. (2015). Theoretical foundations of functional data analysis, with an introduction to linear operators , volume 997. John Wiley & Sons

2015
[23]

Hu, Q. (2024). Change point analysis of functional variance function with stationary error. Journal of Multivariate Analysis , 202:105311

2024
[24]

and Li, J

Hu, Q. and Li, J. (2024). Statistical inference for mean function of longitudinal imaging data over complicated domains. Statistica Sinica , 34:955--982

2024
[25]

H., and Yau, C.-Y

Jiao, S., Chan, N. H., and Yau, C.-Y. (2023). Enhanced change-point detection in functional means. Statistica Sinica . DOI: 10.5705/ss.202022.0312

2023
[26]

and Liebl, D

Kneip, A. and Liebl, D. (2020). On the optimal reconstruction of partially observed functional data . The Annals of Statistics , 48(3):1692--1717

2020
[27]

Li, D., Li, R., and Shang, H. L. (2023). Detection and estimation of structural breaks in high-dimensional functional time series. arXiv preprint arXiv:2304.07003

2023 arXiv
[28]

Li, J., Chen, L., Wang, W., and Wu, W. B. (2024). 2 inference for change points in high-dimensional time series via a two-way mosum. The Annals of Statistics , 52(2):602--627

2024
[29]

and Hsing, T

Li, Y. and Hsing, T. (2010). Uniform convergence rates for nonparametric regression and principal component analysis in functional/longitudinal data . The Annals of Statistics , 38(6):3321--3351

2010
[30]

Liebl, D. (2013). Modeling and forecasting electricity spot prices: A functional data perspective. The Annals of Applied Statistics , 7(3):1562--1592

2013
[31]

Liebl, D. (2019). Nonparametric testing for differences in electricity prices: The case of the Fukushima nuclear accident . The Annals of Applied Statistics , 13(2):1128--1146

2019
[32]

M., Wang, D., Zhao, Z., and Yu, Y

Madrid Padilla, C. M., Wang, D., Zhao, Z., and Yu, Y. (2022). Change-point detection for sparse and dense functional data in general dimensions. Advances in Neural Information Processing Systems , 35:37121--37133

2022
[33]

and Steland, A

Mies, F. and Steland, A. (2023). Sequential gaussian approximation for nonstationary time series in high dimensions. Bernoulli , 29(4):3114--3140

2023
[34]

and Zhang, C

Rice, G. and Zhang, C. (2022). Consistency of binary segmentation for multiple change-point estimation with functional data. Statistics & Probability Letters , 180:109228

2022
[35]

Schumaker, L. (2007). Spline functions: basic theory . Cambridge university press

2007
[36]

and Hosseini-Nasab, S

Sharghi Ghale-Joogh, H. and Hosseini-Nasab, S. M. E. (2021). On mean derivative estimation of longitudinal and functional data: from sparse to dense. Statistical Papers , 62(4):2047--2066

2021
[37]

Trenberth, K. E. and Hoar, T. J. (1996). The 1990-1995 el ni \ n o-southern oscillation event: Longest on record. Geophysical research letters , 23(1):57--60

1996
[38]

Trenberth, K. E. and Hoar, T. J. (1997). El ni \ n o and climate change. Geophysical Research Letters , 24(23):3057--3060

1997
[39]

Wang, M., Harris, T., and Li, B. (2023). Asynchronous changepoint estimation for spatially correlated functional time series. Journal of Agricultural, Biological and Environmental Statistics , 28(1):157--176

2023
[40]

Xue, G., Xu, H., and Yu, Y. (2024). Change point localisation and inference in fragmented functional data. arXiv preprint arXiv:2405.05730

2024 arXiv
[41]

and Li, Y

Zhang, H. and Li, Y. (2022). Unified principal component analysis for sparse and dense functional data under spatial dependency. Journal of Business & Economic Statistics , 40(4):1523--1537

2022
[42]

and Wang, J.-l

Zhang, X. and Wang, J.-l. (2016). From sparse to dense functional data and beyond. The Annals of Statistics , 44:2281--2321

2016
[43]

and Dette, H

Zhou, Z. and Dette, H. (2023). Statistical inference for high-dimensional panel functional time series. Journal of the Royal Statistical Society Series B: Statistical Methodology , 85(2):523--549

2023

Pith tools

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

[1] [1]

Aue, A., Gabrys, R., Horv \'a th, L., and Kokoszka, P. (2009). Estimation of a change-point in the mean function of functional data. Journal of Multivariate Analysis , 100(10):2254--2269

work page 2009

[2] [2]

Aue, A., Rice, G., and S \"o nmez, O. (2018). Detecting and dating structural breaks in functional data without dimension reduction. Journal of the Royal Statistical Society Series B: Statistical Methodology , 80(3):509--529

work page 2018

[3] [3]

Baek, C., Kokoszka, P., and Meng, X. (2024). Test of change point versus long-range dependence in functional time series. Journal of Time Series Analysis , 45(4):497--512

work page 2024

[4] [4]

Bai, L., Hu, Q., and Wu, W. (2024+). Inference for structural changes in nonstationary functional time series with partial measurement error. Manuscript

work page 2024

[5] [5]

Bastian, P., Basu, R., and Dette, H. (2023). Multiple change point detection in functional data with applications to biomechanical fatigue data. arXiv preprint arXiv:2312.11108

work page Pith review arXiv 2023

[6] [6]

Berger, M., Hermann, P., and Holzmann, H. (2023). From dense to sparse design: Optimal rates under the supremum norm for estimating the mean function in functional data analysis. arXiv preprint arXiv:2306.04550

arXiv 2023

[7] [7]

Berkes, I., Gabrys, R., Horv \'a th, L., and Kokoszka, P. (2009). Detecting changes in the mean of functional observations. Journal of the Royal Statistical Society Series B: Statistical Methodology , 71(5):927--946

work page 2009

[8] [8]

Bosq, D. (2000). Linear processes in function spaces: theory and applications , volume 149. Springer Science & Business Media

work page 2000