← Back to Research & Projects

Time-Series Analysis of 20+ Year Treasury Bond ETF

Abstract

This paper analyzes the behavior of daily returns of the iShares 20+ Year Treasury Bond ETF (TLT), focusing on the impact of macroeconomic and financial market variables, including market volatility index, S&P500 returns, changes in Fed rates, Treasury bond market yield, and inflation. Models utilized daily data from January 2004 to July 2024, the analysis incorporates advanced time-series models to capture the dynamics of returns in TLT. An ARIMAX model with state-dependent dynamics using a Threshold Autoregressive (TAR) framework reveals the significant role of Federal Funds Rate changes in regime-dependent effects on returns. Additionally, a GARCH-ARIMAX model was employed to address volatility clustering and persistence, highlighting the impact of recent shocks and long-term volatility persistence. Both models demonstrate the critical influence of market volatility and equity returns on the daily performance of TLT. The findings provide a comprehensive understanding of the drivers of long-term Treasury bond returns and underscore the importance of accounting for regime changes and time-varying volatility in modeling.

Data

  • Asset: TLT daily returns (R_TLT)
  • Sample: Jan 5, 2004 – Jul 1, 2024 (5,156 trading-day observations)
  • Sources: Yahoo Finance (prices/returns) + FRED (macro/market series aligned to daily frequency)
  • Core drivers: VIX, S&P 500 returns (R_GSPC), inflation (INFL), 20-year Treasury yield (DGS20), changes in Fed funds rate (dFEDRATE)
  • Regime definition: Monetary-policy state based on dFEDRATE (rate-cut regimes vs stable/increasing-rate regimes)

Methods

  • Baseline benchmark: ARIMA(2,0,2) selected by AIC
    AIC: 13,765.91
    Ex-post RMSFE: 1.0598
  • Preferred mean-dynamics extension (regime switching): ARIMAX(2,0,2) + TAR
    Regimes defined by dFEDRATE and implemented using regime-specific lag interactions (TLT lag effects vary by policy state)
    AIC: 13,148.81
    Ex-post RMSFE: 0.9097872
  • Volatility modeling (best overall in the writeup): GARCH(1,2) variance layered on an ARIMA(1,0,3)
    Captures volatility clustering and persistent variance dynamics
    ARCH effects: ARCH LM χ² = 725.141 (p = 0.0000)
    AIC: 12,191.62
    Ex-post RMSFE: 0.90935889

Key Results

What the models agree on:

  • Flight-to-safety is the dominant daily mechanism. Across specifications, VIX is positively associated with TLT returns and S&P 500 returns are negatively associated, consistent with risk-off periods supporting long-duration Treasuries.
  • Monetary policy regimes materially change return dynamics. The ARIMAX+TAR results show that the impact of lagged TLT returns differs between rate-cut and stable/increasing-rate regimes (dFEDRATE-based state interactions).
  • Volatility clustering is large and persistent. Strong ARCH effects and GARCH dynamics indicate that volatility shocks propagate forward, making constant-variance models inadequate for this series.
  • Regime + volatility structure improves forecast performance vs baseline. Relative to the ARIMA benchmark (RMSFE 1.0598), both the TAR and GARCH specifications reduce forecast error to approximately 0.91, with the GARCH-ARIMAX model also delivering the strongest overall fit (lowest AIC).

Conclusion

Across 2004–2024, TLT's daily returns are best explained by risk-off conditions (higher VIX and weaker equities) and monetary-policy regime shifts captured by changes in the Fed funds rate. Incorporating regime dependence and time-varying volatility improves fit and reduces forecast error relative to a baseline ARIMA model, and the GARCH-ARIMAX specification provides the most realistic description of volatility behavior in TLT returns.

Key Model Outputs

GARCH(1,2) with ARIMAX(1,0,3) Model Results

Variable Coefficient Std. Error t-value p-value Sig
X-Var (Explanatory Variables)
VIX (CBOE Volatility Index)0.0300.0093.240.001***
VIX L (VIX Lag 1)-0.0290.009-3.150.002***
R_GSPC (S&P 500 Returns)-0.2230.015-14.920.000***
R_GSPC L (S&P 500 Returns Lag 1)0.0110.0091.290.198
INFL (Inflation)-0.0170.006-2.680.007***
DGS20 (20-Year Treasury Constant Maturity Rate)0.0170.0062.760.006***
dFEDRATE (Change in Federal Funds Rate)-0.6900.327-2.110.035**
AR (Autoregressive)
L (AR Lag 1)-0.9040.077-11.690.000***
MA (Moving Average)
L (MA Lag 1)0.8860.07811.360.000***
L2 (MA Lag 2)-0.0480.019-2.490.013**
L3 (MA Lag 3)-0.0440.014-3.040.002***
ARCH (Autoregressive Conditional Heteroskedasticity)
L (ARCH Lag 1)0.0790.00710.930.000***
GARCH (Generalized Autoregressive Conditional Heteroskedasticity)
L (GARCH Lag 1)0.3890.1023.800.000***
L2 (GARCH Lag 2)0.5210.0985.290.000***
Constant
Constant0.0090.0025.100.000***

*** p<0.01, ** p<0.05, * p<0.1. Number of obs: 5155. Mean dependent var: 0.018. SD dependent var: 0.921. Chi-square: 1220.283 (Prob > chi2: 0.000). Akaike crit. (AIC): 12191.621.

GARCH ARIMAX Forecast Plot

ARIMAX(2,0,2) with TAR Model Results

Variable Coefficient Std. Error t-value p-value Sig
TAR (Threshold Autoregressive)
state_neg (Rate Cut Regime Intercept)0.0230.0360.630.527
state_neg_TLT (Rate Cut Regime: TLT Lag 1)0.0730.0125.860.000***
state_neg_TLT2 (Rate Cut Regime: TLT Lag 2)-0.0820.012-7.010.000***
state_pos_same (Stable/Increasing Rate Regime Intercept)0.0370.0311.210.228
state_pos_same_T (Stable/Increasing Rate Regime: TLT Lag 1)-0.0470.014-3.300.001***
state_pos_same_TLT2 (Stable/Increasing Rate Regime: TLT Lag 2)-0.0300.015-2.000.045**
X-Var (Explanatory Variables)
VIX (CBOE Volatility Index)0.0260.0073.860.000***
VIX L (VIX Lag 1)-0.0240.006-3.730.000***
R_GSPC (S&P 500 Returns)-0.2210.010-21.900.000***
R_GSPC L (S&P 500 Returns Lag 1)0.0310.0084.100.000***
INFL (Inflation)-0.0140.005-2.800.005***
AR (Autoregressive)
L (AR Lag 1)1.9330.006299.610.000***
L2 (AR Lag 2)-0.9840.007-151.140.000***
MA (Moving Average)
L (MA Lag 1)-1.9410.005-415.110.000***
L2 (MA Lag 2)0.9920.005211.050.000***
Constant
Constant0.8640.006144.910.000***

*** p<0.01, ** p<0.05, * p<0.1. Number of obs: 5155. Mean dependent var: 0.018. SD dependent var: 0.921. Chi-square: 1268245.504 (Prob > chi2: 0.000). Akaike crit. (AIC): 13148.809.

ARIMAX TAR Forecast Plot