Consistency of Extended Nelson-Siegel Curve Families with the Ho-Lee and Hull and White Short Rate Models

Abstract

Nelson and Siegel curves are widely used to fit the observed term structure of interest rates in a particular date. By the other hand, several interest rate models have been developed such their initial forward rate curve can be adjusted to any observed data, as the Ho-Lee and the Hull and White one factor models. In this work we study the evolution of the forward curve process for each of these models assuming that the initial curve is of Nelson-Siegel type. We conclude that the forward curve process produces curves belonging to a parametric family of curves that can be seen as extended Nelson and Siegel curves. We show that the forward rate curve evolution has a linear or an exponential growth, depending on the particular short rate interest model. We applied the results to Argentinian short and forward rates obtained from the Lebac’s bills yields using the Hull and White short rate model, showing a good estimation of the observed forward rate curve for near dates when the initial forward curve is adjusted with a Nelson and Siegel one.

Share and Cite:

Kisbye, P. and Meier, K. (2017) Consistency of Extended Nelson-Siegel Curve Families with the Ho-Lee and Hull and White Short Rate Models. Journal of Mathematical Finance, 7, 919-930. doi: 10.4236/jmf.2017.74050.

1. Introduction

A standard procedure when dealing with concrete interest rate models is to calibrate the initial forward curve with the market observed data. That is the case of the Ho and Lee, and Hull and White models, where every curve can be perfectly fitted by adjusting the model parameters. By the other hand, some parametric curves are extendedly used to fit daily data, as is the case of the Nelson-Siegel curves [1] .

${f}_{NS}\left(\tau \right)={z}_{1}+{z}_{2}{\text{e}}^{-\lambda \tau }+{z}_{3}\tau {\text{e}}^{-\lambda \tau },\text{\hspace{0.17em}}\text{\hspace{0.17em}}\tau \ge 0$ (1)

with ${z}_{1},{z}_{2},{z}_{3}$ and $\lambda$ being specified parameters. So this means that it is possible to choose the Ho-Lee and Hull-White models parameters in such a way that the initial forward rate curve fits with a specific Nelson and Siegel curve. In this work we show that in this particular case, the following forward rate curves move on a manifold generated by specific parametric forward curves that can be written as a sum of a Nelson and Siegel curve and a linear or an exponential function, depending upon the short rate model. A Nelson-Siegel curve can be decomposed in three factors: 1, ${\text{e}}^{-\lambda \tau }$ and $\tau {\text{e}}^{-\lambda \tau }$ . The constant factor is related with the long term interest rate level. The exponential decay is the second factor, with an upward slope if ${z}_{2}>0$ or downward if ${z}_{2}<0$ . The third factor gives a hump or o a trough, depending on ${z}_{3}$ . Finally, λ is called the shape parameter, and it determines the critical point of the third factor and the steepness of the hump/trough. (See [2] ).

We prove that the forward rate curves produced by the Ho-Lee model and the Hull and White models when starting with a Nelson and Siegel curve is decomposed in four factors. Three of them are the same as in the Nelson and Siegel curve, and the fourth is a linear function (τ) in the Ho-Lee model or an exponential function $\left({c}_{1}{\text{e}}^{-a\tau }+{c}_{2}{\text{e}}^{-2a\tau }\right)$ in the Hull and White model, where a is a model parameter. This result extends part of Bjork and Christensen (see [3] ) paper results, where they proved that the above two models are inconsistent with a strictly Nelson-Siegel manifold.

In sections §3 and §4, we present the Ho-Lee and Hull-White models, and derive the formula for the corresponding forward curve. In each case, we choose a Nelson and Siegel curve as the initial forward curve and then prove that the following ones are extended Nelson-Siegel curves in the sense that they can be written as (1) plus a linear function or an exponential function.

In particular we also prove that each of these two short rate models is consistent with a forward curve manifold ${\mathcal{G}}^{\lambda }$ , for each λ > 0.

1.1. Notation and Facts

In this section we assume a probability space $\left(\Omega ,\mathcal{F},\mathcal{Q}\right)$ . Let $W\left(t\right),t\ge 0$ be a Wiener process, and ${\left\{{\mathcal{F}}_{t}\right\}}_{t\ge 0}$ be the filtration generated by $W\left(t\right)$ . Let $\left(\Omega ,\mathcal{F},\mathcal{Q},{\left\{{\mathcal{F}}_{t}\right\}}_{t\ge 0}\right)$ denote the filtered probability space. A stochastic process $\alpha$ is called an adapted process if $\alpha \left(t\right)$ is ${\mathcal{F}}_{t}$ -measurable for every $t\ge 0$ . An Ito process with drift $\mu$ and volatility $\sigma$ is a stochastic process $X\left(t\right)$ such that

$X\left(t\right)=X\left(0\right)+{\int }_{0}^{t}\mu \left(s,X\left(s\right)\right)\text{d}s+{\int }_{0}^{t}\sigma \left(s,X\left(s\right)\right)\text{d}W\left(s\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}t\ge 0,$ (2)

where $\mu$ and $\sigma$ are adapted process and the second integral at the right hand side is an Ito integral [4] . Equation (2) is usually written in terms of a stochastic differential equation

$\text{d}X\left(t\right)=\mu \left(t,X\left(t\right)\right)\text{d}t+\sigma \left(t,X\left(t\right)\right)\text{d}W\left(t\right)$ (3)

We also introduce the relationship between the Ito and the Stratonovich integral forms. If X is an Ito process as in (3), then its Stratonovich integral form is as follows,

$\text{d}X\left(t\right)=\left(\mu \left(t,X\left(t\right)\right)+\varphi \left(t,X\left(t\right)\right)\right)\text{d}t+\sigma \left(t,X\left(t\right)\right)o\text{d}W\left(t\right)$ (4)

where $o$ denotes the Stratonovich integration an $\varphi \left(t,X\left(t\right)\right)$ is a quadratic co-variance term, (see [4] ). If $\sigma$ in (3) is deterministic, then $\varphi \left(t,X\left(t\right)\right)=0$ .

We assume the existence of zero coupon bond market $\left\{P\left(t,T\right),0\le t\le T\right\}$ , where $P\left(t,T\right)$ denotes the price at time t of a bond with maturity T. We call this bond a T-bond. We assume that for each t, the curve $T\to P\left(t,T\right)$ is differentiable, with positive values and that $\left\{P\left(t,T\right),t\ge 0\right\}$ follows an Ito process for each $T\ge 0$ . The forward rate curve associated to these bonds is given by

$f\left(t,T\right)=-\frac{\partial \mathrm{ln}P\left(t,T\right)}{\partial T},$

and the short interest rate is given by $r\left(t\right)=f\left(t,t\right)$ . We shall suppose an arbitrage free model and we denote with $\mathcal{Q}$ the corresponding martingale measure. Every Ito process shall be described in terms of the $\mathcal{Q}$ -measure. As in the Heath, Jarrow y Morton (HJM) [5] framework, we assume that the forward rate curve dynamics is given by a family of stochastic differential equations, which expression under $\mathcal{Q}$ is given by

$\text{d}f\left(t,T\right)=\alpha \left(t,T\right)\text{d}t+{\sigma }_{0}\left(t,T\right)\text{d}W\left(t\right),$ (5)

with $\alpha$ and ${\sigma }_{0}$ adapted processes. The hypothesis of an arbitrage free market implies the HJM-drift condition on $\alpha$ . More precisely

$\alpha \left(t,T\right)={\sigma }_{0}\left(t,T\right){\int }_{t}^{T}{{\sigma }^{\prime }}_{0}\left(t,s\right)\text{d}s$ (6)

where the superscript in ${\sigma }^{\prime }$ denotes transpose in case of a vectorial process.

Given a T-bond, we denote $\tau =T-t$ the time up to maturity of the bond. The Brace and Musiela parametrization [6] describes the forward curve process in terms of t and $\tau$ as follows,

${f}_{r}\left(t,\tau \right)=f\left(t,t+\tau \right),$ (7)

and so $r\left(t\right)={f}_{r}\left(t,0\right)$ . Under this parametrization, Equation (5) can be written as

$\text{d}{f}_{r}\left(t,\tau \right)=\left(\frac{\partial }{\partial \tau }{f}_{r}\left(t,\tau \right)+\sigma \left(t,\tau \right){\int }_{0}^{\tau }{\sigma }^{\prime }\left(t,s\right)\text{d}s\right)\text{d}t+\sigma \left(t,\tau \right)\text{d}W\left(t\right),$ (8)

where $\sigma \left(t,\tau \right)={\sigma }_{0}\left(t,t+\tau \right)$ .

In particular, the forward rate process (8) can be expressed in terms of the Stratonovich integral form as

$\text{d}{f}_{r}\left(t,\tau \right)=\left(\frac{\partial }{\partial \tau }{f}_{r}\left(t,\tau \right)+\sigma \left(t,\tau \right){\int }_{0}^{\tau }{\sigma }^{\prime }\left(t,s\right)\text{d}s+\varphi \left(t,\tau \right)\right)\text{d}t+\sigma \left(t,\tau \right)o\text{d}W\left(t\right).$ (9)

1.2. Consistency

Consistency between short rate models and forward curves manifolds were stated by Bjork and Christensen in [3] . To make this work more self-contained, we recall some of their definitions and the main theorem. Let $\mathcal{M}$ a given one factor interest rate model specifying a forward rate process ${f}_{r}\left(t,\cdot \right)$ . In terms of the Musiela parametrization, ${f}_{r}$ satisfies a stochastic differential equation:

$\text{d}{f}_{r}\left(t,\tau \right)=\left(\frac{\partial }{\partial \tau }{f}_{r}\left(t,\tau \right)+\alpha \left(t,\tau \right)\right)\text{d}t+\sigma \left(t,\tau \right)\text{d}W\left(t\right),$ (10)

$t\ge 0,\tau \ge 0$ , where $\alpha$ and $\sigma$ are adapted processes. In particular, the no arbitrage Heath, Jarrow and Morton (HJM) drift condition implies that $\alpha \left(t,x\right)=\sigma \left(t,x\right){\int }_{t}^{x}\sigma \left(t,s\right)\text{d}s$ . Also, using the Stratonovich integral form, (10) can be written as

$\text{d}{f}_{r}\left(t,\tau \right)=\left(\frac{\partial }{\partial \tau }{f}_{r}\left(t,\tau \right)+\sigma \left(t,x\right){\int }_{t}^{\tau }\sigma \left(t,s\right)\text{d}s+\varphi \left(t,\tau \right)\right)\text{d}t+\sigma \left(t,\tau \right)o\text{d}W\left(t\right),$ (11)

Bjork and Christensen [3] stated the following definitions of consistency between a short rate model and a parametric family of curves. At first, let $\mathcal{Z}\subseteq {ℝ}^{d}$ be a set of parameters, and let $G:\mathcal{Z}\to C\left[0,\infty \right)$ be a smooth function. The forward curve manifold $\mathcal{G}$ is defined as $\mathcal{G}=Im\left(G\right)$ . That is,

$\mathcal{G}=\left\{G\left(\cdot ;z\right):\left[0,\infty \right)\to ℝ\right\},$

where, with some abuse of notation $G\left(\cdot ;z\right)$ denotes the function $G\left(z\right)$ .

Definition 1. (invariance). Consider a given interest rate model $\mathcal{M}$ , specifying a forward rate process ${f}_{r}\left(t,\cdot \right)$ , and a forward curve manifold $\mathcal{G}$ . We say that $\mathcal{G}$ is invariant under the action of ${f}_{r}$ if, for every fixed initial time $s$ , the condition ${f}_{r}\left(s,\cdot \right)\in \mathcal{G}$ implies that ${f}_{r}\left(t,\cdot \right)\in \mathcal{G}$ , for all $t\ge s$ , a.s.

Bjork and Christensen also stated a more restricted concept of invariance, the ${f}_{r}$ -invariance.

Definition 2. ( ${f}_{r}$ -invariance). Consider a given interest rate model $\mathcal{M}$ , specifying a forward rate process ${f}_{r}\left(t,\cdot \right)$ as in (11), as well as a forward curve manifold $\mathcal{G}$ . We say that $\mathcal{G}$ is -invariant under the action of the forward rate process ${f}_{r}\left(t,\cdot \right)$ if there exists a stochastic process $Z$ with state process $\mathcal{Z}$ and possessing a Stratonovich differential of the form $\text{d}Z\left(t\right)=\gamma \left(t,Z\left(t\right)\right)\text{d}t+\psi \left(t,Z\left(t\right)\right)o\text{d}W\left(t\right)$ , such that, for every fixed choice of initial time $s$ , whenever $y\left(s,\cdot \right)\in \mathcal{G}$ , the stochastic process defined by

$y\left(t,\tau \right)=G\left(\tau ,Z\left(t\right)\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall t\ge s,\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\ge 0,$

satisfies the SDE (11) with initial condition ${f}_{r}\left(s,\cdot \right)=y\left(s,\cdot \right)$ .

In this case, we say that the short rate model $\mathcal{M}$ and the manifold $\mathcal{G}$ are consistent. It is easy to prove that ${f}_{r}$ -invariance implies invariance. Moreover, Bjork and Christensen proved the following theorem.

Theorem 3. The forward curve manifold G is ${f}_{r}$ -invariant for the forward rate process ${f}_{r}\left(t,\cdot \right)$ in $\mathcal{M}$ if and only if

${G}_{\tau }\left(\cdot ,z\right)+\sigma \left(t,\cdot \right){\int }_{0}^{\left(\cdot \right)}{\sigma }^{\prime }\left(t,s\right)\text{d}s+\varphi \left(t,\cdot \right)\in Im\left[{G}_{z}\left(\cdot ;z\right)\right]$ (12)

$\sigma \left(t,\cdot \right)\in Im\left[{G}_{z}\left(\cdot ;z\right)\right]$ (13)

for all $\left(t,z\right)\in \left[0,\infty \right)×\mathcal{Z}$ . ${G}_{\tau }$ and ${G}_{z}$ denote the Frechet derivatives of G with respect to $\tau$ and $z$ , which are assumed to exist.

Definition 4. An interest rate model $\mathcal{M}$ is consistent with the forward rate manifold G if the consistent drift and volatility Conditions (12)-(13) hold.

2. The Ho-Lee Short Rate Model

The short rate model proposed by Ho and Lee [7] , (henceforth HL) has a dynamic given by the r-process

$\text{d}r\left(t\right)=\theta \left(t\right)\text{d}t+\sigma \text{d}W\left(t\right)$ (14)

In (14), $\left\{W\left(t\right),t\ge 0\right\}$ is a Wiener process, $\sigma >0$ and $\theta$ is a deterministic function. The HL model belongs to the family of affine short rate models. That is, if $P\left(t,T\right)$ denotes the price at time t of a zero coupon bond with maturity T, then the term structure of the interest rate is given by:

$P\left(t,T\right)={\text{e}}^{A\left(t,T\right)-r\left(t\right)B\left(t,T\right)},$ (15)

for certain functions A and B. In particular, in the case of the Ho-Lee model, A and B are given by:

$B\left(t,T\right)=T-t$ (16)

$A\left(t,T\right)={\int }_{t}^{T}\theta \left(s\right)\left(s-T\right)\text{d}s+\frac{{\sigma }^{2}{\left(T-t\right)}^{3}}{6}$ (17)

(see for instance [8] ). The forward rate curve is related with the term structure by the equation:

$f\left(t,T\right)=-\frac{\partial \mathrm{ln}\left(P\left(t,T\right)\right)}{\partial T}=-\frac{\partial A\left(t,T\right)}{\partial T}+\frac{\partial B\left(t,T\right)r\left(t\right)}{\partial T}$ (18)

Then, in this case replacing A and B by the expressions in (17) and (16), we get:

$\begin{array}{c}f\left(t,T\right)=\frac{\partial }{\partial T}\left(-{\int }_{t}^{T}\theta \left(s\right)\left(s-T\right)\text{d}s-\frac{{\sigma }^{2}{\left(T-t\right)}^{3}}{6}\right)+r\left(t\right)\frac{\partial }{\partial T}\left(T-t\right)\\ =-\left(\theta \left(T\right)T-{\int }_{t}^{T}\theta \left(s\right)\text{d}s-T\theta \left(T\right)\right)-\frac{{\sigma }^{2}}{2}{\left(T-t\right)}^{2}+r\left(t\right)\\ ={\int }_{t}^{T}\theta \left(s\right)\text{d}s-\frac{{\sigma }^{2}}{2}{\left(T-t\right)}^{2}+r\left(t\right)\end{array}$ (19)

In particular it holds that $f\left(t,t\right)=r\left(t\right)$ . If $T\to {f}^{*}\left(0,T\right)$ is the observed initial forward curve, and $\theta$ is defined as

$\theta \left(t\right)={\sigma }^{2}t+\frac{\partial {f}^{*}}{\partial T}\left(0,t\right),$

then $f\left(0,T\right)={f}^{*}\left(0,T\right)$ . That is, the model parameters can be adjusted such that the initial forward curve fits the observed one.

We now assume that the initial forward curve is given by a Nelson and Siegel parametric curve. That is, we define:

${f}^{*}\left(0,T\right)={z}_{1}+{z}_{2}{\text{e}}^{-\lambda T}+{z}_{3}T{\text{e}}^{-\lambda T},\text{\hspace{0.17em}}\text{\hspace{0.17em}}T\ge 0$ (20)

where ${z}_{1},{z}_{2},{z}_{3}$ and $\lambda$ are fixed real numbers, $\lambda >0$ . We want to study the evolution of this initial curve in the $t$ variable. With this particular choice of ${f}^{*}$ , $\theta$ is given by:

$\theta \left(t\right)={\sigma }^{2}t+\left({z}_{3}-{z}_{2}\lambda \right){\text{e}}^{-\lambda t}-{z}_{3}\lambda t{\text{e}}^{-\lambda t}$

and the solution of the Ho-Lee stochastic differential Equation (14) is given by:

$\begin{array}{c}r\left(t\right)=r\left(0\right)+\frac{{\sigma }^{2}{t}^{2}}{2}-\frac{{z}_{3}-{z}_{2}\lambda }{\lambda }{\text{e}}^{-\lambda t}-\frac{{z}_{3}}{\lambda }\left(1-{\text{e}}^{-\lambda t}\left(\lambda t+1\right)\right)+\sigma W\left(t\right)\\ =r\left(0\right)+\frac{{\sigma }^{2}{t}^{2}}{2}-\frac{{z}_{3}}{\lambda }+\left({z}_{2}+{z}_{3}t\right){\text{e}}^{-\lambda t}+\sigma W\left(t\right)\end{array}$

where $r\left(0\right)$ is the short rate value at time $t=0$ .

We compute the integral term in (19):

$\begin{array}{l}{\int }_{t}^{T}\theta \left(s\right)\text{d}s={\int }_{t}^{T}{\sigma }^{2}s+\left({z}_{3}-{z}_{2}\lambda \right){\text{e}}^{-\lambda s}-{z}_{3}\lambda t{\text{e}}^{-\lambda s}\text{d}s\\ =\frac{{\sigma }^{2}}{2}\left({T}^{2}-{t}^{2}\right)+\left({z}_{2}+{z}_{3}T\right){\text{e}}^{-\lambda T}-\left({z}_{2}+{z}_{3}t\right){\text{e}}^{-\lambda t}\end{array}$

We can now derive an explicit formula for the forward rate curve process:

$\begin{array}{c}f\left(t,T\right)=\frac{{\sigma }^{2}}{2}\left({T}^{2}-{t}^{2}\right)+\left({z}_{2}+{z}_{3}T\right){\text{e}}^{-\lambda T}-\left({z}_{2}+{z}_{3}t\right){\text{e}}^{-\lambda t}-\frac{{\sigma }^{2}}{2}{\left(T-t\right)}^{2}+r\left(t\right)\\ ={\sigma }^{2}t\left(T-t\right)+\left({z}_{2}+{z}_{3}T\right){\text{e}}^{-\lambda T}-\left({z}_{2}+{z}_{3}t\right){\text{e}}^{-\lambda t}+r\left(t\right)\end{array}$ (21)

The above computations allow us to state the following theorem.

Theorem 5. Let r denote the HL short rate with a dynamic as stated in (14). Then, if the initial forward rate curve is a Nelson and Siegel parametric curve as in (20), the corresponding forward rate curve at time t is given by the formula

$f\left(t,T\right)=r\left(t\right)+{\sigma }^{2}t\left(T-t\right)+\left({z}_{2}+{z}_{3}T\right){\text{e}}^{-\lambda T}-\left({z}_{2}+{z}_{3}t\right){\text{e}}^{-\lambda t}$

Let $\tau =T-t$ the time up to maturity. Using the Brace and Musiela parametrization, we denote ${f}_{HL}$ the forward curve given by ${f}_{HL}\left(t,\tau \right)=f\left(t,t+\tau \right)$ . Then the forward rate curve ${f}_{HL}$ has the expression

${f}_{HL}\left(t,\tau \right)={\sigma }^{2}t\tau +{C}_{1}\left(t\right)+{C}_{2}\left(t\right){\text{e}}^{-\lambda \tau }+{C}_{3}\left(t\right)\tau {\text{e}}^{-\lambda \tau },\text{\hspace{0.17em}}\text{\hspace{0.17em}}t\ge 0,\text{\hspace{0.17em}}\tau \ge 0,$ (22)

where ${C}_{1},{C}_{2},{C}_{3}$ are coefficients that depends on t and the Nelson and Siegel parameters:

${C}_{1}\left(t\right)=r\left(t\right)-\left({z}_{2}+{z}_{3}t\right){\text{e}}^{-\lambda t}$

${C}_{2}\left(t\right)=\left({z}_{2}+{z}_{3}t\right){\text{e}}^{-\lambda t}$

${C}_{3}\left(t\right)={z}_{3}{\text{e}}^{-\lambda t}$

Proof. The proof follows arranging terms after replacing $T$ by $\tau +t$ in Equation (21).

The expression of the function $\tau \to {f}_{HL}\left(t,\tau \right)$ given in the formula (22) is a sum of a linear function plus a Nelson and Siegel parametric curve.

Definition 6. Let $\lambda >0$ and $g:\left[0,\infty \right)\to ℝ$ be a function defined as

$g\left(\tau \right)={z}_{0}\tau +{z}_{1}+{z}_{2}{\text{e}}^{-\lambda t}+{z}_{3}{\text{e}}^{-\lambda t},$

where ${z}_{0},{z}_{1}$ and ${z}_{3}$ are constant real numbers. We call g a linearly extended Nelson-Siegel curve.

In particular, in the following subsection we study the consistency of the Ho-Lee model with a family of forward curve manifolds ${\mathcal{G}}^{\lambda }$ generated by linearly extended Nelson-Siegel curves.

Consistency between the HL Model and Forward Curve Manifolds

The forward rate curve f given in (21) satisfies the Heath, Jarrow and Morton drift condition:

$\alpha \left(t,T\right)=\sigma \left(t,T\right){\int }_{t}^{T}\sigma \left(t,s\right)\text{d}s={\sigma }^{2}\left(T-t\right),$

and the stochastic differential equation

$\text{d}f\left(t,T\right)=\alpha \left(t,T\right)\text{d}t+\sigma \left(t,T\right)\text{d}W\left(t\right)={\sigma }^{2}\left(T-t\right)\text{d}t+\sigma \text{d}W\left(t\right),$

with $T\ge 0$ and $0\le t\le T$ . In terms of the Brace and Musiela parametrization, fHL satisfies:

$\text{d}{f}_{HL}\left(t,\tau \right)=\left(\frac{\partial }{\partial \tau }{f}_{HL}\left(t,\tau \right)+{\sigma }^{2}\tau \right)\text{d}t+\sigma \text{d}W\left(t\right),$ (23)

$t\ge 0,\tau \ge 0$ . Starting from Equations (22) and (23), our conjecture is that there exists a forward curve manifold containing linearly extended Nelson-Siegel curves that is consistent with the HL model. In fact, this is stated and proved by the following theorem.

Theorem 7. Let $\mathcal{Z}={ℝ}^{4},\lambda >0$ and

${G}^{\lambda }\left(\tau ;\beta \right)={\beta }_{0}\tau +{\beta }_{1}+{\beta }_{2}{\text{e}}^{-\lambda \tau }+{\beta }_{3}\tau {\text{e}}^{-\lambda \tau },\text{\hspace{0.17em}}\text{\hspace{0.17em}}\tau \ge 0$

Let ${f}_{HL}\left(t,.\right)$ be the Ho-Lee forward rate process. Then, for each $\lambda >0$ the forward curve manifold ${\mathcal{G}}^{\lambda }$ is ${f}_{HL}$ -invariant.

Proof. We shall apply Theorem 3 to see that ${\mathcal{G}}^{\lambda }$ is ${f}_{HL}$ -invariant. Because the volatility term in (23) is deterministic, the standard differential equation is the same for the Ito and the Stratonovich integral formulation. The Frechet derivatives of ${G}^{\lambda }$ are given by:

${G}_{\beta }^{\lambda }\left(\tau ,\beta \right)=\left[\tau ,1,{\text{e}}^{-\lambda \tau },\tau ,{\text{e}}^{-\lambda \tau }\right]$

${G}_{\tau }^{\lambda }\left(\tau ,\beta \right)={\beta }_{0}+\left(-{\beta }_{2}\lambda +{\beta }_{3}\right){\text{e}}^{-\lambda \tau }-{\beta }_{3}\lambda \tau {\text{e}}^{-\lambda \tau }$

In order to prove that ${\mathcal{G}}^{\lambda }$ is ${f}_{HL}$ -invariant, we must check the drift and volatility consistency Conditions (12) and (13). We shall first prove that ${G}_{\tau }^{\lambda }\left(\cdot ,\beta \right)+{\sigma }^{2}\left(\cdot \right)\in Im\left[{G}_{\beta }^{\lambda }\left(\cdot ,\beta \right)\right]$ . This means that there must be real numbers $A,B,C$ and $D$ such that:

${\beta }_{0}+\left(-{\beta }_{2}\lambda +{\beta }_{3}\right){\text{e}}^{-\lambda \tau }-{\beta }_{3}\lambda \tau {\text{e}}^{-\lambda \tau }+{\sigma }^{2}\tau =A\tau +B+C{\text{e}}^{-\lambda \tau }+D\tau {\text{e}}^{-\lambda \tau }$

In fact, this is possible taking

$A={\sigma }^{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}B={\beta }_{0},\text{\hspace{0.17em}}\text{\hspace{0.17em}}C=-{\beta }_{2}\lambda +{\beta }_{3}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}D=-{\beta }_{3}\lambda$

so Condition (12) is satisfied. To prove Condition (13), we must find $A,B,C$ and $D$ such that

$\sigma =A\tau +B+C{\text{e}}^{-\lambda \tau }+D\tau {\text{e}}^{-\lambda \tau }$

and this can be done taking $A=B=D=0$ and $B=\sigma$ .

Theorem 7 implies that the forward curve manifold ${G}^{\lambda }$ is ${f}_{HL}$ -invariant, so the next corollary follows:

Corollary 1. For every $\lambda >0$ , the forward curve manifold ${G}^{\lambda }$ is consistent with the Ho-Lee short rate model.

Corollary 1 implies that, in the particular case that $\theta$ in (14) is chosen such that the initial forward rate curve fits the Nelson and Siegel curve

${f}^{*}\left(0,\tau \right)={z}_{1}+{z}_{2}{\text{e}}^{-\lambda \tau }+{z}_{3}\tau {\text{e}}^{-\lambda \tau },\text{\hspace{0.17em}}\text{\hspace{0.17em}}\tau \ge 0$

then, for each $t\ge 0$ the corresponding forward rate curve ${f}_{HL}\left(t,\tau \right)$ can be written as a linearly extended Nelson and Siegel curve. That is:

${f}_{HL}\left(t,\tau \right)={\beta }_{0}\tau +{\beta }_{1}+{\beta }_{2}{\text{e}}^{-{\beta }_{4}\tau }+{\beta }_{3}\tau {\text{e}}^{-{\beta }_{4}\tau },\text{\hspace{0.17em}}\text{\hspace{0.17em}}\tau \ge 0$

then the drift consistency condition is not satisfied, because it requires ${\beta }_{4}=\lambda$ .

3. The Hull and White Model

The short rate model proposed by Hull and White [9] , (henceforth HW) or extended Vasicek model has the following stochastic differential equation:

$\text{d}r\left(t\right)=\left(\theta \left(t\right)-ar\left(t\right)\right)\text{d}t+\sigma \text{d}W\left(t\right)$ (24)

where $a,\sigma$ are positive real numbers and $\theta$ is a deterministic function. $\theta$ can be chosen in such a way that the initial forward curve $f\left(0,\cdot \right)$ fits with the observed data at $t=0$ . Let ${f}^{*}\left(0,\cdot \right)$ be a particular forward curve observed at $t=0$ . Then $\theta$ is defined as:

$\theta \left(t\right)=\frac{\partial {f}^{\text{*}}}{\partial T}\left(0,t\right)+a{f}^{\text{*}}\left(0,t\right)+\frac{{\sigma }^{2}}{2a}\left(1-{\text{e}}^{-2at}\right)$ (25)

The Hull and White short rate model belongs also to the class of affine models, and in this case the functions A and B in Equation (15) are given by [8] :

$B\left(t,T\right)=\frac{1}{a}\left(1-{\text{e}}^{-a\left(T-t\right)}\right)$ (26)

$A\left(t,T\right)={\int }_{t}^{T}\left(\frac{1}{2}{\sigma }^{2}{B}^{2}\left(s,T\right)-\theta \left(s\right)B\left(s,T\right)\right)\text{d}s$ (27)

The corresponding forward rate curve is given by:

$\begin{array}{c}f\left(t,T\right)=-\frac{\partial A\left(t,T\right)}{\partial T}+\frac{\partial B\left(t,T\right)r\left(t\right)}{\partial T}\\ =-\frac{\partial }{\partial T}\left({\int }_{t}^{T}\left(\frac{1}{2}{\sigma }^{2}{B}^{2}\left(s,T\right)-\theta \left(s\right)B\left(s,T\right)\right)\text{d}s\right)+r\left(t\right){\text{e}}^{-a\left(T-t\right)}\end{array}$ (28)

Now, by Leibniz rule, we have:

$\begin{array}{c}\frac{\partial }{\partial T}\left({\int }_{t}^{T}\frac{1}{2}{\sigma }^{2}{B}^{2}\left(s,T\right)\text{d}s\right)=\frac{1}{2}{\sigma }^{2}{B}^{2}\left(T,T\right)+{\int }_{t}^{T}\frac{1}{2}{\sigma }^{2}\frac{\partial {B}^{2}\left(s,T\right)}{\partial T}\text{d}s\\ ={\sigma }^{2}{\int }_{t}^{T}B\left(s,T\right)\frac{\partial B\left(s,T\right)}{\partial T}\text{d}s\\ ={\sigma }^{2}{\int }_{t}^{T}\frac{1-{\text{e}}^{-a\left(T-s\right)}}{a}{\text{e}}^{-a\left(T-s\right)}\text{d}s\\ =-\frac{1}{2}{\sigma }^{2}{B}^{2}\left(t,T\right)=-\frac{{\sigma }^{2}}{2{a}^{2}}{\left(1-{\text{e}}^{-a\left(T-t\right)}\right)}^{2}\end{array}$ (29)

$\begin{array}{c}\frac{\partial }{\partial T}\left({\int }_{t}^{T}\theta \left(s\right)B\left(s,T\right)\right)=\theta \left(T\right)B\left(T,T\right)+{\int }_{t}^{T}\theta \left(s\right)\frac{\partial B\left(s,T\right)}{\partial T}\text{d}s\\ ={\int }_{t}^{T}\theta \left(s\right){\text{e}}^{-a\left(T-s\right)}\text{d}s\end{array}$

We now assume that the initial forward curve is fitted to a Nelson and Siegel parametric curve,

${f}^{*}\left(0,T\right)={z}_{1}+{z}_{2}{\text{e}}^{-\lambda T}+{z}_{3}T{\text{e}}^{-\lambda T},\text{\hspace{0.17em}}\text{\hspace{0.17em}}T\ge 0$

Then, the solution of (24) and the function $\theta$ are given by the following expressions:

$\begin{array}{c}r\left(t\right)=r\left(0\right){\text{e}}^{-at}+{\int }_{0}^{t}{\text{e}}^{-a\left(t-u\right)}\theta \left(u\right)\text{d}u+\sigma {\int }_{0}^{t}{\text{e}}^{-a\left(t-u\right)}\text{d}W\left(u\right)\\ =r\left(0\right){\text{e}}^{-at}+\alpha \left(t\right)-\alpha \left(0\right){\text{e}}^{-at}+\sigma {\int }_{0}^{t}{\text{e}}^{-a\left(t-u\right)}\text{d}W\left(u\right)\end{array}$

$\theta \left(t\right)=a{z}_{1}+\left({z}_{3}-{z}_{2}\lambda +a{z}_{2}\right){\text{e}}^{-\lambda t}+\left(a{z}_{3}-{z}_{3}\lambda \right)t{\text{e}}^{-\lambda t}+\frac{{\sigma }^{2}}{2a}\left(1-{\text{e}}^{-2at}\right)$

where

$\alpha \left(t\right)={f}^{*}\left(0,t\right)+\frac{{\sigma }^{2}}{2{a}^{2}}{\left(1-{\text{e}}^{-at}\right)}^{2}.$

Now we can compute explicitly

$\begin{array}{c}{\int }_{t}^{T}\theta \left(s\right){\text{e}}^{-a\left(T-s\right)}\text{d}s=\alpha \left(T\right)-\alpha \left(t\right){\text{e}}^{-a\left(T-t\right)}\\ ={f}^{*}\left(0,T\right)+\frac{{\sigma }^{2}}{2{a}^{2}}{\left(1-{\text{e}}^{-aT}\right)}^{2}-\alpha \left(t\right){\text{e}}^{-a\left(T-t\right)}\end{array}$ (30)

Replacing the Expressions (29) and (30) in the forward rate curve formula (28), we get:

$f\left(t,T\right)=-\frac{{\sigma }^{2}}{2{a}^{2}}{\left(1-{\text{e}}^{-a\left(T-t\right)}\right)}^{2}+{f}^{\text{*}}\left(0,T\right)+\frac{{\sigma }^{2}}{2{a}^{2}}{\left(1-{\text{e}}^{-aT}\right)}^{2}-\alpha \left(t\right){\text{e}}^{-a\left(T-t\right)}+r\left(t\right){\text{e}}^{-a\left(T-t\right)}$

With the above computations we arrive to the next theorem.

Theorem 8. Let r denote the Hull and White short rate model with the dynamics stated in (24). Then, if the initial forward rate curve is a Nelson and Siegel parametric curve as in (20), the corresponding forward rate curve at time t is given by the formula

$f\left(t,T\right)=-\frac{{\sigma }^{2}}{2{a}^{2}}{\left(1-{\text{e}}^{-a\left(T-t\right)}\right)}^{2}+{f}^{*}\left(0,T\right)+\frac{{\sigma }^{2}}{2{a}^{2}}{\left(1-{\text{e}}^{-aT}\right)}^{2}-\alpha \left(t\right){\text{e}}^{-a\left(T-t\right)}+r\left(t\right){\text{e}}^{-a\left(T-t\right)}$

Let $\tau =T-t$ the time up to maturity. Using the Brace and Musiela parametrization, we denote ${f}_{HW}$ the forward curve given by ${f}_{HW}\left(t,\tau \right)=f\left(t,t+\tau \right)$ . Then the forward rate curve ${f}_{HW}$ has the expression

${f}_{HW}\left(t,\tau \right)={C}_{1}\left(t\right){\text{e}}^{-a\tau }+{C}_{2}\left(t\right){\text{e}}^{-2a\tau }+{C}_{3}\left(t\right)+{C}_{4}\left(t\right){\text{e}}^{-\lambda \tau }+{C}_{5}\left(t\right)\tau {\text{e}}^{-\lambda \tau }$ (31)

where ${C}_{1},{C}_{2},{C}_{3},{C}_{4}$ and ${C}_{5}$ are coefficients that depends on $t,r\left(t\right)$ and the Nelson and Siegel curve parameters:

${C}_{1}\left(t\right)=\frac{{\sigma }^{2}}{{a}^{2}}\left(1-{\text{e}}^{-at}\right)-\alpha \left(t\right)+r\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{2}\left(t\right)=\frac{{\sigma }^{2}}{2{a}^{2}}\left({\text{e}}^{-2at}-1\right),$

${C}_{3}\left(t\right)={z}_{1},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{4}\left(t\right)={z}_{2}{\text{e}}^{-\lambda t}+{z}_{3}t{\text{e}}^{-\lambda t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{5}\left(t\right)={z}_{3}{\text{e}}^{-\lambda t}\text{.}$

The expression of the function $\tau \to {f}_{HW}\left(t,\tau \right)$ given in the formula (31) is a sum of an exponential function plus a Nelson and Siegel parametric curve.

Definition 9. Let $\lambda >0$ and $g:\left[0,\infty \right)\to ℝ$ be a function defined as

$g\left(\tau \right)={c}_{1}{\text{e}}^{-a\tau }+{c}_{2}{\text{e}}^{-2a\tau }+{z}_{0}\tau +{z}_{1}+{z}_{2}{\text{e}}^{-\lambda \tau }+{z}_{3}{\text{e}}^{-\lambda \tau },$

with ${c}_{1},{c}_{2},{z}_{0},{z}_{1}$ and ${z}_{2}$ constant real numbers. We call g an exponentially extended Nelson-Siegel curve.

In particular, in the following subsection we study the consistency of the Hull and White model with a family of forward curve manifolds ${\mathcal{G}}^{\lambda }$ generated by exponentially extended Nelson-Siegel curves.

Consistency between the Hull and White Model and Forward Curve Manifolds

The dynamics of the forward curve process is given in terms of the Ito integral formulation by

$\text{d}f\left(t,T\right)=\frac{{\sigma }^{2}}{a}{\text{e}}^{-a\left(T-t\right)}\left(1-{\text{e}}^{-a\left(T-t\right)}\right)\text{d}t+\sigma {\text{e}}^{-a\left(T-t\right)}\text{d}W\left(t\right),$

$0\le t , and in the Musiela parametrization $\tau =T-t$ ,

$\text{d}{f}_{HW}\left(t,\tau \right)=\left(\frac{\partial }{\partial \tau }{f}_{HW}\left(t,\tau \right)+\frac{{\sigma }^{2}}{a}{\text{e}}^{-a\tau }\left(1-{\text{e}}^{-a\tau }\right)\right)\text{d}t+\sigma {\text{e}}^{-a\tau }\text{d}W\left(t\right)$ (32)

$\tau \ge 0,t\ge 0$ . Because the volatility term is a deterministic function, the Stratonovich and the Ito integral formulation are the same.

We state the next theorem:

Theorem 10. Let $\mathcal{Z}={ℝ}^{5}$ and for each $\lambda >0$ let ${G}^{\lambda }:\mathcal{Z}\to C\left(0,\infty \right)$ be defined as

${G}^{\lambda }\left(\tau ;\beta \right)={\beta }_{1}{\text{e}}^{-a\tau }+{\beta }_{2}{\text{e}}^{-2a\tau }+{\beta }_{3}+{\beta }_{4}{\text{e}}^{-\lambda \tau }+{\beta }_{5}\tau {\text{e}}^{-\lambda \tau }.$ (33)

Let ${f}_{HW}\left(t,\cdot \right)$ be the Hull and White forward rate process given by the standard differential Equation (32). Then, the forward curve manifold ${\mathcal{G}}^{\lambda }$ is ${f}_{HW}$ - invariant.

Proof. We note that the Frechet derivatives ${G}_{\beta }^{\lambda }$ and ${G}_{\tau }^{\lambda }$ are given by:

${G}_{\beta }^{\lambda }\left(\tau ;\beta \right)=\left[{\text{e}}^{-a\tau },{\text{e}}^{-2a\tau },1,{\text{e}}^{-\lambda \tau },\tau {\text{e}}^{-\lambda \tau }\right]$

${G}_{\tau }^{\lambda }\left(\tau ;\beta \right)={\beta }_{1}\left(-a\right){\text{e}}^{-a\tau }+{\beta }_{2}\left(-2a\right){\text{e}}^{-2a\tau }+\left(-{\beta }_{4}\lambda +{\beta }_{5}\right){\text{e}}^{-\lambda \tau }-{\beta }_{5}\lambda \tau {\text{e}}^{-\lambda \tau }.$

So, first we prove that ${G}_{\tau }^{\lambda }\left(\cdot ;\beta \right)+\frac{{\sigma }^{2}}{a}{\text{e}}^{-a\left(\cdot \right)}\left(1-{\text{e}}^{-a\left(\cdot \right)}\right)\in Im\left[{G}_{\beta }^{\lambda }\left(\cdot ;\beta \right)\right]$ . So we

look for real numbers A, B, C, D and E such that

$\begin{array}{l}{\beta }_{1}\left(-a\right){\text{e}}^{-a\tau }+{\beta }_{2}\left(-2a\right){\text{e}}^{-2a\tau }+\left(-{\beta }_{4}\lambda +{\beta }_{5}\right){\text{e}}^{-\lambda \tau }-{\beta }_{5}\lambda \tau {\text{e}}^{-\lambda \tau }+\frac{{\sigma }^{2}}{a}{\text{e}}^{-a\tau }\left(1-{\text{e}}^{-a\tau }\right)\\ =A{\text{e}}^{-a\tau }+B{\text{e}}^{-2a\tau }+C+D{\text{e}}^{-\lambda \tau }+E\tau {\text{e}}^{-\lambda \tau }.\end{array}$

This is true setting

$A=-\alpha {\beta }_{1}+\frac{{\sigma }^{2}}{a},\text{\hspace{0.17em}}\text{\hspace{0.17em}}B=-2a{\beta }_{2}-\frac{{\sigma }^{2}}{a},\text{\hspace{0.17em}}\text{\hspace{0.17em}}C=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}D=-{\beta }_{4}\lambda +{\beta }_{5},\text{\hspace{0.17em}}\text{\hspace{0.17em}}E=-{\beta }_{5\lambda }\text{.}$

We next prove that $\sigma {\text{e}}^{-a\left(\cdot \right)}\in Im\left[{G}_{\beta }\left(\cdot ,\beta \right)\right]$ , or equivalently, we look for real numbers such that

$\sigma {\text{e}}^{-a\tau }=A{\text{e}}^{-a\tau }+B{\text{e}}^{-2a\tau }+C+D{\text{e}}^{-\lambda \tau }+E\tau {\text{e}}^{-\lambda \tau }$

Setting $A=\sigma$ and $B=C=D=E=0$ the identity follows.

As a conclusion, we have the following corollary:

Corollary 2. For each λ > 0, the Hull and White short rate model is consistent with the forward curve manifold ${\mathcal{G}}^{\lambda }$ given in (33).

Remark. As in the case of the HL-model, we see that if the initial forward curve is a Nelson-Siegel curve, then the following forward curves are exponentially extended Nelson-Siegel curves, belonging to a particular manifold ${\mathcal{G}}^{\lambda }$ . It must be noted that, even when the forward curve process moves on a wider manifold $\mathcal{H}$ with parameter set $\mathcal{Z}={ℝ}^{6}$ ,

$H\left(\tau ;\beta \right)={\beta }_{1}{\text{e}}^{-a\tau }+{\beta }_{2}{\text{e}}^{-2a\tau }+{\beta }_{3}+{\beta }_{4}{\text{e}}^{-{\beta }_{6}\tau }+{\beta }_{5}\tau {\text{e}}^{-{\beta }_{6}\tau },\text{\hspace{0.17em}}\text{\hspace{0.17em}}\tau \ge 0,$

it is not true that this particular manifold is consistent with the HW model. The parameter ${\beta }_{6}$ must be equal to the parameter λ in (1) to get the drift consistency condition (12).

4. Application to Argentina’s Interest Rate

In this section we apply the Hull-White model to the Argentinian Central Bank bills, named Lebac (Letras del Banco Central). There are seven maturities for the Lebac bills: 35, 63, 98, 120, 147, 203 and 252 days. We have taken the 35d-maturity Lebacs as the short interest rate and we roughly estimated the forward rate curve using the other six maturities:

$f\left(0,{t}_{i}\right)\left({t}_{i+1}-{t}_{i}\right)=\left(\frac{Lebac\left(0,{t}_{i+1}\right)}{Lebac\left(0,{t}_{i}\right)}-1\right),$

where $Lebac\left(0,{t}_{i}\right)$ is the price of the Lebac bill with maturity ${t}_{i}$ at the initial date, ${t}_{1}=35d$ , ${t}_{2}=63d$ , and so on. The data used to calibrate the model parameters are weekly data, from January 6 to December 27, 2016 In order to model the forward rate curves, we set January 4, 2017 as the initial date and calibrated the initial curve to a Nelson and Siegel rate curve. The observed forward rates $f\left(0,{t}_{i}\right)$ for each of the seven maturities are plotted in “Figure 1”, together with a Nelson and Siegel curve with ${z}_{1}=0.19508$ , ${z}_{2}=0.040724$ and ${z}_{3}=0.0061793$ .

The scale parameter λ was estimated from the data: λ = 0.051237 which corresponds to a maximum at thirty five days, or λ = 0.3587 respectively at five weeks.

In “Figure 2” we show the Hull-White forward curves consistent with Nelson and Siegel manifold for specifics days and the corresponding Lebac forward curve observed in the market in the same day, so as to do a back-testing of the models.

The Hull-White model presents a good fit when it is used to make projections in this market for few weeks. The estimated forward curve at the closest date to the initial one presents a good picture comparing it with what actually happened. It follows the same decreasing shape and with less than one basic point of difference, and it continues providing an accurate projection along the time. As t grows, the estimated curve becomes a flat curve, due to the exponential decreasing terms in t that appears in the coefficients of the forward curve formula (31).

Another interesting issue is to test the projection of the Lebac bills in each maturity date. “Figure 3” shows the cross-section of the forward curve for each maturity and the estimated one. Showing the projected evolution for Lebacs at different maturities and the corresponding market data. As in the further case, the decreasing exponential term in t produces a very small estimated volatility and the estimation is acceptable only for some few weeks.

5. Conclusions

In the previous sections we presented two families of parametric curves that are consistent with the Ho-Lee and the Hull and White short rate models, respectively. These families contain curves that are extensions of the Nelson and Siegel

Figure 1. Initial forward rate curve.

Figure 2. Forward rate curve evolution.

classical curves, in the sense that each element can be written as the sum of a Nelson and Siegel curve plus a linear function in τ or plus an exponential function in τ. We proved that each of these short rate models is consistent with a such family of forward curve manifolds ${\mathcal{G}}^{\lambda }$ . Moreover, for each $\lambda$ , ${G}^{\lambda }$ is a set of forward curves driver by four factors. Three of them are the Nelson and Siegel factors 1, ${\text{e}}^{-\lambda \tau }$ and $\tau {\text{e}}^{-\lambda \tau }$ . In the HL model, the fourth factor is a linear function in $\tau$ multiplied by $\sigma$ , the short rate volatility. In the HW model it is an exponential decaying function ${c}_{1}{\text{e}}^{-a\tau }+{c}_{2}{\text{e}}^{-2a\tau }$ , where a is the model parameter. In both cases, the scale parameter $\lambda$ remains the same along the forward rate process.

Nevertheless, it is worth noting that the dynamics of the Nelson and Siegel curve under the hypothesis of an HJM framework does not seem to be a good estimator for the evolution of the observed forward rate curve. In fact, the decreasing exponential terms in the case of the HW-model produce flat curves as t grows, and the linear term $t\tau$ in the HL-model produces too big values that does not represent rate values. A clearly and known reason for this poor estimation is that they are one factor short rate models and so they cannot capture the

Figure 3. Cross-section of the forward curve.

hole forward curve evolution.

A further research work is to study the consistency of multi-factor short rate models with families of Nelson and Siegel curves, so as to incorporate other volatilities term structures and make the forward rate curve evolution more realistic.

Conflicts of Interest

The authors declare no conflicts of interest.

 [1] Nelson, C. and Siegel, A. (1998) Parsimonious Modeling of Yield Curves. The Journal of Business, 60, 473-489. http://www.jstor.org/stable/2352957 [2] Annaert J., Claes, A.G.P., de Ceuster, M.J.K. and Zhang, H. (2013) Estimating the Yield Curve Using the Nelson-Siegel Model: A Ridge Regression Approach. International Review of Economics and Finance, 27, 482-496. https://doi.org/10.2139/ssrn.2054689 [3] Björk, T. and Christensen, B.J. (1999) Interest Rate Dynamics and Consistent Forward Rate Curves. Mathematical Finance, 9, 323-348. https://doi.org/10.1111/1467-9965.00072 [4] Øksendal, B. (1998) Stochastic Differential Equations. Springer Verlag, Berlin. https://doi.org/10.1007/978-3-662-03620-4 [5] Heath, D., Jarrow, R. and Morton, A. (1992) Bond Pricing and the Term Structure of Interest Rates. Econometrica, 60, 77-106. https://econpapers.repec.org/article/ecmemetrp/v_3a60_3ay_3a1992_3ai_3a1_3ap_3a77-105.htm https://doi.org/10.2307/2951677 [6] Brace, A. and Musiela, M. (1994) A Multifactor Gauss Markov Implementation of Heath, Jarrow, and Morton. Mathematical Finance, 4, 259-283. https://doi.org/10.1111/j.1467-9965.1994.tb00095.x [7] Ho, T.S.Y. and Lee, S.B. (1986) Term Structure Movements and Pricing Interest Rate Contingent Claims. The Journal of Finance, 41, 1011-1029. https://doi.org/10.1111/j.1540-6261.1986.tb02528.x [8] Brigo, D. and Mercurio, F. (2006) Interest Rate Models—Theory and Practice. Springer Verlag, Berlin. [9] Hull, J. and White, A. (1990) Pricing Interest-Rate-Derivative Securities. The Review of Financial Studies, 3, 573-592. http://www.jstor.org/stable/2962116 https://doi.org/10.1093/rfs/3.4.573