Nelson-Siegel Model: Factors, Fitting, and Extensions

The Nelson-Siegel model is a formula that describes an entire yield curve using just three factors and one decay parameter. Charles Nelson and Andrew Siegel introduced it in 1987, and it remains the foundation for how most central banks estimate their official yield curves today.1Bank for International Settlements. Zero-Coupon Yield Curves Estimated by Central Banks Its appeal is parsimony: a handful of numbers can capture the level, slope, and curvature of interest rates across every maturity from overnight lending out to 30 years.

The Three Factors

Every Nelson-Siegel curve is built from three structural components, each tied to a different segment of the maturity spectrum. Practitioners label them the level factor, the slope factor, and the curvature factor.

Level

The level factor is a baseline that every maturity eventually converges toward. It corresponds to the parameter labeled Beta 0, and it equals the yield the model predicts for an infinitely long bond. Economically, it captures long-run inflation expectations and the real equilibrium interest rate. When the level factor rises, the whole curve shifts up; when it falls, every maturity drops with it.

Slope

The slope factor, controlled by Beta 1, measures the gap between short-term and long-term rates. At the shortest maturity the model produces a yield equal to Beta 0 plus Beta 1, and at very long maturities Beta 1’s influence decays to zero. A negative Beta 1 means the short end sits below the long end, which is the normal upward-sloping curve most people picture. A positive Beta 1 flips that relationship and produces an inverted curve, the shape often associated with recession expectations.2EFMA. Estimating the Yield Curve Using the Nelson-Siegel Model: A Ridge Regression Approach

Curvature

The curvature factor, governed by Beta 2, creates the hump or trough in the middle of the curve. Its loading function starts at zero for the shortest maturities, rises to a peak at a maturity determined by the decay parameter, then fades back toward zero further out. The location of that peak is not fixed. Diebold and Li, for example, calibrated the decay parameter so the curvature peaked around 2.5 years, arguing that most humps fall between the second and third year.2EFMA. Estimating the Yield Curve Using the Nelson-Siegel Model: A Ridge Regression Approach The curvature factor tends to reflect investor uncertainty about the medium-term path of policy rates, risk premiums for intermediate bonds, or shifts in supply and demand at specific maturities.

How the Parameters Shape the Curve

The full Nelson-Siegel formula produces a yield for any maturity using three linear parameters (Beta 0, Beta 1, Beta 2) and one nonlinear decay parameter, usually labeled Lambda or Tau.

Beta 0 sets the asymptotic level. As time to maturity grows, both the slope and curvature loading functions decay to zero, and the yield converges to Beta 0 alone.2EFMA. Estimating the Yield Curve Using the Nelson-Siegel Model: A Ridge Regression Approach Under normal conditions it should be positive, since it represents a long-run nominal interest rate.

Beta 1 determines how much the short end deviates from that long-run level. The instantaneous yield, at maturity approaching zero, equals Beta 0 plus Beta 1. Because the slope loading decays exponentially, Beta 1’s effect is strongest at short maturities and progressively weaker further out. Central bank rate decisions and overnight funding conditions are the primary economic drivers here.

Beta 2 controls the hump. Its loading is the difference between two exponential decay terms, which is why it starts at zero, rises, and then falls. The sign of Beta 2 dictates whether you get a hump (positive) or a trough (negative) in the middle maturities. The magnitude sets its height or depth.

Lambda governs how quickly the slope and curvature loading functions decay. A small Lambda stretches those effects further along the maturity axis and pushes the curvature peak to longer maturities. A large Lambda concentrates the action at the short end, so the hump peaks earlier. Lambda is the one parameter that makes the model nonlinear, and it is the trickiest to estimate because it interacts with every other parameter at once.

Fitting the Model to Market Data

Estimation means finding the Beta and Lambda values that minimize the gap between the model’s predicted yields and observed market yields. The standard objective is the sum of squared differences. Practitioners split into two camps depending on how they treat Lambda.

Standard practice calls for using zero-coupon (spot) yields rather than par yields. A coupon-bearing bond bundles cash flows at many different maturities into one price, so its yield blurs the information about any single maturity. If only coupon-bearing bond prices are available, a bootstrapping procedure extracts the implied zero-coupon rates first. Bonds trading far from par, callable issues, and illiquid bonds are typically filtered out before estimation, since their prices carry noise unrelated to the pure term structure.

Grid Search With OLS

Nelson and Siegel’s original approach sidesteps the nonlinearity by fixing Lambda at a specific value. With Lambda held constant, the formula becomes linear in Beta 0, Beta 1, and Beta 2, and ordinary least squares produces exact solutions. The analyst repeats this across a grid of Lambda values and picks whichever Lambda yields the lowest overall fitting error.2EFMA. Estimating the Yield Curve Using the Nelson-Siegel Model: A Ridge Regression Approach It is computationally simple and avoids convergence issues, but the grid resolution limits how precisely Lambda can be pinpointed.

Simultaneous Nonlinear Optimization

The alternative is to estimate all four values at once with a nonlinear optimizer. The Nelder-Mead simplex algorithm is a common choice, sometimes as a first stage that generates reasonable starting values before a gradient-based method refines the result. This approach can find a more precise fit than the grid method, but it comes with hazards. The objective function has local minima, and results are sensitive to starting values. Researchers have documented cases where different starting points produce meaningfully different parameter estimates from the same day’s data.2EFMA. Estimating the Yield Curve Using the Nelson-Siegel Model: A Ridge Regression Approach Production systems often use the previous day’s parameters as starting values to keep the daily curve series continuous.

The output in either case is a smooth continuous function that spans from the shortest to the longest maturity. Analysts can interpolate yields for maturities without active trading, price derivatives against it, and compare rate expectations across dates or markets.

Known Limitations

The model’s simplicity creates real constraints, and ignoring them leads to poor pricing or misleading risk estimates.

The three-factor structure can produce at most one hump or one trough. Some market conditions create curves with two distinct humps, and the original framework simply cannot replicate those shapes. Forcing it to try typically distorts the fit at both ends.

Time series of the estimated parameters have been documented as very unstable, particularly when Lambda is estimated freely rather than fixed. The instability is partly driven by collinearity between the slope and curvature loading functions: at certain Lambda values, the two functions become nearly identical, and the optimizer trades off wildly between Beta 1 and Beta 2 without much change in the fitted curve.2EFMA. Estimating the Yield Curve Using the Nelson-Siegel Model: A Ridge Regression Approach

The model does not guarantee the absence of arbitrage opportunities, a limitation formally demonstrated by Björk and Christensen in 1999.3European Central Bank. How Arbitrage-Free Is the Nelson-Siegel Model? The curve fitted today can be inconsistent with the yields it would imply for the future. For descriptive purposes like fitting today’s term structure, this rarely matters. For pricing path-dependent derivatives or running a relative-value strategy, the lack of internal consistency is a genuine risk.

Under certain parameter combinations, the model can produce negative yields at long maturities, which violates basic economic logic in most environments. This typically happens when the estimation wanders into an unstable region of the parameter space.

Major Extensions

Each of those limitations has motivated at least one extension. Three are widely used.

Nelson-Siegel-Svensson

Lars Svensson’s 1994 extension adds a second curvature term with its own Beta (Beta 3) and its own decay parameter (Lambda 2). The second hump factor lets the model capture a broader range of shapes, including curves with two humps or an S-shaped twist the original three factors cannot reproduce.2EFMA. Estimating the Yield Curve Using the Nelson-Siegel Model: A Ridge Regression Approach The trade-off is two extra parameters that amplify the collinearity and stability problems. Even so, the Svensson version is the dominant choice among central banks, adopted by Germany, Switzerland, Norway, Sweden, and the United States, among others.1Bank for International Settlements. Zero-Coupon Yield Curves Estimated by Central Banks The European Central Bank also uses the Svensson extension for its daily euro area yield curve estimates.4European Central Bank. Euro Area Yield Curves

Diebold-Li Dynamic Model

Diebold and Li (2006) reinterpreted the three factors as time-varying latent variables rather than static parameters estimated day by day. They fit the model cross-sectionally at each date with Lambda fixed, then modeled how level, slope, and curvature evolve over time using autoregressive processes. That turns Nelson-Siegel from a curve-fitting tool into a forecasting model. Their results showed the dynamic factor models outperformed standard benchmarks, including random walks, at longer forecast horizons.5IDEAS/RePEc. Forecasting the Term Structure of Government Bond Yields The approach is now a staple in fixed-income research and central bank forecasting.

Arbitrage-Free Nelson-Siegel

Christensen, Diebold, and Rudebusch developed an arbitrage-free version, often abbreviated AFNS, that preserves the Nelson-Siegel factor loadings but adds a yield-adjustment term derived from no-arbitrage conditions. The adjustment depends only on maturity and the volatility structure of the factors, and it ensures that today’s fitted curve is internally consistent with the expected future evolution of yields.6University of Pennsylvania. An Arbitrage-Free Generalized Nelson-Siegel Term Structure Model AFNS is more complex to estimate than the original, but it bridges the gap between the Nelson-Siegel tradition and the rigorous no-arbitrage frameworks used in derivatives pricing.

Software to Run It

Open-source tools have made estimation accessible to anyone with basic programming experience. In Python, the nelson-siegel-svensson package on PyPI implements both the three-factor Nelson-Siegel and the four-factor Svensson variant, using OLS for the Betas and Nelder-Mead optimization from SciPy for the decay parameters.7PyPI. nelson-siegel-svensson 0.5.0 R users have equivalent options in the YieldCurve package on CRAN. Commercial platforms like Bloomberg and Refinitiv embed the fitting directly into their fixed-income analytics, so traders and portfolio managers can generate curves without writing code.

The workflow is the same regardless of the tool: feed in a vector of maturities and a corresponding vector of observed yields, decide whether to fix Lambda or optimize it, and let the solver run. The output is a set of parameters defining a continuous curve, from which any arbitrary maturity’s yield can be read off.