Next Article in Journal
Assessing the Current Situation of the World Wheat Market Leadership: Using the Semi-Parametric Approach
Next Article in Special Issue
Volatility Forecasting for High-Frequency Financial Data Based on Web Search Index and Deep Learning Model
Previous Article in Journal
Between Scylla and Charybdis: The Bermudan Swaptions Pricing Odyssey
Previous Article in Special Issue
GO-GJRSK Model with Application to Higher Order Risk-Based Portfolio
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applying Heath-Jarrow-Morton Model to Forecasting the US Treasury Daily Yield Curve Rates

by
Valerii Maltsev
1,† and
Michael Pokojovy
2,*,†
1
Department of Mathematics and Statistics, University of Konstanz, 78464 Konstanz, Germany
2
Department of Mathematical Sciences, The University of Texas at El Paso, El Paso, TX 79968, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2021, 9(2), 114; https://doi.org/10.3390/math9020114
Submission received: 28 November 2020 / Revised: 21 December 2020 / Accepted: 29 December 2020 / Published: 6 January 2021

Abstract

:
The Heath-Jarrow-Morton (HJM) model is a powerful instrument for describing the stochastic evolution of interest rate curves under no-arbitrage assumption. An important feature of the HJM approach is the fact that the drifts can be expressed as functions of respective volatilities and the underlying correlation structure. Aimed at researchers and practitioners, the purpose of this article is to present a self-contained, but concise review of the abstract HJM framework founded upon the theory of interest and stochastic partial differential equations in infinite dimensions. To illustrate the predictive power of this theory, we apply it to modeling and forecasting the US Treasury daily yield curve rates. We fit a non-parametric model to real data available from the US Department of the Treasury and illustrate its statistical performance in forecasting future yield curve rates.
MSC:
35R60; 60H15; 60H20; 60H35; 62M20; 65C30

1. Introduction

Stocks and bonds play a crucial role in making investment decisions. However, there is a significant difference between stock and bond markets. A stock market is a platform enabling companies to issue their shares and trade them on a stock exchange. Buying these shares and having ownership stake, investors can benefit from the performance of the issuer in the future. In contrast, a bond market is used for issuing and trading debt securities. Usually, an investor and a bond issuer make a contract enabling the former to benefit from the interest paid on a regular basis and full nominal amount received when reaching maturity. The bond issuer may be a government, a corporation, etc.
Besides, there are bonds that do not pay an intermediate interest but are traded at a deep discount, rendering profit at the maturity when they are redeemed for their full face value. These bonds are called zero-coupon bonds. For a maturity T 0 , let P ( t , T ) denote the price of a zero-coupon bond at time t [ 0 , T ] . It is known from the theory of interest rates (viz. Definition 10 in Section 2) that the instantaneous forward rate of a zero-coupon bond is given by
f ( t , T ) = ln P ( t , T ) T   for   all   t [ 0 , T ] ,   T > 0 .
Under appropriate assumptions, Heath, Jarrow and Morton [1] demonstrated that f ( t , T ) can be modeled as a Markovian diffusion that evolves according to the following stochastic differential equation (SDE):
d f ( t , T ) = α ( t , T ) d t + σ ( t , T ) · d W ( t ) , f ( 0 , T ) = f 0 ( T ) ,
where W ( t ) t [ 0 , T ] is a d-dimensional standard Brownian motion, T f 0 ( T ) is the market instantaneous forward curve at time t = 0 , α ( t , T ) t [ 0 , T ] and σ ( t , T ) t [ 0 , T ] are 1- and d-dimensional stochastic processes adapted to the natural filtration of W ( t ) t [ 0 , T ] and σ ( t , T ) · d W ( t ) = k = 1 d σ k ( t , T ) d W k ( t ) is formally defined as the scalar product. Under suitable regularity assumptions, integrating Equation (2) over [ 0 , t ] for t [ 0 , T ] , the instantaneous forward rate f ( t , T ) can be represented as an Itô-process:
f ( t , T ) = f ( 0 , T ) + 0 t α ( s , T ) d s + 0 t σ ( s , T ) · d W ( s ) .
However, the dynamics in Equation (2) is not necessarily arbitrage-free. In the spirit of [2], we understand the presence of arbitrage as the existence of a self-financing strategy which yields a positive income in the future without investing a positive amount today (cf. Section 4). It was further proved in [3] that the existence of an equivalent martingale measure implies the absence of arbitrage opportunities. It was also shown that the unique equivalent martingale measure exists if the drift α depends on the volatility σ (see ([2], Chapter 2) for details). In particular, the requirement that the evolution in Equation (2) is taking place under a risk-neutral measure Q necessitates the following representation:
α ( t , T ) = σ ( t , T ) · t T σ ( t , s ) d s = k = 1 d σ k ( t , T ) t T σ k ( t , s ) d s   for   all   t [ 0 , T ] .
Therefore, for all t [ 0 , T ] , we conclude that f ( t , T ) is fully specified (as a stochastic process) as long as the volatility vector σ is provided:
f ( t , T ) = f ( 0 , T ) + k = 1 d 0 t σ k ( τ , T ) τ T σ k ( τ , s ) d s d τ + k = 1 d 0 t σ k ( s , T ) d W k ( s ) .
Recalling Equation (1), Itô’s lemma furnishes the following SDE for the zero-coupon bond price P ( t , T ) for all t [ 0 , T ] :
d P ( t , T ) = P ( t , T ) r t   d t t T σ ( t , s ) d s d W ( t ) ,
where r t is the instantaneous short-term interest rate at t given by
r t = f ( t , t ) = f ( 0 , t ) + 0 t σ ( τ , t ) τ t σ ( τ , s ) d s d τ + 0 t σ ( s , t ) d W ( s )   = f ( 0 , t ) + k = 1 d 0 t σ k ( τ , t ) τ t σ k ( τ , s ) d s d τ + k = 1 d 0 t σ k ( s , t ) d W k ( s ) .
In his dissertation [4], Filipović proposed to extend this approach to an infinite-dimensional framework and consider the following stochastic partial differential equation (SPDE):
d f t ( x ) = ( d d x f t ( x ) + α t ( x ) ) d t + σ t ( x ) · d W ( t ) , f t ( x ) = f 0 ( x ) ,
with the so-called Musiela’s notation (cf. [5]) f t ( x ) : = f ( t , t + x ) , x [ 0 , T t ] being employed (see also [6], Chapter 6). Here, W ( t ) t [ 0 , T ] is an infinite-dimensional vector of (coutably many) independent real-valued Brownian motions, the forward rate curve x f t ( x ) is a random variable taking its values in a function space H and the drift α is given by the no-arbitrage condition
α t ( x ) : = α ( t , t + x ) = k = 1 σ k ( t , t + x ) 0 x σ k ( t , t + s ) d s .
To facilitate the analytic treatment of Equation (4), we put it into the framework of abstract operator SDEs with stochastic coefficients:
d X t = ( A X t + F ( t , X t ) ) d t + B ( t , X t ) d W ( t ) , X 0 = X 0 ,
where X is an H -valued random process, X 0 is a random H -valued F 0 -measurable initial value, A : H H is a linear (possibly unbounded) operator, W ( t ) t [ 0 , T ] is a standard cylindrical Wiener process over U , F : [ 0 , T ] × Ω × H H , B : [ 0 , T ] × Ω × H L HS ( U , H ) and H , U are separable Hilbert spaces. Here and in the sequel, L HS ( · , · ) denotes the space of Hilbert-Schmidt operators (cf. Definition A8). Letting
A : = d d x ,   X t : = f t ,   F : = α t ,   B : = σ t , · .
we can easily see that Equation (4) is a special case of Equation (5). A straightforward integration of Equation (5) over [ 0 , t ] yields the integral formulation of Equation (5):
X t = X 0 + 0 t A X s + F ( s , X s ) d s + 0 t B ( s , X s ) d W ( s )   for   all   t [ 0 , T ] .
The rest of the paper is structured as follows. Section 2, Section 3 and Section 4 are meant to provide a concise review of the HJM framework and related results. With most of the results being scattered over a number of articles and monographs, our thrust was to give a coherent presentation on this matter with consistent terminology and notation to make the theory better accessible to a broad scientific and professional audience. In Section 2, a brief summary of the interest rate theory is given, in particular, a relation between the price of a zero-coupon bond and its instantaneous forward rate from Equation (1) is established. In Section 3, we introduce the abstract HJM model first presented and studied by [4,6], discuss major advantages of this approach and protocol some known well-posedness results. In Section 4, we summarize central properties of the abstract HJM model, in particular, we give sufficient conditions under which the abstract HJM model defines an arbitrage-free market, describe the long rates associated with the model, and give an example of a state space compatible with the instantaneous forward rate dynamics. Section 5 contains the main contribution of this paper. A (semi-)implicit numerical scheme is developed and applied to Equation (4). Using real data from the US Department of the Treasury [7], the resulting non-parametric model is calibrated using a custom-developed least-squares-like procedure involving principal component analysis (PCA) for parameter estimation and dimension reduction. The simplicity of our approach is one of its main advantages. Indeed, it does not rely on any sophisticated optimization techniques or third-party libraries, thus, making itself attractive and accessible to a wide audience. The model is then deployed to forecast future yield curve rates, estimate the mean future yield curve and compute prediction regions. Using the predictive power of the model as a benchmark measure in lieu of the mere calibration quality is another novelty of our contribution. To the best of our knowledge, this is the first time an infinite-dimensional model is applied to analyze the US Treasury’s daily yield curve rates. We refer [8,9,10] for the finite-dimensional treatise. In Section 6, conclusions are presented. Appendix A, Appendix B and Appendix C contain all seminal probabilistic and functional-analytic concepts and results [6,11,12,13,14,15] used in the article.

2. Interest Rate Theory

Before presenting the abstract Heath-Jarrow-Morton model in Section 3, we give a concise summary of the theory of interest rates based on ([2] pp. 2–18). The goal of this section is twofold. First, we want to give motivation for using the theory of interest rates in this paper. Second, we aim to prove relation (1) between the price and the instantaneous forward rate of a zero-coupon bond in order to define the abstract HJM model afterwards.
As pointed out in ([16], Chapter 25), multivariate (or multidimensional) models have numerous advantages over single-state short rate models, in particular, in describing the intrinsic volatility structure or explaining the inversion of the yield curve, etc. The HJM framework, being at “the far end of this [multivariate] spectrum” (cf. [16], p. 388), is, probably, the most comprehensive one as it employs the forward rate curve as an infinite-dimensional state variable. We refer to the aforecited monograph as well as ([17], pp. 242–245) for further discussions.
Application areas of the HJM framework are very diverse. They include, but are not limited to, modeling of interest futures [18], risk management [19], American option pricing [20], credit risk modeling [21], investment management [22], etc.

2.1. Bank Account and Short Rate

Depositing a certain amount of money in a bank account is typically expected to guarantee its growth at some rate. We need appropriate definitions to express this concept in mathematical terms. Consider a bank account which represents a riskless investment with the profit accruing continuosly.
Definition 1
(Bank account). Let B ( t ) be the value of a bank account at time t 0 . We assume that the bank account evolves according to the differential equation
d B ( t ) = r t B ( t ) d t , B ( 0 ) = B 0 ,
where r t is a positive function of time. Letting B 0 = 1 , we have
B ( t ) = exp 0 t r s d s .
Definition 1 suggests that a unit amount invested at t = 0 has value (8) at t > 0 , implying also that B ( t ) accrues at the short rate r t .
In fact, the first-order Taylor expansion of B ( t + Δ t ) is
B ( t + Δ t ) = B ( t ) ( 1 + r t Δ t ) ,
whence
B ( t + Δ t ) B ( t ) B ( t ) = r t Δ t   as   Δ t 0 .
Definition 2
(Numéraire). A numeraire is any positive non-dividend-paying asset.
Remark 1.
The concept of bank account numéraire is important for relating amounts of currencies available at different times. Assume for simplicity that r t is deterministic. If K units of currency are deposited into a bank account at the time t = 0 , the amount of units of currency at time t 0 is K B ( t ) . Analogously, the amount is K B ( T ) at time T t . If we want to have one unit of currency at time T, i.e., K B ( T ) = 1 , we have to initially invest K = 1 / B ( T ) . Hence, the value of the invested amount K at time t is
K B ( t ) = B ( t ) B ( T ) .
Remark 1 motivates the following definition (cf. [2], Chapter 1).
Definition 3
(Stochastic discount factor). The amount at time t “equivalent” to one unit of currency payable at time T
D ( t , T ) = B ( t ) B ( T ) = exp t T r s d s .
is referred to as stochastic discount factor.

2.2. Zero-Coupon Bonds and Spot Interest Rates

Definition 4
(Zero-coupon bond).A zero-coupon bond is a contract that entitles its holder to a (single) payment of one unit of currency at maturity T. Denoting the contract value at time t T by P ( t , T ) , we observe P ( T , T ) = 1 for all T 0 .
Definition 5
(Time to maturity). The amount of time T t from the present moment of time t to the maturity time T > t is referred to as time to maturity T t .
Obviously, T t is well-defined if t and T are real numbers associated with some time instants. However, if t and T denote two dates D 1 , D 2 , e.g., of format (month/day/year), we can define the amount of time between them in terms of the number of days. Moreover, the number of days between D 1 and D 2 is calculated according to appropriate market convention, which specifies how these days are to be counted.
Definition 6
(Year fraction). Let τ ( t , T ) denote a time measure between t and T, which is usually referred to as year fraction between the dates t and T. If t and T are less than one day apart (typically, occuring when dealing with limiting quantities which involve time to maturities tending to zero), τ ( t , T ) is roughly given by T t . A particular choice is known as the day-count convention (see [2,23] for examples).
Additionally, we have to exploit another feature of the zero-coupon bonds called the compounding type in order to switch from the price of a zero-coupon bond and its instantaneous forward rate and vice versa.
Definition 7
(Continuously-compounded spot interest rate). The (constant) rate
R ( t , T ) : = ln P ( t , T ) τ ( t , T )
at which an investment of amount P ( t , T ) at time t accrues continuously to produce a unit amount of currency at maturity T is referred to as continuously-compounded spot interest rate.
Definition 8
(Simply-compounded spot interest rate). Given an initial capital of P ( t , T ) at time t and assumning accruing takes place proportionally to the investment time, the (constant) rate
L ( t , T ) : = 1 P ( t , T ) τ ( t , T ) P ( t , T )
at which money needs to be invested to yield a single unit of currency at maturity T is called the simply-compounded spot interest.

2.3. Forward Rates

Forward rates are interest rates that can be locked in on a particular day for an investment in a future time period. They are characterized by three time instants: the time t at which the rate is considered, its expiry T and its maturity S, t < T < S .
The forward rate can be also viewed as a typical forward rate agreement represented by a contract involving three time characteristics: the current time t, the expiry time T > t and the maturity S > T . It enables the investor to benefit from the interest paid between T and S. Moreover, it is possible to fix the interest between T and S at rate K with simply-compounded contract rates, implying a payment of τ ( T , S ) L ( T , S ) N and a receipt of τ ( T , S ) K N , with N denoting the contract nominal value. This renders the contract value to be N τ ( T , S ) K L ( T , S ) . Recalling Definition 8, the former quantity rewrites it as
N τ ( T , S ) K 1 P ( T , S ) + 1 .
Let 1 / P ( T , S ) be the amount of currency held at time S. Its value at time T is obtained by multiplying this amount with the zero coupon price P ( T , S ) , i.e.,
P ( T , S ) 1 / P ( T , S ) = 1 ,
so this term is equivalent to holding one unit of currency at the time T. Besides, one unit of currency at the time T is worth P ( t , T ) units of currency at the time t. Therefore, the amount 1 / P ( T , S ) in S is equivalent to an amount of P ( t , T ) in t. Now, consider the remaining terms in Equation (9). The worth of W : = τ ( T , S ) K + 1 at the maturity S at the time t is
P ( t , S ) W = P ( t , S ) τ ( T , S ) K + P ( t , S ) .
Therefore, the total value of the contract at the time t is
N P ( t , S ) τ ( T , S ) K P ( t , T ) + P ( t , S ) .
Remark 2.
The value of K rendering the contract in Equation (10) fair at the time t is unique.
Definition 9
(Simply-compounded forward interest rate). Given expiry T > t and maturity S > T , the simply-compounded forward interest rate at time t reads as
F ( t ; T , S ) : = 1 τ ( T , S ) P ( t , T ) P ( t , S ) 1 .
As the maturity of the forward rate approches its expiry, the concept of instantaneous forward rate emerges. If S is sufficiently close to T, then τ ( T , S ) = S T and we have:
lim S T F ( t ; T , S ) = lim S T 1 P ( t , S ) P ( t , S ) P ( t , T ) S T   = 1 P ( t , T ) P ( t , T ) T   = ln P ( t , T ) T .
This motivates the following definition.
Definition 10
(Instantaneous forward interest rate). The instantaneous forward interest rate prevailing at time t for the maturity T > t is denoted by f ( t , T ) and given as
f ( t , T ) = lim S T f ( t ; T , S ) = ln P ( t , T ) T .
Remark 3.
Using Musiela’s notation (viz. [5]), f t ( x ) : = f ( t , t + x ) for x = T t and t , x 0 . The space of functions f t will be denoted by S and called the state space.

3. The Abstract Heath-Jarrow-Morton Model

Filipović [4] proposed to model the instantaneous forward rate of a zero-coupon bond via an infinite-dimensional SDE driven by the standard cylindrical Wiener process. In contrast to the standard finite-dimensional Heath-Jarrow-Morton (HJM) framework, the underlying bond market captures the entire (unobservable) infinite-dimensional interest rate term structure. Among many other advantages, the infinite-dimensional HJM model mitigates the so-called “maturity specific risk” typically present in finite-dimensional models. For a more detailed discussion of infinite-dimensional models, the reader may be referred to the monograph ([6], Chapter 6.5). Also, it gives a natural way to “extrapolate” outside of the original finite maturity grid. From the empirical point of view, the infinite-dimensional structure of the model not only helps better mimic the observed dynamics by merely increasing the “number” of model parameters, but makes it possible to capture the non-trivial topological aspects of the dynamics, thus, improving the predictive power and allowing for natural dimension reduction or regularization without resorting to ad hoc penalties. Arguably, this justifies the increased price associated with the infinite-dimensional functional-analytic and probabilistic framework.
In this section, we introduce the abstract HJM model following the approach of [4,6]. The resulting model is a semilinear stochastic PDE (SPDE) driven by the standard cylindrical Wiener process (see Appendix C). The well-posedness of these abstract SPDEs in the class of mild solutions is obtained using the abstract results from ([11], Chapter 7) and discussed in Section 3.1. The results are then applied to justify the well-posedness of the abstract HJM model to be introduced in Section 3.2.

3.1. Existence and Uniqueness of Mild Solutions to Semilinear Spdes

We consider the abstract SPDE
d X t = A X t + F ( t , X t ) d t + B ( t , X t ) d W ( t ) , X 0 = X 0
and investigate the existence and uniqueness of mild solutions. Depending on the context, the solution process will be denoted as ( X t ) t [ 0 , T ] or ( X ( t ) ) t [ 0 , T ] . Let Ω , F , { F t } t 0 , P be a filtered probability space and let T > 0 . Let P T and P denote the predictable σ -algebras on [ 0 , T ] × Ω and [ 0 , ) × Ω , respectively. We also assume that U , H are separable Hilbert spaces and W ( t ) t [ 0 , T ] is a standard cylindrical Wiener process over U given in Definition A7. Therefore, using the notation of Appendix C.1, we can let Q = I and U 0 = Q 1 2 U U . The spaces U , H , L HS ( U , H ) induce corresponding Borel σ -algebras B ( U ) , B ( H ) , B L HS ( U , H ) . (Recall that L HS ( · , · ) denotes the space of Hilbert-Schmidt operators.) Lastly, X 0 is a random H -valued F 0 -measurable initial value. See Appendix A and Appendix C for a brief review of probabilistic concepts and results used in this section.
Assumption 1.
Let A, F and B in Equation (11) satisfy:
  • F : [ 0 , T ] × Ω × H H is measurable from [ 0 , T ] × Ω × H , P T × B ( H ) to H , B ( H ) ,
  • B : [ 0 , T ] × Ω × H L HS ( U , H ) is measurable from [ 0 , T ] × Ω × H , P T × B ( H ) to L HS ( U , H ) , B L HS ( U , H ) ,
  • There exists a constant C > 0 such that for all u , v H , t [ 0 , T ] , ω Ω we have:
    F ( t , ω , u ) F ( t , ω , v ) H + B ( t , ω , u ) B ( t , ω , v ) L HS ( U , H ) C u v H , F ( t , ω , u ) H 2 + B ( t , ω , u ) L HS ( U , H ) 2 C 2 1 + u H 2 ,
  • A is the generator of a strongly continuous operator semigroup ( S t ) t 0 on H .
Definition 11
(Mild Solution). An H -valued predictable stochastic process ( X ( t ) ) t [ 0 , T ] with
P 0 T X ( s ) H 2 d s < = 1
is referred to as a mild solution to Equation (11) if it satisfies the integral equation (well-defined on the strength of Proposition A9):
X t = S t X 0 + 0 t S t s F ( s , X s ) d s + 0 t S t s B ( s , X s ) d W ( s )   for   all   t [ 0 , T ]
P -almost surely.
Note that the predictability of ( X ( t ) ) t [ 0 , T ] implies it possesses a continuous modification, i.e., P X C 0 [ 0 , T ] , H = 1 , which (in addition to squared-integrability) is essential for the latter integral to be well-defined.
Theorem 1
([11], Section 7.1). Under Assumption 1, there exists a unique (up to indistinguishability) H -valued mild solution X ( t ) t [ 0 , T ] to Equation (11). Moreover, for any p 2 , there exists a positive constant C, which only depends on p and T, such that
sup t [ 0 , T ] E X ( t ) H p C 1 + E X 0 H p .

3.2. Definition of Abstract HJM Model

Let ( Ω , F , { F t } t 0 , P ) be a filtered probability space. Further, let I either be a compact interval of the form [ 0 , t 0 ] for some positive t 0 or the half-line [ 0 , ) and let ( W ( t ) ) t I be a standard cylindrical Wiener process over a separable Hilbert space U . Let S be a state space and P I be a predictable σ -algebra on I × Ω . We make the following assumption (cf. [6], Chapter 6).
Assumption 2.
Suppose:
  • S is a separable Hilbert space consisting of continuous, real-valued functions f : I R . Additionally, we assume that for every x I the functionals
    δ x ( f ) = f ( x ) ,   J x ( f ) = 0 x f ( s ) d s
    are continuous on S . (The former assumption can also be encountered in statistics, for example, in the context of Reproducing Kernel Hilbert Spaces (RKHS)).
  • The left shift semigroup ( S t ) t 0 defined via
    ( S t f ) ( x ) = f ( t + x )
    is a strongly continuous operator semigroup on S generated by a (linear) operator A.
  • F HJM is a measurable mapping from a non-empty subset of L HS ( U , S ) into S and defined by
    F HJM ( σ ) ( x ) = σ * δ x , σ * J x U
    for all σ L HS ( U , S ) .
The crucial property of the space S is that its elements f S are (locally) integrable functions so that P ( t , T ) = e 0 T t f t ( s ) d s is well-defined. Indeed, the local integrability holds true as the linear functional J x is bounded on S for each x I t since
0 x f ( s ) d s x sup s [ 0 , x ] f ( s ) x sup s [ 0 , x ] δ s S * f S <
where sup s [ 0 , x ] δ s S * is finite on the strength of Banach-Steinhaus theorem.
The financial implication of the first part of Assumption 2 is the identity r t = f t ( 0 ) for the short interest rate. The bank account is then defined via B ( t ) = e 0 t r s d s . The semigroup ( S t ) t 0 given in the second part of Assumption 2 enables us to switch to Musiela’s notation f t ( x ) , where f t ( x ) = f ( t , t + x ) . The third part of Assumption 2 is intimately related to the no-arbitrage principle and is needed to define the drift term in the HJM equation. Further, it is generally false that F HJM ( σ ) is an element of S unless σ is an element of a proper subspace of L HS ( U , S ) .
Assumption 3.
Assume there exists a subspace S of S such that the binary operator defined by
( f g ) ( x ) = f ( x ) 0 x g ( s ) d s ,
mapping from S × S into S is bounded, i.e., there exists C > 0 such that the following bound
f g S C f S g S
holds for all f , g S .
The following estimate for F HJM is known after ([6], p. 167).
Proposition 1.
Let the space S satisfy Assumption 3. Then the mapping F HJM satisfies the estimate
F HJM ( σ 1 ) F HJM ( σ 2 ) S C σ 1 + σ 2 L HS ( U , S ) σ 1 σ 2 L HS ( U , S )
for all Hilbert-Schmidt operators σ 1 , σ 2 L HS ( U , S ) with the range contained in S . In particular, the mapping F HJM is measurable from L HS ( U , S ) to S .
Now we can give a rigorous formulation of the abstract HJM model adopted from ([6], Section 6.3). It is rooted in the standard HJM framework [1], which describes a certain class of arbitrage-free finite-dimensional models for the instantaneous forward rate. We provide a more powerful framework, referred to as the abstract HJM model, which incorporates arbitrage-free infinite-dimensional models.
Definition 12
(The abstract HJM model). The abstract HJM model in a functional space S is a pair of functions ( λ , σ ) such that
  • λ is a measurable function from ( I × Ω × S , P I × B ( S ) ) into ( U , B ( U ) ) ,
  • σ is a measurable function from ( I × Ω × S , P I × B ( S ) ) into ( L HS ( U , S ) , B ( L HS ( U , S ) )
so that there exists a non-empty set of S -valued F 0 -measurable initial conditions Φ 0 , which admits a unique time-continuous mild S -valued solution ( Φ t ) t I to the HJM equation:
d Φ t = ( A Φ t + α ( t , Φ t ) ) d t + σ ( t , Φ t ) d W t , α ( t , Φ t ) = F HJM σ ( t , Φ t ) + σ ( t , Φ t ) λ ( t , Φ t ) , Φ 0 = Φ 0 .
The mild formulation of Equation (12) is given via:
Φ t = S t Φ 0 + 0 t S t s α ( s , Φ s ) d s + 0 t S t s σ ( s , Φ s ) d W s , α ( t , Φ t ) = F HJM σ ( t , Φ t ) + σ ( t , Φ t ) λ ( t , Φ t ) , Φ 0 = Φ 0 .
Imposing typical Lipschitzianity and boundedness assumptions, as a direct corollary of Theorem 1, we obtain (cf. [6], p. 169):
Proposition 2
(Sufficient existence condition). Let S satisfy Assumption 3 with a closed subspace S S such that f g S C f S g S for f , g S . Let the operator σ map into S for every ( t , ω , f ) I × Ω × S . Assume σ is bounded and the Lipschitzian
σ ( t , f ) σ ( t , g ) L HS ( U , S ) C f g S , σ ( t , f ) λ ( t , f ) σ ( t , g ) λ ( t , g ) S C f g S
for some constant C > 0 and all t I , f , g S , then the pair ( λ , σ ) is an abstract HJM model on S . Furthermore, for any initial S -valued forward curve Φ 0 , there exists a unique time-continuous mild solution ( Φ t ) t I to Equation (12) such that
E sup t I Φ t S p <   for   all   finite   p 0 .

4. Absence of Arbitrage

A fundamental economic assumption behind any adequate financial model is the absence of arbitrage opportunities in the associated financial market (cf. [24], pp. 637–654). The absence of arbitrage means impossibility to invest nothing today and receive a positive amount in the future. We restrict ourselves to continuous-time economies (cf. [25], pp. 381–408).
Consider a filtered probability space ( Ω , F , { F t } t [ 0 , T ] , Q 0 ) with a right-continuous filtration { F t } t [ 0 , T ] . Assume ( K + 1 ) non-dividend paying securities are traded from time 0 until time T continuously with their modeled by an adapted ( K + 1 ) -dimensional semimartingale M = { M t   |   0 t T } having non-negative components M 0 , M 1 , , M K . The asset indexed by 0 is a bank account. We refer to Appendix B for a brief review of the martingale theory.
Definition 13
(Trading strategy). A trading strategy is a ( ( K + 1 ) -dimensional) process ϕ = { ϕ t   |   0 t T } with locally bounded, predictable components ϕ 0 , ϕ 1 , , ϕ K . The value process associated with a strategy ϕ is defined by
V t ( ϕ ) = ϕ t · M t = k = 0 K ϕ t k M t k   for   all   t [ 0 , T ] ,
and the gains process associated with a strategy ϕ are given by
G t ( ϕ ) = 0 t ϕ s · d M s = k = 0 K 0 t ϕ s k d M s k   for   all   t [ 0 , T ] .
The quantity ϕ t k corresponds to the number of units of k-th security held by an investor at time t. Predictability of ϕ k suggests that ϕ t k is known directly before t.
Definition 14
(Self-financing trading strategy). We refer to a trading strategy ϕ as self-financing if V T ( ϕ ) 0 and V t ( ϕ ) = V 0 ( ϕ ) + G t ( ϕ ) for t [ 0 , T ) .
Definition 15
(Equivalent martingale measure). A probability measure Q on ( Ω , F ) is referred to as equivaletn martinage measure if it satisfies:
  • Q and Q 0 are equivalent martingale measures, i.e., Q ( A ) = 0 if and only if Q 0 ( A ) = 0 for all A F ,
  • The Radon-Nikodym derivative d Q d Q 0 belongs to L 2 ( Ω , F , Q 0 ) ,
  • The “discounted asset price” process D ( 0 , · ) M is an ( { F t } t [ 0 , T ] , Q ) -martingale, i.e.,
    E Q D ( 0 , t ) M t k | F s = D ( 0 , s ) M s k for all k { 1 , 2 , , K } and all 0 s t T .
Definition 16
(Locally equivalent martingale measure). We say that Q is a locally equivalent martingale measure if for all t > 0 the restriction of Q to F t is equivalent to the restriction of P to F t .
An arbitrage opportunity is defined by a self-financing strategy ϕ such that V 0 ( ϕ ) = 0 but Q 0 ( V T ( ϕ ) > 0 ) > 0 . It is known (cf. ([3], pp. 215–260) and ([26], pp. 17–19)) that the existence of an equivalent martingale measure implies the absence of arbitrage opportunities:
Theorem 2
(Fundamental theorem of asset pricing ([6], Theorem 2.1)). There are no arbitrage strategies if there exists a probability measure Q , equivalent to P , such that the discounted asset prices process P D ( t ) t 0 is a local martingale under Q .
Chapter 3.1.1 of [6] provides sufficient conditions for the abstract HJM model from Section 3 to be free of arbitrage. These results are summarized below. Consider an abstract HJM model ( σ , λ ) with an initial condition Φ 0 and denote by ( Φ t ) t [ 0 , T ] the unique mild solution to the HJM equation. For brevity, denote λ t : = λ ( t , ω , f t ) and σ t : = σ ( t , ω , f t ) .
Theorem 3
(Conditions of no arbitrage for the abstract HJM model). Provided
E exp 1 2 0 t λ s U 2 d s + 0 t λ s d W ( s ) = 1     and   0 t E 0 t σ s * δ s τ U 2 d τ 1 2 d s <
for all t [ 0 , T ] , then the market given by the abstract HJM model ( σ , λ ) admits no arbitrage.
In its current form, the setup of Theorem 3 appears to be too general for practical use. Indeed, it would suffice to show the Novikov condition for Q
E Q exp 1 2 0 T σ s * J T s U 2 d s <
is satisfied in order to prove that ( P D ( t ) ) t 0 is a martingale (cf. [6], Theorem 6.1) since every martingale is a local martingale.

4.1. Long Rates

Modeling the forward rate as a function on a bounded interval I is different from modeling it on the real half-line. Namely, we encounter an unexpected phenomenon that the long rate defined by the limit
l t = lim x f t ( x ) = : f t ( )
never decreases. This phenomenon for long forward and zero-coupon rates was previously studied in [27]. The intrinsic nature of this phenomenon is rooted in the fact that the (appropriately defined) discounted bond prices P ˜ ( t , T ) t [ 0 , T ] are local martingales for all T > 0 (cf. [6], p. 172). The following result is known after ([6], Chapter 6.3.2).
Assumption 4.
Every f S is a function f : R + : = [ 0 , ) R such that the limit lim x f ( x ) exists. The functional δ : f lim x f ( x ) is an element of S * .
Theorem 4.
If 0 s t , then l s l t P -almost surely.

4.2. Example of State Space

An example of a state space S was first proposed in ([4], Chapter 3) and further investigated in ([6], Chapter 3.3). Obviously, this space should fulfill the first and the second parts of Assumption 2 under which the abstract HJM model is well-defined, i.e., it should consist of real-valued continuous functions, for which the functionals δ x , J x should be continuous on S and the strongly continuous semigroup ( S t ) t 0 should exist on S . Parenthetically, we point out that, though the following presentation assumes I = [ 0 , ) , our empirical example in Section 5 will rather utilize a bounded interval [ 0 , x max ] .
Definition 17
(Example of state space). Let w : R + R + be a monotonically increasing function such that
0 1 w ( s ) d s < .
denote by H w the space
H w = { f L loc 1 ( R + ) | f L loc 1 ( R + )   and   f H ω < } ,
where f is the distributional derivative of f and
f H w 2 : = f , f H w   with   f , g H w = f ( 0 ) g ( 0 ) + 0 f ( s ) g ( s ) w ( s ) d s .
If f L loc 1 ( R + ) has a weak derivative f L loc 1 ( R + ) , e.g., some function ϕ such that
0 f ( s ) ϕ ( s ) d s = 0 f ( s ) ϕ ( s ) d s   for   all   ϕ C c 1 ( R + ) ,
then there exists an absolutely continuous representative f, still denoted by f, such that
f ( t 2 ) f ( t 1 ) = t 1 t 2 f ( s ) d s   for   all   t 1 , t 2 R +
(cf. [28], Section VIII.2). In particular, this implies that the trace δ x : f f ( x ) is well-defined.
Proposition 3
([6], Chapter 3.3). If w satisfies inequality (14), then the inner product space H ω is a separable Hilbert space. The evaluation functional δ x and the definite integration functional J x defined by
δ x ( f ) = f ( x ) ,   J x ( f ) = 0 x f ( s ) d s
respectively, are continuous on H w for all positive x. Further, the operator semigroup ( S t f ) ( x ) = f ( t + x ) on H w is strongly continuous.

5. Example: Modeling and Predicting US Treasuries

Over the past two decades, the problem of parameter estimation, also referred to as model calibration or estimation, for the classical finite-dimensional HJM model [1] (viz. Equation (2)) has attracted major attention from various researchers in econometrics, quantitative finance, stochastic analysis, numerical analysis, etc. In the finite-dimensional context, estimation for the HJM model is typically performed using the method of least squares, the method of moments, maximum likelihood estimation and penalized versions thereof, stochastic filtration, kriging, parametric and nonparametric Bayesian inference, etc. See, e.g., [29,30,31,32,33,34] for relevant studies in both parametric and non-parametric settings. We also refer to “classic” papers [8,9,10] as well as recent contributions [35,36] and references therein for closely related results for the Vasicek model [37].
As for the abtsract HJM model (viz. Equation (4)) in the infinite-dimensional context, with the exception of the numerical paper [38], no results are known in the literature. With the primary focus on studying the convergence properties of the proposed numerical splitting schemes, ([38], Section 2.3.1) presents a brief empirical study based on a dataset from the dissertation [39]. As opposed to the finite-dimensional situation, infinite-dimensional parameter estimation or calibration consitutes a (possibly ill-posed) statistical inverse problem. To the best of our knowledge, no rigorous investigations of HJM model calibration in the context of statistical inverse problems are currently available in the literature. A promissing angle of attack is though provided in [40]. Also, no calibration and/or forecasting studies on the infinite-dimensional HJM model in connection with real-world data are also available up to date. In this section, we investigate a non-parametric calibration problem for the linear version of Equation (4). To this end, we discretize the equation with an implicit/explicit “splitting scheme” and propose a PCA- and least-squares based parameter estimation/dimension reduction procedure to fit the discretized model. We illustrate our approach by applying our proposed scheme to a real-world dataset presented below.
Daily US Treasury yield curve rates are available from the US Department of the Treasury [7]. For any given business day (i.e., all days excluding holidays and weekends), the yield curves relate the closing market bid yields on actively traded Treasury securities in the over-the-counter market to their time to maturity. Further details on the methodology used to calculate the yield rates are provided in [7]. The time-to-maturity grid includes 1 mo, 2 mo (starting 16 October 2018), 3 mo(nths), 6 mo, 1 y(ea)r, 2 yr, 3 yr, 5 yr, 7 yr, 10 yr, 20 yr and 30 yr. Figure 1 displays historical yield rates for the year 2018. Since the 2-month rates were missing for the dates between 1 January 2018 and 15 October 2018, linear interpolation was used to estimate these values.
Employing Musiela’s notation [5], let y t ( x ) denote the yield rate of a US treasury traded at time t with the time to maturity x. According to ([4], Definition 3.13), y t ( x ) can be expressed with the aid of the instantaneous forward rate f t ( x ) as
y t ( x ) = 1 x 0 x f t ( x ) d x 1 x I x ( f t ) .
In particular, invoking the mean value theorem, the former equation suggests
y t ( 0 ) = f t ( 0 ) r t ,
where r t stands for the short rate. Multiplying both sides of Equation (15) by x and differentiating in x, we can solve for f t ( x ) , i.e.,
f t ( x ) = x x y t ( x ) = y t ( x ) + x x y t ( x ) .
Unfortunately, no “natural” boundary condition is readily available for f t ( x ) . In contrast, since y t ( x ) can be assumed bounded in x = 0 , a homogeneneous Dirichlet boundary condition emerges for x y t ( x ) . Therefore, in spirit of ([4], Chapter 3.3), we will rather work with the (observable) quantity
Y ( t , x ) = x y t ( x ) ,
referred to as the “integrated forward rate” since
Y ( t , x ) = I x ( f t ) 0 x f t ( η ) d η .
This quantity, Y ( t , x ) , satisfies the homogeneous Dirichlet boundary condition Y ( t , 0 ) = 0 at x = 0 .
Since the shortest time-to-maturity available in our dataset is x min = 1 month, y t ( 0 ) will be approximated/extrapolated as y t ( x min ) . Further, as Figure 1 suggests, the x-grid becomes very coarse for x > 3 years ( = 36 month ) . In particular, this means that the infinitesimal effects of the operator A in the abstract HJM model in Equation (13) will become negligible as the semigroup S t ( · ) will practically “collapse” to the identity operator, essentially reducing the model (for large x) to the classic HJM model. Therefore, the maturities were truncated at 36 months. Lastly, the observed yield rates are interpolated onto [ 0 , 36 ] month using cubic splines. Figure 2 displays the historical Y-curves for the year 2018 over the truncated 36-month x-horizon.

5.1. The Model

Consider the space H w defined in Section 4.2. The “spatio-”temporal evolution of the forward rates will be assumed to follow the linear autonomous HJM model from Equation (4), i.e.,
d f t ( x ) = ( d d x f t ( x ) + α ˜ t ( x ) ) d t + k = 1 σ k ( x ) · d W ˜ k ( t ) , f 0 ( x ) = f 0 ( x ) ,
with some initial datum f 0 and α ˜ t ( x ) determined via the “no arbitrage” condition
α ˜ ( x ) = k = 1 σ k ( x ) 0 x σ k ( η ) d η ,
where the evolution in Equation (19) is taking place under the risk-neutral measure P ˜ . To assure the model is identifiable, similar to [8,10], σ k ( x ) ’s will be assumed orthogonal with respect to an appropriate inner product to be discussed later.
Following the results of Section 4, we can write
W ˜ t k = W t k + 0 t λ s k d s
for a standard cylindircal Wiener process ( W t ) t 0 with respect to the physical measure P , where ( λ t ) t 0 must satisfy the conditions of Theorem 3. See also ([10], Section 2). In the spirit of [8], assuming an “affine specification” for the market price of risk
λ t k 0 x σ k ( η ) d η λ k   with   some   constant   λ k s
the former conditions are trivially satisfied. Thus, Equation (19) can be equivalently expressed under the physical measure P :
d f t ( x ) = ( d d x f t ( x ) + α ( x ) ) d t + k = 1 σ k ( x ) · d W k ( t ) , f 0 ( x ) = f 0 ( x ) ,
with
α ( x ) = k = 1 σ k ( x ) 0 x σ k ( x ) d s + λ t k = k = 1 ( 1 + λ k ) σ k ( x ) 0 x σ k ( x ) d s .
Since the parameter vectors σ k ( x ) k N and ( λ k ) k N are infinite-dimensional, the model (20), though being linear and autonomous, is genuinely non-parametric from the statistical point of view.
Arguing similar to ([4], p. 42), the stochastic process Y ( t , x ) t 0 is expected to solve the SPDE:
d Y ( t , x ) = x Y ( t , x ) r t + 1 2 k = 1 ( 1 + λ k ) I x σ k 2 ( x ) d t     + k = 1 I x σ k ( x ) · d W k ( t )   for   t > 0 , x ( 0 , x max ) , Y ( t , 0 ) = 0   for   t > 0 , Y ( 0 , x ) = x y 0 ( x )   for   x [ 0 , x max ] ,
where y 0 is the initial value for the yield curve y t ( · ) at t = 0 and
r t = f t ( 0 ) y t ( 0 )
denotes the short rate. We would like to draw the reader’s attention to the minus sign in front of the first term on the right-hand side of the first equation in (21). Since the boundary condition is now prescribed at x = 0 , the information propagates from left to right, i.e., the abstract operator is now A = x subject to homogeneneous Dirichlet boundary conditions at x = 0 .
On the strength of Hille & Yosida theorem, the operator A with the maximal domain
D ( A ) = Y H   |   A Y H ,   Y ( 0 ) = 0
(note that the point evaluation is well-defined in view of Sobolev embedding theorem) generates a strongly continuous semigroup ( S t ) t 0 of contractions on
H = L 2 ( 0 , x max ; R )
endowed with the standard scalar product. The semigroup ( S t ) t 0 is now the right-shift semigroup given as
S t Y ( x ) = 0 , 0 x t x max , Y ( x t ) , 0 t x x max .
Assuming r t is a scalar Itô process independent of the cylindircal Wiener process W ( t ) t 0 , the concept of a mild solution can be developed similar to Section 3.1.
To make Equation (21) useful for prediction purposes, we need to specify the evolution of r t . To this end, we utilize the classical Vasiček’s short rate model [37]:
d r t = α ( r ˜ * r t ) d t + β   d W ˜ t 0   for   t > 0 , r 0 = r 0 ,
where α , β , r ˜ * > 0 are parameters, r 0 is a given random variable independent of W ˜ 0 ( t ) t 0 and W ˜ 0 ( t ) t 0 is, in turn, independent of y 0 and W ˜ ( t ) t 0 . Since the evolution in Equation (22) is taking place under the risk-neutral measure P ˜ 0 , we perform a change of measures using Girsanov’s theorem assuming
W ˜ t 0 = W t 0 + 0 t λ 0 d s = W t 0 + λ 0 t .
Thus, under the physical measure P 0 , the evolution in Equation (22) assumes the form
d r t = α ( r * r t ) d t + β   d W t 0   for   t > 0 , r 0 = r 0
with r * = r ˜ * + λ 0 α and a standard Wiener process ( W t 0 ) t 0 . In the following, we will work with Equations (21) and (23) under the physical measures P and P 0 . The key feature of Equations (21) and (23) is that they are autonomous, which is crucial requirement if a model is to be used in forecasting.
Since our dataset is discrete, it is natural to discretize Equations (21) and (23) as well. To this end, we employ (semi-)implicit and explicit Euler & Maruyama schemes for time integration and a two-point finite-difference difference quotient approximation for A. Selecting a time grid { t j   |   j = 0 , , m 1 } and a “space” grid { x i   |   i = 0 , , n 1 } with t j = j ( Δ t ) , x i = i ( Δ x ) , Δ t = T max m 1 , Δ x = x max n 1 , we obtain
Y ( t j + 1 , x i ) Y ( t j , x i ) Δ t = Y ( t j + 1 , x i + 1 ) Y ( t j + 1 , x i ) Δ x r t j + 1 2 k = 1 ( 1 + λ k ) I x i σ k 2 ( x i )     + k = 1 I x σ k ( x i ) W k ( t j + 1 ) W k ( t j ) Δ t       for   j = 0 , , m 2 ,   i = 0 , , n 2 , Y ( t j , 0 ) = 0   for   j = 0 , , m 1 , Y ( t 0 , x i ) = x i y 0 ( x i )   for   i = 0 , , n 1
and
r t j + 1 r t j Δ t = α r * r t j + β W t j + 1 0 W t j 0 Δ t   for   j = 0 , , m 2 , r t 0 = r 0 .
Using standard results from numerical analysis of S(P)DEs, both numerical schemes can be proved to be unconditionally feasible and stable (i.e., for any Δ t , Δ x > 0 ). The consistency also holds true, but is irrelevant for our purposes as Δ t , Δ x remain fixed. A numerical implementation in Matlab® is given in the Supplementary Materials.
The historical yield rate dataset [7] used in this paper contains m = 249 time perios for the year 2018. Restricting the times to maturity to 0 to 12 months, the equidistant x-grid contains n = 769 nodes ranging from 0 to 32 months. Thus, Δ t = 12 m 1 = 0.0484 and Δ x = 32 n 1 = 0.0417 .

5.2. Parameter Estimation

Before the model in Equations (24) and (25) can be used for prediction purposes, the unknown functional parameter vector σ σ k ( · ) j N and the scalar parameters α , β need to be estimated. Since σ is infinite-dimensional, the underlying estimation problem is a genuinely non-parametric and will be regularized through a type of PCA (Principal Component Analysis) truncation ([6], Chapter 1.7.1). In the PDE optimization literature, this approach is commonly known as POD (Proper Orthogonal Decomposition) and used as a dimension reduction technique [41].
We begin with Equation (25). The maximum likelihood (ML) estimators for α , r * and β are given as
α ^ = 1 Δ t log m j = 1 m r t j r t j 1 j = 1 m r t j j = 1 m r t j 1 m j = 1 m R t j 1 2 j = 1 m r t j 1 2 , r ^ * = 1 m 1 exp ( α ^ Δ t ) j = 1 m r t j exp ( α ^ Δ t ) j = 1 m r t j 1 , β ^ 2 = 2 α ^ m 1 exp ( 2 α ^ Δ t ) j = 1 m r t j r t j 1 exp ( α ^ Δ t ) r ^ 1 exp ( α ^ Δ t ) 2
(cf. [42,43,44]). The ML estimates for the 2018 historic data are displayed in Table 1 below.
Since no maximum likelihood estimator is known for σ in Equation (24), we will develop a custom least-squares-like method instead. To this end, we rewrite Equation (24) as
V t j = 1 2 k = 1 ( 1 + λ k ) I x i σ k 2 ( x i ) + k = 1 I x σ k ( x i ) W k ( t j + 1 ) W k ( t j ) Δ t       for   j = 0 , , m 2 ,   i = 0 , , n 2 , Y ( t j , 0 ) = 0   for   j = 0 , , m 1 , Y ( t 0 , x i ) = x i y 0 ( x i )   for   i = 0 , , n 1
with
V t j = Y ( t j + 1 , x i ) Y ( t j , x i ) Δ t + Y ( t j + 1 , x i + 1 ) Y ( t j + 1 , x i ) Δ x r t j .
According to Equation (26), V t j ’s are i.i.d. Gaussian elements of H . In view of the fact that we interested in predicting the yields y t ( x ) given as y t ( x ) = 1 x Y ( t , x ) , for regularization purposes, we will rather employ the stronger topology of D ( A ) .
As mentioned before, to render σ k ( x ) ’s identifiable, we assume the latter are orthogonal in H . This implies I x σ k are orthogonal in D ( A ) since
σ k , σ l H = A I x σ k , A I x σ l H = I x σ k , I x σ l D ( A ) .
Also note that ( I x σ ) 2 H since I x σ k D ( A ) L ( 0 , x max ) for σ k H on the stength of Sobolev’s imbedding theorem and Hölder’s inequality.
We proceed by forming the sample { A Δ x V ( t 0 ) , , A Δ x V ( t m 2 ) } with the discretized negative derivative operator
A Δ x = 1 Δ x 1 0 0 0 0 0 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 1 1 0 0 0 0 0 1 1
and compute the covariance operator in the scalar product of
H Δ x = Y = ( Y t 1 , , Y t m 1   with   Y H Δ x 2 = ( Δ x ) i = 1 m 1 Y 2 ( x i ) .
Let C Δ x denote the resulting (bounded) covariance operator/matrix given as
C Δ x = 1 m 2 j = 0 m 1 A Δ x V ( t j ) m Δ x A Δ x V ( t j ) m Δ x
with the sample mean element/vector
m Δ x = 1 m 1 j = 0 m 2 A Δ x V ( t m ) ,
Performing spectral decomposition
C Δ x = k = 1 p ν k ϕ k ϕ k
for some “principal variances” ν 1 ν 2 ν p 0 and unit ϕ k H Δ x (cf. [6], Chapter 1.7.1), we choose the first p p “principal modes” responsible for, say, 99% of variation, i.e., we choose the smallest p such that k = 1 p ν k k = 1 p ν k 0.99 .
This truncation step, in addition to the D ( A ) -topology employed, renders the procedure numerically and statistically well-posed. Figure 3 displays the six principal modes responsible for 99.32% of variation in A V ( t j ) ’s for the 2018 integrated forward rate data. Summarizing, I x σ k ’s and, thus, σ k can now be estimated via
I x σ ^ k = ( Δ t ) ν k 1 / 2 I Δ x ϕ k , 1 k p , 0 , k > p   and   σ ^ k = ( Δ t ) ν k 1 / 2 ϕ k , 1 k p , 0 , k > p ,
where the discrete negative antiderivative operator
I Δ x : Y ( Δ x ) l = 0 i Y ( t l ) i = 0 , , n 1
comprises the (isometric isomorphic) inverse of A Δ x . Since the covariance operator C Δ x is (in general) non-spherical, we expect that our method is consistent, but does not necessarily need to exhibit optimal efficiency. The latter can be improved with a type of “backfitting” procedure, which is beyond the scope of this article.
Once σ k ( x ) ’s are estimated, we can estimate λ k ’s using the least squares method. Computing the “design matrix”
X Δ x = A Δ x I x σ 1 ^ , , A Δ x I x σ p ^
and recalling the sample mean m Δ x from Equation (28), the Gauss & Markov theorem yields an estimate for the “truncated” λ vector:
λ ^ = 2 ( X Δ x ) X Δ x 1 / 2 ( X Δ x ) m Δ x 1 .
(Note that we first estimated 1 2 ( λ k 1 ) and then solved for λ k .) Applied to the 2018 data, the estimates for the p = 6 “principal” λ k ’s in given in Table 2.

5.3. Prediction

Now, that the estimates α ^ , r ^ * , β ^ and σ ^ k , λ ^ k are at hand, Equations (24) and (25) can be used to predict the dynamics of r t and Y ( t , x ) and, thus, that of y t ( x ) = Y ( t , x ) x . 1 January 2019 was used as t = 0 with the yield curve from 31 December 2018 (multiplied by x) playing the role of the initial value. We generated 10 , 000 sample paths for the standard Wiener increments W k ( t j + 1 ) W k ( t j ) 0 j m 2 , 1 k 6 and W 0 ( t j + 1 ) W k ( t j ) 0 j m 2 . Here, we chose m = 22 which corresponds to predicting the integrated forward rates over a one-month time horizon. Further, 10 , 000 sample paths of the discretized solution processes r Δ t , Δ x , l ,   Y Δ t , Δ x , l , l = 1 , , 10 , 000 were computed by solving the difference Equations (24) and (25). The predicted yield curves were then formed as y t j Δ t , Δ x , l ( x i ) = x i Y Δ t , Δ x , l .
To empirically assess the statistical performance of our method, we chose t = 1 February 2020 and compared our predicted curves with the historically observed one. Figure 4 displays the observed yield rate curve (solid) and the estimated mean yield curve y ¯ t j Δ t , Δ x , l ( x i ) = 1 10 , 000 l = 1 10 , 000 y t j Δ t , Δ x , l ( x i ) (dotted). Further, the dashed curves represent the lower and upper 99% pointwise prediction bounds. Note that the overall confidence level of the resulting tube region is less than 99% since no Bonferroni correction was performed. Figure 4 further contains 5 sample solution paths, which are all (essentially) contained in the prediction region. The estimated mean yield curve is also contained in the prediction region and does not substantially deviate from the observed curve. The prediction region is relatively tight suggesting good predictive power of our model.
The left panel of Figure 5 displays the historical yield curve rates for January 2019 (22 business days), while the right panel gives the estimated mean yield curve rate for the same time window based on 10 , 000 replications. We conjecture that the slight overestimation of the yield curve for large maturities may be related with the non-decreasing behavior of the long rates (see Section 4.1).
In sum, the presented model can be used both for prediction and estimation purposes as well as uncertainty quantification. Further potential applications include stochastic optimization and optimal stopping time problems, etc. The full set of Matlab® codes used to produce the results of Section 5 are available at https://github.com/mpokojovy/abstract.HJM.

6. Discussion and Conclusions

High level of abstraction and numerous prerequisites in various areas of probability theory, stochastic processes and martingale theory, functional analysis and operator semigroup theory, etc., are among the barriers encountered by researchers and practitioners when working with the abstract Heath-Jarrow-Morton (HJM) framework. We attempted to alleviate these difficulties by giving an accessible self-contained review of mathematical results related to the abstract HJM model. Based on real data, we also demonstrated how an infinite-dimensional HJM model can be callibrated and applied to forecasting US Treasury yield curve rates. Our future work will include empirical investigation of the nonlinear models as well as analysis of multiple zero coupon bonds.

Supplementary Materials

The following are available online at https://www.mdpi.com/2227-7390/9/2/114/s1, Matlab® codes implementing the example in Section 5.

Author Contributions

Formal analysis, V.M.; Investigation, V.M.; Methodology, M.P.; Software, M.P. and V.M.; Visualization, M.P.; Writing—original draft, V.M. (Section 2, Section 3 and Section 4 and Appendix A, Appendix B and Appendix C) and M.P. (Section 1 and Section 5); Writing—review & editing, M.P. and V.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The dataset analyzed in this paper is available at https://github.com/mpokojovy/abstract.HJM.

Acknowledgments

The Authors thank the anonymous reviewers for their valuable comments and improvement suggestions.

Conflicts of Interest

The authors declare no conflict of interests.

Appendix A. Q-Wiener Process

Appendix A.1. Representation of Q-Wiener Process

Definition A1
(Q-Wiener process). Let U be a Hilbert space and Q L ( U ) be symmetric, non-negative and with finite trace. A U -valued stochastic process ( W ( t ) ) t [ 0 , T ] on a probability space ( Ω , F , P ) is called a Q-Wiener process if
  • W ( 0 ) = 0 ,
  • W ( t ) has continuous trajectories P -almost surely,
  • W ( t ) has independent increments, e.g., the random variables W ( t 1 ) , W ( t 2 ) W ( t 1 ) , · , W ( t n ) W ( t n 1 ) are independent for all 0 t 1 < t 2 < < t n T and n N ,
  • The increments have Gaussian laws, e.g., W ( t ) W ( s ) N 0 , ( t s ) Q for all 0 s t T .
Proposition A1
(Representation of Q-Wiener process). Let Q L ( U ) be symmetric, non-negative and with finite trace.
Let { e k } k N denote an orthonormal basis of U comprised of eigenvectors of Q associated with eigenvalues ( λ k ) k N . For a stochastic U -valued process ( W ( t ) ) t [ 0 , T ] to be Q-Wiener process, it is necessary and sufficient that
W ( t ) = k N λ k β k ( t ) e k   for   all   t [ 0 , T ] ,
where { β k } k N is a vector of independent real-valued Brownian motions with respect to a probability space ( Ω , F , P ) . The latter series is convergent in the topology of L 2 ( Ω , F , P ) ; C ( [ 0 , T ] , U ) .

Appendix A.2. Q-Wiener Process with Respect to Filtration

Definition A2
(Filtration). Let T [ 0 , ) . A filtration on a probability space ( Ω , F , P ) is a family of σ-algebras { F t } t [ 0 , T ] with the property F s F t F for all 0 s t T .
Definition A3
(Adapted process). A stochastic process ( X ( t ) ) t [ 0 , T ] is said to be adapted to a filtration { F t } t [ 0 , T ] if for every t [ 0 , T ] the random variable X t is F t -measurable.
Definition A4
(Q-Wiener process). A stochastic process ( W ( t ) ) t [ 0 , T ] satisfying
  • ( W ( t ) ) t [ 0 , T ] is adapted to { F t } t [ 0 , T ] ,
  • the increments W ( t ) W ( s ) are independent of F s for all 0 s t T
is referred to as a Q-Wiener process.
Proposition A2.
Let ( Ω , F , P ) be a probability space. Any U -valued Q-Wiener process ( W ( t ) ) t [ 0 , T ] is a Q-Wiener process with respect to the normal filtration { F t } t [ 0 , T ] .

Appendix B. Conditional Expectation and Martingales in Banach Spaces

Proposition A3
(Existence of conditional expectation). Let ( Ω , F , P ) be a probability space, B a separable Banach space, X a Bochner integrable B-valued random variable and Y a σ-algebra contained in F . Then, up to a P -negligible set, there exists a unique integrable B-valued, Y -measurable random variable Z such that
A X d P = A Z d P   for   all   A Y .
The random variable Z will be denoted by E ( X | Y ) and referred to as the conditional expectation of X given Y . Moreover, we have E ( X | Y ) B E X B | Y .
Proposition A3 allows us to extend the definition of a martingale to an F t -martingale.
Definition A5
( F t -martingale). A stochastic process ( X ( t ) ) t 0 on a probability space ( Ω , F , P ) with values in a separable Banach space B is called a martingale with respect to a filtration { F t } t 0 (or an F t -martingale) if
  • E ( X ( t ) B ) < ,
  • X ( t ) is F t -measurable for all t 0 ,
  • E ( X ( t ) | F s ) = X ( s ) P -almost surely for all 0 s t < .
Definition A6
( F t -submartingale). A real-valued integrable adapted stochastic process ( X ( t ) ) t 0 is called a submartingale if E ( X ( t ) | F s ) X ( s ) for all 0 s t < .
Proposition A4.
If ( X ( t ) ) t 0 is a B-valued F t -martingale and p [ 1 , ) , then X B p is a real-valued F t -submartingale.
Proposition A5.
Let ( Ω , F , P ) be a probability space and T > 0 . Then any B-valued Q-Wiener process ( W ( t ) ) t [ 0 , T ] with respect to a normal filtration { F t } t [ 0 , T ] is a continuous square integrable martingale.

Appendix C. Stochastic Integrals Driven by a Standard Cylindrical Wiener Process

Appendix C.1. Cylindrical Q-Wiener Process

Definition A7
(Cylindrical Q-Wiener process). Let Q be a linear bounded operator, but not necessarily with finite trace, on a Hilbert space U . We call ( W ( t ) ) t [ 0 , T ] a cylindrical Q-Wiener process over U if
W ( t ) = k N λ k 1 2 β k ( t ) e k ,
where { ( β k ( t ) ) t [ 0 , T ] } k N is a family of independent Brownian motions, { λ k } k N is a set of eigenvalues of Q with corresponding eigenvectors { e k } k N which form a basis in U . If Q = I , ( W ( t ) ) t [ 0 , T ] is called a standard cylindrical Wiener process.
Since this series may not converge in L 2 ( ( Ω , F , P ) ; U ) if Q has no finite trace, we will use the following propositions.
Proposition A6.
Let Q be a linear bounded operator, but not necessarily with finite trace, on a Hilbert space U . Let U 1 be a Hilbert space such that there exists a Hilbert-Schmidt embedding J : U 0 : = Q 1 2 ( U ) U 1 . If { e k } k N is an orthonormal basis of U 0 and { ( β k ( t ) ) t [ 0 , T ] } k N is a family of independent real-valued Brownian motions, then there exists Q 1 L ( U 1 ) symmetric, non-negative and with finite trace such that the series
W ( t ) = k N β k ( t ) J ( e k )   for   all   t [ 0 , T ] ,
converges in · M T 2 ( U 1 ) = E sup t [ 0 , T ] ( · ) ( t ) U 1 2 1 2 and defines a Q 1 -Wiener process on U 1 .

Appendix C.2. Integration for Predictable Processes

Definition A8
(Hilbert-Schmidt operator). Given an orthonormal basis { e k } k N of U , we refer to an operator A L ( U , H ) as a Hilbert-Schmidt operator if
k N A e k , A e k H < .
The space of all Hilbert-Schmidt operators from U to H will be denoted by L HS ( U , H ) .
Definition A9
(Elementary process). Let T be positive. A stochastic process ( Φ ( t ) ) t [ 0 , T ] is said to be an elementary process on a Hilbert space U if there is a finite family of disjoint intervals
( s k , t k ] [ 0 , T ] , k { 1 , 2 , , n }   f o r   s o m e   n N ,   k = n ( s k , t k ] = ( 0 , T ] ,
and a corresponding finite family of F s k -measurable random variables { Φ k } k { 1 , 2 , , n } taking a finite number of values in the space of Hilbert-Schmidt operators L HS ( U , H ) such that for all t [ 0 , T ]
Φ ( t , ω ) = k = 1 n Φ k ( ω ) · 1 ( s k , t k ] ( t )   P - a l m o s t   s u r e l y .
Definition A10
(Stochastic integral for elementary process). For a given elementary process Φ L 2 [ 0 , T ] × Ω , L HS ( U , H ) the stochastic integral is defined by
0 T Φ ( t ) d W ( t ) = k = 1 n Φ k W ( t k ) W ( s k )   f o r   s o m e   n N ,
where ( W ( t ) ) t [ 0 , T ] is a standard cylindrical U -valued Wiener process.
Definition A11
(Predictable σ -algebra). A σ-algebra over [ 0 , T ] × Ω generated by all subsets of the form ( s , t ] × A , A F s with { 0 } × A 0 for A 0 F 0 is called a predictable σ-algebra and is denoted by P T .
Definition A12
(Predictable process). A stochastic process ( X ( t ) ) t 0 with values in a Hilbert space U is called predictable if it is P T -measurable.
Proposition A7.
The stochastic integral can be uniquely defined for every predictable process Φ L 2 [ 0 , T ] × Ω , L HS ( U , H ) .
Definition A13
( F t -stopping time). An F t -stopping time is a non-negative real-valued random variable τ such that for all times t 0 the event { τ t } F t .
Definition A14
(Simple F t -stopping time). An F t -stopping time τ is said to be simple if it takes only a finite number of values.
Theorem A1
(Itô isometry). If Φ ( t ) t [ 0 , T ] L 2 ( [ 0 , T ] × Ω , L HS ( U , H ) ) is predictable and ( W ( t ) ) t [ 0 , T ] is a standard U -valued cylindrical Wiener process, then we have
E 0 T Φ ( s ) d W ( s ) H 2 = E 0 T Φ ( s ) L HS ( U , H ) 2 d s = E 0 T tr Φ * ( s ) Φ ( s ) d s .
Proposition A8.
Let ( Φ ( t ) ) t [ 0 , T ] L 2 ( [ 0 , T ] × Ω ,   L HS ( U , H ) ) be predictable, ( W ( t ) ) t [ 0 , T ] be a standard cylindrical U -valued Wiener process and τ be an F t -stopping time such that P ( τ T ) = 1 . Then we have P -almost surely
0 t Φ ( s ) 1 [ 0 , τ ] ( s ) d W ( s ) = 0 τ t Φ ( s ) d W ( s ) .
Proposition A9.
The stochastic integral defined by
0 t Φ ( s ) d W ( s ) : = lim n 0 t 1 [ 0 , τ n ] ( s ) Φ ( s ) d W ( s )   for   all   t [ 0 , T ]
does not depend on the choice of the sequence of stopping times { τ n } n N if lim n τ n = T .

References

  1. Heath, D.; Jarrow, R.; Morton, A. Bond Pricing and the Term Structure of Interest Rates: A New Methodology for Contingent Claims Valuation. Econometrica 1992, 60, 77–105. [Google Scholar] [CrossRef]
  2. Brigo, D.; Mercurio, F. Interest Rate Models—Theory and Practice: With Smile, Inflation and Credit, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  3. Harrison, J.; Pliska, S. Martingales and Stochastic Integrals in the Theory of Continuous Trading. Stoch. Proc. Appl. 1981, 11, 215–260. [Google Scholar] [CrossRef] [Green Version]
  4. Filipović, D. Consistency Problems for HJM Interest Rate Models. Ph.D. Thesis, Swiss Federal Institute of Technology, Zurich, Switzerland, 2001. [Google Scholar]
  5. Musiela, M. Stochastic PDEs and Term Structure Models. In Journées Internationales de Finance; IGR-AFFI: La Baule, France, 1993. [Google Scholar]
  6. Carmona, R.; Tehranchi, M. Interest Rate Models: An Infinite Dimensional Stochastic Analysis Perspective, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  7. US Department of the Treasury. Daily Treasury Yield Curve Rates. Available online: https://www.treasury.gov/resource-center/data-chart-center/interest-rates/Pages/TextView.aspx?data=yield (accessed on 14 April 2020).
  8. Dai, Q.; Singleton, K. Specification analysis of affine term structure models. J. Financ. 2000, 55, 1943–1978. [Google Scholar] [CrossRef]
  9. Duffee, G. Term premia and interest rate forecasts in affine models. J. Financ. 2002, 57, 405–443. [Google Scholar] [CrossRef]
  10. Cheridito, P.; Filipović, D.; Kimmel, R. Market price of risk specifications for affine models: Theory and evidence. J. Financ. Econ. 2007, 83, 123–170. [Google Scholar] [CrossRef]
  11. Da Prato, G.; Zabczyk, J. Stochastic Equations in Infinite Dimensions, 2nd ed.; Cambridge University Press: New York, NY, USA, 1992. [Google Scholar]
  12. Reed, M.; Simon, B. Methods of Modern Mathematical Physics, 1st ed.; Academic Press Inc.: Cambridge, MA, USA, 1972. [Google Scholar]
  13. Capinski, M.; Kopp, E. Measure, Integral and Probability, 2nd ed.; Springer Undergraduate Mathematics Series; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  14. Rudin, W. Functional Analysis, 2nd ed.; International Series in Pure and Applied Mathematics; McGraw-Hill: New York, NY, USA, 1991. [Google Scholar]
  15. Knoche, C.; Frieler, K. Solutions of Stochastic Differential Equations in Infinite Dimensional Hilbert Spaces and Their Dependence on Initial Data. Ph.D. Thesis, University of Bielefeld, Bielefeld, Germany, 2001. [Google Scholar]
  16. Björk, T. Arbitrage Theory in Continuous Time, 4th ed.; Oxford Finance Series; Oxford University Press: Oxford, UK, 2020. [Google Scholar]
  17. Cheyette, O. Interest Rates Models. In Investment Management for Insurers; Babbel, D., Fabozzi, F., Eds.; Frank J. Fabozzi Associates: New Hope, PA, USA, 1999. [Google Scholar]
  18. Kamizono, K.; Kariya, T. An implementation of the HJM model with application to Japanese interest futures. Financ. Eng. Jpn. Market. 1996, 3, 151–170. [Google Scholar] [CrossRef]
  19. Sabelli, C.; Pioppi, M.; Sitzia, L.; Bormetti, G. Multi-curve HJM modelling for risk management. Quant. Financ. 2018, 18, 563–590. [Google Scholar] [CrossRef] [Green Version]
  20. Fernando, W. American Options Pricing Using HJM Approach. Ph.D. Thesis, UNC Charlotte, Charlotte, NC, USA, 2017. [Google Scholar]
  21. Schmidt, T. An infinite factor model for credit risk. Int. J. Theor. Appl. Financ. 2006, 9, 43–68. [Google Scholar] [CrossRef]
  22. Babbel, D.; Fabozzi, F. Investment Management for Insurers; Babbel, D., Fabozzi, F., Frank, J., Eds.; Fabozzi Associates: New Hope, PA, USA, 1999. [Google Scholar]
  23. Miron, P.; Swannell, P. Pricing and Hedging Swaps, 1st ed.; Euromoney Publications PLC. Nestor House, Playhouse Yard: London, UK, 1992. [Google Scholar]
  24. Black, F.; Scholes, M. The Pricing of Options and Corporate Liabilities. J. Pol. Econ. 1973, 81, 637–654. [Google Scholar] [CrossRef] [Green Version]
  25. Harrison, J.; Kreps, D. Martingales and arbitrage in multiperiod securities markets. J. Econ. Theor. 1979, 20, 381–408. [Google Scholar] [CrossRef]
  26. Delbaen, F.; Schachermayer, W. The Mathematics of Arbitrage, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  27. Dybvig, P.H.; Ingersoll, J.E.; Ross, S.A. Long forward and zero-coupon rates can never fall. J. Bus. 1996, 59, 1–25. [Google Scholar] [CrossRef]
  28. Brezis, H. Analyse Fonctionnelle, 4th ed.; Thèorie et Applications (French) [Theory and Applications]; Masson: Paris, France, 1993. [Google Scholar]
  29. Bhar, R.; Chiarella, C. The estimation of the Heath-Jarrow-Morton model by use of Kalman filtering techniques. In Computational Approaches to Economic Problems; Springer: Boston, MA, USA, 1997; pp. 113–126. [Google Scholar] [CrossRef] [Green Version]
  30. Pang, K. Calibration of Gaussian Heath, Jarrow and Morton and random field interest rate term structure models. Rev. Derivat. Res. 1998, 2, 315–345. [Google Scholar] [CrossRef]
  31. Jeffrey, A.; Kristensen, D.; Linton, O.; Nguyen, T.; Phillips, P.C. Nonparametric estimation of a multifactor Heath-Jarrow-Morton model: An integrated approach. J. Financ. Econom. 2004, 2, 251–289. [Google Scholar] [CrossRef] [Green Version]
  32. Broszkiewicz-Suwaj, E.; Weron, A. Calibration of the multi-factor HJM model for energy market. Acta Phys. Polonica B 2006, 37, 1455–1466. [Google Scholar]
  33. Bhar, R.; Chiarella, C.; Tô, T.D. A Maximum Likelihood Approach to Estimation of a Class of the Heath-Jarrow-Morton Models. Available online: https://econpapers.repec.org/paper/utsrpaper/80.htm (accessed on 28 October 2020).
  34. Falcò, L.N.; Nave, J. On the calibration of a Gaussian Heath-Jarrow-Morton model using consistent forward rate curves. Quant. Financ. 2011, 11, 495–504. [Google Scholar] [CrossRef]
  35. Li, S.; Dong, Y. Parametric estimation in the Vasicek-type model driven by sub-fractional Brownian motion. Algorithms 2018, 11, 197. [Google Scholar] [CrossRef] [Green Version]
  36. Funahashi, H. An approximate swaption formula in Heath-Jarrow-Morton models. J. Derivat. 2020, 27, 1–21. [Google Scholar] [CrossRef]
  37. Vasicek, O. An equilibrium characterization of the term structure. J. Financ. Econom. 1977, 5, 177–188. [Google Scholar] [CrossRef]
  38. Dörsek, P.; Teichmann, D. Efficient simulation and calibration of general HJM models by splitting schemes. SIAM J. Financ. Math. 2013, 4, 575–598. [Google Scholar] [CrossRef] [Green Version]
  39. Kluge, W. Time-Inhomogeneous Lévy Processes in Interest Rate and Credit Risk Models. Ph.D. Thesis, University of Freiburg, Freiburg, Germany, 2005. [Google Scholar]
  40. Kim, S.; Pokojovy, M.; Wan, X. The taut string approach to statistical inverse problems: Theory and applications. J. Comput. Appl. Math. 2021, 382, 113098. [Google Scholar] [CrossRef]
  41. Sachs, E.W.; Volkwein, S. POD-Galerkin approximations in PDE-constrained optimization. GAMM-Mitteilungen 2010, 33, 194–208. [Google Scholar] [CrossRef]
  42. Tang, C.; Chen, S. Parameter estimation and bias correction for diffusion processes. J. Econom. 2009, 149, 65–81. [Google Scholar] [CrossRef] [Green Version]
  43. Fergusson, K.; Platen, E. Application of maximum likelihood esitmation stochastic shortrate models. Ann. Financ. Econom. 2015, 10, 1–26. [Google Scholar]
  44. Ballestra, L.; Pacelli, G.; Radi, D. A note on Fergusson and Platen: “Application of maximum likelihood estimation stochastic short-rate models”. Ann. Financ. Econom. 2016, 11, 1–7. [Google Scholar] [CrossRef]
Figure 1. Historical yield rate curves for 2018.
Figure 1. Historical yield rate curves for 2018.
Mathematics 09 00114 g001
Figure 2. Historical values for integrated forward rates Y ( t , x ) .
Figure 2. Historical values for integrated forward rates Y ( t , x ) .
Mathematics 09 00114 g002
Figure 3. Six principal modes responsible for over 99% of variation.
Figure 3. Six principal modes responsible for over 99% of variation.
Mathematics 09 00114 g003
Figure 4. Predicting the yield rates on 1 February 2019 based on 2018 historical data.
Figure 4. Predicting the yield rates on 1 February 2019 based on 2018 historical data.
Mathematics 09 00114 g004
Figure 5. Observed yield curve rates (left) and estimated mean yield curve rates (right) for January 2019.
Figure 5. Observed yield curve rates (left) and estimated mean yield curve rates (right) for January 2019.
Mathematics 09 00114 g005
Table 1. ML estimates for the Vasicek model based on 2018 historic short rate data.
Table 1. ML estimates for the Vasicek model based on 2018 historic short rate data.
α ^ r ^ * β ^
0.079473−0.639520.15449
Table 2. Least-squares estimates for the first p = 6 parameteres λ k .
Table 2. Least-squares estimates for the first p = 6 parameteres λ k .
λ ^ 1 λ ^ 2 λ ^ 3 λ ^ 4 λ ^ 5 λ ^ 6
−0.710786.880.48587−0.74969−0.73313−0.84993
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maltsev, V.; Pokojovy, M. Applying Heath-Jarrow-Morton Model to Forecasting the US Treasury Daily Yield Curve Rates. Mathematics 2021, 9, 114. https://doi.org/10.3390/math9020114

AMA Style

Maltsev V, Pokojovy M. Applying Heath-Jarrow-Morton Model to Forecasting the US Treasury Daily Yield Curve Rates. Mathematics. 2021; 9(2):114. https://doi.org/10.3390/math9020114

Chicago/Turabian Style

Maltsev, Valerii, and Michael Pokojovy. 2021. "Applying Heath-Jarrow-Morton Model to Forecasting the US Treasury Daily Yield Curve Rates" Mathematics 9, no. 2: 114. https://doi.org/10.3390/math9020114

APA Style

Maltsev, V., & Pokojovy, M. (2021). Applying Heath-Jarrow-Morton Model to Forecasting the US Treasury Daily Yield Curve Rates. Mathematics, 9(2), 114. https://doi.org/10.3390/math9020114

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop