Skip to content
All library documents

COGARCH Volatility Dynamics and Its Stochastic Differential Form

Article Quant Q&A · Author: KiNest

Summary

This note traces the construction of the continuous-time GARCH model from a discrete GARCH(1,1) recursion. It replaces discrete shocks with increments of a Lévy process, defines an auxiliary spectrally negative process, and expresses volatility using that process. The return process is then driven by the volatility and the same Lévy noise. The note also compares this representation with a stochastic differential equation form for COGARCH volatility.

Its main contribution is to identify a point of confusion: how the volatility representation from the original paper relates to the alternative differential equation. It reports that the original paper gives a proposition connecting its volatility formula to a stochastic differential equation, but the author remains unsure how that matches the alternative form. The document does not resolve the discrepancy, and some displayed formulas or parameter conditions appear inconsistent or mistyped. It is therefore a useful model-definition reference, but not a complete derivation or validation of equivalence.

Key ideas

  • COGARCH extends the discrete GARCH(1,1) volatility recursion to continuous time using Lévy process increments.
  • An auxiliary process provides a compact representation of the volatility path.
  • The return process is defined by integrating volatility against the Lévy process.
  • The note compares two volatility formulations but leaves their relationship unresolved.

Tags

Full text
# COGARCH (continuous GARCH) model definition


# COGARCH (continuous GARCH) model definition












In 2004 in their paper Kluppelberg, Lindner and Maller proposed a continuous vesrion of GARCH(1,1) model (COGARCH). Below is the derivation of model formulae as presented in the original paper:

- Suppose standard GARCH(1,1) model: $$ Y_n=\epsilon_n\sigma_n, \\ \sigma^2_n=\beta+\lambda Y^2_{n-1}+\delta\sigma^2_{n-1}, \\ n\in \mathbb{N} $$

- Take $\delta\ge0$ and $\lambda\ge0$. From previous we obtain: $$ \sigma^2_n=\beta+\lambda Y^2_{n-1}+\delta\sigma^2_{n-1}=\beta+(\delta+\lambda\epsilon^2_{n-1})\sigma^2_{n-1} $$ where $\epsilon_{n-1}$ and $\sigma^2_{n-1}$ are independent.



- The summation in previous can be written as: $$ \beta\sum^{n-1}_{i=0}\prod^{n-1}_{j=i+1}(\delta +\lambda\epsilon^2_j) = \beta\int_0^n\exp\Biggl(\sum^{n-1}_{j=\lfloor s \rfloor+1}\ln(\delta +\lambda\epsilon^2_j)\Biggr)ds $$

- Now switching to continuous time. In order to preserve single source of randomness as in discrete GARCH(1,1) model authors suggest to replace the noise variable $\epsilon_j$ by increments of Levy process $\Delta L=L_t-L_{t-}, \ t\ge0$, where $L$ is (càdlàg) Levy process.

- Now keeping $0<\delta<0$ and $\lambda\ge0$ and with (4) in mind define an auxiliary (càdlàg) process $(X_t)_{t\ge0}$ by: $$ X_t=-t\ln\delta-\sum_{0<s\le t}(\ln(1+(\lambda/\delta)\Delta L_s^2)), \ t\ge0 $$

- This $(X_t)_{t\ge0}$ process will help us to rewrite volatility process from (3) in more compact way as left-continuous process: $$ \sigma^2_t=\Biggl(\beta\int_0^te^{X_s}ds+\sigma_0^2\Biggr)e^{-X_{t-}}, \ t\ge0 $$ $(X_t)_{t\ge0}$ process here is a spectrally negative Levy process of bounded variation with drift $\gamma_{X,0}=-\ln\delta$, gaussian component $\tau^2_X=0$ and Levy measure $\Pi_X((-\infty,-x])=\Pi_L(\{y\in\mathbb{R}:|y|\ge\sqrt{(e^{x}-1)\delta/\lambda}\})$.

- Finally define the intergrated continuous time GARCH (COGARCH) process $(G_t)_{t\ge0}$ as the càdlàg process satisfying: $$ dG_t=\sigma_tdL_t, \ t\ge0, \ G_0=0 $$

Last three equations in above derivation define the COGARCH model (as presented in the paper of Kluppelberg, Lindner and Maller). However, if we take a look at the corresponding Wikipedia page, we will see a slightly different definition of COGARCH model using stochastic differential equations: $$ dG_t=\sigma_{t-}dL_t \\ d\sigma^2_t=(\beta-\lambda\sigma^2_t)dt+\delta\sigma^2_{t-}d[L,L]^d_t $$ where: $$ [L,L]^d_t=\sum_{s\in[0,t]}(\Delta L_t^2) $$ is the purely discontinuous part of the quadratic variation process of $L$.

I have no problem understanding formulas from original paper of Kluppelberg, Lindner and Maller, however I fail to understand how volatility process $\sigma^2_t$ in wikipedia version is derived. I believe I'm missing some basic detail here but I still can't figure it out. Any clarification will be appreciated.

UPDATE: In the original paper I've found the proposition (with proof) that volatility process from (7) above satisfies the stochastic differential equation: $$ d\sigma^2_{t+}=\beta dt+\sigma^2_te^{X_{t-}}d(e^{-X_t}) $$

and we have: $$ d\sigma^2_{t+}=\beta t+\ln\delta\int_0^t\sigma^2_sds+(\lambda/\delta)\sum_{0<s<t}\sigma^2_s(\Delta _sL)^2+\sigma^2_0 $$

This seems similar to wikipedia version of volatility process but still not the same.

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.