Ordinary least squares estimators of variogram parameters in long-memory stochastic volatility are studied in this paper. We use the discrete observations for practical purposes under the assumption that the Hurst parameter is known. Based on the ordinary least squares method, we obtain both the explicit estimators for drift and diffusion by minimizing the distance function between the variogram and the data periodogram. Furthermore, the resulting estimators are shown to be consistent and to have the asymptotic normality. Numerical examples are also presented to illustrate the performance of our method.

1. Introduction

There has been a long history about the research of volatility modelling and estimation from historical data since the stochastic volatility was considered (see Taylor [1]). A series of studies found that the volatility of stock price often had long range dependence (LRD), meaning that the autocorrelations of volatility process decreased at a hyperbolic rate and were nonsummable (see Robinson [2]). Some work about the long-memory stochastic volatility (LMSV) has been done to describe and explain these slowly decaying autocorrelations. There were a wide variety of LMSV models, in terms of discrete and continuous time. Breidt et al. [3] and Harvey [4] independently considered an LMSV case where the log-volatility was modelled as a fractionally integrated ARMA (autoregressive, integrated, and moving average) or FARIMA process. Robinson and Zaffaroni [5] formulated a nonlinear time series model of long-memory stochastic volatility. Comte and Renault [6] and Comte et al. [7] proposed a continuous-time fractionally stochastic volatility (FSV) model which adopted a simpler representation of the fractional Brownian motion. Chronopoulou and Viens [8] did a comparison of how poor many long-memory-parameter estimators were for real data, using various nonparametric methods, and they proposed implied approach.

Correspondingly, the statistical inference problem of a volatility process under discrete and continuous cases is quite important in the implements of LMSV. In terms of the discrete-time models, Geweke and Porter-Hudak [9] proposed a log-periodogram regression method. Furthermore, the expressions for asymptotic bias and variance of the estimators were presented by Deo and Hurvich [10]. Arteche [11] developed parametric and semiparametric methods to estimate the memory parameters. Chen and Deo [12] specified and developed the generalized method of moments (GMM) in a linear long memory process. Brockwell and Davis [13] also showed frequency domain quasimaximum likelihood (FDQML) estimators that were used in long-memory stochastic volatility. However, very few papers developed the parameter estimations for long-memory stochastic volatility in the continuous-time case. Comte and Renault [6] used the log-periodogram regression approach to estimate the long-memory parameter. Their empirical study on CAC40 of the Paris Stock Exchange confirmed the short-memory feature of stock price but illustrated the long-memory feature of the volatility. Based on a first-order approximate version of the volatility process, Casas and Gao [14] estimated the parameters in a special class of LMSV models by Whittle estimation method. Chronopoulou and Viens [8, 15] used real S&P 500 data to compare the numerical performance of the GPH estimator and the implied value of by calibrating it to option price. Their empirical study showed the instability of GPH estimator.

Perhaps, the most popular approach to modelling long memory in continuous time is to use fractional Brownian motion (fBm) as a long-memory driving source in the stochastic volatility. We assume that the stochastic volatility follows a fractional Ornstein-Uhlenbeck process (fOU, see Cheridito et al. [16], Xiao et al. [17], and Hu and Nualart [18]). It is worthwhile to mention that it is usually only possible to observe asset price at points of discrete time. Therefore, statistical inference of long-memory stochastic volatility through discretely observed asset price is of great interest for practical purposes. It would be optimal to estimate the parameters of fractional Ornstein-Uhlenbeck process and the memory-parameter jointly. However, in the majority of the long-memory models, a rigorous way to jointly estimate this vector of parameters is unknown. Rosenbaum [19] provides a strongly consistent estimator of for unobserved stochastic volatility, and he proves that the convergence speed is very slow, even though his estimator reaches optimal speed. His research indicates another reason why one should try to select using calibration instead of estimation. Chronopoulou and Viens [8, 15] suggest to estimate the parameters separately. They propose using a calibration technique to fit models with various Hurst indexes to the options market, picking the model with the best . In this paper, we consider the estimators of the variogram parameters from discrete observations through the ordinary least squares method under the assumption that is known. Here, the main obstacles are the fact that volatility itself is never directly observed and variogram method is limited to stationary process. In this paper, we refer to Viens [20] for the idea of filter in the conditional probability distribution of stochastic volatility through the discrete observations of asset price. In addition, we also refer to Cheridito et al. [16] for the stationary version solution of fractional Ornstein-Uhlenbeck process and the asymptotic form of autocovariance function when the time lag goes to infinity.

This paper contributes to the expositions on theoretical and practical aspects of parameter estimation for long-memory stochastic volatility which follows a fractional Ornstein-Uhlenbeck process. Firstly, we obtain the variance of the stationary version solution to the fractional Ornstein-Uhlenbeck process, which is important in identifying the variogram function. Then, we deduce the approximate form of the ordinary least squares estimators of in long-memory stochastic volatility. Finally, the asymptotic properties are given in this paper.

The remainder of this paper proceeds as follows. In Section 2, we state the estimation procedure. In Section 3, we prove the asymptotic properties of the estimators. In Section 4, we give simulation examples to show the performance of these estimators. Concluding remarks are given in Section 5.

2. The Method of Estimation at Discrete Observation

Let be a complete probability space on which a standard Brownian motion and a fractional Brownian motion with Hurst exponent are well defined, and . We assume that the asset price process is conditionally lognormal and the volatility process is a fractional Ornstein-Uhlenbeck process. and satisfy the following equations: where is the rate of reversion, is the volatility of the volatility process, , and is compact. is a chosen deterministic function. Here, we choose Scott [21] model to make sure that the volatility remains positive; that is, . In this paper, we also assume that is independent of . Since there are no leverage effects in our model, thus the model is identical whether one uses a signed function or its absolute value.

Before discussing the problem in this paper, we need the following results taken from Cheridito et al. [16]. The stationary version solution to (1) admits the following representation: with initial condition which is explicitly contained in the representation (2). The asymptotic autocovariance of is given in Theorem 2.3 in Cheridito et al. [16] as for fixed , , and .

The following lemma, whose proof is given in Appendix, gives the variance of which is important for the identification of the variogram function.

Lemma 1. Let ; then is a Gaussian process with and for every .

Now, suppose that the price process is observed at the discrete time instants and , for some fixed length . Use and to represent and , respectively. Nourdin [22] and Neuenkirch and Nourdin [23] give the error of Euler approximations for the homogeneous one-dimensional SDEs involving fBm and having bounded coefficients. Mishura [24] presents the rate of convergence of Euler approximations of the solutions for SDE with nonstationary coefficient. Inspired by their work, we consider the Euler discretization in the form of where and are independent sequences and is an independent sequence of random variables.

For simplicity, we denote

In what follows, it is necessary to estimate the parameters first. We still denote below where is compact. Then, we will get the estimators of when is known.

Volatility itself cannot be directly observed. In practice, even for the most liquid assets or indexes, one would be restricted to using the values of the asset at discrete times. Therefore, we normalize the fluctuation sequence which is proved to be the discrete observed realization of in Lemmas 2 and 3.

Lemma 2. Let and . Then,

Proof of Lemma 2. Since is independent of , so we have Thus, where is the martingale term of .
According to Belfadli et al. [25], for all , admits a modification with ()-Hölder continuous paths. Replace by at the order . Then, Therefore, and are equivalent at the order . Additionally, Thus, this completes the proof.

Lemma 3. Let , , . Then,

Proof of Lemma 3. Since and are independent, we replace by and by at the order , respectively; it follows that which completed the proof of Lemma 3.

Applying Lemmas 2 and 3, we can get the discrete observations of volatility through asset price.


Define the variogram function as half the variance of the increments of :

For each fixed , we consider an ordinary least squares variogram (see Lahiri et al. [26]) to obtain the estimators of by minimizing a certain distance between the variogram and the data periodogram at a finite number of lags. The ordinary least squares method of estimating variogram parameters used in this paper first needs to get the variogram function.


For fractional Ornstein-Uhlenbeck process, it is not possible to get the exact variogram function. In Lemma 4, we give conditions on the time lag that ensure that is an approximation of variogram function with an approximation error less than any .

Lemma 4. For any , there are some and , such that, for any ,

Proof of Lemma 4. Applying Theorem  2.3 in [16], for any , there exist and , such that, for any , From Lemma 1, it follows that, for any , if , we have This completes the proof.

Accordingly, from (15), we have at the discrete time instants. For notational simplicity, when there is no chance of confusion, we will drop from the subscript and write them as . Let be a finite set of lag vectors. We denote and by and , respectively, where .

In this paper, our results will require the following assumption for time lag in the parameter estimation problems to be satisfied.

Assumption 5. Consider (1) , (2) ,  , where and , and (3) and is known.

Ordinary least squares (OLS) estimate method for variogram parameters fits a parametric model of valid variograms to a pointwise nonparametric variogram estimator by minimizing their certain distance at a finite number of lags. Suppose that is a variogram estimator based on a sample of size and suppose that is a valid variogram model. The most ordinary choice of distance function is given by which, when minimized with respect to and , yields the ordinary least squares (OLS) estimators of . Therefore, the OLS estimation of variogram parameters and is defined as where

In the following, we take as the moment estimator of , which follows

The main result is summarized in the following theorem.

Theorem 6. Let Assumption 5 be satisfied. Then,

Proof of Theorem 6. Let . Since is the moment estimator of , from Lemma 4, it follows that Then, follows that Substituting by its explicit expression, we have Then, (28) can be written as Removing the common term gives Solving for in (30) gives Substituting for in (28) gives Thus, the proof is completed.

Remark 7. According to the estimators of parameters and (6), we can easily get the representation of as .

3. The Asymptotic Properties of the Estimators

Theorem 8. Let Assumption 5 be satisfied; then the estimators converge to with probability 1 as goes to infinity.

Proof of Theorem 8. Theorem 3.1 in Lahiri et al. [26] gives the sufficient conditions for that OLS estimators are consistent as follows.(C1)For any , , .(C2), and has continuous partial derivatives of order () with respect to and .(C3)For any , is consistent for . Therefore, it remains to check the conditions required to complete this proof.
Since is compact, it is obvious that, for any , we have Then, it is not difficult to show that is continuous with respect to and as well as when , for any .
According to Lemma 4, for any , if ,  , we have Since is the moment estimator of , using the consistent results of Hall [28], for any , then we get Therefore, conditions (C1)–(C3) hold for our estimators. It is sufficient to deduce the desired result.

Next, we consider the distribution of the ordinary least squares estimators. We refer to Lahiri et al. [26] for the detailed notations. For any lag , let be an vector of partial derivatives of with the th parameter; that is, Also, set to be a matrix of dimension .

As derived in Theorem 3.2 in Lahiri et al. [26], if the nonparametric variogram estimator is asymptotically normal, then the OLS estimators of are asymptotically normal with the same scaling constant that yield the asymptotic distribution for . Therefore, we need to prove the asymptotically normal property of first. The asymptotic normality of is presented as follows.

Theorem 9. Let Assumption 5 be satisfied; then one has the following results:
(1) as , where the th element of is given by (2) as , where

Proof of Theorem 9. In order to prove (37), by the Cramer-Wold device, it is enough to show that, for any , and , . Fix and let and Assumption 5 hold; then define Based on Cauchy-Schwarz theorem and Lemma  A.1 in Lahiri et al. [26], there exist and satisfying the following inequations: Then, from Lemma 4, we have By Lemma A.1 in Lahiri et al. [26], we could obtain that, as , This implies that converges to in distribution as goes to infinity. Thus, we complete the proof of (37).
Finally, (37) combined with conditions (C1)–(C3) verified in Theorem 8, we conclude the desired result in (39).

4. Numerical Illustration

In this section, we present numerical examples for different values of , , and to illustrate the performance of our estimators. We sum up the estimation procedures by Monte Carlo simulation method as follows.(a)Based on Paxson [27] method, we generate fractional Gaussian noise by fast Fourier transform. Then using the result that a fractional Brownian motion is defined as the partial sums of the fractional Gaussian noise, we get the fractional Brownian motion.(b)We simulate the processes of (4) and of (5), using Euler-Mauyama method for different values of , , and . For simplicity, we set the valuation of the two variables and according to the simulation part in Xiao et al. [17]. For a fixed length , we simulate 100 sample paths on the interval using a regular partition of 2000 intervals.(c)Calculate the term . Here, we take (in order to make sure that is sufficiently small) as an example to show how to calculate in Assumption 5. If , from the formula , we need to choose . Hence, for the estimation procedure, we only use , .(d)We implement these generated data sets to obtain the estimators by (25).

The simulated mean, median value, and standard deviation of these estimators are given in Table 1. As shown in Table 1, both mean and median value of all parameters are close to the true value. The standard deviations are small which also demonstrated a good finite sample behavior of our method. There is not a decrease in quality as increases. The theoretical results in Section 2 are confirmed by these simulation studies. In summary, the simulation results match the chosen parameters exactly, which shows that our estimators perform well.

Finally, in order to illustrate the asymptotical normality for the estimators, we plot the histogram of the empirical distribution that we obtain from the Monte Carlo procedure (see Figures 1, 2, 3, 4, 5, and 6).

5. Conclusion

The process of long memory has evolved into the description of stochastic volatility during the last decades, as researchers in empirical finance have sought to use long-memory stochastic volatility models in terms of discrete time and continuous time for practical application. Among them, fractional Ornstein-Uhlenbeck process has been used by many authors for its brief form and good properties. However, no volatility process can be directly observed from the financial market. We get the approximate logarithm form of in the discrete time, which is available for us to adopt the least squares variogram method. For this reason, obtaining the estimators and would be feasible. This paper establishes an ordinary least squares estimation of the variogram parameters to deal with the case where the volatility process is described by a fractional Ornstein-Uhlenbeck process. We also give the results for the asymptotic properties of the two estimators. The numerical results show that this is an accurate and reliable method to estimate the parameters.

As stated in this paper, we assume the mutual independence between and for showing the estimation procedure clearly. In order to consider leverage effects, we would need to impose some kind of dependent structure on the covariance matrix of these two processes. Such dependent structure would make the current discussion much more complicated and we thus wish to leave such an extension for future research.


Proof of Lemma 1. Note that the fBm is a Gaussian process with for every . Hence, we can easily obtain that and is also a Gaussian process. Based on Theorem 3.1 in Nualart [29], the variance of is
Let and . Then, we obtain that
Let , . Similarly, we get that
Hence, can be written as
Thus, we deduce the desired conclusion.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors are very grateful to many anonymous reviewers for their insightful comments and suggestions. This research is supported by the National Social Science Foundation of China (no. 11&ZD156).