Volatility, Tails, and Financial Risk

Conditional variance, GARCH quasi-likelihood, leverage, heavy tails, and translation into VaR and expected shortfall.

Volatility, Tails, and Financial Risk

Why finance changes the second moment

Daily returns often have weak linear autocorrelation but strong dependence in squared or absolute returns. A useful decomposition is

rt=μt+εt,εt=htzt,E(ztFt1)=0,E(zt2Ft1)=1.r_t=\mu_t+\varepsilon_t,\qquad \varepsilon_t=\sqrt{h_t}z_t,\qquad E(z_t\mid\mathcal F_{t-1})=0,\quad E(z_t^2\mid\mathcal F_{t-1})=1.

The conditional mean μt\mu_t and conditional variance hth_t are separate modelling problems. An unconditional ARMA model can describe E(rtFt1)E(r_t\mid\mathcal F_{t-1}) while missing the risk clustering in hth_t.

From ARCH to GARCH

Engle's ARCH model lets recent squared shocks drive variance. Bollerslev's GARCH model adds lagged variance:

ht=ω+αεt12+βht1,ω>0,α,β0.h_t=\omega+\alpha\varepsilon_{t-1}^2+\beta h_{t-1}, \qquad \omega>0,\quad\alpha,\beta\ge0.

If α+β<1\alpha+\beta<1, the unconditional variance is

E(ht)=ω1αβ.E(h_t)=\frac{\omega}{1-\alpha-\beta}.

Interpretation:

  • α\alpha is the immediate variance response to news;
  • β\beta is persistence in the variance state;
  • α+β\alpha+\beta controls the decay of a volatility shock;
  • ω\omega sets the long-run scale jointly with persistence.

When α+β\alpha+\beta is close to one, volatility forecasts revert slowly. That is not the same as a unit root in the return level.

Conditional covariance as a matrix

Given past information, a block of future shocks has covariance

Ht:t+ht=diag(ht+1t,,ht+ht)H_{t:t+h\mid t} =\operatorname{diag}(h_{t+1\mid t},\ldots,h_{t+h\mid t})

only when future standardised innovations are conditionally uncorrelated. In a multivariate model, HtH_t is a full positive-definite matrix and must simultaneously control variances and correlations.

The one-step GARCH forecast is

ht+1t=ω+αεt2+βht.h_{t+1\mid t}=\omega+\alpha\varepsilon_t^2+\beta h_t.

For k2k\ge2,

ht+kt=ω+(α+β)ht+k1t.h_{t+k\mid t}=\omega+(\alpha+\beta)h_{t+k-1\mid t}.

Gaussian quasi-likelihood

Ignoring constants, the conditional Gaussian objective is

(θ)=12t=1T[loght(θ)+εt2ht(θ)].\ell(\theta)= \frac12\sum_{t=1}^{T} \left[\log h_t(\theta)+\frac{\varepsilon_t^2}{h_t(\theta)}\right].

It is a quasi-likelihood when ztz_t is not Gaussian. Under suitable conditions it can still estimate the variance recursion consistently, but Gaussian tail probabilities and conventional standard errors need not be valid.

Always inspect standardised residuals

z^t=ε^t/h^t.\hat z_t=\hat\varepsilon_t/\sqrt{\hat h_t}.

Remaining dependence in z^t2\hat z_t^2 signals an inadequate variance recursion; heavy tails in z^t\hat z_t signal an inadequate innovation distribution.

R laboratory: estimate GARCH by recursion

R

Simulate and estimate a GARCH(1,1) with base R

Idle

The constraints are built into the parameter transformation. This is safer than accepting negative variances after optimisation. The simulated returns are in percentage-like units, so the reported VaR and ES use the same units.

Heavy tails and leverage

Two extensions answer different diagnostics:

  • Student-tt innovations: large standardised shocks occur more often than under a Normal distribution.
  • Asymmetric variance dynamics: negative and positive shocks of equal magnitude can affect future volatility differently.

For example, a GJR-style recursion is

ht=ω+αεt12+γI(εt1<0)εt12+βht1.h_t=\omega+\alpha\varepsilon_{t-1}^2+ \gamma I(\varepsilon_{t-1}<0)\varepsilon_{t-1}^2+\beta h_{t-1}.

γ>0\gamma>0 represents a larger variance response after negative news. Do not use asymmetry and heavy tails as interchangeable explanations.

From volatility to a risk decision

For lower-tail probability pp, define the return quantile qp,tq_{p,t} by

Pr(rtqp,tFt1)=p.\Pr(r_t\le q_{p,t}\mid\mathcal F_{t-1})=p.

If losses are Lt=rtL_t=-r_t, then VaR1p,t=qp,t\operatorname{VaR}_{1-p,t}=-q_{p,t}. Expected shortfall is the expected loss beyond VaR. State the sign convention; otherwise a “negative VaR” can mean either a gain or a loss.

A volatility model is not validated merely because its likelihood is high:

LayerDiagnostic
variance dynamicsdependence remaining in z^t2\hat z_t^2
distributionquantile and tail calibration
decisionVaR exceptions, expected-shortfall loss, capital or utility impact
stabilityrolling performance through calm and stressed periods
Normal GARCH in the example teaches recursion and likelihood. It is not automatically suitable for regulatory capital, intraday microstructure, options, or a portfolio with changing weights.

Finance versus economics

GARCH is central in finance because the conditional scale of returns is itself a decision input. In macroeconomics it is useful only when time-varying uncertainty is part of the question—for example inflation uncertainty or exchange-rate volatility. For GDP nowcasting, publication timing and latent-state estimation are usually more central than a daily GARCH recursion.

Practice

  1. If α+β=0.98\alpha+\beta=0.98, what fraction of a variance deviation remains after ten forecast steps?
  2. Why can uncorrelated returns still be dependent?
  3. A model has correct average VaR coverage but all exceptions occur in one crisis month. Is it calibrated?
Answers
  1. Approximately 0.98100.8170.98^{10}\approx0.817; mean reversion is slow.
  2. Their conditional variances can be dependent, so squared returns cluster.
  3. Not dynamically: unconditional coverage can hide clustered exceptions.

Next: Unit Roots, Cointegration, and Error Correction.

Copyright © 2026