# GARCH model using high frequency price return

I would like to forecast variance at time length $$k\delta$$ based on a price (return) time series of time step length $$\delta$$. I will apply a GARCH(1,1) model to subsamples at time intervals length $$k\delta$$ on a stock return time series $$\big(r(i\delta,(i+1)\delta)\big)_{i=0}^I$$ each element of which is the return between time $$i\delta$$ and $$(i+1)\delta$$. I take the recursion formula to be $$h(t,t+k\delta) = c+a\,u(t-k\delta,t)^2 +b\,h(t-k\delta,t) \tag1$$ where $$h(t-k\delta,t)$$ is the estimated variance and $$r(t-k\delta,t)$$ is the return for time interval $$(t-k\delta,t)$$. I would like to use the full return time series for Equation (1).

1. Is it correct to use the following variance estimation for time interval $$(t-k\delta,t)$$? $$u(t-k\delta,t)^2 := \sum_{i=1}^k r\big(t-i\delta,r-(i-1)\delta\big)^2.$$ This is then substituted into the maximum likelihood optimizer as the variance for time interval $$(t-k\delta,t)$$ in place of the usual simple estimator $$r(t-k\delta,t)^2$$. $$\big(u(jk\delta,(j+1)k\delta)\big)_{j=0}^{q-1}$$ forms a new time series. Its Gaussian negative logarithmic likelihood $$l(a,b,c):=\sum_{j=0}^{q-1} \bigg( \frac{u(jk\delta,(j+1)k\delta)^2}{h(jk\delta,(j+1)k\delta)}+\ln h(jk\delta,(j+1)k\delta)\bigg).$$

2. Do I have to use something like the realized kernel as constructed in Realised Kernels in Practice: Trades and Quotes, by Ole E. Barndorff-Nielsen, Peter R. Hansen, Asger Lunde and Neil Shephard?

Scaling the length of the discrete time step in GARCH models is, from memory, not at all straightforward. For example, you can't just multiply things by the square root of time, like we do for some other, simpler, processes.

For the case where $$\delta \rightarrow 0$$, the convergence is derived in Nelson (1990) "ARCH Models as Diffusion Approximations". I think there has also been some work done on scaling from one discrete step length to another, but the reference escapes me now (but this would effectively provide the academic answer to your question if you can find it).

Regarding your specific questions, replacing $$u^2$$ with $$\sum_i r_i^2$$ is equivalent to replacing the squared return with the realized variance estimator. The result is not a standard GARCH model (I don't think...), although I suspect it would perform quite well empirically, so just because it doesn't fit nicely into the standard framework is not a reason to abandon it per se.

As to whether you would need to use something like Realized Kernels, this would depend entirely on the the length $$\delta$$. At lower frequencies, e.g. 5-minutes, you can just sum squared returns fairly safely without worrying too much about microstructure noise (although of course this all depends on the liquidity of the underlying asset - the more frequently traded the asset, the higher the frequency you can go to without having to worry about microstructure noise). At higher frequencies, e.g. 5-seconds, you would probably need to use an estimator that corrects for microstructure noise, such as Realized Kernels.

One final point of interest, if $$c = 0$$, and $$a + b = 1$$, and if we're prepared to make some additional heroic assumptions, like the realized variance converging to the true variance over the interval $$(t - k \delta, t)$$, then I think the model as you've written it converges, as $$\delta \rightarrow 0$$, to standard simple exponential smoothing, e.g.: $$\begin{equation} \hat{y}_t = \alpha y_{t-1} + (1 - \alpha) \hat{y}_{t-1} \end{equation}$$ Note, in this equation $$\hat{y}_t$$ is the variance forecast and $$y_{t}$$ is the true variance. The "heroic assumptions" are those necessary for the realized variance to converge to the true variance.

ADDITIONAL: You seem to be very focused on whether your model is "correct" or "legitimate". It is important to understand that those words don't really have any meaning here. All models are wrong. And all estimation techniques are inefficient, except under ideal theoretical assumptions that are rarely (never?) satisfied in the real world. So I'm not really sure how to respond to that aspect of the question. Perhaps this:

What matters is whether the model helps you to reach your end goal. For example, one possible end goal might be find a set of volatility forecasts that are useful predictors for implied volatility (i.e. option prices) or the VIX so you can design a profitable trading rule. In this case, perhaps your metric for "better" and "worse" might be something like: "find a volatility forecast model with smaller out-of-sample squared error than the standard GARCH(1,1)". In my experience, in this context, replacing squared returns with realized variances constructed from higher frequency returns spanning the same interval almost always reduces estimation error and will help you towards that goal.

In regards to methods of estimation, we can apply the same logic. Sure, theory might tell you to construct a maximum likelihood estimator with Gaussian innovations, and as far as I can tell what you've written looks good (but I make no promises here - I'm not going to grind through the maths myself to check your working :-). But that might not give you the best out-of-sample forecast. Indeed, in my experience, robust methods of estimation (i.e. those that reduce/eliminate the impact of observations in the tail of the distribution) tend to provide better estimates when working with financial data.

• 1. I have edited my question to add more details of my parameter estimation procedure. Please see if it matches what you thought I was doing. By constructing a new time series $(u)$, I transform the original subsampling problem into a canonical time series. I then apply the GARCH model together with its maximal likelihood parameter estimation to the latter time series. I can apply more complicated kernel in place of the realized variance estimator if, as you say, the original sampling frequency is higher than 5 minutes. Is it legitimate?
– Hans
Sep 24, 2020 at 19:40
• 2. Is your $\hat y_t$ the price or volatility or variance which is my $h_t$?
– Hans
Sep 24, 2020 at 19:40
• @Hans I've amended my answer, but I'm not sure how happy you'll be with my amendment. The fundamental problem here is that you seem very focused on whether something is "correct" or "legitimate", when those words don't really have much meaning in this context. Sep 25, 2020 at 0:08
• @Hans The answer to 1 is definitely affirmative. In most models using squared returns we are (usually) implicitly using them as a proxy for true variance. Any technique that improves the accuracy of the proxy (such as replacing it with realized variance) is likely to improve the model and is consistent within that general framework. Regarding 2) the answer is also yes. See eg Hansen, Huang, Shek (2011) "Realized GARCH". As a bonus, I think Peter Hansen is one of the best financial econometricians on the planet. Sep 28, 2020 at 3:03
• @Hans I really don't think I could improve on the answer that someone else there has already posted. A great response! Mar 13, 2021 at 9:05