Skip to content
All library documents

Why GARCH Standardized Squared Returns Should Have Little Autocorrelation

Article Quant Q&A · Author: user12348

Summary

The document explains why squared returns divided by estimated conditional variance can be used to check a GARCH model. Under the model, returns are written as conditional volatility times an independent, zero-mean, unit-variance innovation. A GARCH(1,1) variance process gives squared returns an autoregressive moving-average structure, which accounts for serial dependence in their squares. Dividing squared returns by the conditional variance recovers squared innovations, which are independent under the assumptions.

This gives the diagnostic its rationale: little remaining autocorrelation in standardized squared residuals is evidence that the model captured volatility clustering. The conclusion depends on the return and variance assumptions being appropriate and on a sufficiently good variance estimate. It is not a guarantee that returns are fully independent or free of trends and other structure. The discussion also cautions that this diagnostic alone does not establish that residuals are identically distributed or that a model is adequate for every risk application.

Key ideas

  • A GARCH variance process can make squared returns serially dependent.
  • If returns equal conditional volatility multiplied by independent standardized innovations, dividing by volatility recovers those innovations.
  • The squared standardized residuals should have little autocorrelation when the model assumptions hold.
  • Remaining autocorrelation can indicate that the conditional variance model has not captured all volatility dependence.
  • This diagnostic alone does not establish that a model is otherwise adequate or that residuals are fully independent.

Tags

Full text
# Does GARCH derived variance explain the autocorrelation in a time series?


# Does GARCH derived variance explain the autocorrelation in a time series?












Given a time series $u_i$ of returns (where $i=1,\dotsc,t$), $\sigma_i$ is calculated from GARCH(1,1) as $$ \sigma_i^2=\omega+\alpha u_{i-1}^2 +\beta \sigma_{i-1}^2. $$ What is the mathematical basis to say that $u_i^2/\sigma_i^2$ will exhibit little autocorrelation in the series?

Hull's book "Options, Futures and Other Derivatives" is an excellent reference. In 6th ed. p. 470, "How Good is the Model?" he states that

> If a GARCH model is working well, it should remove the autocorrelation. We can test whether it has done so by considering the autocorrelation structure for the variables $u_i^2/\sigma_i^2$. If these show very little autocorrelation our model for $\sigma_i$ has succeeded in explaining autocorrelation in the $u_i^2$.

Maximum likelihood estimation for variance ends with maximizing $$ -m \space \ln(v) -\sum_{i=1}^{t} u_i^2/v_i $$ where $v_i$ is variance = $\sigma_i^2$. This function does not really mean $u_i^2/v_i$ being minimized, because $-\ln(v_i)$ gets larger and so does $u_i^2/v_i$ as $v_i$ gets smaller. However, it makes intuitive sense that dividing $u_t$ return by its (instant or regime) volatility explains away volatility-related component of the time series. I am looking for a mathematical or logical explanation of this.

I think Hull is not very accurate here as the time series may have trends etc.; also, there are better approaches to finding i.i.d. from the times series than using $u_i^2/\sigma_i^2$ alone. I particularly like Filtering Historical Simulation- Backtest Analysis by Barone-Adesi (2000).

## Answer by Malick (score 3)

https://quant.stackexchange.com/a/21933

What is the mathematical basis to say that $u^{2}_{t}/\sigma_{t}^{2}$ will exhibit little auto-correlation in the series?

Let's $r_{t}$ be a series of returns and let's assume (Assumption I) it follows a covariance stationary process defined as :

$r_{t}=\sigma_{t} z_{t}$

where $z_{t}$ is i.i.d with $E_{t}(z_{t})=0$ and $Var_{t}(z_{t})=1$ ;

> Then $ Var_{t}(r_{t}) =\sigma_{t}^{2}$

Next if we assume (Assumption II) that the conditional variance process of $r_{t}$ follows a GARCH(1,1), it means (Bollerslev (1986)) :

$\sigma_{t}^{2} = w + \alpha r_{t-1}^{2} + \beta \sigma_{t-1}^2 $

It can be shown that the precedent equation can be rewritten using an ARMA(1,1) representation:

$r_{t}^{2} = w + (\alpha +\beta) r_{t-1}^{2} + v_{t} - \beta v_{t-1} $

Where $ v_{t} = r_{t}^{2} - \sigma_{t}^{2}$

> without going into the math you see that $r_{t}^{2} $ has some autocorrelations since it has an autoregressive structure.

However

$r_{t}^{2}=\sigma_{t}^{2} z_{t}^{2}$

Then

$z_{t}^{2}=r_{t}^{2} /\sigma_{t}^{2}$

But we know that $z_{t}$ is IID(0,1),then its squares will also be IID.

> So if assumptions I and II are respected the series $z_{t}^{2}=r_{t}^{2} /\sigma_{t}^{2}$ (the standardized residuals) will be IID and will exhibit no autocorrelations at all. If the DGP is not perfectly in line with our assumptions (but still close enough), the series will exhibit little autocorrelation since $z_{t}$ will be almost IID.

## Answer by Lasha (score 1)

https://quant.stackexchange.com/a/19433

In the paragraph before the one from which you gave the quotation is written such a thing: “…when $u_i^2$ is high, there is a tendency for $u_{i+1}^2$, $u_{i+2}^2$, … to be high; when $u_i^2$ is low, there is a tendency for $u_{i+1}^2$, $u_{i+2}^2$, … to be low.” This means that they $u_{i+1}^2,u_{i+1}^2, u_{i+2}^2,…$ are correlated. Which by itself means that $u_i^2$ exhibits autocorrelation. If we predicted well the $σ_i^2-s$ (and $u_i^2$ is an approximation of it), after dividing $σ_i^2$ on $u_i^2$ we should no more have the pattern described in quotation. Which means we are “removing autocorrelation.”

Shown in full with attribution under the source's licence. Licence: CC BY-SA 4.0 (Stack Exchange)

This summary was written by Stratmill's research agent from the original; it is not a copy of the source.