V-Lab
V-Lab

Documentation>Volatility Analysis>Spline-GARCH

Definition

Consider a return time series rt=μ+εt, where μ is the expected return and εt is a zero-mean white noise. Despite of being serially uncorrelated, the series εt does not need to be serially independent. For instance, it can present conditional heteroskedasticity. The Spline-GARCH model assumes a specific parametric form for this conditional heteroskedasticity. More specifically, we say that εt~Spline-GARCH if we can write εt=σt2τtzt, where zt is standard Gaussian:

σt2=ω+αεt-12+βσt-12

is a usual GARCH specification, and:

τt=expi=1kγit-ti2

is the exponential of a quadratic Spline with kknots t1t2... tk. The idea of adding this τt is to capture lower frequency variations on the volatility, like seasonalities and trends. The number of knots, k, is chosen by minimizing the Bayesian Information Criterion (BIC), also known as Schwarz Information Criterion (SIC), while the knots are placed equidistantly.

Estimation

V-Lab estimates all the parameters μωαβ and γ1...γk simultaneously, by maximizing the log likelihood. The assumption that zt is Gaussian does not imply the the returns are Gaussian. Even though their conditional distribution is Gaussian, it can be proved that their unconditional distribution presents excess kurtosis (fat tails). In fact, assuming that the conditional distribution is Gaussian is not as restrictive as it seems: even if the true distribution is different, the so-called Quasi-Maximum Likelihood (QML) estimator is still consistent, under fairly mild regularity conditions.

Besides leptokurtic returns, the Spline-GARCH model, as the GARCH model, captures other stylized facts in financial time series, like volatility clustering. The volatility is more likely to be high at time t if it was also high at time t-1. Another way of seeing this is noting that a shock at time t-1 also impacts the variance at time t.

Usual restrictions on the parameters are ωαβ>0 . There are no restrictions on the parameters of the spline, γ1...γk . The GARCH model is in fact a restricted version of the Spline-GARCH, with γ1=...=γk=0.

Prediction

Let rT be the last observation in the sample, and let ω^, α^, β^ and γ^1... γ^kbe the QML estimators of the parameters ω, α, β and γ1...γ k, respectively.

For out-of-sample forecasts, we need to extrapolate the spline. We have good results when the extrapolation simply uses the last in-sample value. That is, V- Lab assumes that, h0:

τ^t+h=τ^T=expi=1kγ^iT-ti2

Hence, the Spline-GARCH model implies that the forecast of the conditional variance at time T+h is:

Hence, the Spline-GARCH model implies that the forecast of the conditional variance at time T+h is:

σ^T+h2τ^T+h=ω^+α^+β^σ^T+h-12expi=1kγ^iT-ti2

and so, by applying the above formula iteratively, we can forecast the conditional variance for any horizon h. Then, the forecast of the compound volatility at time T+h is:

σ^T+1:T+hτ^T+1:T+h=i=1hσ^T+i2τ^T+i

Notice that, for large h, this forecast of the compound volatility converges to:

hω^1-α^-β^expi=1kγ^iT-ti2

scaling over the forecast horizon with the well known square-root law.

Zero Slope Spline-GARCH

The Zero Slope Spline-GARCH model requires that the low- frequency component (i.e. the exponential of the spline) has zero slope in the end of the sample. That is, the coefficients are estimated by QML with the additional restriction that:

2i=1kγiT-tiexpi=1kγiT-ti2=0

i=1kγiT-ti=0

Note that we can solve the last equation for, say, γ1. That is, from the zero-slope constraint, we can write γ1 as a function of γ2... γk. But then we can substitute this expression wherever γ1 appears in the maximum likelihood estimation. In this way, we can transform a constrained optimization in an unconstrained optimization with one reduced dimension. This leads to a more reliable, fast optimization.

Spline-GARCH(p,q)

The specific model just described can be generalized to account for more lags in the conditional variance. A Spline-GARCHpq model assumes that:

σt2=ω+i=1pαiεt-i2+j=1qβjσt-j2

The best model (p and q) can be chosen, for instance, by Bayesian Information Criterion (BIC), or by Akaike Information Criterion (AIC). The former tends to be more parsimonious than the latter. V-Lab uses p=1 and q=1 though, because this is usually the option that best fits financial time series.

References

Engle, R. F. and J. G. Rangel, 2008. The Spline-GARCH Model for Low-Frequency Volatility and Its Global Macroeconomic Causes. Review of Financial Studies 21(3): 1187-1222. https://www.jstor.org/stable/2961959