Skip to content
All library documents

Choosing Gaussian or Student Errors for ARMA-GARCH Models

Article Quant Q&A · Author: Rodrigo Guinea

Summary

The document considers whether non-normal residuals make a Gaussian ARMA-GARCH model unsuitable when a Student error distribution fits a QQ plot better. Its central guidance is that model quality depends on the intended use: estimating a mean or autoregressive effect, forecasting volatility, calculating risk measures, or simulating volatility over different horizons can favor different specifications. Residual shape alone does not settle the choice.

The discussion compares qualitative implications in the example. The Gaussian fit has volatility persistence above one, while the Student specification has lower persistence but an extremely heavy estimated tail; the answers caution that each may imply implausible volatility behavior. They also raise structural shifts or jumps as possible features that a basic GARCH fit may be absorbing. Forecast evaluation should use a loss function suited to the task, and information criteria may inform selection. These are diagnostic considerations, not a universal ranking; the suggested checks depend on the model’s purpose and the data.

Key ideas

  • Choose an ARMA-GARCH specification according to its intended forecasting, risk, simulation, or inference use.
  • A better-fitting residual distribution does not by itself establish that a model is preferable.
  • Compare volatility persistence and tail estimates for plausible behavior over the forecast horizon.
  • Forecast evaluation depends on the chosen loss function and the quantity being predicted.
  • Structural changes or jumps may explain patterns that a basic GARCH model attributes to volatility dynamics.

Tags

Full text
# Is it too important that my residuals be normal? I am Using an ARMA/GARCH model


# Is it too important that my residuals be normal? I am Using an ARMA/GARCH model












I am trying to fit an ARMA/GARCH model to a time series. I found that the best candidate is an ARMA(1,0) + GARCH(1,1) with gaussian white noise

It has coefficients with p-values near cero and the residuals are white noise. The problem is that the Jarque Bera Test says the residuals are not normal. The QQ normal plot

confirm that. And when I try with several ARMA/GARCH models with t-student white noise, for example:

the QQ t-student plot

fits very well (except for some outliers), but the rest don't (I mean not as good as the first one).

Which one is better?

I have been stuck in this problem for a while.

Thank you very much

Rodrigo

## Answer by Kiwiakos (score 2, accepted)

https://quant.stackexchange.com/a/21916

To get it out the way: you cannot ask 'what model is better' without a reference to what its use is. Do you want to test for the mean or the AR parameter to trade it? Do you want to calculate VaR? Do you want to forecast volatility over one period? Or over 1000 periods? Or higher moments? Do you want to simulate volatility over one period? Or longer?

For some of these your first model can be perfectly adequate, while for other is will not be.

In general I do not take 'statistical significance' at face value. It does not mean much. The (very small) standard errors you get are not reliable.

Having said that, you can look at the qualitative differences that your models reveal.

- They agree on their point estimates of mu and ar1

- The Gaussian model gives very high persistence, alpha1+beta is actually above one. It tries to capture an extremely heavy tail and produces a volatility unit root to achieve that.

- If you were to simulate vols over longer horizons with the Gaussian model, they would be all over the place.

- The Student model gives more stable volatility, the sum is .96. Perhaps it reverts too fast. To achieve that it estimated a very heavy tail, if 'shape' means 'degrees of freedom'. So heave that it barely has a variance.

- Is that mean reverting though? When the error is that fat tailed volatility will look very choppy. A large jump and then exponential decay. Is that how the world looks like?

- What shocks me is the increase in log-likelihood. If you see the second model as an extension of the first by one parameter, then you have a gigantic improvement. Big enough to make me suspicious.

- To me it looks like a process with some Garch perhaps, but with structural changes superimposed. Something very bad happens occasionally, which Garch is trying to mimic. Jump diffusion, regime switching, multi fractal, exponential vol, anything that gives rapid structural changes. If you could give us a plot of the time series (before and after differencing) then we could speculate on that.

But as I said, it all depends on what you want your model to do. Simple Gaussian Garch(1,1) is hard to beat out of sample consistently.

## Answer by Malick (score 1)

https://quant.stackexchange.com/a/21913

A good rule of thumb is to "test" your models by doing forecasts and to choose the best one. Note however that your choice will be based upon the loss function you selected. If you are concerned about outliers you should (for instance) use Median Squared Errors, if you don't you can use Mean Square Errors. In your particular case the Information Criteria Statistics highly advocates for the use of the Student distribution.

So as general answer it depends of your aim, the choice will differs if you’re estimating for 1) forecasting, 2) parameters calibration, 3) hedging 4) research … I don’t think that a formal answer covering all case can be given. The essential is to understand the implication of (not) choosing a model. (ex: if you are particularly concerned with extreme events, the QQ Plot warns you that the Gaussian hypothesis underestimate these events.)

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.