Research Article  Open Access
Parameter Estimation for LongMemory Stochastic Volatility at Discrete Observation
Abstract
Ordinary least squares estimators of variogram parameters in longmemory 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 longmemory 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 logvolatility 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 longmemory stochastic volatility. Comte and Renault [6] and Comte et al. [7] proposed a continuoustime 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 longmemoryparameter 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 discretetime models, Geweke and PorterHudak [9] proposed a logperiodogram 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 longmemory stochastic volatility. However, very few papers developed the parameter estimations for longmemory stochastic volatility in the continuoustime case. Comte and Renault [6] used the logperiodogram regression approach to estimate the longmemory parameter. Their empirical study on CAC40 of the Paris Stock Exchange confirmed the shortmemory feature of stock price but illustrated the longmemory feature of the volatility. Based on a firstorder 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 longmemory driving source in the stochastic volatility. We assume that the stochastic volatility follows a fractional OrnsteinUhlenbeck 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 longmemory stochastic volatility through discretely observed asset price is of great interest for practical purposes. It would be optimal to estimate the parameters of fractional OrnsteinUhlenbeck process and the memoryparameter jointly. However, in the majority of the longmemory 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 OrnsteinUhlenbeck 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 longmemory stochastic volatility which follows a fractional OrnsteinUhlenbeck process. Firstly, we obtain the variance of the stationary version solution to the fractional OrnsteinUhlenbeck process, which is important in identifying the variogram function. Then, we deduce the approximate form of the ordinary least squares estimators of in longmemory 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 OrnsteinUhlenbeck 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 onedimensional 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.
Let
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.
Let
For fractional OrnsteinUhlenbeck 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 CramerWold device, it is enough to show that, for any , and , . Fix and let and Assumption 5 hold; then define
Based on CauchySchwarz 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 EulerMauyama 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).
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
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 longmemory stochastic volatility models in terms of discrete time and continuous time for practical application. Among them, fractional OrnsteinUhlenbeck 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 OrnsteinUhlenbeck 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.
Appendix
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
Set
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.
Acknowledgments
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).
References
 S. J. Taylor, Modeling Financial Time Series, John Wiley & Sons, Chichester, UK, 1986.
 P. M. Robinson, Time Series with Long Memory, Advanced Texts in Econometrics, Oxford University Press, Oxford, UK, 2003. View at: MathSciNet
 F. J. Breidt, N. Crato, and P. de Lima, “The detection and estimation of long memory in stochastic volatility,” Journal of Econometrics, vol. 83, no. 12, pp. 325–348, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. C. Harvey, “Long memory in stochastic volatility,” in Forecasting Volatility in Financial Markets, pp. 307–320, ButterworthHaineman, Oxford, UK, 2002. View at: Google Scholar
 P. M. Robinson and P. Zaffaroni, “Nonlinear time series with long memory: a model for stochastic volatility,” Journal of Statistical Planning and Inference, vol. 68, no. 2, pp. 359–371, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Comte and E. Renault, “Long memory in continuoustime stochastic volatility models,” Mathematical Finance, vol. 8, no. 4, pp. 291–323, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Comte, L. Coutin, and E. Renault, “Affine fractional stochastic volatility models,” Annals of Finance, vol. 8, no. 23, pp. 337–378, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 A. Chronopoulou and F. G. Viens, “Estimation and pricing under longmemory stochastic volatility,” Annals of Finance, vol. 8, no. 23, pp. 379–403, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 J. Geweke and S. PorterHudak, “The estimation and application of long memory time series models,” Journal of Time Series Analysis, vol. 4, no. 4, pp. 221–238, 1983. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. S. Deo and C. M. Hurvich, “On the log periodogram regression estimator of the memory parameter in long memory stochastic volatility models,” Econometric Theory, vol. 17, no. 4, pp. 686–710, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Arteche, “Gaussian semiparametric estimation in long memory in stochastic volatility and signal plus noise models,” Journal of Econometrics, vol. 119, no. 1, pp. 131–154, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 W. W. Chen and R. S. Deo, “Power transformations to induce normality and their applications,” Journal of the Royal Statistical Society. Series B. Statistical Methodology, vol. 66, no. 1, pp. 117–130, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. J. Brockwell and R. A. Davis, Time Series: Theory and Methods, Springer Series in Statistics, Springer, New York, NY, USA, 2nd edition, 2009.
 I. Casas and J. Gao, “Econometric estimation in longrange dependent volatility models: theory and practice,” Journal of Econometrics, vol. 147, no. 1, pp. 72–83, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 A. Chronopoulou and F. G. Viens, “Stochastic volatility and option pricing with longmemory in discrete and continuous time,” Quantitative Finance, vol. 12, no. 4, pp. 635–649, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Cheridito, H. Kawaguchi, and M. Maejima, “Fractional OrnsteinUhlenbeck processes,” Electronic Journal of Probability, vol. 8, no. 3, pp. 1–14, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. Xiao, W. Zhang, and W. Xu, “Parameter estimation for fractional OrnsteinUhlenbeck processes at discrete observation,” Applied Mathematical Modelling, vol. 35, no. 9, pp. 4196–4207, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Hu and D. Nualart, “Parameter estimation for fractional OrnsteinUhlenbeck processes,” Statistics & Probability Letters, vol. 80, no. 1112, pp. 1030–1038, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Rosenbaum, “Estimation of the volatility persistence in a discretely observed diffusion model,” Stochastic Processes and Their Applications, vol. 118, no. 8, pp. 1434–1462, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. G. Viens, “Portfolio optimization under partially observed stochastic volatility,” in Proceedings of the 8th International Conference on Advances in Communication and Control (COMCON '02), W. Wells, Ed., pp. 1–12, Optimization Software, Inc., 2002. View at: Google Scholar
 L. Scott, “Option pricing when the variance changes randomly: theory, estimation, and an application,” Journal of Financial and Quantitative Analysis, vol. 22, no. 3, pp. 419–438, 1987. View at: Google Scholar
 I. Nourdin, “Schémas d'approximation associés à une équation différentielle dirigée par une fonction höldérienne; cas du mouvement brownien fractionnaire,” Comptes Rendus de l'Académie des Sciences. Series I. Mathematics, vol. 340, no. 8, pp. 611–614, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Neuenkirch and I. Nourdin, “Exact rate of convergence of some approximation schemes associated to SDEs driven by a fractional Brownian motion,” Journal of Theoretical Probability, vol. 20, no. 4, pp. 871–899, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. S. Mishura, Stochastic Calculus for Fractional Brownian Motion and Related Processes, vol. 1929 of Lecture Notes in Mathematics, Springer, Berlin, Germany, 2008. View at: Publisher Site  MathSciNet
 R. Belfadli, K. EsSebaiy, and Y. Ouknine, “Parameter estimation for fractional OrnsteinUhlenbeck processes: nonergodic case,” http://arxiv.org/abs/1102.5491. View at: Google Scholar
 S. N. Lahiri, Y. Lee, and N. Cressie, “On asymptotic distribution and asymptotic efficiency of least squares estimators of spatial variogram parameters,” Journal of Statistical Planning and Inference, vol. 103, no. 12, pp. 65–85, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Paxson, “Fast approximate synthesis of fractional Gaussian noise for generating selfsimilar network traffic,” Computer Communication Review, vol. 27, no. 5, pp. 5–18, 1997. View at: Google Scholar
 A. R. Hall, Generalized Method of Moments, Advanced Texts in Econometrics, Oxford University Press, Oxford, UK, 2005. View at: MathSciNet
 D. Nualart, “Fractional Brownian motion: stochastic calculus and applications,” in International Congress of Mathematicians. Vol. III, pp. 1541–1562, European Mathematical Society, Zürich, Switzerland, 2006. View at: Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2014 Xiaohui Wang and Weiguo Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.