EWMA, GARCH Volatilities

Let’s see the few of the volatility models based on past innovations (meaning past observations). There are few models such as moving average which gives equal weights to past observations but we will see most commonly used models in the industry, one such is EWMA meaning, Exponential Weighted Moving Average.

The general formula is first to calculate the weights,

p = \frac{1 - \lambda}{\lambda(1 - \lambda^M)} \lambda^m

where, m = 1, 2.....M

Then applying to predicting volatility,

\sigma_{t+1}  = \sqrt{p^m R^2_{t-m+1}}

Supposing, M is large, we can formulate more generic solution, termed RiskMetrics formerly modelled by JP Morgan using the lambda values close to 1 for volatility persistence and market reactivity.

\sigma_{t+1} = \lambda \sigma_t + (1-\lambda) R_t

Let’s now see for the another volatility model, most widely used and industry adopted solution, GARCH, Generalised AutoRegressive (meaning current return based on past lagged observations) Conditional heteroskedasticity (meaning non constant / non-homoskedastic variance across observations).

The formula to model for GARCH volatility is given below,

\sigma_{t+1} = \omega + \alpha R^2_t + \beta \sigma_t

where, \hat{\sigma}^2 = \frac{\omega}{1 - (\alpha + \beta)}

Estimating GARCH

Estimating GARCH parameters is done in 2 step process,

Estimate the long run volatility,

\hat{\sigma}^2 = \frac{1}{M} \sum R^2_m

Then, estimate GARCH,

E[R_{t+1}^2 | I_t] = \sigma_{t+1}^2 = \lambda \hat{\sigma}^2 + \alpha \sigma^2_t + \beta R^2_t

Now the parameters can be estimated using Maximum Likelihood estimation, given the past innovations are i.i.d.

GARCH Forecasting

With GARCH, you can forecast the volatility using generalised formula,

E[\sigma^2_{t+k}|I_t] = \hat{\sigma^2} +(\alpha + \beta)^{k-1} \sigma^2_{t+1}

Where, \hat{\sigma^2} = \frac{\omega}{1 - (\alpha + \beta)}

Leave a comment