The GJR model is a GARCH variant that includes leverage terms for modeling asymmetric volatility clustering. In the GJR formulation, large negative changes are more likely to be clustered than positive changes. The GJR model is named for Glosten, Jagannathan, and Runkle . Close similarities exist between the GJR model and the threshold GARCH (TGARCH) model—a GJR model is a recursive equation for the variance process, and a TGARCH is the same recursion applied to the standard deviation process.
The GJR(P,Q) model has P GARCH coefficients associated with lagged variances, Q ARCH coefficients associated with lagged squared innovations, and Q leverage coefficients associated with the square of negative lagged innovations. The form of the GJR(P,Q) model in Econometrics Toolbox™ is
The indicator function equals 1 if , and 0 otherwise. Thus, the leverage coefficients are applied to negative innovations, giving negative changes additional weight.
For stationarity and positivity, the GJR model has the following constraints:
The GARCH model is nested in the GJR model. If all leverage coefficients are zero, then the GJR model reduces to the GARCH model. This means you can test a GARCH model against a GJR model using the likelihood ratio test.
 Glosten, L. R., R. Jagannathan, and D. E. Runkle. "On the Relation between the Expected Value and the Volatility of the Nominal Excess Return on Stocks." The Journal of Finance. Vol. 48, No. 5, 1993, pp. 1779–1801.