Home Volatility Modeling with Leverage Effect under Laplace Errors
Article
Licensed
Unlicensed Requires Authentication

Volatility Modeling with Leverage Effect under Laplace Errors

  • Zhengjun Jiang and Weixuan Xia EMAIL logo
Published/Copyright: July 18, 2017
Become an author with De Gruyter Brill

Abstract

This paper discusses four GARCH-type models (A-GARCH, NA-GARCH, GJR-GARCH, and E-GARCH) in representing volatility of financial returns with leverage effect. In these models, errors are assumed to follow a Laplace distribution in order to deal with the typical leptokurtic feature of financial returns. The properties of these models are analyzed theoretically in terms of unconditional variance, kurtosis, autocorrelation function and news impact, and are further examined in the applications to real financial time series. Comparison is made with other choices of error distributions such as normal, Student-5, and Student-7 distributions, respectively. We also conduct residual analyses to justify the choice of error distributions and find that Laplace-E-GARCH model still performs very well. Our main purpose is to study and compare the proposed models’ relative adequacies and underlying limitations.

Appendix

A Derivations of Results in Section 3

In the A-GARCH model under Laplace errors, eq. (7) is deduced from

(38)σr2E[σt2]=ω+αE[rt12]+βE[σt12]+δE[rt1]=ω+(α+β)σr2.

Equation (9) is obtained by noting that

(39)16E[rt4]E[σt4]=E[ω2+α2rt14+β2σt14+δ2rt12+2ωαrt12+2ωβσt12+2ωδrt1+2αβrt12σt12+2αδrt13+2βδrt1σt12]=ω2+(δ2+2ωα+2ωβ)σr2+α2E[rt14]+β2E[σt14]+2αβE[rt12σt12]=ω2+(2ωp+δ2)σr2+(α2+16β2+13αβ)E[rt4]

so that using eq. (7) leads to

(40)κr=E[rt4](σr2)2=6ω2+6(2ωp+δ2)σr216α2β22αβ(1p)2ω2=6(1p)2σr4+6(2p(1p)σr2+δ2)σr216α2β22αβσr4=6(1p)2+12p(1p)+6δ216α2β22αβ=6(1p2+δ2)1p25α2.

Also, since

(41)E[rt2rt+12]=E[rt2(ω+αrt2+βσt2+δrt)]=ωσr2+(α+16β)E[rt4]

then we have

(42)R2(1)=E[rt2rt+12]E[rt2]E[rt+12]Var[rt2]=(α+16β)E[rt4]pσr4E[rt4]σr4=(α+16β)p(1q)6(1p2+δ2)1(1q)6(1p2+δ2)=p5β(1p2+δ2)6(1p2+δ2)(1q)

and hence

(43)C2(2)=p(α+16β)E[rt4]p2σr4=pC2(1),

which leads to the recursive formula (10). Besides, eq. (12) is by induction with

(44)σˆt+12=ω+αE[rt2|Ft1]+βE[σt2|Ft1]+δE[rt|Ft1]=ω+ασt2+βσt2=σr2+p(σt2σr2)
(45)σˆt+22=ω+αE[rt+12|Ft1]+βE[σt+12|Ft1]+δE[rt+1|Ft1]=ω+ασˆt+12+βσˆt+12=σr2+p(σˆt+12σr2)=σr2+p2(σt2σr2).

In the NA-GARCH model under Laplace errors, eq. (14) is due to

(46)σr2E[σt2]=ω+αE[(rt1θσt1)2]+βE[σt12]=ω+(α+αθ2+β)σr2,

while eq. (15) is a result of

(47)16E[rt4]E[σt4]=E[ω2+α2(rt1θσt1)4+β2σt14+2ωα(rt1θσt1)2+2ωβσt12+2αβ(rt1θσt1)2σt12]=ω2+2ωpσr2+(α2+α2θ2+16α2θ4+16β2+13αβ(1+θ2))E[rt4]

and eq. (16) is from

(48)E[rt2rt+12]=E[rt2(ω+α(rtθσt)2+βσt2)]=ωσr2+(α+16αθ2+16β)E[rt4].

In the GJR-GARCH model under Laplace errors, eq. (20) results from

(49)σr2E[σt2]=ω+(α+ϕE[1{ϵt1<0}])E[rt12]+βE[σt12]=ω+(α+β+12ϕ)σr2

and eq. (21) is implied by

(50)16E[rt4]E[σt4]=E[ω2+α2rt14+β2σt14+ϕ21{ϵt1<0}rt14+2ωαrt12+2ωβσt12+2ωϕ1{ϵt1<0}rt12+2αβrt12σt12+2αϕ1{ϵt1<0}rt14+2βϕ1{ϵt1<0}rt12σt12]=ω2+2ωpσr2+(α2+16β2+12ϕ2+13αβ+αϕ+16βϕ)E[rt4].

Similarly, eq. (22) is a result of

(51)E[rt2rt+12]=ωσr2+(α+16β+12ϕ)E[rt4]

and thus

(52)R2(1)=(α+16β+12ϕ)E[rt4]pσr4E[rt4]σr4=p5β(1p2)6(1p2)(1q).

In the E-GARCH model under Laplace errors, eq. (26) is deduced by

(53)σr2E[σt2]=E[exp(ω+βlnσt12+g(ϵt1))]=eωE[eg(ϵt1)]E[exp(ωβ+β2lnσt22+βg(ϵt2))]=eω(1+β)E[eg(ϵt1)]E[eβg(ϵt2)]E[exp(ωβ2+β3lnσt32+β2g(ϵt3))]==exp(ωi=1βi1)i=1E[exp(βi1g(ϵti))]E[exp(limi(βilnσti2))]

so that we require 0β<1 for the limit factor to converge, and due to

(54)mg(k)E[ekg(ϵt)]=22e22kα(0e(k(ψα)+2)xdx+0e(k(ψ+α)2)xdx)

we have that for 0k<1, |ψ|<2α is needed for the MGF of g(ϵt) to converge. Proof of the convergence of the infinite product is given in the next appendix. Likewise, eq. (28) is derived as

(55)κrE[rt4]σr4=6σr4E[exp(2ω+2βlnσt12+2g(ϵt1))]==6σr4exp(2ωi=1βi1)i=1E[exp(2βi1g(ϵti))]E[exp(limi(2βilnσti2))]=6i=1mg(2βi1)(mg(βi1))2,

where |ψ|<2/2α is required to ensure that the MGF converges for 0k<2. Furthermore, for l1, since

(56)E[rt2rt+l2]=E[rt2exp(ω+βlnσt+l12+g(ϵt+l1))]==E[ϵt2exp(ωi=1lβi1+ω(1+βl)i=1βi1+limi((1+βl)βilnσti2)+i=1l(βi1g(ϵt+li))+i=1((1+βl)βi1g(ϵti)))]=exp(2ω1β)E[ϵt2exp(βl1g(ϵt))]i=1l1E[exp(βi1g(ϵt+li))]i=1E[exp((1+βl)βi1g(ϵti))],

we have the expression (29). The forecast eq. (33) is merely calculated from

(57)σˆt+l2=E[exp(ω+βlnσt+l12+g(ϵt+l1))|Ft1].

B Proof of Convergence of Infinite Product in eq. (26)

Although Nelson (1991) verified that convergent MGF of errors automatically leads to convergent unconditional variance in the E-GARCH model, we expressly give a proof in the case of Laplace errors. We attempt to prove that 0β<1 with |ψ|<2α is also sufficient for the infinite product in eq. (26) to converge. To do this, note that a sufficient and necessary condition for an infinite product i=1ai with positive terms to converge to a nonzero finite number is that the infinite series i=1lnai converges to a finite number.

Further, notice that i=1|lnai|< guarantees i=1lnai<. Also, it is necessary that limiai=1 for both i=1|lnai| and i=1|ai1| to converge. According to the comparison test of positive infinite series, if this condition is satisfied,

(58)limi|lnai||ai1|=1,

and then i=1|ai1| converges if and only if i=1|lnai| converges.

In Appendix A we have shown that, with k=βi1, for i=1,2,... and 0β<1, the MGF of g(ϵt) in the E-GARCH model under Laplace errors converges to

(59)mg(βi1)=2αβi12(ψ2α2)β2(i1)+22αβi12e22αβi1>0

provided that |ψ|<2α. Hence, in the convergence problem of i=1mg(βi1), it suffices to show that i=1|mg(βi1)1|<.

Clearly, for a new sequence (βi1),

(60)limi|mg(βi1)1|βi1=limi|2α2β1i(ψ2α2)β2(i1)+22αβi12e22αβi1β1i|=|2limiβ1i2α2limiβ1i|=|22α|<,

which says that the series i=1|mg(βi1)1| has the same convergence rate as i=1βi1. With 0β<1, i=1βi1=1/(1β)< leads to i=1|mg(βi1)1|< as well and therefore i=1mg(βi1)< as desired.

C Selected Moment Formulas under Normal or Student Errors

Recall that the error distribution does not affect the unconditional variance in each model but changes the kurtosis. The following expressions for kurtosis have been used in our analysis of constraints as well as empirical modeling for comparison. They directly follow from generalizing E[ϵt4]=κϵ in Section 3. Recall that normal errors have κϵ=3 while Student-5 errors and Student-7 errors have κϵ=9 and κϵ=5, respectively. Similar versions can also be found in Tsay (2005) .

Assume that p has the same value as that in the corresponding Laplace model. In the A-GARCH model, if p2+(κϵ1)α2<1, then

(61)κr=κϵ(1p2+δ2)1p2(κϵ1)α2;

in the NA-GARCH model, if p2+α2((κϵ1)+(κϵ2)θ2)<1, then

(62)κr=κϵ(1p2)1p2α2((κϵ1)+(κϵ2)θ2);

in the GJR-GARCH model, if p2+(κϵ1)α(α+ϕ)+(2κϵ1)/4ϕ2<1, then

(63)κr=κϵ(1p2)1p2(κϵ1)α(α+ϕ)2κϵ14ϕ2.

For the E-GARCH, it is useful to study the MGF of g(ϵt) in the presence of normal or Student errors. Using the fact that the standard normal density function is n(x)=1/2πex2/2 with xR and that E[|ϵt|]=2/π, the MGF with rate k0 is hence calculated as

(64)mg(k)=E[exp(k(ψϵt+α(|ϵt|2π)))]=e2πkα(012πek(ψα)xx22dx+012πek(ψ+α)xx22dx)=e2πkαe12k2(αψ)2k(αψ)ex222πdx+e2πkαe12k2(α+ψ)2k(α+ψ)ex222πdx=exp(k2(αψ)222πkα)N(k(αψ))+exp(k2(α+ψ)222πkα)N(k(α+ψ)),

where N() is the the standard normal cumulative distribution function. Indeed, under normal errors, mg(k)< as long as k<. On the other hand, when errors have a Student distribution, they have indeterminate MGF (see, e.g., Walck (2007)). As such, the MGF of g(ϵt) cannot exist unless |ψ|<α, in which case the exponential is nonincreasing in ϵt>0 for any given t. This constraint is unrealistically inappropriate by requiring that the returns impact level α must be negative. As illustrated in Section 5, such a constraint is hardly satisfied in practice[14].

References

Aho, K., D. Derryberry, and T. Peterson. 2014. “Model Selection for Ecologists: The Worldviews of AIC and BIC.” Ecology 95 (3): 631–636.10.1890/13-1452.1Search in Google Scholar

Bollerslev, T. 1986. “Generalized Autoregressive Conditional Heteroscedasticity.” Journal of Econometrics 31: 307–327.10.1016/0304-4076(86)90063-1Search in Google Scholar

Bollerslev, T. 1987. “A Conditional Heteroskedastic Time Series Model for Speculative Prices and Rates of Return.” The Review of Economics and Statistics 9 (3): 542–547.10.2307/1925546Search in Google Scholar

Box, G. E. P., and G. M. Jenkins. Holden-Day, 1976 Time Series Analysis: Forecasting and Control Rev Ed. San Francisco, CA, USA.Search in Google Scholar

Campbell, J. Y., and L. Hentschel. 1992. “No News is Good News: An Asymmetric Model of Changing Volatility in Stock Returns.” Journal of Financial Economics 31 (3): 281–318.10.3386/w3742Search in Google Scholar

Chen, C. W. S., R. Gerlach, and D. C. M. Wei. 2009. “Bayesian Causal Effects in Quantiles: Accounting for Heteroscedasticity.” Computational Statistics and Data Analysis 53 (6): 1993–2007.10.1016/j.csda.2008.12.014Search in Google Scholar

Christie, A. A. 1982. “The Stochastic Behavior of Common Stock Variances: Value, Leverage and Interest Rate Effects.” Journal of Financial Economics 10 (4): 407–432.10.1016/0304-405X(82)90018-6Search in Google Scholar

Christofferson, P. F. Elsevier, 2012 Elements of Financial Risk Management 2nd Ed. Waltham, MA, USA.Search in Google Scholar

Engle, R. F. 1982. “Autoregressive Conditional Heteroscedasticity with Estimates of the Variance of United Kingdom Inflation.” Econometrica 50 (4): 987–1008.10.2307/1912773Search in Google Scholar

Engle, R. F. 1990. “Discussion: Stock Market Volatility and the Crash of ‘87.” The Review of Financial Studies 3 (1): 103–106.10.1093/rfs/3.1.103Search in Google Scholar

Engle, R. F., and V. K. Ng. 1993. “Measuring and Testing the Impact of News on Volatility.” The Journal of Finance 48 (5): 1749–1778.10.3386/w3681Search in Google Scholar

French, K. R., G. W. Schwert, and R. F. Stambaugh. 1987. “Expected Stock Returns and Volatility.” Journal of Financial Economics 19: 3–29.10.1016/0304-405X(87)90026-2Search in Google Scholar

Gerlach, R., C. W. S. Chen, and N. Y. C. Chan. 2011. “Bayesian Time-Varying Quantile Forecasting for Value-At-Risk in Financial Markets.” Journal of Business and Economic Statistics 29 (4): 481–492.10.1198/jbes.2010.08203Search in Google Scholar

Glosten, L. R., R. Jaganathan, and D. E. Runkle. 1993. “Relationship Between the Expected Value and the Volatility of the Nominal Excess Return on Stocks.” The Journal of Finance 48 (5): 1779–1801.10.1111/j.1540-6261.1993.tb05128.xSearch in Google Scholar

McKenzie, C. R., and S. Takaoka. 2012. “Software Reviews: Eviews 7.2.” Journal of Applied Econometrics 27 (7): 1205–1210.10.1002/jae.2303Search in Google Scholar

Nelson, D. B. 1991. “Conditional Heteroskedasticity in Asset Returns: A New Approach.” Econometrica 59: 347–370.10.2307/2938260Search in Google Scholar

Tsay, R. S. Wiley-Interscience, 2005 Analysis of Financial Time Series 2nd Ed. NJ, USA.10.1002/0471746193Search in Google Scholar

Walck, C. Fysikum, University of Stockholm: Internal Report SUF-PFY/96–01, 2007 Handbook on Statistical Distributions for Experimentalists.Search in Google Scholar

Wichitaksorn, N., J. J. J. Wang, S. T. B. Choy, and R. Gerlach. 2015. “Analysing Return Asymmetry and Quantiles Through Stochastic Volatility Models Using Asymmetric Laplace Error via Uniform Scale Mixtures.” Applied Stochastic Models in Business and Industry 31 (5): 584–608.10.1002/asmb.2062Search in Google Scholar

Incorporation, Wolfram Research “Mathematica, version 10.3.” URL https://www.wolfram.com/mathematica. 2015.Search in Google Scholar

Published Online: 2017-7-18

© 2018 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 24.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/jtse-2016-0019/html
Scroll to top button