1 Introduction
The (simplest) definition of fractional Brownian motion (fBm), , is that is a continuoustime Gaussian process with continuous paths which is centered, for all , and with covariance function
has stationary Gaussian increments and is selfsimilar. Except the case , when is Brownian motion, increments are not independent and is not a semimartingale. We refer to Nourdin (2012) for more information on properties of fBm. We observe a path of fBm , with a volatility (or scaling) parameter , over the unit interval at equidistant discrete observation times, . Denote , the observed increments called fractional Gaussian noise. We establish limit results under infill asymptotics as , .
The selfsimilarity parameter , called Hurst exponent, determines the persistence and the smoothness regularity of paths of . Modelling data by fBm asks for statistical inference on and . There is a rich literature on estimators for , but local asymptotic normality (LAN) and thus asymptotic efficiency of maximum likelihood (mle) and the whittle estimator under infill asymptotics has been established only recently by Brouste and
Fukasawa (2018)
. Since the computation of the mle is infeasible for too large sample sizes due to the inversion of a covariance matrix, method of moment approaches based on power variations are prominent alternative methods which are simple to implement and fast to compute. We first restrict to the more standard case
, ^{†}^{†}†For , noncentral marginal limit theorems with slower rates are given in Nourdin et al. (2010). We shall use a functional central limit theorem for a sum of squared secondorder increments from Bégyn (2007) for the test in this case.in that the following central limit theorem for the normalized quadratic variation is available:
(1) 
with given in (5), see, for instance, Nourdin et al. (2010). This setup is covered by the wellknown theorem of Breuer and Major (1983). From this convergence, we derive one standard estimator for the Hurst exponent
(2) 
and based on the delta method, when , that . When is unknown, a small modification of this statistic using a ratio of discrete quadratic variations at different sampling frequencies
(3) 
yields a consistent estimator with rate. More refined related estimation methods have been presented by Coeurjolly (2001) and by Bardet and Surgailis (2011)
using increment ratios. Lower bounds for the convergence rates of unbiased estimators of
are known to be when is known and in case of unknown volatility, see Coeurjolly and Istas (2001) or Brouste and Fukasawa (2018). Hence, the estimators using power variations attain optimal rates. The volatility can be estimated when we plugin the estimated Hurst exponent. If , these estimators preserve consistency, but the central limit theorems do no longer apply.While there is a strand of literature that addresses changepoint analysis for changes in the mean of fBm or related timeseries models, see, for instance, Betken (2017) for a recent contribution, there is scant groundwork on changes points of the Hurst exponent . There are a few works in the timeseries literature, see Lavancier et al. (2013) and references therein, which do not include a classical cusum test based on power variations. Although cusum tests are in general very popular because of their appealing asymptotic and finitesample properties, to the best of the author’s knowledge the presented methods have not yet been discussed in the literature.
A necessary prerequisite for a cusum changepoint test is a functional central limit theorem, also referred to as invariance principle. Proving such a result for power variations of fBm is difficult due to the dependence of increments. Nevertheless, exploiting properties of Gaussian processes such results have been established in the last years. Building upon two functional central limit theorems by Bégyn (2007) and Nourdin and Nualart (2018), we construct cusum changepoint tests for a change in the Hurst exponent or the volatility parameter of the fBm. We remark that we use the general result from Nourdin and Nualart (2018) only for one simple specific function, while the authors focus on necessary and sufficient conditions under which functional weak convergence holds true.
2 Statistical application of functional central limit theorems for power variations
The result by Nourdin and Nualart (2018) applied with the Hermite polynomial of rank 2 provides for a functional limit theorem for the normalized discrete quadratic variation
(4) 
as , weakly in the Skorokhod space of cádlág functions with a standard Brownian motion independent of
. The limit process being continuous, convergence holds with respect to the uniform topology. The longrun variance is determined by
(5) 
where for some function , the secondorder increment operator is . The estimate , , , shows the convergence of the series, while the sum of squared increments exhibits longrange dependence for . By continuous mapping, (4) readily implies that
(6) 
that is, the convergence of the cusum process to a Brownian bridge. This functional central limit theorem can be exploited for a changepoint test with the help of the next simple lemma.
Lemma 0.
Let and be realvalued cádlág functions defined on . For any , there exists a , such that
Proof.
If , there exists , such that
and thus . Hence, for and , we obtain the claimed continuity. ∎
Based on Lemma 1 and continuous mapping, we obtain the weak convergence of the supremum of the absolute lefthand side in (6), multiplied with , to the law of the supremum of the absolute value of a standard Brownian bridge, referred to as the KolmogorovSmirnov law. In the vein of Phillips (1987), this weak convergence can be exploited to test for structural breaks in the observed path of . Depending, however, on , and , this statistic is yet infeasible and we thus propose the following modification
(7)  
where we write . is a cusum process of squared increments standardized with the empirical longrun variance, that is, a sum over empirical autocovariances for different lags. It holds that
(8) 
with the KolmogorovSmirnov limit law. The standardization with the empirical longrun variance thus takes out the factor in (4). Moreover, we do not rescale squared increments by in (7). Multiplying numerator and denominator with , the squared denominator consistently estimates the longrun variance as , and it is positive definite by standard results from time series analysis. Slutsky’s lemma in combination with the weak convergence for the infeasible statistic based on (6) thus gives the stated asymptotic distribution of . In the next paragraph, we work out statistical properties of and construct a changepoint test. Let us recall that (8) holds true only in the case . Instead of considering functional noncentral limit theorems when , if they were available at all, and the law of the supremum of a Rosenblatt limit process, a simpler solution to address the general case, , is to use an analogous statistic as (7) with secondorder increments. For observations , set
(9)  
where are secondorder increments. A normalized sum of squared secondorder increments of fBm satisfies a (functional) central limit theorem for any , see Section 3.1 of Bégyn (2007) and we thus readily derive weak convergence of to the KolmogorovSmirnov limit law for any .
3 A cusum test for changes in the Hurst exponent
We consider the statistical hypothesis test
It is standard in the theory of statistics for highfrequency data to address such questions pathwise. This means that and are formulated for the one observed path of , and the statistical decision is based on the discretization of the given path. It is not important here how the interdependence structure evolves under around . One can think for instance of two independent fBms intertwined at the change, but a gradual change of persistence in a vicinity of is possible as well.
Consistency of the test is implied by the following lower bound for under . A standard decomposition in the expectation and a bound for the stochastic deviation by Markov’s inequality yield with an almost surely finite constant and an almost surely positive constant that
where we write and . Hence, when and are fix, with rate . Naturally, the power of the test decreases for closer to one boundary of the observation interval or for smaller . Considering a null sequence , an expansion of yields the speed . We have proved the following main result of this note.
Theorem 3.1.
The sequence of tests with critical regions
(10) 
where denotes the fractile, that is, the (1
)quantile, of the KolmogorovSmirnov law, provides an asymptotic distributionfree test for the null hypothesis
against the alternative hypothesis , which has asymptotic level and asymptotic power 1. The consistency is valid for decreasing sequences of in , as long as .Consistency of the test directly extends to the case of more than one change in . Further, a locally bounded drift term added to will not affect the results, which follows with standard estimates for highfrequency data provided by Jacod and
Protter (2012). In case of additional nuisance jumps, truncation methods can be used to obtain a robust version.
Using (9), we can drop the restriction under .
Corollary 3.2.
The sequence of tests with critical regions , provides an asymptotic distributionfree test for , for any , against with asymptotic level and asymptotic power 1. The consistency is valid for decreasing sequences of in , as long as .
For the lower bound of under , we use that are centered and
(11) 
with , and decreasing in . Analogously as for , if without loss of generality , a lower bound is thus given by
with an almost surely positive constant . This proves Corollary 3.2.
4 A cusum test for changes in the volatility
We prove that reacts as well to a change of and thus provides a statistical hypothesis test for
A related lower bound for under as in Section 3 yields with an almost surely finite constant and an almost surely positive constant that
We obtain the following corollary to Theorem 3.1.
Corollary 4.1.
The sequence of tests with critical regions (10), provides an asymptotic distributionfree test for the null hypothesis against the alternative hypothesis , which has asymptotic level and asymptotic power 1. The consistency is valid for decreasing sequences of in , as long as .
Thus, rejects the null hypothesis under both types of changes. There are opportunities to discriminate the two types of changes based on the different behavior of the ratio of expected squared increments before and after the break at time . We present one in the next paragraph. With (11) we obtain an analogous result for the test based on (9) for any under as well. In view of the rates in the LAN result by Brouste and Fukasawa (2018), we can see that our tests attain asymptotic minimaxoptimal rates. This means that parameter differences smaller than the orders stated in Theorem 3.1 and Corollary 4.1 are impossible to detect and determine the minimax detection boundary.
5 Estimation of the change point and discriminating the type of change
If a changepoint test rejects the null hypothesis of no change, the estimation of the time of a change, referred to as the break date, becomes of interest. We prove consistency of the argmaxestimator associated with our statistic (7) under the alternative hypothesis with one break in at time . The result analogously extends to a break in , and using an iterative algorithm it may be extended to multiple changes.
Proposition 5.0.
Under the alternative hypothesis , the estimator
(12) 
satisfies that .
Proof.
Without loss of generality, we consider observations , with
where and , . Generalizations of the proof to , and when one increment is affected by and are obvious. Define , a piecewise constant function with
and . We obtain for , that
We see that is nonnegative and increasing for , and decreasing for , with a unique maximum at . is the expectation of the numerator of the cusum process in (7). Thus, uniformly in . For , it holds that
, with a probability tending to 1, and the leading term is larger than the
term with a probability tending to 1. Hence, the convergence rate is implied by the relationas long as . More precisely, we apply Lemma 1 from page 45 of Bibinger and Madensoy (2018) to the function , which yields for that with a probability converging to 1 we have that
An analogous application of the same lemma to the function , yields that with a probability converging to 1 it holds that
∎
One opportunity to distinguish which type of change occurs when the test rejects the null hypothesis uses the estimator .
Corollary 5.1.
Proof.
We have by definition and since , almost surely for sufficiently large. By Proposition 2, , as . Thus, we obtain that
Analogous results can be proved for statistics based on (9) along the same lines.
6 Simulations and data example
We simulate discretizations of using the code by Coeurjolly (2000) with the Cholesky method. First, the plot in the middle of Figure 1 compares the performance of three related estimators for the Hurst exponent:

from (2).

from (3).

The increments ratio estimator from Bardet and Surgailis (2011) using secondorder increments.
The parameter configuration is and . The boxplots based on 10,000 Monte Carlo iterations reveal that has a much smaller finite sample variance compared to the two other estimators for . The robustification for general from (3) has a much larger variance. The more sophisticated increments ratio estimator has a better performance, but still cannot keep up with . According to Brouste and Fukasawa (2018) the mle, which can be computed for moderate sample sizes, allows to attain a smaller variance for general .
Next, we illustrate the empirical size and power of the test by plotting the relative amount of realizations of from (7) smaller or equal to the percentiles of the KolmogorovSmirnov law against those percentiles. The results are given in the left and right plots of Figure 1 under the null hypothesis and alternative hypotheses, respectively.
For and constant and , the left plot in Figure 1 confirms a highly accurate fit under by the KolmogorovSmirnov limit law based on 10,000 Monte Carlo iterations. For other values of , we obtain as well a high accuracy of the fit. Figure 2 shows histograms of 10,000 Monte Carlo iterations under , for , and , and different alternative hypotheses, when for and , and when for and . As expected, the power increases with increasing and with larger sample sizes. The left histogram under closely tracks the asymptotic KolmogorovSmirnov law. The factors between medians (or quantiles close to the center) for the same changes , and for the two different sample sizes, are close to the expected factors determined by the rate we have found in Section 3. For the large sample sizes, the histograms under show fat tails which, however, barely reduce the power of the test in these cases. The right plot in Figure 1 depicts the percentage type II error rates. The testing levels of interest are located on the right of the xaxis, that is, the xaxis gives the percentage values of when is the testing level. For sample size , we basically have a power of one at all reasonable testing levels in both considered cases. The power curves for in the two different scenarios are hence more informative. The ordering of the different curves is clear as we have larger power for larger and larger (compare to the histograms). Using instead of , we obtain analogous empirical distributions under the alternative hypotheses while the size under the null hypothesis is a bit better for . However, when under the null hypothesis, the left plot of Figure 3 shows that the empirical distribution of is not close to the KolmogorovSmirnov law any more, while the empirical distribution of attains a good size. Under alternative hypotheses with a change of large Hurst exponents, both statistics show similar empirical distributions again. Overall, a very good power for moderate finite sample sizes of our tests is confirmed. This is in line with classical findings in changepoint theory for parametric cusum tests.
For a real data example, we use the daily total sunspot number^{‡}^{‡}‡Source: WDCSILSO, Royal Observatory of Belgium, Brussels, http://www.sidc.be/silso/datafiles, accessed on April 3, 2019. from 1848/12/23 to 2019/03/31 having daily observations. Time series of sunspots are often analyzed using fBm models, see, for instance, Shaikh et al. (2008). Estimator (3) and the secondorder increments ratio estimator from Bardet and
Surgailis (2011) applied to all data both yield . However, , such that we clearly reject the hypothesis that is constant. A classical R/S estimate from the increments yields a different value of ca. 0.32 for . This method is often used in the applied literature. We do not rely on this R/S method, however, since Bardet (2018) notes that “a convincing asymptotic study of such an estimator” does not exist, and Taqqu
et al. (1995) have demonstrated that the obtained estimates are in general not accurate. Sunspots have a periodic behavior with at least one cycle of about 11 years, compare Figure 3. Especially for nonstationary time series, we obtained inaccurate R/S Hurst exponent estimates also in simulation experiments.
The rightplot of Figure 3 shows 21 point estimates for the Hurst exponent on time blocks with 3,000 days. Since more than 10,000 daily increments and around 2,000 secondorder increments are zero, we need to adjust the increments ratio estimator and we see some relevant differences between both estimators. Nevertheless, the empirical findings allow to conclude that a multifractional Brownian motion model is better suited and that the Hurst exponent is not constant and has larger values in a period after 1950 than before. Looking at subsamples of the time series, we find that for 3years block length the test does not reject the hypothesis of a constant Hurst exponent in these blocks at 10%level in 7 of the 63 blocks with a minimum value of . For 1year block length it is not rejected for 81 from 170 years. This indicates that within smaller time intervals fBm may be used as a suitable model.
References
 Bardet (2018) Bardet, J.M. (2018). Theoretical and numerical comparisons of the parameter estimator of the fractional Brownian motion. In Mathematical structures and applications, STEAMH: Sci. Technol. Eng. Agric. Math. Health, pp. 153–173. Springer, Cham.
 Bardet and Surgailis (2011) Bardet, J.M. and D. Surgailis (2011). Measuring the roughness of random paths by increment ratios. Bernoulli 17(2), 749–780.
 Bégyn (2007) Bégyn, A. (2007). Functional limit theorems for generalized quadratic variations of gaussian processes. Stochastic Process. Appl., 117(12), 1848 – 1869.
 Betken (2017) Betken, A. (2017). Change point estimation based on wilcoxon tests in the presence of longrange dependence. Electron. J. Statist. 11(2), 3633–3672.
 Bibinger and Madensoy (2018) Bibinger, M. and M. Madensoy (2018). Changepoint inference on volatility in noisy Itô semimartingales. Stochastic Process. Appl., forthcoming.
 Breuer and Major (1983) Breuer, P. and P. Major (1983). Central limit theorems for nonlinear functionals of Gaussian fields. J. Multivariate Anal. 13(3), 425–441.
 Brouste and Fukasawa (2018) Brouste, A. and M. Fukasawa (2018). Local asymptotic normality property for fractional Gaussian noise under highfrequency observations. Ann. Statist. 46(5), 2045–2061.
 Coeurjolly (2000) Coeurjolly, J.F. (2000). Simulation and identification of the fractional Brownian motion: a bibliographical and comparative study. J. Stat. Softw. 50(7), 1–53.
 Coeurjolly (2001) Coeurjolly, J.F. (2001). Estimating the parameters of a fractional Brownian motion by discrete variations of its sample paths. Stat. Inference Stoch. Process. 4(2), 199–227.
 Coeurjolly and Istas (2001) Coeurjolly, J.F. and J. Istas (2001). CramérRao bounds for fractional Brownian motions. Statist. Probab. Lett. 53(4), 435–447.
 Jacod and Protter (2012) Jacod, J. and P. Protter (2012). Discretization of processes. Springer.
 Lavancier et al. (2013) Lavancier, F., R. Leipus, A. Philippe, and D. Surgailis (2013). Detection of nonconstant long memory parameter. Econometric Theory 29(5), 1009–1056.
 Nourdin (2012) Nourdin, I. (2012). Selected aspects of fractional Brownian motion. Milano: Springer.
 Nourdin and Nualart (2018) Nourdin, I. and D. Nualart (2018). The functional BreuerMajor theorem. Preprint, arXive: 1808.02378.
 Nourdin et al. (2010) Nourdin, I., D. Nualart, and C. A. Tudor (2010). Central and noncentral limit theorems for weighted power variations of fractional Brownian motion. Ann. Inst. Henri Poincaré Probab. Stat. 46(4), 1055–1079.
 Phillips (1987) Phillips, P. C. B. (1987). Time series regression with a unit root. Econometrica 55(2), 277–301.
 Shaikh et al. (2008) Shaikh, Y. H., A. R. Khan, M. I. Iqbal, S. H. Behere, and S. P. Bagare (2008). Sunspots data analysis using time series. Fractals 16(03), 259–265.
 Taqqu et al. (1995) Taqqu, M. S., V. Teverovsky, and W. Willinger (1995). Estimators for longrange dependence: An empirical study. Fractals 03(04), 785–798.