---
title: Factor Augmented Dynamic Nelson–Siegel Model
url: https://www.emergentmind.com/topics/factor-augmented-dynamic-nelson-siegel-model
type: topic
---

# Factor Augmented Dynamic Nelson–Siegel Model

The Factor Augmented Dynamic Nelson–Siegel (FADNS) model with Wishart stochastic volatility is a Bayesian state-space framework for modeling and forecasting the term structure of commodity futures, extending the classic Nelson–Siegel (NS) representation to incorporate time-varying multivariate volatility and multiple latent factors governing the dynamics of forward curves. The model builds on the dynamic 3-factor Nelson–Siegel model and its 4-factor Svensson extension, employing a Gaussian vector autoregression structure for the latent factors and a Wishart process to capture stochastic volatility. Its development enables robust estimation and uncertainty quantification for large cross-sections and long time series, yielding competitive performance for in-sample fit, out-of-sample prediction, and risk management in empirical applications such as WTI crude oil futures [1908.07798].

## 1. Model Structure: Measurement Equation

The observed data consist of daily log-prices \( y_{i,t} \), where \( i=1,\dots,N \) enumerates futures contracts with times-to-maturity \( \tau_{i,t} \) at time \( t \). The core observation model follows the dynamic Nelson–Siegel (NS) form:
\[
y_{i,t} = \beta_{1,t} + \beta_{2,t} \frac{1-e^{-\lambda_1 \tau_{i,t}}}{\lambda_1 \tau_{i,t}}
+ \beta_{3,t} \left( \frac{1-e^{-\lambda_1 \tau_{i,t}}}{\lambda_1 \tau_{i,t}} - e^{-\lambda_1 \tau_{i,t}} \right) + \varepsilon_{i,t}, \quad \varepsilon_{i,t} \sim N(0,\sigma_{y,i}^2)
\]
where \(\beta_{1,t}\), \(\beta_{2,t}\), \(\beta_{3,t}\) denote the time-varying level, slope, and curvature factors.

The 4-factor Svensson extension augments this with a second curvature component:
\[
y_{i,t} = \beta_{1,t} + \beta_{2,t}\frac{1-e^{-\lambda_1\tau_{i,t}}}{\lambda_1\tau_{i,t}}
+ \beta_{3,t}\left( \frac{1-e^{-\lambda_1\tau_{i,t}}}{\lambda_1\tau_{i,t}} - e^{-\lambda_1\tau_{i,t}} \right)
+ \beta_{4,t}\left( \frac{1-e^{-\lambda_2\tau_{i,t}}}{\lambda_2\tau_{i,t}} - e^{-\lambda_2\tau_{i,t}} \right) + \varepsilon_{i,t}
\]

Stacked notation gives:
\[
y_t = Z_t \beta_t + \varepsilon_t, \qquad \varepsilon_t \sim N(0, \Sigma_y)
\]
where \(Z_t\) is the design matrix of NS or Svensson loadings.

## 2. Latent Factor Dynamics: Vector Autoregression

The evolution of latent factors \(\beta_t \in \mathbb{R}^m\) for \(m=3\) (NS) or \(m=4\) (Svensson) is governed by a first-order vector autoregression (VAR), typically set to a random walk:
\[
\beta_t = \alpha + \Phi \beta_{t-1} + \eta_t, \qquad \eta_t \sim N(0, H_t^{-1})
\]
where \(\alpha \in \mathbb{R}^m\), \(\Phi \in \mathbb{R}^{m \times m}\) (commonly \(\Phi=I_m\)), and \(H_t\) is the time-varying precision matrix encoding the stochastic volatility. The generic transition density is
\[
p(\beta_t | \beta_{t-1}, H_t) = \mathcal{N}(\alpha+\Phi \beta_{t-1}, H_t^{-1})
\]
which allows for parsimonious yet flexible modeling of factor persistence and interdependence.

## 3. Wishart Stochastic Volatility Specification

Multivariate time-varying volatility is introduced via a Wishart–Beta process (Uhlig 1994, 1997; Windle & Carvalho 2014), driving the evolution of the factor innovation precision matrix \(H_t\). Its transition equation in scaled-Beta form is
\[
H_t = \frac{1}{\gamma} H_{t-1}^{1/2} \Psi_t H_{t-1}^{1/2}, \qquad \Psi_t \sim \mathcal{B}_m \left( \tfrac{\nu}{2}, \tfrac{1}{2} \right)
\]
where \(\nu > m-1\) (degrees of freedom) and \(\gamma > 0\) (scale) with the identification constraint \(1/\gamma = 1 + \tfrac{1}{\nu-m-1}\) ensuring \(E(H_t|H_{t-1})=H_{t-1}\).

Initial precision \(H_1\) follows
\[
H_1 \sim \mathcal{W}_m\left(\nu, \Sigma_0^{-1}/\gamma\right)
\]
with user-selected positive-definite matrix \(\Sigma_0\).

## 4. Bayesian Posterior Inference and MCMC Estimation

The model is estimated in a fully Bayesian framework. The joint posterior for all parameters (\(\theta = (\lambda, \sigma_y, \alpha, \nu, \beta_0, \Sigma_0)\)), latent factors \(\beta_{1:T}\), and stochastic volatilities \(H_{1:T}\) is proportional to the product of Gaussian densities (likelihood and factor transition), the Wishart/Beta prior, and priors on remaining parameters:
\[
\pi(\beta_{0:T}, H_{1:T}, \theta | y_{1:T}) \propto p(y_{1:T}| \beta_{1:T}, \sigma_y) \times p(\beta_{1:T}| H_{1:T}, \alpha) \times p(H_{1:T}|\nu,\Sigma_0) \times p(\theta)
\]

Key conjugacies allow:
- Given \(H_{1:T}\), \(\beta_{0:T}\) form a linear-Gaussian state-space model with tractable Gaussian updates.
- Given \(\beta_{1:T}\), the Wishart specification produces a closed-form backward sampler for \(H_{1:T}\).
- Parameters \(\alpha\) (Normal) and \(\sigma_y^2\) (inverse gamma) have conjugate updates.

A collapsed-Gibbs sampler cycles through:
1. Jointly sampling \((\beta_{0:T}, \lambda)\): update \(\lambda\) via random-walk Metropolis (marginalizing \(\beta\)), update \(\beta_{0:T}\) with a sparse Gaussian precision sampler (as in Chan & Jeliazkov 2009).
2. Jointly sampling \((H_{1:T}, \nu)\): update \(\nu\) via marginal likelihood, then backward-sample \(H_{1:T}\) from the shifted singular Wishart.
3. Standard updates for \(\alpha\) and \(\sigma_y^2\).

Posterior quantities are obtained by averaging across post-burn-in MCMC cycles.

## 5. Empirical Performance in Crude Oil Futures

An empirical application to 24 monthly WTI crude-oil log-futures (Jan 1996–May 2016, \(T=5118\)) compares several model variants: 3F (homoscedastic NS), 3F-SV (NS + Wishart SV), 4F (homoscedastic Svensson), and 4F-SV (Svensson + Wishart SV). Weakly-informative priors are specified.

Estimation and forecasting results include:
- Posterior means (4F-SV, second estimation window): \(\lambda_1 \approx 0.0036\), \(\lambda_2 \approx 0.0158\), \(\sigma_y \approx 0.00316\), \(\nu \approx 23.97\) (implying \(\gamma \approx 0.958\)). The intercept \(\alpha\) is close to zero, indicating factor processes are approximately unit-root random walks.
- Effective sample sizes for all parameters exceed 200.
- Deviance Information Criterion (DIC) ranks: 4F-SV best, followed by 4F, 3F-SV, then 3F.
- Out-of-sample log-predictive likelihoods (2008 window, 24 contracts, 576 dates): 3F–24,853; 3F-SV–25,019; 4F–26,034; 4F-SV–26,195. Thus, 4F-SV displays superior predictive density performance.
- One-day RMSE (2008 window): 3F–0.0290; 3F-SV–0.0290; 4F–0.0289; 4F-SV–0.0289 versus the random-walk benchmark of 0.0286.
- In VaR forecasting, 4F-SV produces the most accurate hit rates at 1%, 5%, and 10%, and passes unconditional/conditional coverage tests more consistently than alternatives.

Posterior smoothed volatility and cross-correlation estimates closely track rolling-window realized measures. Both curvature factors are empirically distinct: one captures longer-horizon term-structure curvature, while the second extracts a short-horizon “bump.”

## 6. Model Comparison and Interpretation

Model comparisons indicate that the 4-factor Svensson with Wishart stochastic volatility (4F-SV) outperforms lower-order or homoscedastic alternatives in both in-sample goodness-of-fit and out-of-sample density/point forecasting. The parsimonious stochastic volatility specification captures temporal variation and cross-sectional dependence among risk factors with high persistence (\(\nu \approx 24\), \(\gamma \approx 0.96\)), recovering both the level and complex curvature dynamics of the commodity forward curve.

A plausible implication is that incorporating Wishart SV enhances risk management and density forecasting capability compared to homoscedastic or classic factor models, notably improving Value-at-Risk assessment for portfolios including long–short (bull spread) positions.

| Model Variant | Factors         | SV Process | Best Use            |
|---------------|----------------|------------|---------------------|
| 3F            | NS (level, slope, 1 curvature) | No         | Baseline comparison |
| 3F-SV         | NS             | Yes (Wishart)   | Volatility modeling |
| 4F            | Svensson (2 curvature) | No         | Enhanced fit        |
| 4F-SV         | Svensson       | Yes (Wishart)   | Forecast/risk       |

## 7. Conclusions and Practical Value

The FADNS model with Wishart stochastic volatility, as formulated in Kleppe et al. [1908.07798], provides a flexible, computationally efficient, and Bayesian-coherent approach to modeling the term structure of commodity futures. The fully conjugate state-space and stochastic volatility structure enables robust inference even for large panels (\(N=24\), \(T\approx5000\)), with jointly estimated latent factors and time-varying covariance, and competitive or superior performance versus common benchmarks such as linear Gaussian state-space models and random walk processes. The empirical evidence suggests that the model is particularly effective for forecasting time-varying forward curves, volatility, cross-factor covariances, and portfolio risk in commodity markets.

Source: https://www.emergentmind.com/topics/factor-augmented-dynamic-nelson-siegel-model