Abstract
The paper illustrates a new procedure for estimating asymmetric stochastic volatility models. These models shape the asymmetric effect of negative and positive financial returns on the expected volatility, behaviour often observed in the stock prices, and known as “leverage effect”. The procedure is based on the iterative application of the quasi-maximum likelihood (QML) method and is proposed as an alternative to the procedure presented by Harvey and Shephard in 1996 and based on the application of the QML method on a modified auxiliary model. The estimation results generally converge to constant values after a few iterations. The volatility predictor provided by the new method is conceptually similar to the EGARCH predictor and different from the predictor of the other procedure. A simulation study shows that the iterative QML method provides parameter estimators with RMSEs decreasing as series length increases. The distribution of the estimates is approximately normal, and the approximation improves as the series size increases. Empirical applications of the method provide results similar to ones of the method known in literature. However, the two methods provide two different predictors and smoothers of volatility, which should be compared on a case-by-case basis.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
According with a very common pattern for financial returns, \(r_t\), the term volatility refers to the coefficient \(\sigma _t\) into the equation:
where \(\varepsilon _t\) is a unit-variance innovation, Gaussian distributed in the simplest case; when \(\sigma _t\) does not depend on \(\epsilon _t\), \(\mu _t\) is the expected value of \(r_t\) given all information known at time \(t-1\) (Taylor 1986).
Two main classes of volatility models are known in literature:
-
Models for conditional heteroskedasticity;
-
Stochastic Volatility models.
In the models for conditional heteroskedasticity, generally identified as GARCH models (Bollerslev 1986), volatility is a function of the previous information, \(I_{t-1}\), and \(\sigma ^2_t\) corresponds to the conditional variance of the returns: \(\sigma _t^2 = Var(r_t \mid I_{t-1})\).
An alternative approach in volatility modelling consists in considering \(\sigma ^2_t\) as a latent stochastic process, whose logarithm, \(h_t\), is usually represented as an autoregressive process, in most cases, of order one:
Models in the form (1)–(2) are reported in literature with the general name of Stochastic Volatility (SV) models (Melino and Turnbull 1990; Taylor 1994) and are discrete approximations to various diffusion processes proposed in the asset-pricing theory (Hull and White 1987; Wiggins 1987).
Unlike GARCH family models, which are easily estimated by the Maximum Likelihood (ML) method, the ML estimation of SV models is a tough challenge since the likelihood is defined by a T-dimensional integral, which is hard to manage (Harvey and Shephard 1996).
In accord with Sandmann and Koopman (1998), methods to face SV models estimation can be subdivided into two groups: (i) methods oriented to rebuild the exact likelihood of the SV model or a related model; (ii) methods based on more workable, but sub-optimal methods. Among all the methods, the Quasi Maximum Likelihood (QML) method (Ruiz 1994; Harvey et al. 1994) appears to be a good compromise between simplicity and efficiency: it maximizes a function that is not the actual exact likelihood function, but it provides an optimal linear estimator (and predictor) of \(h_t\), consistent and asymptotically Gaussian, according with the results of Dunsmuir (1979).
The method is based on the Kalman filter (Harvey 1990), which filters past and present volatilities and predicts future ones; Ruiz (1994) suggest that the method works well for sample sizes usually used in financial economics.
On the other hand, the application of the basicFootnote 1 QML faces some issues when the volatility innovation, \(\eta _t\), and the return innovation, \(\varepsilon _t\), are correlated, as explained in Sect. 2. This correlation is fundamental to model the asymmetric effect of negative and positive returns on the expected volatility, behaviour often observed in stock prices, and known as leverage effect (Christie 1982; Engle and Ng 1993).
This paper presents a new procedure for applying the QML method to asymmetric SV models, which produces a volatility predictor different from the one known in the literature, preferable to the latter in some cases. A brief presentation of the known procedure is reported in Sect. 2, while the new procedure is illustrated in Sect. 3. In Sect. 4, a simulation study evaluates the properties of the estimators obtained with the new proposal. Section 5 illustrates the application of the old and new procedure on three financial series. The achieved goals are summarized in Conclusions.
2 Leverage effect and stochastic volatility
In financial series increases of volatility are often observed after negative returns, and generally the greater the extent of the loss the greater the increase. This evidence is called “leverage effect” as price drop in a stock decreases the value of the firm equity, and increases the leverage-ratio. The increased leverage-ratio will involve higher risk on the equity which will be more volatile during next period (Black 1976).
A way to include the leverage effect in (2) consists in assuming the volatility innovation \(\eta _t\) correlated with the return innovation \(\varepsilon _t\), as in the following general Gaussian SV model:
where: \(x_t=r_t-\mu _t\) is the mean-adjusted return on an asset, simply “return” in the following; \(h_t=\ln \sigma ^2_t\) is the log-square volatility; \(\varepsilon _t\) and \(\eta _t\) are zero-mean Gaussian innovations, serially independent, with \(E(\eta _t\varepsilon _{t-k})=\gamma\) if \(k=0\), and zero otherwise; then, the correlation between \(\varepsilon _t\) and \(\eta _t\) is \(\rho =\gamma /\sigma _{\eta }\).
In Model (3), the leverage effect on \(h_{t+1}\) corresponds to:
which is proportional to the size of \(\varepsilon _t\) and of opposite sign if \(\gamma <0\).
The basic QML estimation of a Gaussian SV model would consists in the estimation, via Kalman Filter, of the auxiliary state space model:
where: \(y_t=\ln x^2_t\), and \(\xi _t = \ln \varepsilon ^2_t + 1.27\); \(-\) 1.27 and 4.93 are the values of \(E(\ln \varepsilon ^2_t)\) and \(Var(\ln \varepsilon ^2_t)\) respectively, when \(\varepsilon _t\) is Gaussian (Zelen and Severo 1972).
Unfortunately, the basic QML does not allow to estimate the parameter \(\gamma\) (or \(\rho\)) directly since the disturbances \(\xi _t\) and \(\eta _t\) are uncorrelated when the distribution of \(\varepsilon _t\) and \(\eta _t\) is symmetric, e.g. Gaussian or Student’s t, so that the original covariance (correlation) between \(\varepsilon _t\) and \(\eta _t\) is definitely lost in (5) (Harvey et al. 1994).
Harvey and Shephard (1996) proposed a method, here briefly named HS-QML, to estimate the parameters of a SV model with correlated disturbances. According with this method, the parameters of (3) can be estimated applying the QML method to the time-varying linear state space modelFootnote 2:
where \(s_t\) indicates the sign of \(x_t\): it is the equal to 1 (\(-\,1\)) when \(x_t\), i.e. \(\varepsilon _t\), is positive (negative).
The one step ahead predictor of \(h_t\) provided by the Kalman Filter (see Appendix A) combines the mechanisms of the EGARCH (Nelson 1991) and the Threshold-ARCH predictors (Glosten et al. 1993; Zakoian 1994):
where \(\kappa _t\) is the gain of the Kalman filter (Harvey 1990); \({\hat{y}}_{t\mid t-1}=-1.27+{\hat{h}}_{t\mid t-1}\) is the one step ahead prediction of \(y_t\); \(P_{t\mid t-1}=E[({\hat{h}}_{t\mid t-1}-h_t)^2]\) is the MSE of \({\hat{h}}_{t\mid t-1}\).
If \(\gamma <0\), Formula (8) entails that: \((\kappa _t\mid P_{t\mid t-1}, s_t=-1)>(\kappa _t\mid P_{t\mid t-1}, s_t=+1)\). Therefore the leverage effect predicted by (7) can be broken down into two parts: (i) a part, \(\gamma s_t\), depending on the sign but not on the size of \(\varepsilon _t\) (as in the Threshold-ARCH model); (ii) a part, \(\kappa _t (y_t-{\hat{y}}_{t\mid t-1})\), depending on the sign and the size of \(\varepsilon _t\).Footnote 3 This forecast only partially reflects the characteristics of the leverage effect as expressed in (4).
3 Iterative QML for asymmetric SV models
Having assumed that \(\eta _t,\) and \(\varepsilon _t\) are bivariate normal with \(E(\eta _t\varepsilon _t)=\gamma\), the auxiliary model (5) can be rewritten as:
where the disturbance \(\eta ^+_t\) represents the exogenous innovation on \(h_t\).
As \(\varepsilon _t=x_t/\exp (h_{t}/2)\), Model (9) is not linear, then it cannot be estimated by means of basic QML. Nevertheless, we can perform the following procedure:
-
1.
\(\varepsilon _t\) is initially set to \(x_t/s_X\), being \(s_X\) the sample standard deviation of the series of returns \((x_1, x_2,\ldots ,x_T)\);
-
2.
Model (9), now linear, is estimated by basic QML;
-
3.
\(\varepsilon _t\) is updated by \(x_t/\exp {({\tilde{h}}_{t}/2)}\), where \({\tilde{h}}_{t}\) are the smoothed \(h_t\), provided by the Kalman smoother;
-
4.
steps 2 and 3 are repeated successively according with a pre-set stopping ruleFootnote 4.
Empirical evidence shows that the parameter estimates converge to realistic values after few steps. The trick of the iterative procedure, IQML in the following, consists in smoothing \(\varepsilon _t\) and \(h_t\) not conjointly: treating separately \(\varepsilon _t\) and \(h_t\) does not compromise the linearity of the second equation into (9).
The IQML can be viewed like a variant of the EM approach (Dempster et al. 1977) proposed by Shumway and Stoffer (1982) for smoothing and forecasting time series, but some differences should be highlighted. In the case of Model (9), the approach of Shumway and Stoffer (1982) would consist in: (i) setting initial values of the model parameters; (ii) smoothing \(h_t\) in order to build a likelihood function (expectation step); (iii) estimating the model parameters maximizing the likelihood function (maximization step); (iv) iterating steps (ii)–(iii) until the estimates and the likelihood function are stable. Nevertheless the direct smoothing of \(h_t\) (step ii) is complicated by the non-linearity of the model. This problem is bypassed in IQML because \(\varepsilon _t\) and \(h_t\) are smoothed separately: \(\varepsilon _t\) is smoothed before the maximisation (estimation) step, whereas \(h_t\) is smoothed after the maximization step in order to update \(\varepsilon _t\). As a result the maximisation step differs between the algorithms: in IQML it consists in maximizing a pseudo (quasi) log-likelihood given \(y_t\) and \({\tilde{\varepsilon }}_t\); in EM it consists in the ML estimation of a regression model, the second equation in (9), given \({\tilde{h}}_t\).Footnote 5 The IQML procedure is more viable than the EM approach described above, but involves an inevitable loss of efficiency due to smoothing \(\varepsilon _t\) and \(h_t\) separately.
It is interesting to note the form of the predictor of \(h_t\) provided by the Kalman Filter:
where \(\kappa _{t}\) is the gain of the Kalman Filter. Since Model (9) satisfies the steady-state conditions,Footnote 6\(\kappa _{t}\) converges to a constant that we name \(\alpha\). Therefore, the steady-state predictor of \(h_t\) can be formalized as:
where the function \(\xi (| {\tilde{\varepsilon }}_{t}| )=2\ln | {\tilde{\varepsilon }}_{t}| +1.27\) is a monotonically increasing, mean-corrected, function of the (approximated) magnitude of \(\varepsilon _t\). Predictor (11) shows a clear similarity with the EGARCH predictor in which \(\xi (| \varepsilon _{t}| )=|\varepsilon _{t}| -\sqrt{2/\pi }\) is the mean-corrected magnitude of \(\varepsilon _t\) when \(\varepsilon _t\sim N(0,1)\). Predictor (11) may be also viewed as a Log-GARCH predictor (Geweke 1986; Pantula 1986) with leverage effect.
The leverage effect predicted by (11) tries to replicate the form of the leverage effect Model (4): it is proportional to the amplitude of the estimated return innovation, \(\tilde{\varepsilon _t}\), and of opposite sign (if \(\gamma <0\)). Formulas (7) and (11) show that the HS-QML and IQML methods involve two different predictors of \(h_t\), whose performances are evaluated in Sects. 4 and 5.
3.1 Student’s \(\textbf{t}\) return innovations
The IQML method can be generalized to the case where \(\varepsilon _t\) has a scaled Student’s t-distribution with v degrees of freedom, scaled in order to have unit variance, i.e. \(\varepsilon _t\sim t_v\sqrt{(v-2)/v}\). In this case (see “Appendix B”):
being \(\psi _0\) and \(\psi _1\) the digamma and trigamma function, respectively (Davis 1972).
If the parameter v is assumed known, the model estimation procedure is the IQML described above, with the values \(-\)1.27 and 4.93 now replaced by the values of \(g_1(v)\) and \(g_2(v)\), respectively. Alternatively, \(-\)1.27 and 4.93 are replaced by the parametric formulas of \(g_1(v)\) and \(g_2(v)\), and v is treated as an additional unknown parameter.
4 Finite sample properties of the IQML estimator
By means of a simulation study, Harvey and Shephard (1996) provided empirical results in favour of the consistency of the estimates obtained with the HS-QML method. In the study, series of different lengths, T, were simulated 1000 times (\(n=1000\)) from Model (3) with “empirically reasonable” parameters values. On each series, the model parameters were estimated using the HS-QML method, under some practical constraints.Footnote 7
Table 1 reports some estimation results of that simulation study: the average and the root mean square error, RMSE, (figures in brackets) of the estimates of \(\beta\), \(\ln \sigma ^2_{\eta }\) and \(\rho\); \(\ln \sigma ^2_{\eta }\) was preferred to \(\sigma ^2_{\eta }\) because the estimates of the first are closer to be normally distributed then ones of the second. We can note that the RMSEs decrease as the series length increases and, ceteris paribus, they decrease if the absolute value of \(\rho\) (i.e. \(\gamma\)) increases.
The same approach, based on simulations, is followed here to asses the finite sample properties of the IQML estimators.Footnote 8 With the same settings adopted by Harvey and Shephard, 1000 series were simulated from Model (3), and on each series the model was estimated with the IQML method (Table 2). As in HS-QML, the RMSEs of the IQML estimators decrease as the strength of the correlation increases and as the length of the series increases. The RMSEs of the new method are slightly lower than ones of the other method.
The normal distribution of the IQML estimates has been tested using the Shapiro-Wilks and the Jarque-Bera tests. Table 3 reports the p-values of the tests on the estimates of \(\omega\), \(\gamma\), \(\beta\), and \(\ln \sigma ^2_{\eta }\), both with \(T=3000\) and \(T=6000\). With series of length \(T = 6000\), the estimates of \(\omega\), \(\gamma\) and \(\ln \sigma ^2_{\eta }\) can be already considered normally distributed with a significance level of 0.03. On the other hand, the sampling distribution of \({\hat{\beta }}\) cannot yet be approximated by the normal, this is because the simulated value, 0.975, is very close to the upper limit in case of stationary volatility (\(\beta =1\)). As a result, the finite sample distribution of \({\hat{\beta }}\) has a slightly longer left tail than the right one, although the skewness decreases the more negative the correlation \(\rho\) (Fig. 1). In this case the Gaussian approximation should be possible with series length greater than those considered.
4.1 Goodness of the IQML filtered and smoothed volatilities
A simulation was conducted to evaluate the goodness of the filtered and smoothed volatilities provided by the IQML and HS-QML methods. To this end: (i) 100 series of length \(T=1000\) were simulated from model (3) with \(\beta =0.97\), \(\rho =-0.90\), \(\ln (\sigma ^2_{\eta })=-4\) (i.e. \(\gamma =-0.122\) and \(\sigma ^2_{\eta }=0.0183\)); (ii) on each series, the IQML and HS-QML methods were applied to filter and smooth the log-square volatility \(h_t\) using the estimated parameters provided by each method, then filtered and smoothed \(\sigma ^2_t\)s were obtained by exponential transformation of the corresponding log-square volatilities; (iii) the closeness of the filtered and smoothed \(\sigma ^2_t\)s to the simulated ones was measured on each series using the (average) loss functions MSE, MAE and Qlike (Patton 2011; Hansen and Lunde 2005):
The loss function (16) corresponds to the variant of Qlike proposed by Patton (2011, p. 252) in order to make the function homogeneous. This variant is formally the relative difference of \(\sigma ^2_t\) from \(\exp {h}_t\) minus the corresponding log-difference. Given the relationship between the two differences, \(\mathrm {Qlike^*}\) grows as the gap between \(\sigma ^2_t\) and \(\exp h_t\) increases, and grows more when the gap is negative, i.e. \(\sigma ^2_t\) is underestimated (Fig. 2). On the other hand, MSE and MAE present symmetric effects of overestimation and underestimation, but MSE is more sensitive than MAE to large gaps between \(\sigma ^2_t\) and \(\exp h_t\) (due to sharp changes in volatility or outliers).
Table 4 reports the overall average of the loss functions for the methods,Footnote 9 and the the number of series, \(n^*\), where IQML determines the lowest loss. We can see that the values of all three loss functions are lower when the IQML method is used, both in case of filtered and in case of smoothed \(\sigma ^2_t\)s. These results seem to suggest that the IQML method provides filtered and smoothed \(\sigma ^2_t\)s that fit better the actual ones. In particular, the IQML method appears to better limit: (i) large gaps between actual volatility and filtered (smoothed) volatility; (ii) underestimation of actual volatility.
Finally, the \(n^*\) counter shows that the IQML method outperforms the HS-QML in a large percentage of simulated series.
5 Empirical Applications
The IQML method was applied on several financial series for the estimation of the asymmetric SV model (3). This section illustrates the application of the method on three series of financial indices:
-
NASDAQ Composite (IXIC)
-
DAX index (GDAXI)
-
CAC40 index (FCHI)
in the period from 04-01-2021 to 30-12-2022. The returns we consider are percentage log differences of the daily index closing values. The estimation results are compared with those obtained with the HS-QML method.Footnote 10
Table 5 reports the estimates of \(\omega , \gamma , \beta\), and the standard deviation of the exogenous innovation on \(h_t\) (i.e. \(\sigma _{\eta ^+}\)). The standard errors of the estimates are calculated on the basis of the “Outer Product of the Gradient” (OPG) method.Footnote 11 The quasi log-likelihood, lly, of models (6) and (9) is also reported.
The two methods provide fairly close estimates, also standard errors are pretty close. We note that the IQML standard errors are moderately smaller than the corresponding HS-QML standard errors in the IXIC and FCHI series, but not in GDAXI. The quasi-log likelihood values, lly, of the two methods are also very close in each series: the lly of IQML is slightly higher in IXIC and FCHI, and slightly lower in GDAXI.
The goodness of the filtered and smoothed \(\sigma ^2_t\)s provided by the methods is assessed by the MSE (14) and \(\mathrm {Qlike^*}\) (16) loss functions computed using the square return \(x^2_t\) as proxy of \(\sigma ^2_t\) (see Table 6). As stated by Patton (2011), these loss functions are robust to noise when the volatilities proxy are the square returns. Based on these results, neither method appears clearly better than the other: the IQML seems little better for the IXIC filtered \(\sigma ^2_t\)s and the GDAXI smoothed \(\sigma ^2_t\)s; the HS-QML seems little better for the IXIC smoothed \(\sigma ^2_t\)s; for the FCHI series, IQML seems preferable if the criterion is MSE, but HS-QML could be better if Qlike is the criterion. Nevertheless, all differences in the loss functions are too small to highlight a clear superiority of one method over the other. In the case of financial series, the choice of the most appropriate method should take place on a case-by-case basis, considering more than one criterion.
6 Conclusion
The IQML method consists in iterating the basic QML method over an asymmetric SV model. The procedure is made possible using a proxy of the return innovation \(\varepsilon _t\). This simple procedure allows the user to estimate the parameters of an SV model in which the return innovation and volatility innovation are correlated. This goal can also be achieved with the modified QML method proposed by Harvey and Shephard (1996) (HS-QML), but the two methods provide different volatility predictors. The IQML predictor is conceptually similar to the EGARCH predictor, whereas the HS-QML predictor is more similar to the Threshold-ARCH predictor. A simulation study shows that the IQML filtered and smoothed square volatilities generally fit simulated \(\sigma ^2_t\)s better than those of the other model do. On the other hand, empirical applications suggest comparing the two methods, using loss functions, to identify the most suitable for the series under study.
The simulation study also shows that IQML estimators exhibit decreasing RMSEs as the series length increases and finite sample distributions that can be approximated by the Gaussian distribution; the approximation improves as the sample size increases.
The IQML method can be viewed as a variant of the Expectation-Maximization (EM) algorithm proposed by Shumway and Stoffer (1982), although the algorithms differ in some methodological aspects as specified in Sect. 3.
Finally, the method can be applied to the case with Student’s t return innovations in order to treat returns with high kurtosis.
Data availibility
The data used in the work can be requested directly from the author.
Notes
See “Appendix A”.
\((y_t-{\hat{y}}_{t\mid t-1})\) can be viewed as an approximation of \(\xi _t=\ln \varepsilon ^2_t+1.27=2\ln | \varepsilon _t| +1.27\).
Given \({\tilde{h}}_t\) (and \({\tilde{\varepsilon }}_t=x_t/\exp ({\tilde{h}}_{t}/2)\)), the model estimation concerns only the second equation in (9).
See Harvey and Shephard (1996) for more details.
The procedure was performed using the econometric software Gretl, version 2022a.
Results are reported as fraction of the corresponding HS-QML loss. Since the functions are homogeneous, the use of relative values is not affected by the scale of measurement of the returns (percentages or decimals).
The HS-QML method is applied using the equivalent model form (A.5) explained in Appendix (A). This form puts in evidence the standard deviation of the exogenous innovation on \(h_t\).
See Gretl User’s Guide, version of January 2021, p. 222 (Cottrell and Lucchetti 2021).
References
Black F (1976) Studies of stock market volatility changes. In: Proceedings of the American statistical association, business & economic statistics section
Bollerslev T (1986) Generalized autoregressive conditional heteroskedasticity. J Econ 31:307–327
Christie AA (1982) The stochastich behaviour of common stock variances: value, leverage, and interest rate effects. J Finan Econ 10:407–432
Cottrell A, Lucchetti R (2021) Gretl user’s guide. Distributed with the Gretl library
Davis P (1972) Chapter 6. Gamma function and related functions. In: Abramovitz M, Stegun IA (eds) Hanbook of mathematical functions. Dover, New York
Dempster AP, Laird NM, Rubin DB (1977) Maximum likelihood from incomplete data via the EM algorithm. J R Stat Soc Ser B (Methodological) 39(1):1–22
Dunsmuir W (1979) A central limit theorem for parameter estimation in stationary vector time series and its application to models for a signal observed with noise. Ann Stat 490–506
Engle RF, Ng VK (1993) Measuring and testing the impact of news on volatility. J Financ 48(5):1749–1778
Geweke J (1986) Modelling the persistence of conditional variances: a comment. Econ Rev 5:57–61
Glosten L, Jagannathan R, Runkle D (1993) On the relation between the expected value and volatility of nominal excess returns on stocks. J Financ 46:1779–1801
Hamilton J (1994) Time series analysis. Princeton University Press, Princeton, NJ
Hansen PR, Lunde A (2005) A forecast comparison of volatility models: does anything beat a Garch(1,1)? J Appl Econ 20(7):873–889
Harvey A, Shephard N (1996) Estimation of an asymmetric stochastich volatility model for asset returns. J Bus Econ Stat 14(4):429–434
Harvey A, Ruiz E, Shepherd N (1994) Multivariate stochastic volatility methods. Rev Econ Stud 61:247–264
Harvey AC (1990) Forecasting, structural time series models and the Kalman filter. Cambridge University Press, Cambridge
Hull J, White A (1987) The pricing of options on assets with stochastic volatility. J Financ 42:281–300
Melino A, Turnbull SM (1990) Pricing foreign currency options with stochastic volatility. J Econom 45(1–2):239–265
Nelson D (1991) Conditional heteroskedasticity in asset returns: a new approach. Econometrica 59:347–370
Pantula S (1986) Modelling the persistence of conditional variances: a comment. Econom Rev 5:71–73
Patton AJ (2011) Volatility forecast comparison using imperfect volatility proxies. J Econom 160(1):246–256
Ruiz E (1994) Quasi-maximum likelihood estimation of stochastic volatility models. J Econom 63(1):289–306
Sandmann G, Koopman SJ (1998) Estimation of stochastic volatility models via Monte Carlo maximum likelihood. J Econom 87(2):271–301
Shumway RH, Stoffer DS (1982) An approach to time series smoothing and forecasting using the EM algorithm. J Time Ser Anal 3(4):253–264
Taylor S (1986) Modelling financial time series. Wiley, Chichester
Taylor SJ (1994) Modeling stochastic volatility: a review and comparative study. Math Financ 4(2):183–204
Wiggins JB (1987) Option values under stochastic volatility: theory and empirical estimates. J Financ Econom 19(2):351–372
Zakoian JM (1994) Threshold heteroskedastic models. J Econom Dyn Control 18(5):931–955
Zelen M, Severo N (1972) Chapter 26. Probability function. In: Abramovitz M, Stegun IA (eds) Hanbook of mathematical functions. Dover, New York
Funding
Open access funding provided by Università degli Studi del Piemonte Orientale Amedeo Avogrado within the CRUI-CARE Agreement. This work was supported by funding from the Università del Piemonte Orientale.
Author information
Authors and Affiliations
Contributions
Paolo Chirico is the sole author of the work.
Corresponding author
Ethics declarations
Conflict of interest
Author declare that there is no conflict of interest.
Consent to participate
The work was carried out in compliance with the ethical code of the Università del Piemonte Orientale.
Code availability
The gretl script for IQML can be requested directly from the author.
Consent for publication
The author agrees for publication.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendix A Proofs of Model (6)
The original method of Harvey and Shephard (1996) concerns the estimation of the following asymmetric SV model:
where \(\varepsilon _t\) is a process of independent, identically distributed random disturbances, symmetrically distributed around 0. The authors demonstrated that model (A.1) can be estimated applying the QML method to the time varying state space model:
where:
-
\(y_t=\ln x^2_t\);
-
\(\omega ^*=\ln \sigma ^2+E[\ln \varepsilon ^2_t]\);
-
\(\xi _t=\ln \varepsilon ^2_t-E[\ln \varepsilon ^2_t]\);
-
\(s_t= +1 (-1)\) when \(x_t\), i.e. \(\varepsilon _t\), is positive (negative);
-
\(\mu ^*\) and \(\gamma ^*\) correspond to \(E(\eta _t\mid s_t=1)\) and \(Cov(\eta _t, \varepsilon _t\mid s_t=1)\) respectively.
Harvey and Shephard (1996) proved that the relations of \(\mu ^*\) and \(\gamma ^*\) with \(\gamma\) depend on the distribution of \((\eta _t, \varepsilon _t)\). When \((\eta _t, \varepsilon _t)\) is bivariate normal, it results (approximately) that:
Using the restrictions (A.3) and the transformation \(h_t=h^*_t+1.27+\omega ^*\), Model (A.2) assumes, after some algebra, the form (6):
where \(\omega =(1-\beta )(1.27+\omega ^*)\).
Model (A.4) can also be represented in the alternative form:
where \(\mathbf {c'}=[\sqrt{4.93}, \hspace{0.2cm} 0]\), \(\mathbf {b'_t}=[1.11\gamma s_t/\sqrt{4.93}, \hspace{0.2cm} \sigma _{\eta ^+}]\); \(\sigma _{\eta ^+}\) is the standard deviation of the exogenous innovation on \(h_t\).
Taking into account Formulas (3.2.3c), (3.2.4a) and (3.2.22) in the book of Harvey (1990), the one step ahead predictor of \(h_t\) is:
where: \({\hat{y}}_{t\mid t-1}=-1.27+{\hat{h}}_{t\mid t-1}\) is the one step ahead prediction of \(y_t\); \(P_{t\mid t-1}=E[({\hat{h}}_{t\mid t-1}-h_t)^2]\) is the MSE of \({\hat{h}}_{t\mid t-1}\).
Appendix B Proofs of Formulas (12) and (13)
A scaled Student’s t-distribution, with v degrees of freedom, is the distribution of a random variable, \(\zeta _v\), defined as:
where Z is a standard normal r.v., and \(\chi ^2_v\) is a chi-square r.v. with v degrees of freedom, independent of Z.
It is easy to derive that: \(E(\zeta _v)=0\) and \(Var(\zeta _v)=1\).
Given \(\xi _v=\ln \zeta _v^2\), it follows:
From Zelen and Severo (1972) we know that:
where \(\psi _0\) and \(\psi _1\) are the digamma and trigamma function, respectively (Davis 1972).
Since \(\psi _0(1/2)\approx -1.963\), \(\psi _1(1/2)\approx 4.93\), and \(\ln (1/2)\approx -0.693\), it follows approximately that:
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chirico, P. Iterative QML estimation for asymmetric stochastic volatility models. Stat Methods Appl 33, 885–900 (2024). https://doi.org/10.1007/s10260-024-00747-z
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10260-024-00747-z