MIDAS Regression

midasdistributed-laghigh-frequencyvolatility-forecastingbeta-polynomialmixed-frequencyforecasting

Definition

MIDAS (Mixed Data Sampling) regression is a flexible distributed-lag model that relates a low-frequency dependent variable to a weighted sum of high-frequency predictors, with the lag weights parameterized by a small number of free parameters. The framework accommodates any mismatch between the frequency of the outcome and the frequency of the predictors without requiring aggregation of the predictors to the lower frequency first.

Key Ideas

How It Works

MIDAS Specification for Volatility Forecasting (Ghysels-Santa-Clara-Valkanov [GSV] 2006)

The baseline MIDAS regression (eq. 1.1 in GSV 2006) is:

Vt+H,t(Hm)=μH+ϕHk=0kmaxbH(k;θ)X~tk,tk1(m)+εHtV_{t+H,t}^{(Hm)} = \mu_H + \phi_H \sum_{k=0}^{k^{\max}} b_H(k;\theta)\,\tilde{X}_{t-k,t-k-1}^{(m)} + \varepsilon_{Ht}

where Vt+H,t(Hm)V_{t+H,t}^{(Hm)} is realized variance over the next HH days, X~(m)\tilde{X}^{(m)} is a daily volatility measure computed from mm-period intraday returns, and bH(k;θ)b_H(k;\theta) are Beta polynomial weights with kbk=1\sum_k b_k = 1.

Five Predictor Specifications

GSV (2006) evaluate five choices of X~(m)\tilde{X}^{(m)}:

Predictor Formula Notes
Realized variance Q~t(m)\tilde{Q}^{(m)}_t j=1mrt(j1)/m2\sum_{j=1}^m r_{t-(j-1)/m}^2 Converges to IVtIV_t (Integrated Variance) as mm\to\infty; contaminated by jumps
Squared daily return rt2r_t^2 (ptpt1)2(p_t - p_{t-1})^2 Noisiest proxy; daily frequency (m=1m=1)
Absolute daily return rt|r_t| ptpt1|p_t - p_{t-1}| Less noise than rt2r_t^2
Daily range [hi-lo]t[hi\text{-}lo]_t (pthiptlo)2/(4ln2)(p_t^{hi} - p_t^{lo})^2/(4\ln 2) Garman-Klass estimator
Realized power P~t(m)\tilde{P}^{(m)}_t j=1mrt(j1)/m\sum_{j=1}^m |r_{t-(j-1)/m}| Barndorff-Nielsen–Shephard 2003/2004; excludes jump component

Andersen-Bollerslev-Diebold-Labys (ABDL) autoregressive fractionally integrated (ARFI)(5,d) Benchmark (eq. 1.10)

The long-memory benchmark against which MIDAS is evaluated is:

(1k=15bA(k)Lk)(1L)dlogQ~t+1,t(m)=μ+εt,d0.401(1-\textstyle\sum_{k=1}^5 b_A(k)L^k)(1-L)^d\,\log\tilde{Q}_{t+1,t}^{(m)} = \mu + \varepsilon_t, \quad d\approx 0.401

where the bA(k)b_A(k) are unconstrained autoregressive (AR)(5) coefficients applied to fractionally-differenced log realized variance. This model captures long memory explicitly via the fractional difference operator.

Empirical Results (GSV 2006)

Data: Dow Jones index + 6 stocks (DIS, GE, JPM, XOM, MCD, HON); April 1993 – October 2003; 2,669 trading days; m=80m=80 five-minute intervals per day.

Performance metric: mean squared error (MSE) ratio MSEMIDAS/MSEABDLMSE^{MIDAS}/MSE^{ABDL} — values below 1 mean MIDAS beats the long-memory benchmark.

Predictor Horizon HH In-sample MSE/ABDL (level) Out-of-sample
Realized power P~(m)\tilde{P}^{(m)} 1–20 days 0.606–0.912 0.714–0.897
Daily range [hi-lo][hi\text{-}lo] 1–20 days 0.7\approx 0.71.051.05 0.8\approx 0.81.01.0
Realized variance Q~(m)\tilde{Q}^{(m)} 1–20 days 0.85\approx 0.851.201.20 0.9\approx 0.91.11.1

Predictor ranking (consistent across assets and horizons): P~(m)[hi-lo]Q~(m)rr2\tilde{P}^{(m)} \succ [hi\text{-}lo] \succ \tilde{Q}^{(m)} \succ |r| \succ r^2

Why realized power wins: The Beta MIDAS with realized power places weight b(1)=38.3%b(1) = 38.3\% on yesterday's power at H=1H=1, declining to 8.6%8.6\% at H=20H=20 days (Table 2). The remaining weight spreads over ~30 lags. The ABDL ARFI model effectively requires >50>50 lags to represent its long-memory tail, making it parameter-heavier despite the long-memory shortcut.

No gain from direct intraday use: Replacing daily P~(m)\tilde{P}^{(m)} with raw 5-minute absolute returns as MIDAS regressors (eqs. 3.1–3.2) yields no improvement in MSE. The information in 5-min data is fully captured by daily realized power.

Why Realized Power Has Higher Persistence and Forecasting Power

Under the Barndorff-Nielsen–Shephard (2004) framework, realized power converges to the continuous component of quadratic variation (integrated variance), while realized variance Q~(m)\tilde{Q}^{(m)} includes jumps. Since jumps are transient, their inclusion in the predictor reduces persistence and degrades multi-step forecast accuracy. Empirically, logP~(m)\log\tilde{P}^{(m)} has a higher AR(1) coefficient and smoother autocorrelation function (ACF) decay than logQ~(m)\log\tilde{Q}^{(m)}.

Why It Matters

Open Questions

Related