# Autoregressive conditional heteroskedasticity

In econometrics, the **autoregressive conditional heteroskedasticity** (**ARCH**) model is a statistical model for time series data that describes the variance of the current error term or innovation as a function of the actual sizes of the previous time periods' error terms;^{[1]} often the variance is related to the squares of the previous innovations. The ARCH model is appropriate when the error variance in a time series follows an autoregressive (AR) model; if an autoregressive moving average model (ARMA) model is assumed for the error variance, the model is a **generalized autoregressive conditional heteroskedasticity** (**GARCH**) model.^{[2]} For forecasting, combining ARIMA and ARCH models could be considered. For instance, a hybrid ARIMA-ARCH model was examined for shipping freight rate forecast.^{[3]}

ARCH models are commonly employed in modeling financial time series that exhibit time-varying volatility and volatility clustering, i.e. periods of swings interspersed with periods of relative calm. ARCH-type models are sometimes considered to be in the family of stochastic volatility models, although this is strictly incorrect since at time *t* the volatility is completely pre-determined (deterministic) given previous values.^{[4]}

## ARCH(*q*) model specification

To model a time series using an ARCH process, let denote the error terms (return residuals, with respect to a mean process), i.e. the series terms. These are split into a stochastic piece and a time-dependent standard deviation characterizing the typical size of the terms so that

The random variable is a strong white noise process. The series is modeled by

- ,
- where and .

An ARCH(*q*) model can be estimated using ordinary least squares. A methodology to test for the lag length of ARCH errors using the Lagrange multiplier test was proposed by Engle (1982). This procedure is as follows:

- Estimate the best fitting autoregressive model AR(
*q*) . - Obtain the squares of the error
and regress them on a constant and
*q*lagged values:- where
*q*is the length of ARCH lags.

- The null hypothesis is that, in the absence of ARCH components, we have
for all
. The alternative hypothesis is that, in the presence of ARCH components, at least one of the estimated
coefficients must be significant. In a sample of
*T*residuals under the null hypothesis of no ARCH errors, the test statistic*T'R²*follows distribution with*q*degrees of freedom, where is the number of equations in the model which fits the residuals vs the lags (i.e. ). If*T'R²*is greater than the Chi-square table value, we*reject*the null hypothesis and conclude there is an ARCH effect in the ARMA model. If*T'R²*is smaller than the Chi-square table value, we do not reject the null hypothesis.

## GARCH

If an autoregressive moving average model (ARMA) model is assumed for the error variance, the model is a generalized autoregressive conditional heteroskedasticity (GARCH) model.^{[2]}

In that case, the GARCH (*p*, *q*) model (where *p* is the order of the GARCH terms
and *q* is the order of the ARCH terms
), following the notation of the original paper, is given by

Generally, when testing for heteroskedasticity in econometric models, the best test is the White test. However, when dealing with time series data, this means to test for ARCH and GARCH errors.

Exponentially weighted moving average (EWMA) is an alternative model in a separate class of exponential smoothing models. As an alternative to GARCH modelling it has some attractive properties such as a greater weight upon more recent observations, but also drawbacks such as an arbitrary decay factor that introduces subjectivity into the estimation.

### GARCH(*p*, *q*) model specification

The lag length *p* of a GARCH(*p*, *q*) process is established in three steps:

- Estimate the best fitting AR(
*q*) model- .

- Compute and plot the autocorrelations of
by
- The asymptotic, that is for large samples, standard deviation of
is
. Individual values that are larger than this indicate GARCH errors. To estimate the total number of lags, use the Ljung-Box test until the value of these are less than, say, 10% significant. The Ljung-Box Q-statistic follows
distribution with
*n*degrees of freedom if the squared residuals are uncorrelated. It is recommended to consider up to T/4 values of*n*. The null hypothesis states that there are no ARCH or GARCH errors. Rejecting the null thus means that such errors exist in the conditional variance.

### NGARCH

#### NAGARCH

**Nonlinear Asymmetric GARCH(1,1)** (**NAGARCH**) is a model with the specification:^{[7]}^{[8]}

- ,
- where and , which ensures the non-negativity and stationarity of the variance process.

For stock returns, parameter
is usually estimated to be positive; in this case, it reflects a phenomenon commonly referred to as the "leverage effect", signifying that negative returns increase future volatility by a larger amount than positive returns of the same magnitude.^{[7]}^{[8]}

This model should not be confused with the NARCH model, together with the NGARCH extension, introduced by Higgins and Bera in 1992.^{[9]}

### IGARCH

Integrated Generalized Autoregressive Conditional heteroskedasticity (IGARCH) is a restricted version of the GARCH model, where the persistent parameters sum up to one, and imports a unit root in the GARCH process. The condition for this is

.

### EGARCH

The exponential generalized autoregressive conditional heteroskedastic (EGARCH) model by Nelson & Cao (1991) is another form of the GARCH model. Formally, an EGARCH(p,q):

where
,
is the conditional variance,
,
,
,
and
are coefficients.
may be a standard normal variable or come from a generalized error distribution. The formulation for
allows the sign and the magnitude of
to have separate effects on the volatility. This is particularly useful in an asset pricing context.^{[10]}^{[11]}

Since may be negative, there are no sign restrictions for the parameters.

### GARCH-M

The GARCH-in-mean (GARCH-M) model adds a heteroskedasticity term into the mean equation. It has the specification:

The residual is defined as:

### QGARCH

The Quadratic GARCH (QGARCH) model by Sentana (1995) is used to model asymmetric effects of positive and negative shocks.

In the example of a GARCH(1,1) model, the residual process is

where is i.i.d. and

### GJR-GARCH

Similar to QGARCH, The Glosten-Jagannathan-Runkle GARCH (GJR-GARCH) model by Glosten, Jagannathan and Runkle (1993) also models asymmetry in the ARCH process. The suggestion is to model where is i.i.d., and

where if , and if .

### TGARCH model

The Threshold GARCH (TGARCH) model by Zakoian (1994) is similar to GJR GARCH. The specification is one on conditional standard deviation instead of conditional variance:

where if , and if . Likewise, if , and if .

### fGARCH

Hentschel's **fGARCH** model,^{[12]} also known as **Family GARCH**, is an omnibus model that nests a variety of other popular symmetric and asymmetric GARCH models including APARCH, GJR, AVGARCH, NGARCH, etc.

### COGARCH

In 2004, Claudia Klüppelberg, Alexander Lindner and Ross Maller proposed a continuous-time generalization of the discrete-time GARCH(1,1) process. The idea is to start with the GARCH(1,1) model equations

and then to replace the strong white noise process by the infinitesimal increments of a Lévy process , and the squared noise process by the increments , where

is the purely discontinuous part of the quadratic variation process of . The result is the following system of stochastic differential equations:

where the positive parameters
,
and
are determined by
,
and
. Now given some initial condition
, the system above has a pathwise unique solution
which is then called the continuous-time GARCH (**COGARCH**) model.^{[13]}

### ZD-GARCH

Unlike GARCH model, the Zero-Drift GARCH (ZD-GARCH) model by Li, Zhang, Zhu and Ling (2018) ^{[14]} lets the drift term
in the first order GARCH model. The ZD-GARCH model is to model
, where
is i.i.d., and

The ZD-GARCH model does not require , and hence it nests the Exponentially weighted moving average (EWMA) model in "RiskMetrics". Since the drift term , the ZD-GARCH model is always non-stationary, and its statistical inference methods are quite different from those for the classical GARCH model. Based on the historical data, the parameters and can be estimated by the generalized QMLE method.

## References

- ↑ Engle, Robert F. (1982). "Autoregressive Conditional Heteroscedasticity with Estimates of the Variance of United Kingdom Inflation".
*Econometrica*.**50**(4): 987–1007. JSTOR 1912773. - 1 2 Bollerslev, Tim (1986). "Generalized Autoregressive Conditional Heteroskedasticity".
*Journal of Econometrics*.**31**(3): 307–327. doi:10.1016/0304-4076(86)90063-1. - ↑ Munim, Ziaul Haque; Schramm, Hans-Joachim (2017). "Forecasting container shipping freight rates for the Far East – Northern Europe trade lane".
*Maritime Economics & Logistics*.**19**(1): 106–125. doi:10.1057/s41278-016-0051-7. - ↑ Brooks, Chris (2014).
*Introductory Econometrics for Finance*(3rd ed.). Cambridge: Cambridge University Press. p. 461. ISBN 9781107661455. - ↑ Lanne, Markku; Saikkonen, Pentti (July 2005). "Non-linear GARCH models for highly persistent volatility".
*The Econometrics Journal*.**8**(2): 251–276. doi:10.1111/j.1368-423X.2005.00163.x. JSTOR 23113641. - ↑ Bollerslev, Tim; Russell, Jeffrey; Watson, Mark (May 2010). "Chapter 8: Glossary to ARCH (GARCH)" (PDF).
*Volatility and Time Series Econometrics: Essays in Honor of Robert Engle*(1st ed.). Oxford: Oxford University Press. pp. 137–163. ISBN 9780199549498. Retrieved 27 October 2017. - 1 2 Engle, Robert F.; Ng, Victor K. (1993). "Measuring and testing the impact of news on volatility" (PDF).
*Journal of Finance*.**48**(5): 1749–1778. doi:10.1111/j.1540-6261.1993.tb05127.x. SSRN 262096.It is not yet clear in the finance literature that the asymmetric properties of variances are due to changing leverage. The name "leverage effect" is used simply because it is popular among researchers when referring to such a phenomenon.

- 1 2 Posedel, Petra (2006). "Analysis Of The Exchange Rate And Pricing Foreign Currency Options On The Croatian Market: The Ngarch Model As An Alternative To The Black Scholes Model" (PDF).
*Financial Theory and Practice*.**30**(4): 347–368.Special attention to the model is given by the parameter of asymmetry [theta (θ)] which describes the correlation between returns and variance.

^{6}...^{6}In the case of analyzing stock returns, the positive value of [theta] reflects the empirically well known leverage effect indicating that a downward movement in the price of a stock causes more of an increase in variance more than a same value downward movement in the price of a stock, meaning that returns and variance are negatively correlated - ↑ Higgins, M.L; Bera, A.K (1992). "A Class of Nonlinear Arch Models".
*International Economic Review*.**33**(1): 137–158. JSTOR 2526988. - ↑ St. Pierre, Eilleen F. (1998). "Estimating EGARCH-M Models: Science or Art".
*The Quarterly Review of Economics and Finance*.**38**(2): 167–180. doi:10.1016/S1062-9769(99)80110-0. - ↑ Chatterjee, Swarn; Hubble, Amy (2016). "Day-Of-The-Week Effect In Us Biotechnology Stocks—Do Policy Changes And Economic Cycles Matter?".
*Annals of Financial Economics*.**11**(2): 1–17. doi:10.1142/S2010495216500081. - ↑ Hentschel, Ludger (1995). "All in the family Nesting symmetric and asymmetric GARCH models".
*Journal of Financial Economics*.**39**(1): 71–104. doi:10.1016/0304-405X(94)00821-H. - ↑ Klüppelberg, C.; Lindner, A.; Maller, R. (2004). "A continuous-time GARCH process driven by a Lévy process: stationarity and second-order behaviour".
*Journal of Applied Probability*.**41**(3): 601–622. doi:10.1239/jap/1091543413. - ↑ Li, D.; Zhang, X.; Zhu, K.; Ling, S. (2018). "The ZD-GARCH model: A new way to study heteroscedasticity".
*Journal of Econometrics*.**202**(1): 1–17. doi:10.1016/j.jeconom.2017.09.003.

## Further reading

- Bollerslev, Tim; Russell, Jeffrey; Watson, Mark (May 2010). "Chapter 8: Glossary to ARCH (GARCH)" (PDF).
*Volatility and Time Series Econometrics: Essays in Honor of Robert Engle*(1st ed.). Oxford: Oxford University Press. pp. 137–163. ISBN 9780199549498. - Enders, W. (2004). "Modelling Volatility".
*Applied Econometrics Time Series*(Second ed.). John-Wiley & Sons. pp. 108–155. ISBN 0-471-45173-8. - Engle, Robert F. (1982). "Autoregressive Conditional Heteroscedasticity with Estimates of Variance of United Kingdom Inflation".
*Econometrica*.**50**(4): 987–1008. JSTOR 1912773.*(the paper which sparked the general interest in ARCH models)* - Engle, Robert F. (1995).
*ARCH: selected readings*. Oxford University Press. ISBN 0-19-877432-X. - Engle, Robert F. (2001). "GARCH 101: The Use of ARCH/GARCH Models in Applied Econometrics".
*Journal of Economic Perspectives*.**15**(4): 157–168. doi:10.1257/jep.15.4.157. JSTOR 2696523.*(a short, readable introduction)* - Gujarati, D. N. (2003).
*Basic Econometrics*. pp. 856–862. - Hacker, R. S.; Hatemi-J, A. (2005). "A Test for Multivariate ARCH Effects".
*Applied Economics Letters*.**12**(7): 411–417. doi:10.1080/13504850500092129. - Nelson, D. B. (1991). "Conditional Heteroskedasticity in Asset Returns: A New Approach".
*Econometrica*.**59**(2): 347–370. JSTOR 2938260.