#### Abstract

This paper deals with the problem of estimating the Hurst parameter in the fractional Brownian motion when the Hurst index is greater than one half. The estimation procedure is built upon the marriage of the autocorrelation approach and the maximum likelihood approach. The asymptotic properties of the estimators are presented. Using the Monte Carlo experiments, we compare the performance of our method to existing ones, namely, R/S method, variations estimators, and wavelet method. These comparative results demonstrate that the proposed approach is effective and efficient.

#### 1. Introduction

It is well-known that many time series, in diverse fields of application, may exhibit the phenomenon of long-memory or long-range dependence. As a result, time series with long- memory are currently used as stochastic models in various applications including telecommunications, hydrodynamics, economics, and environment. Moreover, applications of long-memory are found in areas as diverse as energy market analysis [1], neurosciences and other biological applications [2], and more traditional financial analysis and statistical theory [3–6]. We refer to the monographs Beran [7] and Rao [8] for complete expositions on theoretical and practical aspects of long-memory processes.

In the literature, there exist many stochastic processes which exhibit long-memory property. The most popular self-similar stochastic process that exhibits long-range dependence is of course the fractional Brownian motion (fBm), which is a well-known Gaussian self-similar stochastic process with stationary increments. In fact, up to a multiplicative constant, the fBm is the only Gaussian process with these two properties. Obviously, if some phenomenon can be modeled by fBm, the estimation of the Hurst parameter in fBm is an important problem. In other words, to avoid using an arbitrary value of the unknown Hurst parameter of fBm, we should estimate the Hurst exponent. This leads to a demand for rigorous estimation procedures for fBm, which is the aim of this paper.

Actually, there are several estimation procedures for obtaining the Hurst parameter for fBm (see, e.g., [9]). One of the estimators worth mentioning is the celebrated rescaled range analysis (R/S analysis) since it is very popular among researchers until today. In fact, there are many other methods available in the literature for estimating the Hurst exponent. For example, variance-time analysis, Higuchi’s approach, correlogram method, periodogram method, Whittle estimator, wavelet method, and Detrended fluctuation analysis method. Recently, Chronopoulou and Viens [10] proposed a new estimator of the Hurst parameter based on the discrete variations. Actually, some comparisons of these estimation methods have also been investigated in that paper. Taqqu et al. [11] stated that the Whittle method is better than the classical methods using the empirical study while Abry and Veitch [12] showed that the wavelet method is better than the Whittle method. For a general comparison of Fourier and wavelet approach, see Faÿ et al. [9].

As far as the estimation is concerned, there is not yet a perfect method that is agreed by all researchers. Each method has its own drawbacks and cannot be used as a sole estimator in all cases. For more detailed discussions, see the work by Beran [7]. In this paper, for estimating the Hurst coefficient in fBm, we present a new estimation procedure, which is built upon the marriage of the autocorrelation approach and the maximum likelihood approach. Specifically, when , we present the asymptotic normality of the estimator borrowing the idea of Hosking [13]. Moreover, when , we propose the asymptotic normality of the maximum likelihood estimator under strictly weaker conditions than those employed by Lai [14]. We also describe the numerical implementation based on our method and compare the performance of our method with the other known approaches, namely, R/S method, variations estimators, and wavelet method.

The remainder of this paper is organized as follows. In Section 2, we give a brief description of the key result and present the asymptotic distributions of our estimators. In Section 3, the performance of the proposed estimator is illustrated by some numerical experiments. Finally, Section 4 draws the concluding remarks. All technical details are relegated to the appendix.

#### 2. The Estimators Based on Discrete Observations

With the dramatic increase of the importance of application of long-memory in time series, the fBM model has been successfully applied in many fields of economics, finance, physics, chemistry, medicine, and environmental studies. Indeed, when fBm is used to describe some phenomena, a crucial problem is how to identify the Hurst parameter. Thus, the parameter estimation problem for the Hurst parameter was of great interest and became a challenging theoretical problem in the past decade. In what follows, we consider the problem of estimating the Hurst index in fBm.

Now, let us recall that a fractional Brownian motion with the Hurst coefficient , defined on a complete probability space , is a centered Gaussian process. Its law is thus characterized by its covariance function, which is given by

where and and denotes the long-memory parameter. Moreover, let be the column vector of this observed time series at instances , where the superscript denotes the transpose of a vector.

##### 2.1. The Estimator When

In this subsection, we develop the estimator of the Hurst parameter using the autocorrelation function.

From (1), we can see that the covariance can be written as for , and the process is an ordinary Brownian motion. In this case the increments of the process in disjoint intervals are independent. However, for , the increments are not independent. Set . Then is a Gaussian stationary sequence with unit variance and autocovariance function

Moreover, straightforward computations show that

as .

Therefore, the autocovariance function is nonsummable for . This phenomenon is called long-range dependence, indicating (relatively) a slow decay of the covariance function. Since is exactly second-order self-similar, we have from (2) that

Thus, we can solve for to get

Now, given data , let

denote the sample mean, the sample covariance, the sample variance, and the sample autocorrelation, respectively. Based on (4), we can write the estimator for the Hurst parameter in fBm as

To assess the performance of the proposed estimate, we appeal to the following result due to Hosking [13].

Theorem 1. *Let be an exactly second-order self-similar Gaussian process, that is, a fractional Gaussian noise. Assume in (2) that . Then for a large sample size , is approximately , where
**
when and
**
when .*

*Proof. *This is a special case of Hosking [13].

##### 2.2. The Estimator When

In what follows, we deal with the problem of estimating the Hurst index in fBm when . The technique we employed here is the maximum likelihood method. The reason for choosing this approach is that this technique has been applied efficiently in a large set. Since the fBm is Gaussian the log-likelihood for the discrete observations may be explicitly computed. Borrowing the idea of Lai [14], we can obtain the exact maximum likelihood estimator of the Hurst parameter for the fBm. However, since our assumptions are slightly different from those in Lai [14], we present the detailed derivations and make some comparisons.

Note that for any , is Gaussian. Thus the joint probability density function of is where

As a consequence, the log-likelihood function of is given as

Moreover, let be the maximum likelihood estimator of . Then, based on Taylor expansion we have where is a random point between and ; the single prime and the double prime denote the first and the second derivatives with respect to , respectively. From Lai [14], we have the following result.

Lemma 2. *For , the expectation of the first derivative of the log-likelihood function, , is zero and the variance is .*

Below we will establish the asymptotic distribution for . First we list some technical conditions.

*Assumption 3. *Let denote the largest eigenvalue of the matrix . Then one assumes that .

Actually, we should verify the feasibility of Assumption 3.

Proposition 4. *Assumption 3 is reasonable when *

*Proof. *First, noting that the lower bound for the largest eigenvalue of can be obtained from the result of Walker and Mieghem [15], we have
where are the coefficients and is the minimum element of the matrix .

On the other hand, we obtain the higher bound for the largest eigenvalue of by the Gerschgorin Circle Theorem (see [16]: Theorem , P395)
with being a positive constant.

Consequently, we obtain

From (16), it is obvious that Assumption 3 is practicable when .

We are now ready to state the key result, whose proof is postponed to the appendix.

Theorem 5. *Suppose that Assumption 3 is satisfied and . Then the maximum likelihood estimator of the parameter is approximately normally distributed such that
**
as tends to infinity, where denotes the convergence in distribution.*

#### 3. Simulation Studies

In this section, we conduct Monte Carlo simulations to study the performance of our estimator and compare our estimator with existing methods (the simulation in this paper has been implemented in the Matlab language. Readers should contact the authors if they are interested in obtaining the code of this study). Although there are many approaches available for estimating the Hurst parameter, here we compare some of these methods using the simulated time series. In particular, we consider the R/S, the variations estimators, the wavelet method, and our estimators. First, for a fixed time-step , we generate the fBm for different values of the parameters and for different sample sizes: = 50 and 100. Moreover, for each case we consider independent realizations. Thus, for a given estimation method we obtain as the estimated values for . We calculate where are the estimated values for a single realization. We calculate the standard deviation (S.D.)

where is the mean over all the realizations and the mean- square error (MSE) where is the value of the parameter that we have used to generate the model.

Using this way we can obtain approximate values for the standard deviation and the mean-squares error of the estimators. These results are important because they are necessary to investigate the performances of the estimators. Simulation results are summarized in Tables 1 and 2 with the mean, the standard deviation (S.D.), and the mean- square error (MSE) of these estimators. The CPU times (in seconds), which are the average, are also presented in these tables (all the procedures were programmed using a PC with 2.4 GHz Intel Duo Core CPU and 2-GB RAM).

From these numerical computations, we can conclude that the simulated mean converges to the true value rapidly and the bias tends quickly to zero when the sample size increases. Indeed, by carefully observing Tables 1 and 2, we can see that the mean values obtained by our method and wavelet method are closer to the true values than those obtained based on R/S method and variation approach. Moreover, both the S.D. and the MSE obtained by our method and the wavelet method are smaller than those obtained based on R/S method and variation approach. Furthermore, it is interesting to see that our approach and wavelet method have the same accuracy and bias when the number of sample paths is large enough. Hence both methods are efficient. However, the most important finding is that the computation time costed by our approach is lower than that by the wavelet method. Therefore, our method is effective and efficient. In summary, the Monte Carlo simulations verify our theory and indicate that our estimators perform reasonably well in finite samples.

We next investigate the asymptotic distribution of and when is not so large. Here, the chosen parameters are , 0.70, 0.80, 0.90 and we take = 120 and = 1/12 (i.e., ). We just perform 100 Monte Carlo simulations of the sample paths generated by the process of fBm. The results are presented in Figure 1.

**(a)**

**(b)**

**(c)**

**(d)**

From Figure 1, we can see that the normal approximations of the distributions of the Hurst parameter based on are reasonable even when is not so large. This confirms our theoretical analysis: the convergence of the distributions of these two estimators is fast. In conclusion, our simulation results show that our estimators perform well since the estimating results match the chosen parameters exactly.

#### 4. Conclusion

Research on estimating the Hurst parameters of fBm has been ongoing in the econometric and statistical literatures for more than two decades. But the subject has received its greatest attention in the last decade, as researchers in empirical finance have sought to use fBm to capture the long-range dependency of the prices of financial assets. This paper considered the inference problem for fBm based on the marriage of the autocorrelation approach and the maximum likelihood approach. The main contribution of this paper was in the establishment of asymptotic normality when the Hurst parameter satisfied . This paper therefore extended the Lai [14] seminal contribution, which was done under strong conditions. Finally, we also showed the performance of our method to existing ones and illustrated that the proposed approach is effective and efficient. Certainly, for a future study, to improve the methodology it is required to use different schemes of estimation with a higher order of convergence. The field is therefore of growing importance for both theorists and practitioners.

#### Appendix

*Proof of Theorem 5. *Let . Then we obtain
where denotes the Euclidean norm. The calculation of is shown by the following computation:
where are the coordinates of the matrix . The last asymptotic term follows since .

The Gerschgorin Circle Theorem (see [16]: Theorem , P395), combined with (A.2), implies that
where is a positive constant. Consequently, we have as .

Now let . Then a simple computation shows that . To compute the variance of we first calculate the variance of :
where are the coefficients, is the eigenvalues, is the largest eigenvalue, and is the spectral norm. The last term of the asymptotic relationship follows from the Gerschgorin Circle Theorem and the positive definiteness of the matrix .

Indeed, by applying (A.2) and using (A.4), we immediately obtain that
where is a positive constant. Consequently, we obtain as goes to infinity when . Now, we obtain , which is the condition (C1) of Sweeting [17]. Moreover, the continuity condition (C2) in Sweeting [17] holds trivially since and are continuous functions of .

From Corollary 1 in Sweeting [17], we obtain the desired result.

#### Conflict of Interests

The authors declare that they have no financial and personal relationships with other people or organizations that can inappropriately influence their work; there is no professional or other personal interest of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in or the review of this paper.

#### Acknowledgments

The authors thank the anonymous referee for very helpful suggestions and comments which led to the improvement of the original paper. This research was supported by the National Natural Science Foundation of China (no. 71171086) and the Natural Science Foundation of Guangdong Province, China (no. S2013010016270).