Ensuring Positivity in GJR-GARCH Models with Exogenous Variables
Summary
This note explains why a negative coefficient on an exogenous variable in a GJR-GARCH variance equation does not by itself prove that the conditional variance becomes negative. The sign of the full variance recursion matters: standard positivity restrictions on the intercept and lag terms must be combined with a bound on the exogenous contribution. In the formulation described, requiring the intercept plus the exogenous term to remain nonnegative at each observation provides a sufficient condition, alongside the usual GJR-GARCH restrictions.
The discussion also connects that condition to nonnegative unconditional variance under covariance stationarity, and notes that constraints depend on the behavior and specification of the exogenous variable. For a strictly positive regressor, a nonnegative coefficient is one common choice. The note offers parameter restrictions and model-design guidance, but does not provide a general procedure for choosing the regressor model or establish that plotting in-sample variance alone is sufficient.
Key ideas
- A negative exogenous-variable coefficient does not alone determine whether conditional variance is negative.
- GJR-GARCH positivity also depends on the intercept, ARCH and GARCH terms, and asymmetric contribution.
- A sufficient restriction is to keep the intercept plus the exogenous contribution nonnegative at every time step.
- The suitable coefficient restrictions depend on the range and dynamics assumed for the exogenous variable.
Tags
Full text
# Volatility Modelling negative GJR-GARCH-X coefficient
# Volatility Modelling negative GJR-GARCH-X coefficient
I have estimated GARCH and GJR-GARCH with several exogenous variables. Some of the exogenous variables have negative coefficients that are statistically significant. For instance, I can write my GJR-GARCH estimate as:
$h_t = 0.213 + 0.011 u_{t-1}^2 + 0.847 h_{t-1} + 0.196 u^2_{t-1}I_{t-1<0} - 0.026 X_{t-1} $
Where $u_{t-1}^2$ is the lagged ARCH term, $I_t$ is the dummy variable that models the asymmetric leverage effect and $X_{t-1}$ is the exogenous variable.
My main concern is that the coefficient of $X_{t-1}$ is negative. How can I verify that my volatility is not negative? Is just a plot of the conditional variance over the in-sample period is enough?
I have tried to estimate the GJR-GARCH-X model, however, I obtained some negative coefficients. I am not sure how to verify that these negative coefficients will cause the volatility to be negative.
## Answer by Pleb (score 2)
https://quant.stackexchange.com/a/72166
#### Positivity of GARCH-X models:
When it comes to GARCH models with exogenous regressors, it is more "tricky" to ensure positivity of the model, due to the unspecified model-dynamics of your exogenous regressor (moreover exogenous regressors can have drastically different characteristics). Without any model-specification on the exogenous regressor $X_{t-1}$, it is common to restrict the parameter-space of the model in order to ensure positivity.
For simplicity, let us vaguely define the GJR-GARCH(1,1)-X model with demeaned returns, $r_t$: \begin{align*} r_t \vert \mathcal{F}_{t-1} &= \varepsilon_t\\ \varepsilon_t &= \sigma_t \cdot z_t\\ \sigma^2_t &= \omega + \alpha \varepsilon_{t-1}^2 + \beta \sigma_{t-1}^2 + \gamma I_{t-1} \varepsilon_{t-1}^2 + \nu X_{t-1}, \end{align*} where $z_t \overset{iid}{\sim} D(0,1)$ is a standardized distribution and
$$I_{t-1} =\begin{cases} 1 & \text{if } \varepsilon_{t-1} < 0 \\ 0 & \text{if } \varepsilon_{t-1} \geq 0 \end{cases}.$$
I have detailed some of my observations:
- When working under the GJR-GARCH(1,1) positivity is satisfied when we impose $\omega, \beta,\alpha > 0$ and $\alpha + \gamma > 0$. The latter condition is a broader statement than imposing $\alpha, \gamma >0$, since we can allow one of the parameters to become negative (in this case, $\gamma$). I have made an in-depth answer detailing the GJR-GARCH(1,1) model, positivity, covariance stationarity and economical interpretations of the parameter estimates.
- When working under the GJR-GARCH(1,1)-X we can further ensure positivity by additionally restricting the exogenous regressor such that $\omega + \nu X_{t-1} \geq 0$. Here, we allow $\nu$ to vary freely (since $\omega > 0$) as long as the above condition is satisfied. The primary motivation for the extra restriction, comes from observing the unconditional variance (calculated under assumed covariance stationarity) of the GJR-GARCH-X model:
\begin{equation} \mathbb{V}ar(r_t) = \frac{\omega + \nu \mathbb{E}\left[X_{t-1}\right]}{1 - \alpha - \beta - \kappa \gamma}. \end{equation} In order to ensure non-negative unconditional variance of the return process, we specifically need $\omega + \nu \mathbb{E}\left[X_{t-1}\right] \geq 0$, which is satisfied when imposing $\omega + \nu X_{t-1} \geq 0$ for all $t$.
In conclusion, imposing $\omega, \beta, \alpha > 0$, $\alpha + \gamma > 0$ and $\omega + \nu X_{t-1}>0$, ensures that you obtain non-negative volatility estimates.
When $X_{t-1}$ is strictly positive, it is common to let $\nu \geq 0$. This is also emphasized in the article of Han, H. (2015) that investigates asymptotic results of the GARCH-X model when $X_{t-1}$ follows a fractionally integrated process. In general, it is common in academia to assume a functional form on the exogenous regressor, (see for instance (C1 - C4) in this paper p.699). This is also done in the Realized GARCH model that incorporates intraday data to procure better forecasts (I have detailed this model here, if you're interested). I hope this helps.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.