(Lecture Notes in Economics
(Lecture Notes in Economics
Founding Editors:
M. Beckmann
H. P. Kiinzi
Managing Editors:
Prof. Dr. G. Fandel
Fachbereich Wirtschaftswissenschaften
Femuniversitat Hagen
Feithstr. 140/AVZ II, 58084 Hagen, Germany
Prof. Dr. W. Trockel
Institut fUr Mathematische Wirtschaftsforschung (IMW)
Universitat Bielefeld
Universitatsstr. 25, 33615 Bielefeld, Germany
Co-Editors:
C. D. Aliprantis, Dan Kovenock
Editorial Board:
P. Bardsley, A. Basile, M.R. Baye, T. Cason, R. Deneckere, A. Drexl,
G. Feichtinger, M. Florenzano, W Giith, K. Inderfurth, M. Kaneko, P. Korhonen,
W. Kiirsten, M. Li Calzi, P. K. Monteiro, Ch. Noussair, G. Philips, U. Schittko,
P. Schonfeld, R. Selten, G. Sorger, R. Steuer, F. Vega-Redondo, A. P. Villamil,
M. Wooders
Springer-Verlag Berlin Heidelberg GmbH
B. Philipp Kellerhals
Springer
Author
Dr. B. Philipp Kellerhals
Deutscher Investment-Trust
Gesellschaft fUr Wertpapieranlagen mbH
Mainzer LandstraBe 11-13
60329 Frankfurt am Main, Germany
ISSN 0075-8450
ISBN 978-3-540-42364-5 ISBN 978-3-662-21901-0 (eBook)
DOI 10.1007/978-3-662-21901-0
This work is subject to copyright. All rights are reserved, whether the whole or part
of the material is concerned, specifically the rights of translation, reprinting, re-use
of illustrations, recitation, broadcasting, reproduction on microfilms or in any other
way, and storage in data banks. Duplication of this publication or parts thereof is
permitted only under the provisions of the German Copyright Law of September 9,
1965, in its current version, and permission for use must always be obtained from
Springer-Verlag. Violations are liable for prosecution under the German Copyright
Law.
http://www.springer.de
© Springer-Verlag Berlin Heidelberg 2001
Originally published by Springer-Verlag Berlin Heidelberg New York in 2001.
The use of general descriptive names, registered names, trademarks, etc. in this
publication does not imply, even in the absence of a specific statement, that such
names are exempt from the relevant protective laws and regulations and therefore free
for general use.
Typesetting: Camera ready by author/editors
Cover design: design & production, Heidelberg
Printed on acid-free paper SPIN: 10844260 55/3142/du 5432 I 0
To my parents
Foreword
Straight after its invention in the early sixties, the Kalman filter approach
became part of the astronautical guidance system of the Apollo project and
therefore received immediate acceptance in the field of electrical engineer-
ing. This sounds similar to the well known success story of the Black-Scholes
model in finance, which has been implemented by the Chicago Board of Op-
tions Exchange (CBOE) within a few month after its publication in 1973.
Recently, the Kalman filter approach has been discovered as a comfortable
estimation tool in continuous time finance, bringing together seemingly un-
related methods from different fields.
Dr. B. Philipp Kellerhals contributes to this topic in several respects.
Specialized versions of the Kalman filter are developed and implemented
for three different continuous time pricing models: A pricing model for
closed-end funds, taking advantage from the fact, that the net asset value is
observable, a term structure model, where the market price of risk itself is a
stochastic variable, and a model for electricity forwards, where the volatility
of the price process is stochastic. Beside the fact that these three models can
be treated independently, the book as a whole gives the interested reader a
comprehensive account of the requirements and capabilities of the Kalman
filter applied to finance models. While the first model uses a linear version
of the filter, the second model using LIBOR and swap market data requires
an extended Kalman filter. Finally, the third model leads to a non-linear
transition equation of the filter algorithm.
Having some share in the design of the first two models, I am rather
impressed by the potential and results of the filtering approach presented
by Dr. Kellerhals and I hope, that this monograph will play its part in
making Kalman filtering even more popular in finance.
2 Stochastic Environment 7
4 Filtering Algorithms 15
4.1 Linear Filtering 16
4.2 MMSE . . . . . . 17
4.3 MMSLE . . . . . 21
4.4 Filter Recursions 22
4.5 Extended Kalman Filtering 24
5 Parameter Estimation 29
II Pricing Equities 35
6 Introduction 37
6.1 Opening Remarks . . . . . . . . 37
6.2 The Case of Closed-End Funds 38
7 Valuation Model 43
7.1 Characteristics of Closed-End Funds 43
7.2 Economic Foundation . . . . . . . 47
7.3 Pricing Closed-End Fund Shares .. . 50
XII Contents
Bibliography 231
1 Overview of the Study
The research objective of this study can be grasped by the question: Are
Kalman filter algorithms applicable to efficiently estimate financial models
of contingent claim pricing? This formulation asks for two more clarifica-
tions: First, which financial pricing models for which markets are selected
and second, which type of Kalman filters are applied? In answering these
questions, we hope to generate both theoretical and econometric contribu-
tions to the scientific field of financial economics.
The Kalman filter framework was originally developed by Kalman (1960)
and Kalman and Bucy (1961) for applications in aerospace engineering and
has successfully been used in electrical engineering. Only in the last decade,
marked by Harvey (1989), Kalman filters have found ground as economet-
ric tool for economic and financial estimation problems. The technique of
Kalman filtering generally leads to more efficient parameter estimates since
it imposes greater theoretical modeling restrictions on the data than other
commonly used methods. Given a hypothesized financial pricing model, the
estimation procedure of Kalman filtering imposes cross-sectional as well as
time-series restrictions, making it a true maximum likelihood method. We
use different types of Kalman filter algorithms to estimate current pricing
problems in the field of financial economics. Three major financial markets
are to be examined for which we select the equity market, the bond market,
and the electricity market. In each market we derive new contingent claim
valuation models to price selected financial instruments in continuous time.
The decision criterium for choosing a continuous time modeling framework is
the richness of the stochastic theory available for continuous time processes
with Merton's pioneering contributions to financial economics, collected in
Merton (1992). The continuous time framework, reviewed and assessed
by Sundaresan (2000), allows us to obtain analytical pricing formulae that
2 1. Overview of the Study
2. For the bond market we present a new term structure model in part III
which belongs to the affine term structure models. Motivated by the
incomplete market situation of fixed income instruments we propose a
pricing model that shifts from modeling static tastes to beliefs of the
dynamic behavior of tastes. Besides a stochastic short interest rate
the proposed two-factor term structure model allows for a stochastic
behavior of the market price of interest rate risk. Given the theoretical
1. Overview of the Study 3
Typical financial pricing models in continuous time are built upon hypoth-
esized stochastic processes to assess the models' dynamics over time. A
common and convenient method to describe the evolution of chosen state
variables over time is to specify appropriate stochastic differential equations.
By these equations we try to capture possibly all behavior of the state vari-
able during time evolution or equivalently set up a hypothesis on how the
state variables could evolve over time with the final objective of finding the
real data generating process underlying the analyzed sample of data. To
set up such stochastic differential equations we first take a look at stochas-
tic processes in general defined on a probability space which represents the
uncertainty of the intertemporal stochastic economy.
Definition 2.0.1 (Probability Space) The underlying set-up for stochas-
tic processes consists of a complete probability space (n, F, JPl), equipped with
a filtmtion, i.e. a non-decreasing family IF = {Ft}tET of sub-a-fields of F :
Fs ~ F t ~ F for 0 ~ s < t < 00 where T = [0,00). Here, F t represents
the information available at time t, and the filtmtion {Ft }tET represents
the information flow evolving over time and accruing to all agents in the
economy.
Definition 2.0.2 (Stochastic Process) A scalar (n-vector) stochastic
process, denoted by {XthET' is a family of mndom variables (n-vectors)
indexed by the pammeter set T, where the pammeter t will refer to time
in our applications. The process is defined on a filtered probability space
(n, F, JPl, IF) with values in lR,n. We say the process is adapted if X t E F t
for each time t, i.e. X t is Ft-measumble; thus, X t is known when :ft
is known. Further, if a filtmtion is genemted by a stochastic process,
i.e. F t = a (Xs; 0 ~ s ~ t), we call F t the natuml filtmtion of the process
{XthET' Thus, a process is always adapted to its natuml filtmtion.
8 2. Stochastic Environment
Remark 2.0.4 Note that the process is measurable with respect to two dis-
tinct spaces: For every t E T the process X t is a random variable with
w ~ Xdw);w E n; if we fix the value for wEn we obtain the mapping
t ~ X t (w); t E T, a so-called path of X t . Thereby, we need to look at a
stochastic process as a function of two distinct variables (t, w) ~ X (t, w)
from T x n into the space JR 1 •
lThe earliest modeling is by Bachelier (1900), who used this type of stochastic pro-
cesses to describe stock price movements. To cover jumps in the dynamics of a security,
one can add a Poisson process component; see, for example, Merton (1976).
2. Stochastic Environment 9
E [etl = 0, and
for s = t
E [ete~l = { H t~'l/J)
otherwise
with the transition matrix <Pt ("p), an additive component Ct ("p) and a Gaus-
sian noise term 'fit ("p) with
E [l1tl = 0, and
for s = t
E [l1tl1~l = { Qt ~"p)
otherwise.
of the likelihood must be guaranteed for all possible parameter vectors, be-
cause the true parameter structure is not known. For an unidentified model
example in the case of maximum likelihood estimation using the Kalman
filter see Hamilton (1994b, p. 387 f.).
4 Filtering Algorithms
According to the assumed state space form of equations (3.1) and (3.2),
we specify stochastic processes of explaining factors eo
which describe the
state of the system. Instead of being able to observe the factors directly, we
can only observe some noisy function Yt of et.
The problem of determining
the state of the system from noisy measurements Yt is called estimation.
The special estimation problem of filtering has the object of obtaining an
et
expression for the optimal estimate of given the observations up to time t.
The most successful result of this kind is that obtained for linear systems by
Kalman (1960), Kalman and Bucy (1961), and Kalman (1963) which we are
further dealing with. In general, the estimation problem can be classified
upon the available and processed information into three different problems
according to the following definition.
results in the minimum mean square estimator (MMSE) as the best or opti-
mal estimator with respect to any quadratic function of the estimation error.
i. e. in this case the mean square error matrix is equal to the variance-
covariance matrix.
We will find that the Kalman filter provides an optimal solution to the
presented problems of prediction and updating. In deriving the Kalman
filter we can choose among several approaches available in the literature
that are linked to different interpretations of the filter. Besides the original
derivation in the work by Kalman (1960), where he used the idea of orthog-
onal projection (see, for example, Aoki (1990) and Brockwell and Davies
(1987)), we can derive the Kalman filter from the properties of a multi-
variate normal distribution as shown, for example, in Harvey (1989), while
Burridge and Wallis (1988) and Hamilton (1994b), for example, deduce the
Kalman filter to be the minimum mean square linear estimator. Moreover,
alternative derivations are found in Jazwinski (1970). In the following we
present the elementary derivation of the Kalman filter under the normal-
ity assumption which yields the interpretation of the Kalman filter as the
optimal filter in the sense of a minimum mean square estimator (MMSE).
Also, we obtain the Kalman filter for the general case of non-normality by
exploiting the given linear relationships of the observations and dynamics
stated respectively in the state space equations (3.1) and (3.2); in this case
we speak of the Kalman filter as the optimal filter in the sense of a minimum
mean square linear estimator (MMSLE).
4.2 MMSE
In order for the Kalman filter to yield a MMSE, we assume that the ad-
ditive error terms et and "1t are independently and normally distributed.
FUrthermore, the error terms are treated as independent of the initial state
18 4. Filtering Algorithms
In the last step, the updating step, we update the inference on etlt-l
by including the newly available information at time t, which results in
the filtered estimate etit.
Hereby, we can realize the so-called Kalman gain
K t . In order to obtain the updating equations, we first consider the joint
distribution of et and Yt. Since both variables are normally distributed, the
joint distribution given the information F t - 1 is
E [(et - etlt-l)
(et - etlt-l)' 1Ft-I] B~ = :Etlt-lB~
E [(Yt - Ytlt-I) (et - etlt-l)' 1F t-I] = Bt:Etlt-l.
For the further derivation, we take advantage of the following result for
normal variables.
~tlt-l - ~tlt-1B~F~LlBt~tlt-l) .
Comparing this distribution with the result for the conditional moments of
et given the information F t ,
(4.5)
4.3 MMSLE
Dealing with a non-normally distributed state vector but confining ourselves
at the same time to the class of linear estimators, we derive the Kalman filter
as the MMSLE without requiring the normality assumption for the error
terms. The prediction equations (4.1) and (4.2) are obtained similarly to
the previous derivation by taking the conditional expectation and variance
given the information up to time t - l.
For obtaining the updating equations, we begin with stating a general
linear relationship of the updated estimator etit
and the information of the
present sample (yt, at, Ct) and the past information F t - 1 . We will assume
the general linear form
et et - etit
= Ct + <pet-l + '11t - K t (at + B t (Ct + <ptet-l + '11t) + et)
-Lat - MCt - N (et-l - et-l)
(<pt - N - KtBt<p) et-l - (L + K t ) at + (I - M - KtB t ) Ct
+Net-l + (I - KtB t ) '11t - Ktet. (4.7)
For etitbeing the MMSLE the estimation error et needs to be (i) uncondi-
tionally unbiased, i.e. the estimation error has zero expectation, which gives
us the conditions
0,
0, and
for equation (4.7). This yields, with the substitution N = (I - KtB t ) <Pt ,
the estimation error
22 4. Filtering Algorithms
which is equivalent to the expression derived in equation (4.6) for the case
of the Kalman filter resulting in the MMSE.
Initialization of
~and;
at Time t=O
The system is initialized
with specific (prior) values
for the state variables and
their covariance matrix.
t=T?
The recursion algorithm is
performed until we reach the
yes last observations.
Choose
other
Parameter The conditional likelihood
Values function is evaluated at the
current parameter values using
'V the results of the Kalman
recursions.
L----no'----<
Maximization
Criterium met?
\. The Kalman algorithm is run
. until we reach the predefined
abortion Criterium.
yes
8
Figure 4.1: Flowchart of the Kalman Filter and MLE
24 4. Filtering Algorithms
Vt = Yt - at - Btetlt-i, and
Ftlt- i = Bt~tlt-iB~ + H t .
Given these intermediate results we then update the a priori estimates yield-
ing the filtered estimates
~t
kxl
= h t (~t-l' "'t ('ljJ) , 'ljJ) , (4.9)
with g (~t' E:t ('ljJ) , 'ljJ) and h (~t-l' "'t ('ljJ) , 'ljJ) denoting the possible non-linear
functional relationships. The error terms E:t ('ljJ) and"'t ('ljJ) are assumed to
follow the properties
'" ar "'t('ljJ)
The overview of the Kalman filter estimation procedure in figure 4.1 shows
the likelihood function and the maximization procedure as the last rele-
vant steps. We start to examine the appropriate log-likelihood function
within our framework. Let 'I/J Ewe ]Rn be the vector of the n unknown
parameters, which various system matrices of the state space model for-
mulated in equations (3.1) and (3.2), or (4.10) and (4.11), as well as the
variances of the measurement errors depend on. The likelihood function of
the state space model is given by the joint density of the observational data
Y = (YT, YT-l, ... ,Yl)
which reflects how likely it would have been to have observed the data if 'I/J
were the true values for the parameters. Using the definition of conditional
probability, we can split the likelihood up into conditional densities using
Bayes's theorem recursively and write the joint density as the product of
conditional densities
where we approximate the initial density function P (Yl; 'I/J) by P (YIlyo; 'I/J).
Moreover, since we deal with a model with a Markovian structure according
to the state space model, i.e. future values of Yu, with u > t, only depend
on (Yt, Yt-l,··· ,Yl) through the current value Yt, the expression in equa-
tion (5.1) is equivalent to conditioning on only the last observed vector of
observations, i.e. it reduces to
30 5. Parameter Estimation
For our purpose of estimating the parameter vector 1/; given the data
Y and the structural form of the specified state space model, we use the
approach of Schweppe (1965) known as the prediction error decomposition
of the likelihood function to be maximized with respect to 1/;. That is, we
express the likelihood function in terms of the prediction errors; since the
variance of the prediction error Vt = Yt - lE [YtIFt-1l is the same as the
conditional variance of Yt, i.e.
we are able to state the density function. The general density function
p (YtIYt-l; 1/;) is given by the Gaussian distribution with conditional mean
lE [YtIFt-ll = a + Betlt-l and the conditional variance-covariance matrix
Cov [YtIFt-d = Ftlt- 1, which takes the form:
(5.3)
1 T
L (y; 1/;) = -"2 L kIn (271") + In IFtlt-11 + v~Ftlt-lVt (5.4)
t=l
instead of the product in equation (5.1).
In order to estimate the unknown parameters from the log-likelihood
function in equation (5.4) we need to choose an appropriate optimization
algorithm to minimize L (y; 1/;) with respect to 1/;. Given the simplicity of
the innovation form of the likelihood function, we can find the values for
the conditional expectation lE [YtIFt-1l and conditional variance-covariance
matrix Cov [YtIFt-1l for every given parameter vector 1/; using the Kalman
filter algorithm. Thus, we are able to compute the value of the conditional
5. Parameter Estimation 31
where \II c ]Rn denotes the parameter space. In the case of implementing
extended Kalman filter algorithms where we use approximations to the con-
ditional normal distribution of equation (5.3) the log-likelihood of equation
(5.5) can be used to yield quasi maximum likelihood parameter estimates. 18
However, the general statistical properties of such estimates must be studied
separately and are assessed in detail by White (1982), Gallant and White
(1988), and White (1994).
In the numerical implementation to obtain the optimal parameter es-
timates we actually minimize the log-likelihood function of equation (5.5)
including its opposite sign. In infinite precision the necessary conditions
for a local minimum :¢ML of L(y; 1/J) are defined by the conditions for the
gradient:
(5.6)
(5.7)
is needed.
In order to estimate the model parameters 1/J M L we decided to choose a
quasi-Newton method, which is similar to the methods used by Lund (1997)
and Nunes and Clewlow (1999) based on Dennis and Schnabel (1996). The
iteration rule for Newton's method for unconstrained minimization is given
by the parameter estimates at the ith iteration step of
18For an exploration of quasi (or pseudo) maximum likelihood estimation methods see,
for example, Gourieroux, Monfort, Renault, and Thognon (1984).
32 5. Parameter Estimation
For the construction of \7 L (y; 1/Ji) and \7 2 L (y; 1/Ji) we use numerical ap-
proximations. We found that the available closed-form solutions19 do not
lead to superior results in that they accumulate numerical errors in their re-
quired extensive calculations while asking for multiples of computer time. In
cases where \7 2 L (y; 1/Ji) is not positive definite Dennis and Schnabel (1996,
sec. 5.5) suggest a modification of the modified Cholesky factorization.
This algorithm changes the model when necessary so that it has a unique
minimizer and uses this minimizer to define the Newton step. Specifically,
at each iteration we take steps of the form
where 1/Ji are the parameter estimates at the ith iteration step.
Moreover, looking at the bottom part of the flowchart in figure 4.1, we
need to discuss how to terminate our optimization algorithm in finite preci-
sion, i.e. we further have to examine the conditions stated in equations (5.6)
and (5.7).20 Although \7 L(y; 1/J) = 0 can also occur at a maximum or sad-
dle point, our globalizing strategy and our method of perturbing the model
Hessian to be positive definite make convergence impossible to maxima and
19 As given, for example, in Harvey (1989, ch. 3.4).
20See Dennis and Schnabel (1996, ch. 7.2).
5. Parameter Estimation 33
IV'L(y;'ljJ)1 ~ c, (5.8)
(5.9)
which is independent of any change in the units of L(y; 'ljJ) and'ljJ. However
the idea of relative change in 'l/Ji and L(y; 'ljJ) breaks down if'ljJi or L(y; 'ljJ)
happen to be close to zero. This problem can be fixed by altering the test
of equation (5.9) to
(5.10)
Pricing Equities
6 Introduction
Thompson (1978).
28 Especially, see the studies by Bonser-Neal, Brauer, Neal, and Wheatley (1990) and
Diwan, Errunza, and Senbet (1995).
29See Lee, Shleifer, and Thaler (1991).
3 0 The hypothesis is based on the noise trader model of De Long, Shleifer, Summers,
and Waldmann (1990) who argue that this concept can explain both the persistent and
variable premia on closed-end funds as well as the creation of such funds.
6.2 The Case of Closed-End Funds 41
Definition 7.1.1 The net asset value (NAV) or the 'book value of share-
holder's equity' for time t is given by
20 0.8
. . . . . .. net asset value
16 - - market price 0.6
12 0.4
.~
Co 8 0.2
0.0 E
~
o+_--.--_,_-----,.-----.----r---.------r-____,.------.----+ -0.2
E 30% 0 52 104 156 208 260
.~ 20%
!!! 10%
~ 0%~~4r~~-_4H4~--_r~~------~~.-------------
.g -10%
'~-20%
m_~%+_--.--_,_--____,.------.-----r----.---____r--____,.------~~~
o 52 104 156 208 260
35For the detailed description of the data set see section 8.1 on page 55.
7.1 Characteristics of Closed-End Funds 45
20 1.00
....... net asset value
16 - - market price 0.75
Ul 12
Q)
0
.~ 0.50
8
4 0.25
0
0.00 Ul
-4 E
~
-8 -0.25
36Further information on the empirical premia of our sample is provided in table 8.2
on page 58.
46 7. Valuation Model
28 0.6
24
net asset value
o 0 0 0 0 0 °
0.4
20 - - market price
..........
'" 16
.g 0.2
Q. 12
0.0
4 '"
E
~
o -0.2
5°1< 0 52 104 156 208 260
§Oo~
°Ee -5%10~V' ,tAA~
~ -10%
.g -15%
.~ -20%
m -25%+--~-~-~~-~-~--r-~-~--~--'
o 52 104 156 208 260
take a look at the values of the premia across all closed-end funds in our
sample (see figure 7.4 on page 47) over the examined period of time, we can
confirm this view based on the whole sample. For the average premia we
observe mainly positive values for the first half of the sample period with a
peak of 13 percent at the end of the first year, which tends to continually fall
to discount values of about 14 percent in the fourth and fifth year. Looking
at the maximum positive and negative values, we notice a higher volatility
for the positive premia whereas the discounts oscillate in a narrow band of
minus 10 to 30 percent.
Moreover, for Sharpe and Sosin (1975) the fluctuations of the empir-
ical premia appear to be mean-reverting. For their sample, they find a
slight tendency for the premia to revert to a long-term mean discount value
of approximately seven percent. Further, in an examination of the rela-
tion between closed-end fund premia and returns, Pontiff (1995) finds no
economically motivated explanation for discounted closed-end funds having
higher expected returns than non-discounted closed-end funds, but rather
accounts a mean-reverting nature of the empirical premia for this effect. For
the market of closed-end funds in the UK, the study of Minio-Paluello (1998)
7.2. Economic Foundation 47
60%
50%
40%
Maximum Values
30%
20%
10%
0%
-10%
-20%
-30%
-40%
0 52 104 156 208 260
shows the empirical tendency of premia to revert to their mean. These re-
sults from the literature give rise to capturing this pricing phenomenon in
a stochastic framework of closed-end fund share valuation which we layout
in the next section.
(ii) the secondary market for the shares of the closed-end funds, i. e. the
US host market.
Assumption 7.2.2 (Net Asset Value) For the primary market we model
the price movements for the investment portfolio of the closed-end funds
measured by the reported net asset value as geometric Brownian motion38
as in Ingersoll (1976).
Corollary 7.2.3 For the assumed price behavior of the net asset value we
conveniently work with the stochastic dynamics of the log-transformed vari-
able X t = In N A Vt given by
where dXt is the continuous return on the net asset value portfolio.
(7.3)
which is a type of Ornstein- Uhlenbeck process. In this special case the con-
tinuous premium 7rt is modeled with a volatility of (J1r and reverts with a
speed of /'i, to its long-run mean (). Further, the derived valuation model will
allow for a correlation of the two model factors X t and 7rt as the standardized
Brownian motions are set up as dWx,tdW1r,t = pdt.
Now, we link the price of the closed-end fund shares that are traded on
the secondary market both to the values of their underlying net asset value
and the dynamic premium. We model the market value of the closed-end
fund share prices, denoted by P (X, 7r, t; 'l/J M) with parameter vector 't/J M, as
a derivative security dependent on the values of the N A Vi and the premium
7rt. In analogy to the standard technique known from contingent claim
pricing40, we derive the following stochastic dynamics for P (X, 7r, t; 't/J M)
by applying Ito's formula
dP
(7.4)
P
with J-Lp
'f/x =
39 Specifically,
see the description on page 61 for this relationship.
40See, for example, Merton (1992).
50 7. Valuation Model
model
+ (r _ 12(JX2) ap(X,1r,t;t/J
ax
M) + ( (() _ ) _
K 11"
\) ap(X,1r,t;t/J M) + ap(X,1r,t;t/J
(J7r/\7r a7r at
M)
(7.5)
which describes how the closed-end funds' share price P varies for differ-
ent values of the state variables X and 11" as well as for changing trading
time t. Additionally, to fully reflect the dynamic behavior of the closed-end
funds' share price P (X, 11", t; '¢M) we further need to specify the appropri-
ate boundary condition for this security. To defer the relevant condition, we
exploit the economic behavior of the closed-end funds at the fourth stage
of their life cycle at time T.43 According to Brauer (1984) and Brickley
and Schallheim (1985) the market values of the closed-end funds equal their
net asset values with empirical regularity at the time of the funds liqui-
dation or open-ending. From this fact we obtain the boundary condition
P (X, 11", T; 'l/JM) = eXT for our valuation model.
The next step in deriving the price of closed-end fund shares in our
model is to solve the partial differential equation of equation (7.5). For the
solution to this equation we consider the exponential affine form 44
~/. ) = eA (t)X+B(t)7r+C(t) ,
P (X , 11" , t·,'I'M (7.6)
with model parameters '¢M and factor loadings A (t), B (t), and C (t). Using
the partial derivatives of the proposed solution
P .= ap(x,7r,t;t/JM) = AP
x· ax '
P7r .-
.-
ap(x,7r,t;t/J M) -
a7r -
BP,
PX7r .=
•
a 2p(x,7r,t;t/JM)
axa>.
= ABP' n._ ap(x,7r,t;t/JM) -
.rt·- at -
(A t X t + Bt1l"t + C)t P
43 According to the nomenclature of Lee, Shleifer, and Thaler (1991).
44This type of solution is known from interest rate modeling; see, for example, Duffie
and Kan (1996).
52 7. Valuation Model
(7.7)
A(t) = 1, and
B (t) -.!. (1 - e-K(T-t») .
'"
The last function of interest is C (t). Solving the integral
J
T
+-1 [e - -(77r/\7r
1 \ +1 (7~] (1 - e- K(T - t»)
-P(7x(77r - -
(72
'" '" '" ",2
+~ (1 - e- 2K (T-t») .
4",3
To further derive the market price of the closed-end fund shares, we now
use the calculated functions A (t), B (t), and C (t) along with the suggested
solution in equation (7.6).
7.3 Pricing Closed-End Fund Shares 53
The derived valuation formula states that the log-market price yt should
basically be equal to the log-value X t of the underlying assets in that the
market prices carry the price risk with a factor loading of A (t) = 1. Second,
the closed-end funds' market prices are affected by changing values of the
dynamic premium 7rt. The realizations of the premia influence the market
prices by a negative loading B (t), i.e. closed-end funds with positive values
for the dynamic premium trade at a discount. Finally, the valuation formula
contains a time varying term C (t) which depends on the maturity of the
closed-end funds and the model parameters 'l/JM. In the following chapter
we motivate a possible implementation of the valuation model upon market
data of closed-end funds.
8 First Empirical Results
Notes:
a) Net proceeds are in million $.
b) P-returns denote returns on the market prices of the closed-end fund shares.
c) Returns are annualized log-returns, calculated on a weekly basis, and given in percent.
d) Denotes the correlation coefficient between the P- and NAV-returns.
8.1 Sample Data 57
are not as volatile. The standard deviation of the net asset value returns
averages on 25 percent and is about one-half lower than the 33 percent on
the market returns. 48 Looking at the correlation between the two returns,
we see a strong positive comovement with a correlation coefficient of 0.57
on average.
In table 8.2 we show the descriptive statistics for the empirical premia
of each closed-end fund. Across all closed-end funds in the sample, we find
an average discount of 3.25 percent which is in line with previous research
on emerging market closed-end funds. 49 The range of premia is given by
the spanning from a 20 percent discount to positive premia of 15 percent
on average. However, only eight out of the 33 closed-end funds trade at
positive premia on average. Especially, the four closed-end funds including
the symbols IGF, IF, KF and EMF show premia that seem to be persistently
on a level above zero over time. For these funds we observe positive values
for at least 75 percent of the sample.
To finally characterize the sample data, we further analyze the time-
series properties of the empirical premia. According to the described liter-
ature in section 7.1 the empirical premia show some evidence to behave in
a mean-reverting manner. Therefore, we model the premia P REMt ad hoc
by an Ornstein-Uhlenbeck process with a specification of the form
Notes:
a) The premia are given in percent, measured by In(PtfNAVd. UQ and LQ denote
the upper and lower quartile, respectively. The median is abbreviated by MED.
b) We also report the skewness (SK) and kurtosis (KU) of the premia distribution.
8.1 Sample Data 59
Note:
Statistically significant parameter estimates at the 1%-, 5%- and lO%-levels
are denoted by ***, **, and *, respectively.
60 8. First Empirical Results
for the long-run mean exhibit positive as well as negative values related to
whether the closed-end fund trades mainly with a premium or a discount
over the sample period. The range is from plus 15 percent to discounts of 21
percent which is close to the descriptive results for the mean values of table
8.2. When comparing the estimation results with the three closed-end funds
shown in figures 7.1, 7.3 and 7.2, we can confirm the high mean-reversion
tendency of 5.9 for the GT Global Eastern Europe (GTF). We can also
infer the high volatility for the Indonesia Fund (IF) from the figure. The
long-run means for the closed-end funds IF and GTF are estimated in line
with values of plus 15 percent and minus 12 percent.
Thus, we actually reach a steady state except for the last expression
which remains from the first term of the function C (t). Therefore, we
constrain our model to the case where the long-run mean of the continuous
50There are very rare exceptions which show this feature and no closed-end funds
included in our sample.
8.3. State Space Form 61
premium is given by
(8.2)
which reduces our parameter space by one parameter. Using equation (8.2)
we obtain a functional mapping of the two state variables X t and 7rt to the
value of the observable asset Y (X, 7r, t; 1/JM ). The pricing relationship for
the implementation of the closed-end fund market prices according to the
valuation model is given by
(8.3)
where P (X, 7r, t; 1/JM) depends on the model parameters 1/JM = {O'x, /'i" 0'1r, p,
>'1r}. In further examining the valuation formula of equation (8.3) we are
finally able to show how the second factor we choose in the valuation model
is related to the empirical premium as defined in equation (7.1). By rewrit-
ing the valuation formula in terms of logarithms we obtain the connecting
formula
1
P REMt = --7rt
/'i,
+ const (8.4)
for the two premia, with a constant term denoted by const = f (1/J M ). From
this linear relationship between the empirical and the modeled continuous
premia we expect to observe positive dynamic premia on closed-end funds
with corresponding empirical discounts, and vice versa, in the empirical
analysis.
(8.5)
for a time interval of 6..t. In order to specify the corresponding state space
form to equations (7.2) and (7.3) with the exact discrete time equivalents,
i.e. the discrete approximation yields exactly the same distributions as in
the continuous case, we need to match the following moments. 52
Remark 8.3.1 (Moments of the State Variables) For the first factor
el,t = X t we use equation {7.2} which can be integrated over the time period
6..t = [t - 6..t, t] to yield the almost explicit solution
J
t
J
t
(8.7)
E ['1Jtl = [0,0]' ,
(8.8)
Combining the transition and the measurement equations, the full state
space model depends on the time homogeneous parameter vector tPSSM =
{K;, A7r , 0"7r' hy , /-l, o"x, hx, p}. These parameters are estimated in a statistical
inference using a recursion algorithm based on linear Kalman filtering as
described in chapter 4. With the specified state space form of equations
(8.5) and (8.8), we are able toperform a maximum likelihood estimation of
54The function Oi,j stands for the Kronecker symbol which takes value 1 if i = j, and
o else.
55The small errors et are introduced to capture imperfections in observed price quota-
tions.
56 Further, we assume independent error terms et and 'TJt.
64 8. First Empirical Results
the parameters
Note:
Statistically significant parameter estimates at the 1%-, 5%- and 10%-levels
are denoted by ***, **, and *, respectively.
66 8. First Empirical Results
Note:
Statistically significant parameter estimates at the 1%-, 5%- and lO%-levels are denoted
by ***, **, and *, respectively.
8.4 Closed-End Fund Analysis 67
Fund 1f~)
Name ea ) MEAN STD MAX UQ MED LQ MIN SK c) KU c)
AF 4.335 -0.20 8.60 19.30 6.83 -2.33 -7.17 -15.97 0.34 2.15
APB 1.624 -0.45 10.24 16.68 9.89 -1.77 -9.19 -20.96 0.Dl 1.67
BZF 7.489 4.15 8.13 20.71 12.91 2.50 -2.65 -10.40 0.33 1.86
BZL 9.468 6.51 14.58 37.56 21.50 5.53 -5.48 -26.77 -0.05 2.05
CH 29.341 26.13 23.28 70.94 43.26 29.07 12.91 -31.60 -0.48 2.42
CHN 9.233 0.54 14.87 30.41 12.41 -1.49 -11.51 -31.02 0.27 1.95
CLM 20.841 20.67 5.85 30.11 25.83 21.61 15.81 7.51 -0.33 1.91
MEF 5.114 3.56 12.98 23.45 15.96 2.71 -2.94 -29.03 -0.55 2.64
ETF 3.601 2.67 5.02 9.56 7.49 2.95 -1.12 -8.84 -0.47 2.03
ISL 14.874 8.26 16.71 36.87 20.94 14.03 -1.66 -30.75 -0.74 2.42
FPF 85.265 83.45 27.75 135.50 101.18 87.94 75.33 -26.02 -1.58 6.15
GCH 12.752 8.35 13.15 34.58 21.72 8.44 -3.31 -17.61 0.02 1.78
GTF 28.684 25.65 12.06 45.39 34.32 28.93 18.88 -7.04 -0.82 2.83
IGF -2.854 -14.10 16.09 25.88 -2.50 -12.60 -23.52 -68.83 -0.37 3.40
IF -33.910 -37.88 29.96 30.00 -28.14 -43.34 -60.28 -94.02 0.81 2.82
JGF -6.969 -6.49 18.08 29.98 5.82 -9.40 -20.19 -46.33 0.32 2.30
JFC 11.124 7.13 15.75 39.80 22.58 8.32 -6.70 -31.14 -0.10 1.85
KF -25.960 -15.21 11.84 21.17 -7.27 -13.09 -23.52 -44.24 -0.30 2.77
KIF -5.262 -4.84 19.24 29.31 10.20 1.02 -18.54 -65.34 -0.71 2.69
LAQ 0.779 7.72 10.18 27.64 16.38 7.68 1.25 -16.25 -0.32 2.26
LAM 1.522 7.64 10.73 27.50 16.69 8.53 0.47 -22.35 -0.55 2.68
LDF 2.164 8.40 5.45 19.55 11.82 8.91 4.76 -6.66 -0.31 2.82
MSF 0.029 3.37 7.69 21.03 8.86 2.72 -1.70 -16.41 0.12 2.51
MF -7.817 5.18 19.95 32.59 18.63 10.60 -2.25 -76.84 -1.81 6.79
MXE 2.364 2.65 11.25 20.36 14.41 1.73 -5.20 -26.62 -0.34 2.26
MXF 8.367 8.49 10.60 29.24 18.03 7.31 2.88 -19.34 -0.36 2.72
ROC 1.917 3.12 20.42 51.26 18.29 1.42 -12.36 -41.07 0.26 2.24
SAF 4.356 3.90 12.96 25.66 15.91 2.45 -5.56 -23.49 0.01 1.94
SGF -10.817 -12.15 33.02 56.99 13.36 -14.12 -34.29 -113.76 -0.22 3.01
TWN 1.319 -2.42 13.14 25.98 8.09 -5.24 -12.22 -35.95 0.18 2.31
EMF -58.447 -59.01 26.67 -8.21 -35.08 -58.70 -77.13 -139.25 -0.30 2.36
TC -3.932 -0.42 12.90 16.08 6.96 3.68 -1.98 -42.19 -1.58 4.81
TTF -1.375 -0.17 14.96 22.96 9.15 5.37 -4.14 -49.40 -1.43 4.38
MEAN 3.310 2.85 14.97 32.00 14.44 3.37 -6.74 -37.21 -0.34 2.75
STD 21.975 21.38 6.97 23.24 21.33 22.58 23.95 31.53 0.60 1.17
MAX 85.265 83.45 33.02 135.50 101.18 87.94 75.33 7.51 0.81 6.79
MIN -58.447 -59.01 5.02 -8.21 -35.08 -58.70 -77.13 -139.25 -1.81 1.67
Notes:
a) The long-run means of the dynamic premia (in percent) are given by equation (8.2).
b) In the statistics for the dynamic premia (in percent), UQ and LQ denote the upper
and lower quartile respectively; the median is abbreviated by MED.
c) For each closed-end fund we also report the skewness (SK) and kurtosis (KU) of the
dynamic premium distribution.
68 8. First Empirical Results
the net asset value volatility ax closely resemble the values obtained empir-
ically from the data (compare table 8.1). The differences take a maximum
of 4.6 percentage points with an average of 1.1%. Further, the estimated
model error h y of 0.005 across all funds translates to average pricing devia-
tions of 0.54% with the observable prices on the net asset values. 58 Lastly,
the values for the correlation coefficient p primarily show positive values up
to +0.5 with a mean of +0.13, i.e. changes in the net asset value and the
dynamic premium are related.
160%
40%
-40%
-80%
-120%
58S ee, for example, Aitchison and Brown (1957, ch. 2.2) for the log-normal transfor-
mation.
8.4 Closed-End Fund Analysis 69
40%
20%
.~
E 0%
e!
a. .~,~yv~~~
-20%
-40% ROC Taiwan Fund (ROC) I ',' ••••••' / .....,'.,
80%
60%
ttl
'E 40%
e! 20%
a.
0%
-20%
10%
'E -10%
ttl
0%
;.;;;t'''f.
e! -20%
a. -30% .. / ... /./' :~
:····:······/":.·· ....t···· .. ·. . . . . :..... ··..... r·:· ~" . :'...:. . :\. . ::..":.;. .: . ./
-40% GT Global Eastern Europe (GTF) '.' .......... :.. . .
-50%
o 52 104 156 208 260
- - empirical premium ....... (inverted) dynamic premium
deviations of the AF, APB, TC, and TTF funds. However, the sign needs
to be put into perspective when incorporating the reported standard devia-
tions. Looking at the realized dynamic premia, we find high fluctuations as
indicated by the estimated values for the speed of mean-reversion and the
volatility from above. With 135.5 percent FPF marks the maximum and
EMF shows the minimum of minus 139.3 percent.
In figure 8.1 we graph the evolution of the dynamic premia with inverted
sign over the sample period, i.e. we visualize the numerical results of table
8.6. The figure contains the maximum and minimum values across all funds,
along with the mean values and the 90 percent interval. In general we see a
higher fluctuation in the positive premia from the maximum values as well
as the confidence interval. The negative premia seem to be estimated in a
much closer range. Overall, the visual impression resembles the results for
the empirical premia as shown in figure 7.4 on page 47 except for a scaling
factor. This fact also holds when we look at the three selected funds from
section 7.1 which we graph in figure 8.2. In each of the graphs we plot the
time-series of the empirical as well as the (inverted) dynamic premia for the
three closed-end funds ROC, IF, and GTF. We see that for the funds IF
and GTF the inverted dynamic premia emphasizes the positive and negative
values.
9 Implications for Investment
Strategies
setting. In order to get current values for the state variables at initiation of
the forecast we filter the data based on the estimated four year parameters
for the period of the first week in the sample to the X-th week.
The finite sample forecasts are obtained from the same state space form
as used with the Kalman filter. The expected values for the s-period-ahead
predictions of the state variables and the observable prices are the condi-
tional expectations
yt - m = a (yt-~t - m) + Et + bEt_~t
with b..t = 1/52. This process specification of the closed-end funds market
prices translates with the matrix specifications
Notes:
a) We first report the sum of squared errors (SSE).
b) The two measures to evaluate the goodness of fit are the mean absolute deviation
(MAD) and the root mean square error (RMSE).
c) We further show the results on the maximum and the minimum of the prediction
errors (RESMAX and RES MIN ).
d) The differences are calculated as the result from our model minus the result from
the ARMA(l,l) model, i.e. negative values signal better performance than the reference
model.
9.2. Implementing Trading Rules 75
values of the statistics. We find the realized maximum MAD to lie between
$1.13 and $2.37 which is throughout lower than the results for the reference
model. The highest individual prediction errors take values from $4.26 up
to $5.32 for one month. The corresponding minimum realizations of the
forecasting errors are found to range from -$5.08 to -$12.04. This indi-
cates that the pricing model tends to overestimate the forthcoming price
movements and provides evidence that the forecasting power is better in
times of rising market prices. Negative price jumps are only captured with
higher forecast errors which, however, are lower than those predicted by the
reference model.
The different aspects covered by the implemented forecasting study re-
veal a good prediction quality of the pricing model in terms of the closed-end
fund market prices in our sample. This motivates us to examine the second
application on how the information behind the pricing model could be used
to create abnormal return generating trading portfolios.
Several studies show that abnormal returns can be earned by following trad-
ing strategies based on the information content revealed in the empirical
premia. 62 In the second application of the derived pricing model we follow
the idea of exploiting the information content using the modeled dynamic
premia. For the implemented trading rule study we split the available five
year sample into a four year parameter estimation period and a one year
out of sample implementation period similar to the setup of the forecasting
study.
The filter strategies are based on the information content of the residual
premium difference defined by RESt,i = 7rt,i - (}i for a closed-end fund ion
day t. For this residual value the theoretical mean-reversion property of the
62For the US market of closed-end funds see, for example, Thompson (1978), Anderson
(1986), and Pontiff (1995).
76 9. Implications for Investment Strategies
ARt = ;;:
t i=l
+1 if filter strategy indicates buying,
with lREst,i = { -1 if filter strategy indicates shortselling, and
o otherwise
across the number nt of selected funds contained in the trading portfolios
on week t. The average abnormal returns are then cumulated over the fifth
year and the reported t-tests on a given trading strategy are computed as
in Campbell, Lo, and MacKinlay (1997).
Analogously to the described trading strategies based on the information
content of the dynamic premia we implement the same trading rules using
the empirical values for the premia as defined in equation (7.1). For this
benchmark test we use the residual of the empirical premium RESt,i =
PREMt,i - PREM i , where PREM i serves us as proxy for the long-run
mean value of the empirical premia and is calculated as the sample mean
over the first four years. We then apply the selection criterion of equation
(9.1) with opposite signs, since empirical discounts correspond to positive
dynamic premia. The final procedure to calculate cumulated abnormal
returns is done as described above.
Notes:
For each type of strategy we first report the average number of CEFs (Obs.) contained
in the trading portfolios and the number of trades (nT) needed to realize the trading
strategy. Statistically significant cumulated abnormal returns (CARs) at the 1%-, 5%-
and lO%-levels are denoted by ***, **, and *, respectively.
9.2 Implementing Trading Rules 79
Notes:
For each type of strategy we first report the average number of CEFs (Obs.) contained
in the trading portfolios and the number of trades (nT) needed to realize the trading
strategy. Statistically significant cumulated abnormal returns (CARs) at the 1%-, 5%-
and lO%-levels are denoted by ***, **, and *, respectively.
80 9. Implications for Investment Strategies
return of 21.8 percent which marginally falls to 20.0 percent for an imple-
mentation lag of one week. Second, from the average number of observations
(Obs.) for a given portfolio we see that the long strategies contain two to
six times more closed-end funds in their trading portfolios on average than
do the short strategies. This indicates that the dynamic premia lie above
their long-run equilibrium () in most cases of the fifth sample year which cor-
responds to observable closed-end fund discounts of the empirical premia.
Third, we notice in general that the short strategies earn higher cumu-
lated abnormal returns than do the buying strategies. With an immediate
implementation of a long strategy the cumulated abnormal return results
start from 13.6% and then gradually fall to negative values for higher filter
values. However, the cumulated abnormal returns of around 13% are re-
mained for filter values up to 0.1 which signal that the optimal filter height
for a trading strategy lies between zero and ten percent. The comparable
short strategy selects a daily average of 5.7 funds and yields a cumulated
abnormal return of 58.3% with no filter implemented. Here, the highest
cumulated abnormal return with 65.2% is realized for an intermediate fil-
ter height of 0.125. Comparing the overall results for the immediate and
the lagged implementation of the trading strategies show that the main re-
sults are preserved and the information content persists over one week. The
realizable values for the cumulated abnormal returns, however, are not as
high any more. Fourth, in implementing the portfolio strategies with the
aim of earning high cumulated abnormal returns one needs to consider the
transactions costs associated with these strategies. The reference number
of trades (nT) is 66 which stems from building a buy-and-hold portfolio
containing all closed-end funds and its liquidation at the end of the year. 65
In order to administer the long and the combined strategies the number of
trades lies above 66 for attractive cumulated abnormal returns. The highly
profitable short strategies can be followed with twenty to forty transactions.
Compared to the short strategies we further find the long portfolio holdings
more stable over time in that the relation of the number of trades to the
included closed-end funds is generally lower. With the objective to lower
transactions costs it seems optimal to pursue an intermediate filter strat-
egy of around 0.1 which only marginally sacrifices the returns. Finally, in
order to evaluate the outcomes of the trading test we compare the model
based results to those obtained from the information content of the empir-
ical premia which we report in table 9.3. The facts already described for
the model based results principally hold for the benchmark results as well.
When comparing the realizable cumulated abnormal returns our model per-
forms better for filter values close to zero in terms of long and combined
strategies. Overall, we see the similar pattern that the cumulated abnor-
mal returns remain stable up to certain filter values while the number of
closed-end funds and the number of transactions decrease.
10 Summary and Conclusions
11.1 Overview
In part III of the study we explore a new approach in modeling the term
structure of interest rates and the pricing of fixed income instruments. The
pricing model we present is a type within the affine class and is estimated
using interest rate panel data. This first chapter considers the peculiarities
of modeling and pricing assets in an incomplete market such as the market
for fixed income instruments. We clear the basic definitional terminology
when dealing with bond prices and interest rates and motivate our choice of
model type. In chapter 12 we present our dynamic term structure model and
derive a closed-form solution for the price of discount bonds as the security
of primary interest. In the next chapter we perform a comparative statistic
in order to evaluate the characteristic theoretical properties of our interest
rate model. Here, we especially examine the distinct behavior of the term
structure of interest rates as well as the term structure of volatilities which
are relevant in valuing derivatives. The following chapter 14 is dedicated
to the management of interest rate risk. Starting with a clarification of the
types of risk involved, we show how to implement the duration technique for
our model, and finally show how to price term structure derivatives. Among
the most relevant fixed income instruments in risk management we special-
ize in valuing bond options, swap contracts, and interest rate floor and
cap agreements. Chapter 15 considers the calibration of the proposed term
structure model to standard fixed income instruments. First, we describe
the different econometric methods used in the literature to then specialize
our implemented version that resembles the theoretical model properties at
best. For the different interest rate markets we choose data from US govern-
ment issues, LIBOR rates and swap rates. Upon these data we report our
88 11. Introduction
inferences on the model parameters and analyze the filtered state variables.
Finally, we conclude in the last chapter with a summary of our findings.
Before we describe our proposed term structure model in the next chap-
ter, we first take a look at the incomplete market for fixed income instru-
ments. There, we unfortunately realize that interest rates are not directly
tradeable. Instead different financial instruments are traded on the markets
which depend on interest rates themselves. Guided by this fact, we initially
consider the basic financial instruments in the bond market and how they
are related to each other. Second, we closer examine the peculiarity of the
risk premium concept that generally need to be taken into account in inter-
est rate modeling. Last but not least, we theoretically motivate our choice
of a special interest rate model with a risk premium following stochastic
dynamics.
-P- = lim
t (r)n
1+- ;
Pt - 1 n-+oo n
this result can be stated in the form of the Taylor series expansion
. (
1~m
noon
r)
1+ -
n
= 1 + -1'1. 1 2 1 3
r +.
2' r +
.
3' r + 0
( 3)
r = er .
66S ee, for example, Campbell, Lo, and MacKinlay (1997).
11.2 Bond Prices and Interest Rates 89
where c (u) denotes the cash flows (coupon and notional) received from the
coupon bond holder at time U. 67
The functional mapping T 1----+ Y (t, T), for a given time t, is called the yield
curve. The spot rate y (t, T) collapses at the short end of the term structure
to
which equivalently is the instantaneous short and forward rate. The general
forward rate for [T1' T2J contracted at t is defined as
oR (t T)
f(t,T)=y(t,T)+(T-t) cJ,
which shows that the forward curve superficially resembles the spot curve
dependent on the slope of the spot curve. The forward curve is higher
11.2 Bond Prices and Interest Rates 91
than the spot curve if the yield curve is increasing, and vice versa. More-
over, using the definitions on interest rates we gain the ordinary but crucial
relationships
P(t, u) = e-(u-t).r(t,u), and
pc (t,T) = L CF(u)· e-(u-t).r(t,u)
uE(tjT]
for discount and coupon bond prices, respectively. By restating the price of
a discount bond as
Jf
T
B(O,t)=exp jr(()d(,
o
i.e. the account starts with an initial value of B (0,0) = 1 and accumulates
interest with the dynamics dB (0, t) = r (t) B (0, t) dt over time.
Finally, to value financial instruments such as interest rate swaps, cap
or floor agreements we introduce a special interest rate, the LIBOR rate.
This interest rate uses simple instead of continuous compounding as we see
from its definition.
Definition 11.2.4 (LmOR Rate) The simple spot rate or LIBOR rate
is defined as
P(t,T) -1
L (t, T) = - (T _ t) P (t, T) (11.3)
The relationship between interest rates and forward rates shows that
interest rates of different maturities are heavily correlated. 72 Therefore and
to keep an interest rate model arbitrage-free and analytically tractable, the
stochastic processes governing r (t) or f (t, T) are assumed to be functions
of a fixed number of Brownian motions.
The considered relationships between bond prices and different types
of interest rates along with the fact that there are several interest rate
dependent securities traded on the financial markets and interest rates on
the other hand are non-tradable call for a special treatment as we will discuss
in the next section.
attributed with the level of interest rates. In fact, this is the classical ap-
proach to interest rate theory adopted in the field of financial economics
with the most prominent representatives being the arbitrage-free Vasicek
(1977) model and the utility based methodology of Cox, Ingersoll, and Ross
(1985b). The models make assumptions about the stochastic evolution of
the short rate in a continuous time model. The specification of such dynam-
ics are chosen with the aim of incorporating the fluctuations of the short
rate during time at best. The peculiarity about these types of models is
that the short rate as chosen factor is a non-traded underlying object in the
interest rate market. Therefore, we deal with the bond market as an incom-
plete market in comparison to, for example, the well known stock market
model dating back to Black and Scholes (1973). Since the bond market
is incomplete by assumption, we will further see that it is not possible to
hedge a generic derivative contract as well as we will not be able to derive
a unique price on such a contract. We begin with the assumption of the
relevant model for the bond market.
Assumption 11.3.1 (Bond Market) For the state variable in our model
we assume the a priori specified belief of the short-rote dynamics
In order to find the arbitrage-free price of the discount bonds, we use the
idea of an investor seeking for arbitrage opportunities. He has the possibility
to hold three different assets: The money market account and two discount
bonds P (t,~) of different maturities i E [1,2]. The discount bonds follow
94 11. Introduction
'T/i =
1 [0- aP (t, Ii)]
P (t, Ii) or
found from the above assumption using Ito's formula. The portfolio struc-
ture of the arbitraguer is given by the sum of the investments in the three
tradeable assets in the market:
dV (t) = [no (t) r (t) B (0, t) + nl (t) /-LIP (t, T I ) + n2 (t) /-L2P (t, T2)] dt
+ [nl (t) 'T/IP (t, T I ) + n2 (t) 'T/2P (t, T2)] dW.
i. e. the arbitraguer should not be able to realize any profits on his portfolio.
Using the abbreviations Wo (t) = no (t) B (0, t) and Wi (t) = ni (t) P (t, Ii)
for the relative portfolio weights the conditions of no-arbitrage reduce to
1 1 1] [wo (t) ]
[ o "11 "'2 W1 (t) = o. (11.6)
r (t) ILl IL2 W2 (t)
Therein, solutions for the portfolio weights only exist if the market satisfies
the internal consistency relation, i.e. the matrix in equation (11.6) needs to
be singular. With one degree of freedom we choose, for example, the last
row being a linear combination of the remaining rows. Using>.. for the free
coefficient not determined in the model for the bond market yields
These equations reduce to the internal consistency result for the bond mar-
ket:
\ ILl - r IL2 - r
A=--=--. (11.7)
"'1 "'2
More specifically, the quantity>.. is called the market price of interest rate
risk, since it gives the extra increase in expected instantaneous rate of return
74 For a valuation of derivative contracts based on no-arbitrage see, for example, Schobel
J.Li - r on a bond per an additional unit of risk 'fJi' The result of equation
(11.7) now states that all bonds in an arbitrage-free market have the same
market price of risk regardless of the time to maturity. This result entails
even more information on the bond dynamics. We substitute the drift and
volatility coefficient of equation (11.4) into one of equations (11.7) to obtain
the so-called term structure equation
~
2a
2 1j2 P
ar2
(t, Ii)
+
(-.A)
J.L a
ap ar
(t, Ii) ap (t, Ii)
+ at
- P ( 1',.)
- r t,..
This equation governs the price of discount bonds in an arbitrage-free mar-
ket when we additionally incorporate the appropriate boundary condition
P (Ii, Ii) = 1 for discount bonds. We also note that the price of a particular
bond P (t, Ii) is not uniquely given by the bond market model in that the
market price of risk .A remains undetermined in the model. It constitutes
an exogenous parameter. The reason for this fact is that arbitrage pricing
is always a case of pricing a derivative (in our case P (t, Ii)) in terms of
some underlying assets (in our case r). But, with the bond market choosing
a value for .A to satisfy the internal consistency relation of equation (11.7)
arbitrage possibilities are avoided. At this point we see the possibility to
establish a new dimension in term structure modeling in that we build a
term structure model that captures not a time homogeneous but a time
varying market price of risk .A (t).
12 Term Structure Model
criteria that need to be met by the model. For Rogers (1995) a term struc-
ture model should be:
A class of term structure models that broadly meet these criteria and
has recently been attended considerable focus is the affine class. 77 This
class is defined by factor models in which the drift and volatility coefficients
of the state variable processes are affine functions of the underlying state
vector. 78 Our approach is to work within this class of factor models and
extend the single-factor term structure models, such as the model used in
the previous section, to include the market price of risk A (t) as a second
state variable. Thereby, we establish a two-factor model that enables us to
analyze the investors' beliefs with respect to the interest rate risk inherent
in the observable term structures on the markets. We aim at building a two-
factor model that accommodates the above mentioned criteria. Especially,
as we will see in the following section, we try to build a maximal model that
does not suffer from implicitly imposing potentially strong over-identifying
restrictions. 79
Looking at the extant financial literature on this subject we already find
several propositions of two-factor model specifications. The main idea is
that single-factor models based on specifying a dynamic structure of the
short interest rate exhibit some shortcomings that are rather unrealistic. 8o
The most problematic feature of such models is that they endogenously
result in a perfect correlation among bond prices of different maturities.
76Rogers (1995, p. 93).
77For an overview of different classes in term structure modeling see, for example,
Uhrig (1996, ch. 1).
78See, for example, Duffie and Kan (1996).
79See, for example, the specification analysis of Dai and Singleton (2000).
80See, for example, Chen (1996b, ch. 2).
12.1 Motivation for a Stochastic Risk Premium 99
Several studies have shown that such single-factor models are not suitable
in modeling the term structure of interest rates. 8l
81The empirical studies on this subject include Chen and Scott (1992) and Pearson
and Sun (1994) who use a maximum likelihood estimation, the works by Heston (1992)
and Gibbons and Ramaswamy (1993) based on the generalized methods of moments,
and the factor analysis of Litterman and Scheinkman (1991).
100 12. Term Structure Model
82Chen (1996a), for example, develops a model using the short rate, the short-term
mean, and the short-term volatility as explaining factors.
12.2 Economic Model 101
Under the drift transformations the mean-reversion property and the corre-
lation d~dW~t = pdt of the two processes are retained.
Remark 12.2.3 (Pull System Matrix) Note that the system of stochas-
tic differential equations (12.2) incorporates a full system matrix TQ, i. e. we
model a two-fold interdependence: The two modeled state variables r (t) and
..\ (t) may both depend on their own value as well as on the other factor's
value. Furthermore, the state variables are related through the correlation
between the Brownian increments.
(12.4)
881.e. the stochastic part in the models corresponding stochastic differential equation
drops out. See, for example, Merton (1973).
89Note that the vice versa argument does not hold; freedom of arbitrage does not imply
an equilibrium.
13 Initial Characteristic
Results
P, .-
r>. . -
82 P(t,T,x) p,.=
8r8>. -- ABP, t .
8P(t,T,x)
8t
= (-A t r - B t ).. - C)
t,
P
Looking closer at the system, we fortunately see that function C only ap-
pears in the last equation (13.5) and only in its first derivative with respect
to t. Thus, we can start solving the system by examining the first two equa-
tions (13.3) and (13.4) in order to derive a solution for the functions A (t)
and B (t). In matrix notation we have
with the system matrix [TQ]' and the driving term d = [-1 0]'. The gen-
eral solution to this non-homogeneous subsystem is the sum of a particular
solution to the non-homogeneous equation (denoted by AP, and BP) and the
general solution to the homogeneous equation (complementary functions AC,
and BC):
A (t) AP + A C , and
B (t) BP + B C •
For the particular solution we can gain the constant solution by solving the
following system:
I ~1 ~A
(h
I I ~T
-(h = -aT
AP= = , and BP -aT
~A I ~T(h + ~AaT
;~ I
~T ~T
I -aT BA I -aT
for the particular solution. Now, we look at the general solutions AC and
BC to the homogeneous equations
(13.6)
Putting these with the corresponding time derivatives into equation (13.6),
we get
and
+ A) . a + K),. . /3
(Kr o (13.7)
-ar . a + (0),. + A) . /3 O. (13.8)
We note that equations (13.7) and (13.8) only have non-trivial solutions, if
the characteristic polynomial is equal to zero, i.e.
I Kr-a+ A K),.
0),. +A
I= 0'
r
or
(13.9)
This results in Ai = -!
(Kr + (h) ± h/(Kr + 0),.)2 - 4 (Krfh + K),.ar ) for the
eigenvalues, with i E {1, 2}.
Now, we can resolve the values for a and /3 from, for example, equation
(13.8) to aj/3 = (0),. + A) jar. Choosing arbitrary constants ki' with i E
{1,2}, we get ai = (0),. + Ai) ki' and /3i = arki . Thus we get for the com-
plementary functions
For an intermediate result, we now combine the particular and the comple-
mentary solution for the functions A (t) and B (t). Making use of Vieta's
Theorem, by which the product of the eigenvalues Ai, with i E {I, 2}, equals
the constant term "'rfh. + "'>/l'r in equation (13.9)91, we obtain:
A (t) (13.10)
B (t) (13.11)
In order to derive the constants k1 and k2' we now make use of the first two
boundary conditions A (T) = 0 and B (T) = O. For equations (13.10) and
(13.11) we get in matrix form:
...!!..l......
A1A2
(h +A2
1 1
A1A2 -1
k1e A1T
(h. + Al 0>. +A2
= Al (Al - A2 )
1 1
-1
k1 = A1 (Al - A2) e-A1T , and
1
(h +A1 (h +A2
1 1
1 -A2T
k2 = A2 (A1 _ A2) e .
This finally results in the values for the functions A (t, T) and B (t, T):92
Having derived a solution for A (t), we then resolve equation (13.4) for the
function B (t). With the boundary conditions A(T) = 0 and B(T) = 0, we
finally get the same results for the functions A (t), and B (t).93
122 122 )
Ct = 20'rA - JLrA + 20'),.B - (JL),. - 0'),.A2 B + PO'rO'),.AB,
where we now enter our results of equations (13.12) and (13.13) for A (t)
and B (t). A straightforward calculation gives us in a nice rearrangement
92 For a better understanding of the independent variables of the functions A (t) and
B (t) we add the parameter T to the final expressions, i.e. denote the functions as A (t, T)
and B(t,T).
93This alternative derivation is of interest in terms of studying the equilibrium behav-
ior of the system, because it nicely relates to the nonhomogeneous partial differential
equation which is known from physics to model oscillation. However, our approach is
more symmetric in its formulation and derivation of the solutions to A (t) and B (t).
13.1 Valuing Discount Bonds 111
with eo '122
2 urao - J.trao + '2122
u ),bo - (J.t), - U),A2 ) bo + paru),ao b0
CI = (u~ao - J.tr + paru),bo) al + (u~bo - (J.t), - U),A2) + parU),ao) bl
C2 = (u~ao - J.tr + puru),bo) a2 + (u~bo - (J.t), - U),A2) + puru),ao) ~
1 22 1 22
C3 = '2 Ural + paru),a1bl + '2 U),bl
1 22 1 22
C4 = '2 Ura2 + purU),a2 b2 + '2 U),b2
Cs U~ala2 + parU)' (al~ + a2bl) + U~bl~'
Finally, in view of the boundary condition C (T) = 0, we integrate
equation (13.14) to derive the result for the function C (t)
-C (t) = Jeo
T
+ cle-A1(T-() + C2e-A2(T-() + C3 e- 2A1 (T-()
t
+C4 e - 2A2 (T-() + cse-(A1+A2)(T-()d(,
which yields: 94
C (t, T) = -eo(T - t) + ~ (e- A1 (T-t) - 1) + ~ (e- A2 (T-t) - 1)
Al A2
+~ (e- 2A1 (T-t) _ 1) + ~ (e- 2A2 (T-t) - 1)
2AI 2A2
+ Cs (e-(Al +A2)(T-t) _ 1) , (13.15)
Al +A2
where the constants Ci, with i E {O, 1, ... , 5}, are given in equation (13.14).
Solution 13.1.4 (Discount Bond Price) Putting our results together, we
are able to state the price of a discount bond in the closed-form expression
P (t, T, x) = e-[A(t,T),B(t,T)]x-C(t,T) (13.16)
where we use the formulas for the functions A (t, T), B (t, T) and G (t, T)
from equations (13.12), (13.13), and (13.15), respectively.
94We also add the parameter T to the final expression on the function C (t) for a better
understanding, i.e. we denote the function as C (t, T).
112 13. Initial Characteristic Results
Moreover, with the explicit formulas of the functions A (t, T), B (t, T)
and C (t, T) we can readily calculate the risk neutral dynamics of the dis-
count bond price. We begin with the P-dynamics under the martingale
measure Q of equation (12.3) which generates the stochastic form:
dP
[( -At - ~rA - ~).B) r + (-Bt - (hB + O'rA) A
P
122 - j.£rA + 1
+20'rA 20'),.B2- 2 ()
j.£),. - 0'),.A2 B + PO'rO'),.AB - Ct] dt
-O'rAdWZ - O'),.BdW~t·
In conjunction with the results for A (t, T), B (t, T), and C (t, T) the stochas-
tic dynamics for the discount bond simplifies to
Equation (13.17) shows explicitly that the modeled discount bond P (t, T, x)
actually has a drift rate of r in the risk neutralized world.
1
Y (t, T) = - T [[A (t, T), B (t, T)] x + C (t, T)]; (13.18)
-t
the factor loadings A (t, T) and B (t, T) are given in equations (13.12) and
(13.13), and the intercept term C (t, T) is from equation (13.15). Especially
this affine structure of the yields in the driving factors - the short rate r
and the market price of risk A(t) - will be of special importance for the
empirical implementation and analysis which we perform in chapter 15.
Further, it is of interest how the forward rate curve looks in this model.
Starting from equation (11.2) we can see that forward rates are related to
spot rates by
Thus, the term structure of forward rates can immediately be derived from
the zero bond yield curve given in equation (13.18).
dy (t, T)
with Jl. y
114 13. Initial Characteristic Results
Therein, the partial derivatives Yr, Y).., Yt, Yrr, Yr).., and Y)..).. are easily found
from the result of equation (13.18). Using these derivatives leads us to
T)
Y (t + .6.t, T + .6.t) Y (t, + (A(t,T) B(t,T)).6.
ILl T _ t + IL2 T _ t t
+ J
Hll.t
A (t,T)dW Q
O"r T _ t r,(+
J
Hll.t
where the two stochastic integrals follow a normal distribution. Thus, the
spot rate Y (t, T) follows a bivariate Gaussian distribution. From this dy-
namics we obtain the variance per unit of time, i.e . .6.t = 1, as:
different maturities:
[ O"~A (t, T1) A (t, T2) + O"~B (t, T1) B (t, T2)] dt
(T1 - t) (T2 - t)
+ [A (t, T1) B (t, T2) + A (t, T2) B (t, T1)]
(T1 - t) (T2 - t)
to reveal our final result for the correlation coefficient of the spot rates with
maturities T1 - t and T2 - t:
(13.23)
J
T T
J ~dAc
Ac = -Kr d(
t t
-lnAC(t) -Kr (t - T) +k
A (t)
C k' eK.r(T-t) .
°
Now, we add the particular and complementary solution, and incorporate
the boundary condition A (T) = to find the general solution to be
1
A (t) = -- (1 - eK.r(T-t)) • (13.24)
Kr
°
For the function B (t), we get from the second homogeneous equation in
(13.23) by using B (T) = the following result:
- CT r
Kr
J(1 -
T
eK.r(T-()) d(
t
-B (t) = _ CTr
Kr
((T _t) + ~ (1 _ Kr
eK.r(T-t)))
In order to derive the function C (t) with its boundary condition C (T) = 0,
we use the last equation of equations (13.23):
-C(t)
C (t)
118 13. Initial Characteristic Results
To see the equivalence between the bond price formula from Vasicek
(1977) and our solution
dr = a b - r) dt + dz.
With the substitutions -Kr = a and -J-lr/Kr = 'Y we are easily able to
transfer our results into that found by Vasicek (1977).
lim B (t) =
TltT - t
+~: (I-A2(T-t)+o(T-t)-I)
+ 2~l (1- AdT - t) + o(T - t) -1)
+ 2~2 (1 - 2A2 (T - t) + 0 (T - t) - 1)
By having calculated these limiting values for the functions A (t), B (t), and
C (t) we indirectly verified our pricing formula of equation (13.16) in that
we obtain
Ys (t)
T!t - tT 1 - [[A (t), B (t)] x + C (t)] = r (t)
= limy),t, T) = l i m
T!t - t
for the short end of the term structure of interest rates. Thus the short end
coincides with the first state variable of our term structure model.
Now, we analyze the asymptotic behavior of the term structure of in-
terest rates at its long end. For the long end we need to investigate the
limiting value
Again, we examine the three functions A (t), B (t) , and C (t) separately:
Thus, we obtain the constant for the long end of the term structure
Yl (t)
for the stochastic processes for the state variables including a correlation
coefficient of p = 0.7102 between the two Brownian motions. For the cur-
rent values of the state variables we use rt = 0.0699 and At = -0.2624,
respectively, which reflect the empirical mean values of the realized state
variables. 97
How the term structure of spot interest rates reacts to different levels
of the state variables is shown in figure 13.1. The left graph is based on
a market price of risk of At = -0.2624 and models the impact of short
rates in the range rt = {0.02, 0.04, ... ,0.18}. The different values for the
short rate can be traced back to the crossing points of the yield curves
with the ordinate as demonstrated in section 13.3.2. All yield curves are
a little upward sloping since the market price of risk (At = -0.2624) lies
below its risk neutral value of At = O. The changes in the short rate have
the impact of an approximately parallel shift of the yield curve. The effect is
slightly higher on the short end of the yield curve than for longer maturities.
Therefore, the first state variable can be attributed as a level factor.
95See equation (13.21).
96The empirical parameter estimates will be considered in depth in chapter 15 on the
calibration of the model to US interest rate data.
97The empirical mean values are those obtained for the US Treasury data set shown
in table 15.7 on page 168.
122 13. Initial Characteristic Results
0.11
0.18 r- values: A. - values: ·1.0
0.16 0.10
0.14 -0.5
0.14 0.09
0.12
0.10
0.10
0.07
0.08 0.06
0.06 0.06
0.02
0.04
0.05
0.02
0.04
0 2 4 6 8 10 0 2 4 6 8 10
from time to time with a close fit to the then prevailing term structures
without reestimating the whole range of model parameters 1/J.
The next subject of analysis is the term structure of volatilities. Since
the volatility curve is independent of the levels of the state variables, we
now move to the influences of choosing particular model parameters on the
term structures.
--0.1 0.08
.(l.21 .(l.31 -0.4
0.07 -.(l.5
0.068
2 4 6 8 10 0 2 4 6 8 10 0 2 4 6 8 10
yield curve throughout the different scenarios. The parameter /-Lr gives the
short rate higher presumed values in the future in thus it has an increasing
effect on the yield curve. The additive terms in the process of the market
price of risk are /-LA and A2. Numerically positive changes in /-LA result in
lower yields, but within the empirically estimated values for the parameter
we only observe neglect able effects. The opposite influence is attained by
A2 since it enters the stochastic process with an opposite sign. The values
for these parameters both make the expected value of the market price of
risk either rise or fall. This further gives a feedback on the short rate which
is lower under the measure Q. Thus, higher values for the market price of
risk lower the term structure of interest rates as we demonstrated in the
previous section.
The next parameters we analyze are the direct factor loadings "-r, "-A,
and 0>.. in the stochastic processes of equations (13.27) and (13.28). With
respect to the short rate higher values of "-r assume higher expectations of
the investors in the future in that it has an increasing effect on the yield
curve. Only in the case of high values of A we observe a contrary reaction.
13.4 Possible Shapes of the Term Structures 125
Ii>.. li I il I il il il
a>.. il
fh. il il 1liil I li li
A2 li
p li
Effects on Yield Curve: li: increasing il: decreasing ::::::l: neglect able
The market price of risk is directly influenced by the parameters Ii>.. and
0>...With positive changes in Ii>.. we generally lower the expectations in the
market price of risk which makes the investors require higher returns in the
future. This results in an increase of the yield curve. The impact of the
level of the market price on itself is modeled by 0>... Its influence on the yield
curve varies in dependence of the sign of At, i.e. the current risk attitude of
investors. In risk taking scenarios we attain a negative influence and with
risk averse investors we obtain rising yield curves. For risk neutral investors
with the current short rate being at its long-run level the short rates are
pushed up whereas the long rates are lowered.
Moreover, the state variables are modeled with a possible correlation
factor p. With high values for this parameter it is more likely to observe
high values of the short rate whenever the risk attitude is high. It only
influences the term structure of interest rates by the term C (t, T) according
to equation (13.18). Through all different scenarios of the state variables we
observe increasing yields with positive changes in the correlation coefficient.
Finally, we examine the influences of the volatility parameters a r and a>...
126 13. Initial Characteristic Results
0.Q15 0.015
0.010 K, - values:
---2
0.005 0.005 -1.5/-1/-0.5
--0
0.0000
2 4 6 8 100 .0000 2 4 6 8 100 .0000 2 4 6 8 10
0.030 0.030 0.030
0.025 0.025
0.020
0.015
The volatility of the short rate does not directly impact the expected value
of the short rate. However, with higher values for a r risk averse investors
need to be compensated with higher returns. Thus, the yield curve increases
with higher ar-values. Risk taking investors on the other side love to have
more volatility in that they are satisfied with lower returns. The volatility
a>. on the risk premium models the frequency of changing risk attitudes of
the investors through time. Thereby, high fluctuations result in lower levels
of the yield curve.
Now, we move on to consider the comparative analysis of the term struc-
ture of spot rate volatilities. The possible cases given by the standard pa-
rameter values for the possible cases of volatility structures are shown in
figure 13.3. From equation (13.21) we see that the exogenous parameters
that influence the volatility structure are given by the reduced number of
't/J = {KT) aT) K>., a>., (h, p}. Thus, aside from the state variables the level
parameters J.Lr' J.L>., and '\2 do not impact the volatility structure as it is
expected from finance theory.
The volatility structures show us how the short rate volatility a r is
13.4 Possible Shapes of the Term Structures 127
101 Compare the results presented on the empirical estimations in chapter 15.
14 Risk Management and
Derivatives Pricing
l02See, for example, Macaulay (1938) and Bierwag, Kaufman, and Toevs (1983).
l030ther recent concepts include the functional duration from Klaflky, Ma, and Nozari
(1992) and the partial duration by Waldman (1992).
130 14. Risk Management and Derivatives Pricing
type managers already use interest rate derivative securities such as swap
contracts, or cap and floor agreements. Lastly, the credit risk is the newest
risk exposure fixed income managers start to deal with.104
The risks captured by our proposed term structure model include the
market and the yield curve risk whereas it is not capable to incorporate
any risk associated with default risky securities. As we analyzed in section
13.4 the first model factor can be seen as level factor and thus attributed to
describe general market movements of interest rates. With the market price
of risk factor we additionally consider the steepness of the term structure of
interest rates and thereby capture the yield curve risk. The two factors can
successfully be applied to the concept of duration using factor durations to
manage the market and yield curve risk.
for a factor F.
With respect to the two factors in our term structure model we get the
factor durations
A (t, T), and D>. = B (t, T), where the functions A (t, T) and B (t, T) are
from equations (13.12) and (13.13).
They duration measures of equations (14.1) and (14.2) can directly be
derived from the stochastic dynamics of coupon bonds. Starting with the
infinitesimal time changes of a discount bond P (t, T, x) = exp[-y (t, T)
(T - t)] based on equation (12.1) using 1M's formula yields:
dP(t, T,x)
(14.3)
P(t,T, x)
with J.£p =
1 [1
P(t,T, x) 20'r
2{P P (t, T, x)
or2 +PO'rO'>.
{P P (t, T, x)
ora)'
1 202p(t,T,x) ( ) oP(t,T,x)
+20'>. 0),2 + J.£r + Krr Or
(
+ J.£>. + K>.r
) oP (t, T, x) oP (t, T,
A), + at
X)]
' an
d
1 [OP(t,T,X) OP(t,T,X)]'
"1 = P(t, T, x) O'r Or 0'>. A),
(14.4)
_
II (t,T,x) -E
Q
[
e
- {r()d(
T II(T,T,x) Ft 1
, (14.5)
I&Q [ exp - I
r (() d(1 FJ
14.3. Bond Options 133
by applying the right numeraire asset which makes the contingent claim a
martingale. In evaluating the prices of equity options one conveniently
works with the money market account B (t, T) as numeraire. 110 However,
in pricing term structure derivatives we further use the forward martingale
measure which we describe and use in depth further on. In the coming
sections we demonstrate how to value interest rate derivatives based on the
martingale approach.
(14.6)
110This results in discounted asset prices being martingales with respect to Q; see
Harrison and Pliska (1981) and Geman, Karoui, and Rochet (1995).
111 From now onwards we drop the notation for the x-dependence of discount bonds,
i.e. we denote their price by P (t, T).
134 14. Risk Management and Derivatives Pricing
Now, we use the fact that a call option is a piecewise linear instrument,
which gives us the possibility of using indicator functions from probability
theory within our specification in equation (14.6). This yields
lffiQ [ e
- J r(()d(
tT P (T, U) l{p(T,u)2:K} F t 1
- KlffiQ [ e
- J r(()d(
tT l{p(T,u)2:K} 1
Ft . (14.7)
However, we cannot split the terms in the first expectation since they are
not independent under Q. To further analyze the expectations, we use
the powerful idea of choosing the optimal numeraire in the following when
pricing and hedging a given contingent claim as described in Geman, Karoui,
and Rochet (1995).112 They change measure from the accumulator measure
in use to a new equivalent forward measure by introducing a discount bond
as a practical numeraire.
J
n
x(w)dQT(w) = J
n
x(w)(TdQ(w).
°
The function (T, often denoted by ~, is the Radon-Nikodym derivative of
Q 'l11ith respect to QT defined on F t for ~ t ~ T. The Radon-Nikodym
derivative ~ is given by the ratio of the numeraires N in use, i. e. (T =
NQT / NQ. For the old numeraire we use the money market account
NQ = B(O,t)
B (0, 0)
'l11ith B (0,0) = 1, whereas the numeraire for the measure QT is given by the
discount bond maturing at date T
NQT = P(t,T)
P(O, T)'
112 The technique of changing numeraire was developed by Geman (1989) for the general
case of stochastic interest rates. Jamshidian (1989) implicitly uses this technique in the
Gaussian interest rate framework of Vasicek (1977).
14.3 Bond Options 135
lffi
Q [d QU I ] _ B (0,Pt)(t,PU)(0, U) , and
dQ F t -
or equivalently
lffiQl [Ycll F,]
lffiQ2[YIF,j= t =rlffiQ1 [yr- l lF,]
t lffiQl [cll F t ] ':.t ':. t
Q
lffil =
U QU [ e - J
(t lffi t
r()d( P (T, U) (TU- 1 l{P(T,U)~K} F
t
1
= P(t,U) lffiQU [ -Ir()d(p(T U)
B (0, t) P (0, U) e ,
B(O,T)P(O,U)
P (T, U) l{P(T,U)~K} F t
I]
= P(t,U)P(O,U)lffiQU [ -lr()d(B(O,T)l F,]
P (0, U) e B (0, t) {P(T,U)~K} t
P(t,U)lffiQU [l{P(T,U)~K}IFt].
113See, for example, Dothan (1990, p. 288).
136 14. llisk Management and Derivatives Pricing
IQ
lE2 = (tT lEIQT [ e - f r()d( (TT - 1 l{p(T,U)2:K}
t
]
Ft
When we insert the results for the expectations lE? and lE~ back into equa-
tion (14.7), we get
Ct = I I
P (t, U) lElQu [l{p(T,U)2: K} F t ] - P(t, T) KlElQT [l{p(T,U)2: K} F t]
= P (t, U) QU (P (T, U) ;:::: K) - P (t, T) KQT (P (T, U) ;:::: K) .(14.8)
Using the explicit formula for discount bonds in equation (13.1) we calculate
the forward bond price ~T as
~T = e[A(t,U)-A(t,T)]r+[B(t,U)-B(t,T)].HC(t,U)-C(t,T).
14.3 Bond Options 137
Applying It6's formula based on the r (t) - and A(t) -dynamics of equation
(12.2) we get the dynamics of Y;T under the measure Ql:
O"y,rdW~U + O"y,,\dW~u
-"21 (O"Y,r
2
+ 2{XTy,rO"y,,\ + O"y,,\
2 )
dt. (14.11)
We now integrate equation (14.1,1) over the time period [t, Tl:
JO"y,rdW~U + JO"y,,\d~u
T T
In yj = In y;T +
t t
-"21 J T
2 + 2{XTY,rO"Y,'\ + O"y,,\d(.
O"Y,r 2
Therein, the stochastic integrals are Gaussian random variables with zero
mean. Since It6 integrals have an expectation of zero, the sum of the two
random variables possesses a zero mean. The variance of the two integrals
can be calculated with 1M's isometry.
138 14. Risk Management and Derivatives Pricing
Theorem 14.3.4 (The Ito Isometry) 114 For a function f (t, w) we have
J
T
J0";
T
QU J
(In Y,u - uy"d~ - J Uy"tiWf
1 ~ In ~)
t t
-~ u}d(
QU (J uy"dw,<'U + J uy"dJ¥1lU
2)
t t
1 1
::; In KY,.u + '2Vy .
114See, for example, Oksendal (1995).
14.3 Bond Options 139
Using the fact that the Ito integrals in this equation is Gaussian with
N (0, v~), we finally get
Comparing the result of equation (14.14) with equation (14.10), we note the
equivalence of the volatility terms just by a reversion of the sign, i.e. {7Z,T =
-(7y,).. and {7Z,T = -{7y,)... Further, we obtain
J{7z'TdW~u + J{7z,)..dW~u
T T
In Zj: In Zr +
t t
- 2"1 J T
2
(7 Z,T + 2p{7 Z,T{7 Z,).. + (7z,)..
2 d(.
Thereby, we result in the same variance for the forward bond price P (T, U)
as in the previous case of the relative price ~u, i.e. we get the identity
relationship
2
V = V 2z = v 2y .
140 14. Risk Management and Derivatives Pricing
Restating our results, we have shown that the value of a European call
option can be calculated within our term structure framework using the
closed-form expression
In order to derive the arbitrage price of the European Put option belonging
to the same option series, we conveniently use the put-call parity relation-
ship.
The pricing formulas of equations (14.15) and (14.16) based on our two-
factor term structure model are actually in line with the Gaussian option
pricing formula of Geman, Karoui, and Rochet (1995, p. 456). Having gone
through the calculations for the two fair option prices we need to remark
the advantage of having the price of discount bonds available in closed-form
expressions. This especially eases the work with changing numeraires.
n
SWPt = 8N'L)L(0,Tj- b Tj ) - k)P(t,Tj) ,
j=l
i. e. the difference of the present values of the floating and the fixed payments.
Since financial markets give quotes on swap rates, i.e. the fixed rates at
which financial institutions offer their clients interest rate swap contracts of
different maturities, we further derive a valuation formula for the fair spot
swap rate. In arbitrage-free markets the fair swap value is given by:
n
oN L (L (0, Tj- b Tj ) - k) P (t, Tj ) ~ 0,
j=1
i.e. a fair swap contract is worth zero today for both counterparties; Sub-
stituting the implicit forward rate L (0, Tj - 1 , Tj ) from equation (11.3) we
obtain
8N~(P(t,1j-d-P(t,Tj)
~ 8P(t,Tj)
-k)P( T)
t, J
°
°
n
L P (t, Tj-d - P (t, Tj ) - 8kP (t, Tj )
j=l
n n
LP(t,Tj-d -P(t,Tj ) k8LP(t,1j).
j=l j=l
k = 1 - P (t, Tn)
n (14.17)
8~P(t,Tj)
j=l
14.5. Interest Rate Caps and Floors 143
Definition 14.5.1 (Interest Rate Cap and Floor) An interest rate cap
(floor) with tenor 8 = [Tj- 1 , Tj] is a series of caplets (floorlets) which re-
semble European call (put) options that are settled in arrears at maturity
dates Tj , j = 1, ... ,n, on the LIBOR interest rate L (Tj- b Tj). The payoff
characteristic of a so-called caplet (floorlet) is defined, respectively, as
Now, we discount the payoff at time Tj to calculate the value of the caplet
at time t. We conveniently work with the substitution K'C = 1 + 8Kc:
NE
Q [ - Tjtr()d(
e t
(1
(7j-l, 7j) -
P
*)+ Kc P (Tj - 1, Tj ) F t
]
NEQ [e - Y,«)d( (1- KcP(T;-" T ))+ .r;] j
* Q [ - Tj.t r()d( ( 1 ) + ]
Kc NE e t K'C - P (7j-l, T j ) Ft. (14.18)
In P(t,Tj)Kc ±
P(t'Tj~)
.!v2
2 and v2 = J
t
T
(O"P(,Tj_J) - O"P(,Tj)) 2 d(.
14.5 Interest Rate Caps and Floors 145
119For an overview see, for example, Campbell, Lo, and MacKinlay (1997, ch. 11).
148 15. Calibration to Standard Instruments
tion, we need to obtain parameter estimates for the risk adjusted processes.
However, the parameters denoting the risk premia cannot be estimated sep-
arately from neither a pure time-series approach nor a pure cross-section
analysis only. ScMbel (1995a), for example, uses a two step procedure
where he first estimates parameters associated with the dynamic behavior
of the factors from time-series data and identifies the risk premium in a
second step from a cross-section sample of the term structure.
The described difficulties with estimating term structure models have
created a branch of econometric literature devoted to the analysis of panel
data. This combination of time-series and cross-section data is considered
very suitable for the interest rate market, since on every trading day a
whole spectrum of interest rates, i.e. a term structure, is being provided
from the markets of fixed income instruments. An estimation based on
interest rate panel data takes into account the whole information embedded
in the term structure and should therefore result in more efficient parameter
estimates. The approach we implement is based on a convenient state space
representation of the term structure model and treats the underlying state
variables correctly as unobservables. In conjunction with the customized
state space model we then exploit the econometric methods of Kalman
filtering and maximum likelihood estimation.
This estimation procedure is considered an appropriate method for both
extracting time-series of the unobservable state variables and estimating the
relevant model parameters. 123 The literature along this line can be traced
back to the work by Pennacchi (1991) who estimates an equilibrium inflation
term structure model with a number of four Treasury bills with different
maturities. Duan and Simonato (1995) and Ball and Torous (1996) provide
encouraging simulation results on Kalman filters working with single-factor
models of the Vasicek (1977) and Cox, Ingersoll, and Ross (1985b) type.
Especially, the mean-reversion parameters are estimated accurately when
compared to other estimation procedures. 124 The work by Lund (1997)
uses extended Kalman filters in order to use other marketable securities than
pure yield curve data like, for example, coupon bond prices and swap rates.
123For further advantages see, for example, Geyer and Pichler (1997) and De Jong
(1997).
124For an exploration of different estimation methods in the case of mean-reverting
interest rate dynamics that are sufficiently close to a non-stationary process with a unit
root see, for example, Ball and Torous (1996).
150 15. Calibration to Standard Instruments
FUrther, De Jong (1997) and Babbs and Nowman (1999) demonstrate how
to implement affine multi-factor term structure models on panel interest rate
data. Nunes and Clewlow (1999) present estimation results on simulations
and real market data for extended Kalman filter algorithms using interest
rate caps and swaptions.
In the following sections we calibrate our term structure model to stan-
dard fixed income instruments. For this purpose we customize the general
state space model as presented in chapter 3 to reflect the time-series and
cross-sectional properties of our theoretical term structure model at best
and setup the Kalman filters along the lines of chapter 4. We implement
two special types of Kalman filter algorithms, one for linear and one for
non-linear functional relationships in the measurement equation, in order
to calibrate our term structure model to the liquid US markets of Treasury
Securities, LIB OR rates, and swaps.
(15.1)
for the underlying factors of our term structure model. This stochastic
representation of the evolution of the state variables can be integrated to
yield the solution125
Xt eTP(t-s)xs + J
t
eTP(t-()p,JI>dt+ Jt
eTI'(t-()udWJI>
s s
eTP(t-()udWJI>
for 8 < t. From this almost explicit integral solution we obtain for the dis-
crete time interval of!:1t = [8, t] the following exact discrete time equivalent:
with 'lJt = Jt
eTP(t-()udWIP .
t-Llt
i.e. the constant is given by m =- (TIP) -1 ,."IP, we obtain for equation (15.2):
Thus, we derived our setting for the transition equation of the state space
model as
(15.3)
with Ct ('ljJ)
<I>t ('ljJ)
[0,0]' ,
[ (T..a,)
E ['11tIFt-lltl = 0, and
Rq,q r
Var [11tIFt-lltl = [ a
2 para>. ] flt
par~>. a >.2
for the matrix and error term specifications. The obtained transition equa-
tion can be used in our empirical models independently of the interest data
we calibrate our term structure model to. However, the functional rela-
tionship of the measurement equation is conditional to the fixed income
instruments we use as sample data. Based on the sample data we use for
calibration we need to differentiate between implementations of the linear
and the extended Kalman filter algorithm as we show in the following two
sections.
10y
Figure 15.1: Estimated Zero Bond Yield Curve from 1970 to 1998
In terms of maturities used for our state space model, i.e. for the ob-
Maturity 3M 6M lY 2Y 3Y 5Y 7Y 10 Y
MEAN 0.0684 0.0709 0.0731 0.0761 0.0777 0.0802 0.0816 0.0828
STD 0.0271 0.0272 0.0263 0.0248 0.0237 0.0226 0.0220 0.0209
MAX 0.1619 0.1640 0.1621 0.1586 0.1573 0.1535 0.1507 0.1477
MED 0.0613 0.0644 0.0679 0.0702 0.0724 0.0749 0.0763 0.0776
MIN 0.0278 0.0288 0.0309 0.0380 0.0420 0.0435 0.0430 0.0451
SK 1.1801 1.1135 1.0447 1.0292 1.0448 1.0425 1.0247 1.0133
KU 4.3848 4.1363 3.9071 3.7443 3.6703 3.4985 3.3949 3.3956
Note:
For each maturity we report the sample mean (MEAN), standard deviation
(STD), maximum values (MAX), median (MED), minimum values (MIN),
skewness (SK), and kurtosis (KU).
From table 15.1 we infer that the longer interest rates are less volatile
than the shorter rates. This calls for a closer representation of the shorter
part of the term structure. The 10-years rate, for example, is only 77
percent as volatile as the 3-months rate. As we have seen in section 13.2.2
our interest rate model is capable to cover such volatility structures. Thus,
156 15. Calibration to Standard Instruments
we select the four maturities containing three months, one year, five years,
and ten years to govern the measurement equation (15.5).
Maturity 3M 6M lY 2Y 3Y 5Y 7Y lOY
3M 1
6M 0.9963 1
lY 0.9849 0.9944 1
2Y 0.9584 0.9713 0.9890 1
3Y 0.9355 0.9500 0.9732 0.9959 1
5Y 0.8973 0.9128 0.9419 0.9789 0.9928 1
7Y 0.8713 0.8874 0.9196 0.9639 0.9826 0.9973 1
10 Y 0.8484 0.8644 0.8985 0.9480 0.9702 0.9905 0.9972 1
1.00
~
0.98
3Y
0.96
2Y
0.94
g
~ 0.92
~ 3M 6M
8 0.90
0.88
0.86
0.84 +--...--,.--.----,r-----.----,--.........----r---.----,
o 2 4 6 8 10
Maturity
Figure 15.2: Cross-Correlations among Interest Rates
133For the shortcomings of a single-factor model see, for example, Chen (1996b).
13 4 This multivariate analysis technique is described in detail, for example, in Theil
(1971, sec. 1.9) and Jolliffe (1986). For testing cointegration with principal component
methods see, for example, Phillips and Ouliaris (1988).
135For a discussion of cointegration in the case of the term structure of interest rates see,
for example, Campbell and Shiller (1987). The formal development of the key concepts
of cointegration is given in Engle and Granger (1987).
158 15. Calibration to Standard Instruments
Based on our data set we further examine whether we are able to dras-
tically reduce the number of independent factors with little loss of informa-
tion. This will further support our proposed term structure model in that it
only uses two factors, the instantaneous interest rate and the market price
of risk. The results of the principal components analysis are presented in
tables 15.3 and 15.4. The original variables are the returns of the different
interest rates from the yield curve. 136
Table 15.3 shows the resulting values for the explained fractions of vari-
136For the principal component analysis we use all eight interest rates with maturities
given in table 15.1.
15.3 US Treasury Securities 159
0.12
0.10
0.08
0.06
(/)
C> 2. Factor
/
c:
'6 0.04
<II 4. Factor
0
...J
0.02
----- -.-.-.-.-.-.-.-.-
~
u.
0.00
.- -
-0.02 ~3.Factor
-0.04
·0.06
0 2 4 6 8 10
Time to Maturity
ance. The first line consists of the separate numbers of the accounted vari-
ance. Therein, the first principal component already contains 87.1 percent
of the total variance incorporated in the original variables. The second fac-
tor still has an explaining power of 9.8 percent. Thus, as can be seen from
the second line, the first two factors add up to 96.9 percent of cumulated
accounted variance. An additional third factor can raise the possible vari-
ance by another 1.8 percent, whereas the fourth and fifth factor only reveal
another 0.7 and 0.3 percent, respectively. Thus, by implementing our two
factor interest rate model we at best only loose 3.1 percent of the original
variance inherent in the data sample.
Further, we present the factor loadings of each of the first five princi-
pal components in table 15.4. These weights represent the sensitivity of
the original variables to the common factors. Thus, we can interpret the
loadings as the impact of the factors on the different yields of the term struc-
ture. The results on the weights of the first three principal components are
graphed in figure 15.3. For the first factor weights we obtain positive load-
ings in a fairly constant range of 0.1126 to 0.0570. This means that the yield
changes caused by the first factor essentially represent a parallel change of
the yield curve with the short maturities being somewhat more sensitive to
the factor movement. This result gives evidence to the interpretation of the
first principal component as being the average level of the yield curve. The
second principal component is made up by continually increasing weights
160 15. Calibration to Standard Instruments
across the term structure. The factor loadings are of opposite signs at the
short and long end of the maturity spectrum with a minimum at the 3-
months yield of -0.0521 and a maximum of 0.0320 for the lO-years rate.
Thereby, the second factor lowers the interest rates with shorter maturities
and increases the longer maturity rates which leads to refer to this factor as
steepness factor. Besides the first two principal components the remaining
higher order factors exhibit significant lower weights. The highest weight
of 0.0209 we obtain with the third factor. For this factor the loadings are
positive at the very short end and for maturities above five years. In the
intermediate maturity range we observe negative factor loadings. Thus, the
third principal component makes the yield curve twist in the way that it
raises it at the ends and at the same time lowers interest rates with interme-
diate maturities. The subsequent factors show even lower weights, at most
half of the weights of the third factor, and do not exhibit distinct patterns
in that we see theoretical interpretable facts. We conclude that for our data
set we obtain generally comparable results to those presented in the study
of Litterman and Scheinkman (1991).
In the light of these findings, the reduction of modeling the yield curve
with only specifying the dynamics of two factors, as with our theoretical
term structure model, is further given ground from a statistical analysis of
the historical interest rate data. Next, we present the estimation results on
the parameter values of our term structure model before we turn to analyze
the behavior of the underlying state variables.
1
y(t,T) = - T [[A(t,T),B(t,T)]x+C(t,T)]
-t
from section 13.2.1. Restating this affine structure for a vector of 9 zero
bond yields y (t, T) and the driving factors, the short rate r (t) and the
15.3 US Treasury Securities 161
(15.5)
where the modified factor loadings are given by 0* = [0 (t, Tt) / (Tl - t),
o (t, T2) / (T2 - t), ... ,0 (t, Tg) / (Tg - t)l' with 0* E {A*, B*, C*}. In the
measurement equation we further add a normally distributed measurement
error term €t (1/J) for that the observable yields Yt are only measurable with
noise. For the moment specification of the error term €t (1/J) we assume
that the errors have zero mean, are serially uncorrelated, and have a time
invariant covariance matrix. Notationally, we have the assumptions
IE [€tIFt_~t] 0, and
Var [€tIFt_~tl 17; • Igxg
which adds with 17e one more parameter to our parameter space 1/J.137
Thus, with the functional mapping of the state variables et
on the ob-
servable yield vector Yt from the bond market, we establish our specific
measurement equation. Based on the state space model defined by the tran-
sition equation (15.4) and the measurement equation (15.5) we are able to
run a linear Kalman filter algorithm. Applying the filter algorithm we esti-
mate the parameter values 1/J as laid out in chapter 5 by means of maximum
likelihood. Based on the term structure data set of US Treasury Security
issues we estimate the parameter values for the full available period ranging
from 1970 to 1998.
In table 15.5 we present the results for the entire sample period of our
two-factor model as well as of the reduced model which we derived in section
13.3.1. For the parameters common to both models we come to comparable
estimates besides for parameters /-LA and 17f:. The parameter /-LA actually
corresponds to the constant market price of risk A in the reduced model.
13 7 Notethat this homoskedastic error term specification leaves space to incorporate the
full cross-correlation structure among interest rates of different maturities as given in
equation (13.22) in order to obtain a close mapping with the empirical values reported
in table 15.2.
162 15. Calibration to Standard Instruments
A2 -0.116206 0.172479
P 0.710261*** 0.033481
(7c; 0.001400**· 0.000036 0.006011*** 0.000125
Notes:
Statistically significant parameter estimates at the 1%-, 5%- and
lO%-levels are denoted by ***, **, and *, respectively. The para-
meter J.L),. corresponds to the constant A in the reduced model.
With a value of -0.2552 the constant value actually lies close to the average
value of the stochastic market price of risk in our model which is given in
table 15.7 with a mean value of -0.2624. Looking at the values of the
measurement error term (7c; we see a significant improvement by using the
two-factor model. Adding one more factor lowers the error term from 60
b.p. to 14 b.p. Overall, the parameter estimates are statistically significant
besides the two level parameters of the market price of risk process, J.L)"
and A2. This, however, is not surprising since drift parameters generally
show wide confidence intervals. The mean-reversion parameters /'i,r and /'i,),.
imply mean half-lives, i.e. the expected time for the state variables to return
halfway to their long-run means, of 9.3 years and 1.1 years, respectively.
These results correspond to the estimates on the volatility parameters in
that (7r = 2.25% is much lower than the volatility of the market price of
risk with (7),. = 38.6%. The point estimates on the parameters J.Lr' /'i,r, (7r,
and (7c;, which are comparable to the one and two-factor specifications of
15.3 US 'freasury Securities 163
term structure models tested by Duan and Simonato (1995) and Babbs and
Nowman (1999), are found to be in line with their reported findings. For the
correlation factor p in our model we obtain a point estimate of 0.71. This
demonstrates that in times of high short rates the market price of short
rate risk A is also high which means that market participants then parallel
become more risk taking. 138
0.8 0.025
0.6
0.4 j_---:.;;---_ _ _ _ _ *""---- 0.020
0.2 CJ}../
Figure 15.4: Estimated Values for the Parameters t/J over Time
138For the impact of the state variables on the term structure compare section 13.4.1.
164 15. Calibration to Standard Instruments
f.-lr 0.0100 0.0017 0.0121 0.0115 0.0104 0.0088 0.0069 -0.48 1.90
/'i,r -0.1577 0.0396 -0.0744 0.0097 0.0090 0.0083 -0.2222 0.37 2.45
ar 0.0159 0.0070 0.0238 0.0085 0.0100 0.0118 0.0070 -0.02 1.14
f.-l),. 0.0008 0.0001 0.0013 0.0120 0.0087 0.0117 0.0007 2.10 7.67
/'i,),. -0.6801 0.7048 -0.0109 0.0119 0.0089 0.0096 -1.7787 -0.37 1.50
a),. 0.4049 0.0766 0.6373 0.0097 0.0076 0.0096 0.3238 1.27 4.63
0),. -0.2608 0.1275 -0.0918 0.0113 0.0100 0.0109 -0.4434 -0.01 1.29
>'2 -0.0442 0.1473 0.2139 0.0105 0.0089 0.0096 -0.2379 0.21 1.43
Notes:
a) In the statistics for the parameter estimates, UQ and LQ denote the upper and lower
quartile respectively; the median is abbreviated by MED.
b) For each parameter we also report the skewness (SK) and kurtosis (KU).
the literature that the shift in the monetary policy during the years of 1979
to 1982 caused a structural break in the interest rate process. 139 This can
especially be seen from the volatility of the short rate a r which starts at
225 b.p. fairly constant until 1980 where it falls monotonically to a value
of 7 b. p. for the minimum sample length. Therein, the impact of the high
volatilities during these years becomes evident. Also, the correlation factor
between the two state variables changes from positive values of 0.7 to -0.3
which marks a shift in the market participants changes of risk aversion
attitudes. Further, during time the mean-reversion tendency of the market
price of risk with a maximum of -1.78 does not become significant any
more after 1982 whereas the central tendency parameter of the short rate
/'i,r rather stays constant with maximum and minimum values of -0.07 and
139See, for example, Duan and Simonato (1995) and Hansen (1998).
15.3 US Treasury Securities 165
(15.6)
for the state variables x under measure M E {lP, Q}. Choosing the market
price of risk as the dependent state variable we divide the given equations
by each other which results in the differential equations
These equations define the phase space of the dependent state variable )....
The solutions to these first order differential equations can be plotted in a
vector field as graphed in figure 15.5. 140 The direction, the angle, and the
length of the vector line segments represent, respectively, the movement,
the way, and the speed of the state variables from the starting values at
that point on an orbital trajectory.
In order to algebraically analyze the stability 141 of our state space system
we first calculate the eigenvalues of the system as
eVp =
{ [ 0] [-0.5115]}
1 ; 0.8593 ' and eVQ =
{ [ 0.0557] [-0.5456]}
0.9984; 0.8380 .
140 Thevector field is based on the estimated parameter values reported in table 15.5.
141 Fora systematic classification of the possible equilibria of a two dimensional linear
system see, for example, Brock and Malliaris (1989, ch. 3.6).
166 15. Calibration to St8Jldard Instruments
I I I I 1 ""'~. • • j j! j I
I I I I 8l
f
,
~
f
.~.
,
I ,
• '! j
I I I '"
~ 0.0
...... t
, , , , .,..... ,~ ; • j j j
I I
, ... ~ I • .: • ~ • t " , .. ···t: • ; 1 1 j
~~"··"··l+ttt t "".(-"';11
-0.1 ' , , , ...•. I~ , , ! • j j -0.1 1 I I 'i I ". , • ;
, , f , j • , , ,
, ....• j j tIl ~II"""
" " , , [ ' - - ·. • • '!I
ttl t , , • ~"'~
I I , t t l : • ' .. ~ , • • ; t tt il 'I~
-0.2 I I t t , I: I • ....• , , -0.2 h-t..,.....,~t,...",t~..,.....,,..;..;,...,,,1~"Tt~'.....'
3% 6% 9% 12% 15% 18% 3% 6% 9% 12% 15% 18%
r" - values rQ - values
i.e. they are distinct and lie both below zero. In this case we obtain a
stable tangent node at xl!' = [0.0950, -0.1311]' and xlQ = [0.1106, -0.0517]',
respectively, as indicated by the point of intersection of the two dotted
lines. 143 The dotted lines represent the eigenlines of the system which are
lines through the asymptotic point of the system. They can be derived from
the general solutions to the partial differential equation of (15.6)
0.0950]
[ -0.1311 [ 0] -0.4434 [ -0.5115] -0.0744 d
+ Cl 1 e + C2 0.8593 e , an
0.8
0.6
0.4
Eigenline2
I
0.2
en 0.0 .......
(I)
~
(ij
> -0.2
I
a..
<-< -0.4 V/-'_'--01/82
-0.6
-0.8
01/83
-1.0
0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18
rP - values
with arbitrary constants Gi, with i E {I, 2, ... ,4}. The eigenlines are found
by consecutively setting the constants Ci equal to zero and solving the para-
metric representation of the functional relationships between r and A for the
state variable A. Generated by the eigenvectors [0,1]' and [-0.5115,0.8593],
for the dynamic system under the measure IP' we obtain
Notes:
a) In the statistics for the state variables x, UQ and LQ denote the upper and
lower quartile respectively; the median is abbreviated by MED.
b) For each state variable we also report the skewness (SK) and kurtosis (KU).
From the numerical results we first see the short rate in a range of
realized values from a minimum of 2.3% to a maximum of 16.5% with a
mean value of 7.0%. These values closely resemble the statistics of the
short rate (3-months rate) included in the data sample as can be seen from
the reported values in table 15.1. Also, the positive skewness stemming
from the 1979-1982 period (see the short rate values in graph 15.6) and
the excess kurtosis are retained. Second, with our filtering technique we
can extract the market prices of risk over time. The values we obtain for
At take a minimum of -0.98, a mean value of -0.26, and a maximum of
0.81. With an upper quartile value of close to zero we can infer that the
market participants were risk averse over 75 percent of the sample period
from 01/1970 - 12/1998. Finally, we see from graph 15.6 how the two state
variables evolve over time. Therein, the years of 1979 to 1982 mark a period
of high short rates along with risk taking values for the market price of risk.
We further note that the path of the state variables actually crosses the
asymptotic lines when we include the historical randomness. The central
tendency is especially effective when the realized values are far away from
their long-run means.
15.4. Other Liquid Markets 169
1 - P (t, T)
L (t, T) = (T - t) P (t, T)' and (15.7)
1 - P (t, T)
k (t, T) = n (15.8)
OLP(t,Tj )
j=l
as defined in equations (11.3) and (14.17). From these formulae we see that
the LIBOR curve represents the natural extension of the swap rate curve
144Liquidity in the US interest rate swap market is generally agreed to be due to the
high demand on swaps and the ease of hedging swaps with the US Treasury and futures
markets. Also standard terms introduced by the ISDA and BBA in 1985 have assisted
the growth in the swap markets. See, for example, Das (1994).
170 15. Calibration to Standard Instruments
for shorter maturities for which no swap contracts are traded. 145 Thus, in
order to cover the whole spectrum of the term structure we choose LIBOR
and swap rates as observables 0 (t, Ii) E {L (t, Ii), k (t, Ii)} with which we
specify the non-linear measurement equation as
in the case of the extended Kalman filter. With the discount bonds P (t, T)
and P (t, 1j) being a function of the state variables the derivatives for the
LIBOR and the swap rate can be calculated from equations (15.7) and (15.8)
in closed-form. The moment assumptions regarding the error term et ("p)
are identical with those in the linear case. With equation (15.9) we obtain
a linearized version of the non-linear functional relationship between the
unobservable state variables et and the LIBOR and swap rates Yt which are
observable in the financial markets.
In analogy to the linear case, we are now able to run the extended
Kalman filter algorithm of section 4.5 based on the state space model defined
by the transition equation (15.4) and the measurement equation (15.9). The
empirical implementation and the estimation results are presented in the
next section.
in the previous section we are able to implement our term structure model
on LIBOR and swap rates. The calibration procedure is based on interest
rate data sampled monthly over the longest available period from 04/1987
to 01/2000. We obtain 154 observation dates and at each date we use four
distinct maturities as in the case of the Treasury Security sample. In both,
the LIB OR and the swap market, we select the maturities with the aim of
capturing a maximum spectrum of the available term to maturities. For the
LIB OR rates we choose to include the 1-, 3-, 6-, and 12-months rates and for
the swap rates the maturities of 2, 3, 5, and 10 years. The closer maturity
choices at the short end are intended to better capture the information of
this more volatile part of the term structure.
Note:
For each data series we report the sample mean (MEAN), standard
deviation (STD), maximum values (MAX), median (MED), minimum
values (MIN), skewness (SK), and kurtosis (KU).
In table 15.8 we report the summary statistics separately for each time-
series. For the short-term market of LIB OR rates the values of the statistics
for the different maturities show very similar results. The term structure
rises on average with mean sample values of 5.98% for the one-month rate
up to 6.39% for the 12-months rate. The volatility of the rates slightly
decreases with higher maturities. For the criteria of skewness and kurtosis
we see modest positively skewed rates that do not show excess kurtosis.
172 15. Calibration to Standard Instruments
Note:
Statistically significant parameter estimates at the 1%-, 5%- and lO%-levels
are denoted by ***, **, and *, respectively.
These sample facts mostly resemble those for the swap rates. Due to the
broader maturity spectrum we only see the term structure of swap rates
increasing up to a 10-years rate of 7.71%.
The estimation results on the parameter values 1/J of our term structure
model are given in table 15.9. We calibrate the term structure model to
three different panel data: The four LIBOR rates, the four swap rates,
and a combination of the 1- and 12-months LIBOR rates along with the 5-
and 10-years swap rates. The combined data set is motivated by the close
relation of the two types of interest rate markets as we notice from equations
(15.7) and (15.8) with the aim of obtaining a smooth term structure pattern
over all available maturities. First, we compare the model fit to the three
different samples in terms of the measurement error. The LIBOR rates
exhibit a standard deviation of 7 b.p. which falls to 3 b.p. for the swap
sample but increases to 20 b. p. for the combined data. This indicates that
the swap market is the most homogeneous sample despite its high spanning
of maturities from 2 to 10 years. Second, the results on the point estimates
15.4 Other Liquid Markets 173
of J.Lr' ""r, and a r are comparable among the three samples in use. Third,
the parameters corresponding to the stochastic specification of the market
price of short rate risk are estimated with plausible values for the swap and
the combined data. However, based on the LIB OR data we obtain diverging
estimates for the parameters J.L)..., "")..., a)..., (h and ).2. This provides evidence
that a sample including only a small one year window of the term structure
is not a rich enough sample in order to get reliable estimates of the market
price of risk parameters. Fourth, the mean-reversion parameters ""r and "")...
estimated for the swap and combined data at around -0.045 and -0.65
imply mean half-lives of 15.4 years and 1.1 years, respectively. Finally, the
correlation coefficient is estimated at significantly positive values of around
+0.6. Overall, all parameters are estimated at statistically significant levels
besides the estimates for J.L)..., "")..., ).2, and p based on the LIBOR data and
).2 for the combined data.
16 Summary and Conclusions
In part III of this study we present a new continuous time term struc-
ture model and provide its implications on risk management and the pric-
ing of different fixed income derivatives based on a model calibration to
standard fixed income instruments. The theoretical and empirical analysis
contributes to the term structure literature in that it shifts the theoret-
ical attention of modeling static investor tastes towards capturing beliefs
of investor risk attitudes towards interest rate risk using an affine model
structure.
In the introductory chapter we consider the specific financial instru-
ments most relevant in modeling fixed income markets and describe the
peculiarities of an incomplete market such as the market for fixed income
instruments. In the following chapter 12 we motivate our desire to build
a new term structure model that takes into account time varying market
prices of interest rate risk. The model we suggest chooses the short interest
rate and its market price of risk as the two explaining state variables where-
upon we establish an affine term structure model. In chapter 13 we present
initial characteristic results and implications of the proposed term structure
model. First, we derive a closed-form solution for arbitrage-free prices of dis-
count bonds. Thereupon, we are able to state and analyze the implied term
structures of interest rates and volatilities. In a comparative analysis we
examine limiting cases of the term structures and the influence of the state
variables and the model parameters on the shape of the term structures.
We identify the first model factor as a level factor and the market price of
risk as steepness factor that captures at the same time the investors' risk
attitudes. Chapter 14 shows how to implement risk management techniques
based on the proposed model. With the duration concept we illustrate the
implementation of a simple model for interest rate risk management. More
176 16. Summary and Conclusions
17.1 Overview
In part IV of the study we develop a continuous time pricing model for one
of the latest cash markets to be transformed by derivative securities - the
electricity market. Based on available data on electricity derivatives we fo-
cus our attention on the valuation of short-term electricity forwards. First,
we investigate the more mature futures markets contracted on established
commodity underlyings and describe proposed pricing models for traditional
commodity futures contracts from the literature. Thereby, we gather rele-
vant information about the commodity electricity in order to build an ap-
propriate valuation model to price electricity forwards in chapter 18. From
given model assumptions we derive a closed-form valuation model for elec-
tricity forwards using risk neutral pricing techniques. The suggested model
especially captures the high and time varying volatility seen in electricity
prices. In chapter 19 we present an empirical adaptation of the theoretical
pricing model in state space form. Using maximum likelihood estimation
based on extended Kalman filtering we report empirical results on electric-
ity data from the largely deregulated Californian electricity market. In the
last chapter we conclude.
, ,
[.Com
·Coff..
'Cotton
t·
• Uvecata.
Pork Bellies
• Lean Hogs
r-~r
• AlUminum
'Lead
r-G~
' Silver
-Platinum
[. Crude Oil
'HeetingOiI
'Natural Gas
[ • Electricity
146For historical, institutional and regulatory issues of futures exchanges see, for exam-
ple, Markham (1987).
147See, for example, Duffie (1989).
148See UN (1996).
17.2 Commodity Futures Markets 181
Mercantile Exchange (CME), and the Coffee, Sugar, and Cocoa Exchange
(CSCE). Futures contracts are traded on vegetable goods such as grains
(corn, oats, wheat, and rice), foodstuffs (coffee, cocoa, orange juice, and
sugar), textiles (cotton), oil and meal (soybeans, -meal, -oil, sunflower seed
and oil), timber products (lumber and plywood) and on animal goods such
as live hogs, live cattle, lean hogs, and pork bellies. For many of these
commodities various futures contracts are available due to the necessary
standardization of futures contracts. The specification of contracts for dif-
ferent grades and types of the commodity ensures a homogeneous deliverable
commodity. In the metallurgical futures market the underlying commodities
are classified by either being industrial or precious metals. The spectrum
ranges from copper, aluminum, palladium, lead, zinc, nickel, and tin up
to the other end with silver, platinum, and gold. Major metal futures ex-
changes are the Commodity Exchange (COMEX), the New York Mercantile
Exchange (NYMEX), and the London Metal Exchange (LME).
Among the various commodities the energy market is the most recent
market to be transformed by derivative securities and risk management. A
major reason for this late development is considered to be the postponed
deregulation of energy markets giving the spot markets ultimately the nec-
essary liquidity and price volatility. During the 1960s and 1970s prices of oil
as one of the most crucial commodities in the world were tightly controlled
by the Organization of Petroleum Exporting Countries (OPEC) together
with large globally integrated oil companies. The trading of oil futures at
the NYMEX began with heating oil futures in 1978, followed by crude oil fu-
tures introduced in 1983, and unleaded gasoline futures introduced in 1985
establishing the NYMEX as the world's leading energy futures exchange. A
further major energy commodity besides oil is natural gas which accounts
for almost a quarter of total US energy consumption. The gas industry has
been deregulated since the enactment of the Natural Gas Policy Act of 1978
changing it to an industry that largely operates as a free market today. An
accompanying futures contract to the natural gas spot market centers was
launched in 1990 traded on the NYMEX.
The transformation of the electricity market into a competitive mar-
ket has been essential for the development of electricity futures markets.
Since electricity is considered a regional commodity not only deregulation
in electricity markets but also in the transmission system has been crucial in
182 17. Introduction
149For a detailed record of deregulation history in the US power market at the federal
and the state levels see, for example, Bell and Lilyestrom (1997).
150 See CalPX (2000a).
151See CalPX (2000b).
152See Johnson and Sogomonian (1997).
17.3. Pricing Commodity Futures 183
(17.2)
price including carry costs and returns under the martingale measure Q,
investors follow a cash-and-carry arbitrage. The arbitraguer borrows funds
to buy the spot asset, takes a short position in a futures contract, and
carries the deliverable asset until the futures delivery date. Second, if the
futures price is below the discounted expected spot price a reverse cash-and-
carry arbitrage can be followed. The arbitraguer sells the spot asset short,
lends the proceeds, and takes a long position in a futures contract. Together
both arbitrage possibilities ensure that the equality of equation (17.2) holds.
However, in view of implementing the two arbitrage arguments there is a
marked difference. The cash-and-carry arbitrage can be set up and pursued
easily with financial underlyings as well as storable commodities in that it
only asks for carrying the deliverable asset until the futures delivery date.
Thus, the cost-of-carry model can be used to determine an upper bound on
futures prices in that the futures price F (t, T) can not exceed the cost of
stored underlying assets for delivery. If in this case the future-spot price
relationship is mainly characterized by the carry costs the futures market is
said to be in contango. This is the standard market situation where carry
costs are higher than returns from carrying the asset and the futures prices
lie above the spot prices. However, in the case of reverse cash-and-carry ar-
bitrage we need to sell the deliverable asset short in order to realize risk-less
profits. The implementation of this arbitrage can either be complicated or
impossible when there are no or little amounts of deliverable assets avail-
able in storage. 156 Thus, the cost-of-carry model can not provide a lower
bound for futures pricing. The modified cost-oj-carry Jormula for non-carry
commodities resolves this pricing problem by including a convenience yield
CY:157
The convenience yield is the return that users of the commodity realizes
for carrying inventory of the spot commodity. This scenario of positive
156The reverse cash-and-carry arbitrage also breaks down in cases of storable commodi-
ties for which there exist no leasing market as, for example, the natural gas market. See,
for example, Leong (1997).
157The concept of a convenience yield was first described by Kaldor (1939) and Working
(1949).
17.3 Pricing Commodity Futures 185
returns or negative storage costs arises when the users will not be willing
to sell their inventory in the spot market and buy futures contracts that
will replenish their supplies at a later date. A commodity's convenience
yield can be different for various users and can vary over time. The return
will be at highest when there are spot shortages of the spot commodity. In
instances where the convenience yield plays a major role for users the spot
prices will lie above futures prices, Le. a market situation which is known as
backwardation.
15 8 This
modeling approach can be traced back to Ho and Lee (1986) and Heath, Jarrow,
and Morton (1992).
17.3 Pricing Commodity Futures 187
also followed by Cortazar and Schwartz (1994) and Amin, Ng, and Pirrong
(1995).
The category of equilibrium models assumes the underlying commodity
to be a non-traded asset. Thereby, the no-arbitrage hedging argument does
not apply and the value of a futures contract is determined by the expected
spot price at maturity of the contract. One-factor models based on a dy-
namic specification of the spot price behavior are presented by Ross (1997)
and Schwartz (1997, Modell). Ross (1997) considers an underlying asset
which cannot be stored and held over time except at high cost making the
pricing of derivative contracts impossible with the usual risk neutral hedg-
ing argument. The quite general volatility specification of Ross (1997) is
put in concrete terms by Schwartz (1997, Modell) who assumes a volatility
that is linear in the spot price. In the two-factor model of Schwartz and
Smith (1997) the log-spot price behavior is decomposed into a process to
capture the equilibrium level ~t and a dynamic specification of transitory
disturbances Xt from the long-run leveP59
FUrthermore, two regime models are proposed in the literature which
are hybrid approaches grasping both arbitrage and equilibrium pricing sit-
uations of futures contracts. The regime switching approach of Schobel
(1992) models the spot price P via the dynamics of the level of inventories
I. Dependent on the availability of the spot asset the model captures a com-
plete and an incomplete market situation which are given for a spot asset
that is, respectively, available or too scarce for arbitrage purposes. Based
on Deaton and Laroque (1992) and Deaton and Laroque (1996), Routledge,
Seppi, and Spatt (2000) determine spot and forward prices endogenously
from an immediate net-demand process and the resulting dynamic structure
of equilibrium inventory.160 They additionally discuss a tractable two-factor
augmentation of their model which leads to better calibration results of the
volatility term structures. In the paper by Buhler, Korn, and Schobel (2000)
a unifying approach a la Schobel (1992) is presented. Conditional on the
levels of current spot prices and inventories the derived futures prices span a
wide price range which is equivalently obtainable from using a cost-of-carry
model up to using a pure equilibrium model.
159Schwartz and Smith (1997) show that their model is equivalent to the convenience
yield model of Gibson and Schwartz (1990).
160 However, it is not a continuous time pricing model but formulated in discrete time.
18 Electricity Pricing Model
(i) The non-storability feature means that due to a lack of any power
inventories electricity must be produced at exactly the same time as
it is consumed. This creates a load-matching problem in that the
utility industry needs to discover the value of lost load and decide on
the optimum amount of reserve margin to provide. 161
(iv) Due to the load-matching problem electricity prices show highly volatile
seasonal, weekday and even intra-day patterns. These patterns per-
sist since there are no arbitrage possibilities due to the non-storability
feature of electricity. The hourly changing price behavior needs also
to be considered in specifying appropriate electricity derivative con-
tracts which are useful for pricing, hedging, and risk management
applications. 164
_200
..s:::
~
:iE
-
......
~ 150
Q)
u
.;::
a.
~ 100
'0
.;::
~
iIi 50
04-----r------r~~---,--------~----__,
with the drift of the variance rate process specified as J-t v = K (8 - Vt). We
further assume that dWE,tdW:,t = pdt to cover possible correlation between
the two state variables.
of assets that cannot be stored and held over time except at high cost. Given
electricity as an underlying commodity that is non-storable, we face a market
situation that is incomplete. In this case both state variables - the spot price
of electricity St and its variance rate Vt - cannot be hedged. 172 In thus, the
state variables must be risk adjusted by incorporating market prices of risk in
order to develop arbitrage-free pricing formulas for electricity derivatives. 173
Using the market prices of risk, we are able to derive the risk neutral
dynamics of the state variables in order to price electricity derivatives under
the corresponding martingale measure in an incomplete market. In order to
obtain the corresponding martingale process for the electricity spot price,
we use the Girsanov transformation dW~t = dW%,t + A* yVidt with a time
invariant market price of risk A*. Based on equation (18.1) the stochas-
tic development of the transformed spot price X t = In St leads under the
Girsanov transformation to the stochastic behavior
(18.3)
for the first model factor using the substitution A = ~ + A*. In the case of
the spot price variance rate Vt as the second state variable, we use the same
adjustment as suggested in Heston (1993) . This is given by
(18.4)
with the market price of risk being Av (St, Vt, t) = AvVt where Av is a con-
stant. Based on the risk neutralized processes (18.3) and (18.4) of the two
state variables we are able to price forward contracts arbitrage-free as we
demonstrate in the following section, using martingale methods developed
by Geman, Karoui, and Rochet (1995) and Scott (1997).
(18.5)
J + JVV;d~t
T T
XT = X t + /-l (T - t) - A Vt dt (18.6)
t t
for the relevant time horizon [t, TJ. Therein, we can further get an expression
for the time integral from the variance rate specification of equation (18.4).
Working analogously with an integral solution we get
F (t, T) =]EQ e
Xt+J1.{T-t)- ;>. J J
[t<:IJ{T-t)-{VT-Vt)+". yVidw:!t] + yVidw2 t
[ 0< v t ' t •
F (t, T) = ]E1Qi [ e
Xt+JL(T-t)-
I<
;>. v
[1t1J(T-t)-(VT-vt)+<TJ
t '
ytitd~t]
(18.8)
in which we can work with Ito's isometry on the last integral since the
Brownian motion WtlQi is independent of the variance rate variable. Thereby,
we obtain for the last term
J J
T T
Putting the results together and bringing the constant terms with re-
spect to the information available in F t outside the conditional expectation
we get the expression
F (t, T) = eXt-;Vt+(JL-t;!p)(T-t)
With equation (18.9) we have almost derived the price of the forward
-1 V,dS)
contract F (t, T) besides solving for the remaining expectation
with kl ~ +.Av
.A - -(J-P 1( 2)
- 2 1 - P ,and
p
k2 = (J
18.3 Valuation of Electricity Forwards 197
y(t,x)=Ex [ e
- Jq(X.)ds
0t 1
f(Xt).
holds with the necessary initial condition being y (0, x) = f (x); A denotes
the differential operator. 176
~ a 2y(t,Vt)
2
2
a Vt a 2
vt
+ (()
/'i,
_ (
/'i, + \ )Vt) ay(t,Vt)
Av a
Vt
+ ay(t,Vt)
a
t
_ k ( )
- I VtY t, Vt ,
(18.11)
with the boundary condition for the value at maturity y (T, Vt) = exp (k 1VT).
For the solution for the differential equation we propose the exponential
affine structure
(18.13)
0, (18.14)
where we applied the separation of variables technique. 177 The first equation
(18.13) is known as a lliccati type differential equation 178 for which we get
the solution
-"'(} J
T
"'(}
(J [('" + Av) (T - t)
B(t,T') = 2"
177 See, for example, Ingersoll (1987, p. 397) for an application with interest rate con-
tingent claims.
178 This type of differential equation commonly arises with the squared Gaussian model
of Cox, Ingersoll, and Ross (1985b); see, for example, Rogers (1995, p. 99 f.).
179For a similar solution to equation (18.13) see, for example, Rogers (1995, p. 100)
and Stein and Stein (1991, p. 730).
18.3 Valuation of Electricity Forwards 199
for the second function in equation (18.12). Thus, we have finally calculated
an expression for the remaining expectation y(t, Vt) of equation (18.9).
F (t , T , S t, v·t, ~/.)
0/
= S t e(A(t,T)-;)vt+B(t,T)+(JL-~p)(T-t) , (18.17)
where the functions A (t, T) and B (t, T) are given by the expressions of
equations (18.15) and (18.16), respectively, with 1/J = {/-L, K" (), a, p, A, Av}.
J J J
s s s
s s s
for a period oftime T = [t, s]. These equations fully contain the information
to derive the conditional moments of the state variables. In order to calcu-
late the moments we employ a technique of putting up differential equations
that contain the conditional moments of interest as functions. With the rel-
evant boundary conditions we then solve these differential equations for the
moments. 180
For the conditional first moments we take the expectations on equations
(19.1) and (19.2). To calculate the conditional mean of Vs we take the
derivative with respect to the upper bound s of the integrals which yields
olElP[vsl.Ttl_ 0 IT]
oS - K + KJrJlU'1!'[ vs.rt· (19.3)
Exploiting the initial condition of IEI!' [Vtl .Tt ] = Vt we solve the differential
equation (19.3) as
(19.4)
Using the expression for lEI!> [vsl.Ttl from equation (19.4) we can calculate
the solution
(19.5)
180 This method is, for example, used by De Jong and Santa-Clara (1999) in the context
of a factor based interest rate model a la Heath, Jarrow, and Morton (1992).
181Compare Cox, Ingersoll, and Ross (1985b, p. 392).
19.1 Estimation Model 203
Given the conditional means EIP' [Xa IFtl and EIP' [Va IFtl we further need
to calculate the first terms in equation (19.6) to get solutions for the three
conditional second moments. Therefore, we first obtain the stochastic dy-
namics for the functions it = Xl, h = XtVt, and fa = v; using Ito's
formula. This yields the following dynamics:
(19.7)
with solving the first equation. Finally, we use these results along with
the statistical relationship stated in equation (19.6) to derive the following
expressions for the second moments. First, we obtain the solution to the
conditional variance of the second state variable from the third equation of
(19.7) as:
which matches the formula stated in Cox, Ingersoll, and Ross (1985b, eq.
19). Given the result of equation (19.8) we are further able to calculate the
covariance term of the state variables as
(19.9)
+7Jt(et-~t,¢) (19.11)
for the transition equation. However, within this equation the error term
7Jt (et-~t, 'I/J) still contains the state variable et-~t which is not observ-
able. Thus, we approximate the transition equation around (et-~t, 7Jt)
(et-~tlt-~t' 0) which results in:
206 19. Empirical Inference
= [~ (19.12)
Stlt-At =
(et-at,TJt,1/J )=(et-atlt-at,o,1/J )
[
VarF!' [ Xt 1Ft - Ad 1F
CovF!' [XtVt t - At ] ]
(19.13)
CovF!' [ XtVt 1Ft - Ad 1F
VarF!' [ Vt t - At ]
where we use the expressions of equations (19.10), (19.9), and (19.8) and
substitute values Vt-Atlt-At for Vt-At. Thus, we can finally state our transi-
tion equation of the state space model as
where the formulae of h t (et-Atlt-At, 0, 1/J), <Ptl t - At , and Stlt-At are given in
equations (19.11), (19.12), and (19.13).
Next, we work out the specification of the measurement equation where
we use the forward pricing formula derived in section 18.3. The idea is to
comprehensively use market information contained in both the spot electric-
ity prices as well as the information from the forward prices of electricity.
Bringing both market informations together, we formulate the measurement
equation as an identity between the observable market prices and the the-
oretical prices. However, the two pairs of prices do not need to fully match
in that we indude a measurement error to allow for some discrepancies be-
tween the market and model prices. We obtain the measurement equation
by exploiting the functional relationship of the state variables et with the
forward price from the original pricing formula given in equation (18.17)
(19.15)
o ].
2
ac:,F
(19.16)
ity wholesale and retail market has two core institutions, the California
Power Exchange (CalPX) and the California Independent System Operator
(CaISO). We briefly describe their operations on the Californian electricity
market, before we come to the data analysis and present our estimation
results. 183
The CalPX is a non-profit, public benefit corporation which is regulated
by the FERC. As a power commodity exchange it provides an efficient,
competitive marketplace by conducting an open, non-discriminatory trad-
ing process for qualified electricity suppliers and purchasers. The CalPX
determines the price of electricity in the hour-ahead and day-ahead mar-
kets of forward contracts to deliver a given quantity of electricity over a
one hour period. For each delivery hour the exchange collects demand and
generation bids in prices and quantities to settle for unconstrained market
clearing prices. Since its opening on March 31, 1998 the CalPX's day-ahead
market has traded between 80 and 90 percent of California's electricity mar-
ket. This day-ahead market is open from 6 a.m. to 1 p.m. when sellers and
buyers of electricity can submit their portfolio supply and demand bids for
24 hourly auction periods for electricity delivery on the next day. There-
upon, the market clearing prices are determined separately for each auction
hour. Besides providing marketplaces for electricity forwards, the CalPX
serves as scheduling coordinator in that it submits balanced supply and
demand schedules along with information on the generating units and the
locations for delivery to the CaISO.
The CaISO serves as the control area operator for the Californian elec-
tricity market where it provides open access to the transmission and man-
ages the real time operations of the transmission grid. The objective is to
provide reliable system operations, primarily keeping the right frequency
and providing sufficient generation, which are maintained by buying and
providing ancillary services. The CaISO matches the power output of the
generating units with the power demand within the controlled electric power
system including the scheduled interchange with adjacent control areas.
This real time balancing of load and generation is based on the coordina-
tion of the hour-ahead and day-ahead schedules dispatched by the CalPX.
On the basis of the real time operations, the CaISO determines and pub-
lishes real time spot electricity prices by hour.
183The Californian power market is explained in detail in, for example, CalPX (1999).
210 19. Empirical Inference
184The major factors for the price patterns are seen in the fuel price, load uncertainty,
variations in hydroelectricity production, generation uncertainty, and transmission con-
gestion; see, for example, CalPX (1998).
19.2 Data Analysis and Estimation Results 211
Notes:
Prices are stated in $/MWh. For each distribution of hourly prices
we report the mean, the standard deviation, the maximum, the
median, the minimum, the skewness, and the kurtosis. Summary
statistics across all time-series are given at the bottom of the table.
212 19. Empirical Inference
Notes:
Prices are stated in $/MWh. For each distribution of hourly prices
we report the mean, the standard deviation, the maximum, the
median, the minimum, the skewness, and the kurtosis. Summary
statistics across all time-series are given at the bottom of the table.
19.2 Data Analysis and Estimation Results 213
Notes:
The price differences are stated in $/MWh and their distribution is characterized by
values for the mean, the standard deviation, the maximum, the median, the minimum,
the skewness, and the kurtosis. Further, we report results on a ordinary least squares
regression of forward prices on spot prices. At the bottom of the table we show summary
statistics across all time-series.
214 19. Empirical Inference
prices. The null hypothesis for unbiasedness requires the intercept term
and the slope coefficient to be statistically not different from zero and one,
respectively. Looking at the results indicates that the null hypothesis can be
rejected at conventional significance levels for almost all 24 time-series. The
intercept term lies above zero at an average of 9.67 with a corresponding
standard error of 1.91. The slope coefficients show values that are almost
all significantly below one with an average of 0.71 at a standard deviation of
0.05. Thus, even for a time period of only one-day ahead the forward prices
are not an unbiased predictor of realized spot prices which demonstrates
the need for sound theoretical valuation models.
In order to capture the spot price behavior and to price the day-ahead
forward contracts we select to calibrate two theoretical pricing models to
the given electricity data. First, we choose to implement the one-factor
model as parametrized in Schwartz (1997, Model 1) which originally is in-
tended to capture the price processes on oil spot and futures markets. We,
however, use this equilibrium pricing model to describe the electricity spot
price behavior by the continuous time process
er 2
with a* J..l - 2", - A (19.18)
under the martingale measure Q including the time invariant market price of
risk parameter A. Based on equations (19.17) and (19.18) as transition and
measurement equation, respectively, we estimate the model using the stan-
dard Kalman filter setup. The results on the parameter estimates including
the two measurement errors erg,S and erg,F are shown in tables 19.4 and 19.5.
The coefficient of mean-reversion is estimated by 0.6 on average across all
24 time-series. The values translate to mean half-lives of the electricity spot
price with a minimum of 0.6 years for hour 16 and a maximum period of
3.8 years for the 22nd hour where the average lies at 1.1 years. The rather
not statistically significant estimates for the long-run mean J..l reach values
around 2.8. Related to the coefficient of mean-reversion is the er-parameter
which exhibits how volatile the spot prices fluctuate around their long-run
19.2 Data Analysis and Estimation Results 215
Note:
Statistically significant parameter estimates at the 10/0-, 5%- and lO%-levels
are denoted by ***, **, and *, respectively.
216 19. Empirical Inference
Note:
Statistically significant parameter estimates at the 1%-, 5%- and lO%-levels
are denoted by ***, **, and *, respectively.
19.2 Data Analysis and Estimation Results 217
JL A /'i, (J
Note:
Statistically significant parameter estimates at the 1%-, 5%- and lO%-levels
are denoted by ***, **, and *, respectively.
218 19. Empirical Inference
Table 19.7: Inference Results for the Stochastic Volatility Model (continued)
Note:
Statistically significant parameter estimates at the 1%-, 5%- and lO%-levels
are denoted by ***, **, and *, respectively.
19.2 Data Analysis and Estimation Results 219
the other parameter estimates. The expected rate of return J-t of the spot
prices is estimated at levels ranging from -0.98 up to 7.2. AB expected
for a level parameter, the inferences on the mean rate are not statistically
different from zero. The values for the market price of spot price risk >.
are deduced from the data at a negative spectrum starting from -1.1 and
ending at a value of -13.0. With the stochastic volatility pricing model the
market prices of risk are mostly estimated at significant negative levels. The
inferred negative values are in line with the results from the estimation of
Schwartz (1997, Modell). The highly negative parameter values for the
market price of risk further demonstrate that it is crucial to calibrate an
equilibrium forward pricing model in that the non-storable commodity of
electricity creates unhedgeable risk. This risk need to be compensated for on
the short sellers' side of the forward contracts. The remaining parameters
K, (J, and p describe the time-series behavior of the second state variable,
i.e. the variance rate. The variance rate is estimated with high tendencies
to revert to its long-run equilibrium value; the inferences on the coefficients
of mean-reversion K translate to mean half-lives with a minimum of 1.1
months and a maximum of 36.3 months. The volatility of the variance
rate is inferred at generally highly significant values with an average of 320
percent across all 24 time series. Finally, the correlation coefficient p, which
shows whether the two state variables generally move in the same or in the
opposite direction, is estimated at a negative average value of -0.15. Only
37.5 percent or 9 out of the 24 time-series exhibit correlation coefficients
above zero.
20 Summary and Conclusions
The symbols and notation used throughout this study are introduced at
their first appearance. The choices of symbols and notation are adopted
from common uses in the literature. In the following we enumerate and
explain the most frequently used symbols and notation separately for each
part of the study.
Part I
Cov[···I···] conditional covariance operator
lBJ [ •. ·1 ... ] conditional expectation operator
et vector of measurement errors
"It vector of transition equation noise terms
Ft variance-covariance matrix of the prediction error
F( ... ) distribution function
Ft information set up to time t
gt general non-linear measurement function
ht general non-linear transition function
Kt Kalman gain matrix
L ( ... ) log-likelihood function
MMSE minimum mean square error
MMSLE minimum mean square linear error
MSE (... ) mean square error
N( .. . ) normal distribution function
p ( ... ) density function
'ljJ vector of parameters
Et variance-covariance matrix of the state variables
t variable for calendar time
Xt general stochastic process
224 List of Symbols and Notation
Part II
A,B,C factor loadings
ARt,i abnormal return of security i on date t
BMt benchmark return
CAR cumulated abnormal return
dt infinitesimal time increment
dWt infinitesimal Brownian increment
E[ .. ·I···] conditional expectation operator
HPt,i holding period return of security i on date t
/'i, speed of mean-reversion
KU kurtosis
L( ... ) log-likelihood function
MAD mean absolute deviation
JL expected mean return of the net asset value
NAVt net asset value of closed-end funds
Pt market price of closed-end funds
1ft instantaneous dynamic premium
PREMt empirical premium
1/1 vector of parameters
r riskless interest rate
RESt residual premium
p coefficient of correlation
RMSE root mean square error
(Jx volatility of the net asset value
volatility of the dynamic premium
SK skewness
SSE sum of squared errors
STD standard deviation
f) long-run mean value
VAR[···I···] conditional variance operator
Xt logarithm of the net asset value
et vector of state variables
Yt logarithm of the market price
List of Symbols and Notation 225
Part III
A,B,C factor loadings
B (0, t) money market account
Cov [... 1... ] conditional covariance operator
dt infinitesimal time increment
dWt infinitesimal Brownian increment
1& [... 1... ] conditional expectation operator
f (t, T) forward interest rate
Ft information set up to time t
k (t, T) swap rate
KU kurtosis
L (t, T) LIBOR rate
A market price of risk
P (t, T) price of a discount bond
lP empirical probability measure
r instantaneous interest rate
p coefficient of correlation
Q risk neutral probability measure
QT forward martingale measure
SK skewness
STD standard deviation
Var [... 1... ] conditional variance operator
x vector of theoretical state variables
~t vector of empirical state variables
Yt vector of measurable observations
y (t, T) spot interest rate
r7t forward price of a discount bond
Part IV
Cov[···I···] conditional covariance operator
dt infinitesimal time increment
dWt infinitesimal Brownian increment
1&[···1···] conditional expectation operator
F (t, T) forward or futures price
Ft information set up to time t
K, speed of mean-reversion
KU kurtosis
226 List of Symbols and Notation
Ball, C. A., and W. N. Torous (1996): "Unit Roots and the Estimation of
Interest Rate Dynamics," Journal of Empirical Finance, 3(2), 215-238.
Black, F., and M. Scholes (1973): "The Pricing of Options and Corporate
Liabilities," Journal of Political Economy, 81, 637-654.
Bliss, R R, and E.!. Ronn (1998): "Callable U.S. Treasury Bonds: Optimal
Calls, Anomalies, and Implied Volatilities," Journal of Business, 71(2),
211-251.
Brockwell, P. A., and R. A. Davies (1987): Time Series: Theory and Models.
New York: Springer.
Broze, L., O. Scaillet, and J.-M. Zakoian (1995): "Testing for Continuous-
Time Models of the Short-Term Interest Rate," Journal of Empirical
Finance, 2, 199-223.
BUhler, W., O. Korn, and R. Schobel (2000): "Pricing and Hedging of Oil
Futures - A Unifying Approach," Discussion paper, University of Tubin-
gen.
Chen, R-R (1996b): Understanding and Managing Interest Rate Risk. Sin-
gapore: World Scientific.
Chen, R-R, and L. Scott (1992): "Pricing Interest Rate Options in a Two-
Factor Cox-Ingersoll-Ross Model of the Term Structure," Review of Fi-
nancial Studies, 5(4), 613-636.
Cox, J. C., J. E. Ingersoll, and S. A. Ross (1985b): "A Theory of the Term
Structure of Interest Rates," Econometrica, 53(2), 385-407.
Duffie, D., and R. Kan (1996): "A Yield-Factor Model of Interest Rates,"
Mathematical Finance, 6(4), 379-406.
Duffie, D., and K. J. Singleton (1997): "An Econometric Model of the Term
Structure ofInterest-Rate Swap Yields," Journal of Finance, 52(4), 1287-
1321.
238 Bibliography
Eydeland, A., and H. Geman (1998): "Pricing Power Derivatives," Risk, 10,
71-73.
Gibbons, M. R., and K. Ramaswamy (1993): "A Test of the Cox, Ingersoll,
and Ross Model of the Term Structure," Review of Financial Studies,
6(3), 619-658.
Gourieroux, C., and A. Monfort (1997): Time Series and Dynamic Models.
Cambridge: Cambridge University Press.
Heath, D., R. Jarrow, and A. Morton (1992): "Bond Pricing and the Term
Structure of Interest Rates: A New Methodology for Contingent Claims
Valuation," Econometrica, 60(1), 77-105.
240 Bibliography
Ho, T. S., and S.-B. Lee (1986): "Term Structure Movements and Pricing
Interest Rate Contingent Claims," Journal of Finance, 41(5), 1011-1029.
Hull, J., and A. White (1987): "The Pricing of Options on Assets with
Stochastic Volatilities," The Journal of Finance, 42(2), 281-300.
Love, D. A., and S. F. McShea (1997): "The U.S. Treasury STRlPS: Valua-
tion, Liability, and Asset Allocation Frontier," in Yield Curve Dynamics,
ed. by R J. Ryan, pp. 9-32. Chicago: Glenlake.
Marshall, J. F., and K. R Kapner (1993): The Swaps Market. Miami: Kolb,
2 edn.
McCulloch, J. H., and H.-C. Kwon (1993): "U.S. Term Structure Data.
1947-1991," Working paper, Ohio State University.
Moler, C., and C. van Loan (1978): "Nineteen Dubious Ways to Compute
the Exponential of a Matrix," SIAM Review, 20(4), 801-836.
Pirrong, C., and M. Jermakyan (1999): "The Price of Power: The Valua-
tion of Power and Weather Derivatives," Working paper, Olin School of
Business, Washington University.
Uhrig, M., and U. Walter (1996): "A New Numerical Approach for Fitting
the Initial Yield Curve," Journal of Fixed Income, (3), 82-90.
Wu, X. (1996): "A New Stochastic Duration Measure by the Vasicek and
CIR Term Structure Theories," Working paper, City University of Hong
Kong.
Vol. 405: S. Koml6si, T. Rapcsak, S. Schaible (Eds.), Vol. 431: A. Aulin, Causal and Stochastic Elements in Busi-
Generalized Convexity. Proceedings, 1992. VIII, 404 pages. ness Cycles. XI, 116 pages. 1996.
1994. Vol. 432: M. Tamiz (Ed.), Multi-Objective Programming
Vol. 406: N. M. Hung, N. V. Quyen, Dynamic Timing and Goal Programming. VI, 359 pages. 1996.
Decisions Under Uncertainty. X, 194 pages. 1994. Vol. 433: J. Menon, Exchange Rates and Prices. XIV, 313
Vol. 407: M. Ooms, Empirical Vector Autoregressive pages. 1996.
Modeling. XIII, 380 pages. 1994. Vol. 434: M. W. J. Blok, Dynamic Models ofthe Firm. VII,
Vol. 408: K. Haase, Lotsizing and Scheduling for Production 193 pages. 1996.
Planning. vm, 118 pages. 1994. Vol. 435: L. Chen, Interest Rate Dynamics, Derivatives
Vol. 409: A. Sprecher, Resource-Constrained Project Pricing, and Risk Management. XII, 149 pages. 1996.
Scheduling. XII, 142 pages. 1994. Vol. 436: M. Klemisch-Ahlert, Bargaining in Economic and
Vol. 4\0: R. Winkelmann, Count Data Models. XI, 213 Ethical Environments. IX, 155 pages. 1996.
pages. 1994. Vol. 437: C. Jordan, Batching and Scheduling. IX, 178 pages.
Vol. 411: S. Dauz~re-Pl!res, J.-B. Lasserre, An Integrated 1996.
Approach in Production Planning and Scheduling. XVI, 137 Vol. 438: A. Villar, General Equilibrium with Increasing
pages. 1994. Returns. XIII, 164 pages. 1996.
Vol. 412: B. Kuon, Two-Person Bargaining Experiments Vol. 439: M. Zenner, Learning to Become Rational. VII,
with Incomplete Information. IX, 293 pages. 1994. 201 pages. 1996.
Vol. 4\3: R. Fiorito (Ed.), Inventory, Business Cycles and Vol. 440: W. Ryll, Litigation and Settlement in a Game with
Monetary Transmission. VI, 287 pages. 1994. Incomplete Information. VITI, 174 pages. J996.
Vol. 414: Y. Crama, A. Oerlemans, F. Spieksma, Production Vol. 441: H. Dawid, Adaptive Learning by Genetic
Planning in Automated Manufacturing. X, 210 pages. 1994. Algorithms. IX, 166 pages. 1996.
Vol. 415: P. C. Nicola, Imperfect General Equilibrium. XI, Vol. 442: L. Corch6n, Theories ofImperfectly Competitive
167 pages. 1994. Markets. XIII, 163 pages. 1996.
Vol. 416: H. S. J. Cesar, Control and Game Models of the Vol. 443: G. Lang, On Overlapping Generations Models with
Greenhouse Effect. XI, 225 pages. 1994. Productive Capital. X, 98 pages. 1996.
Vol. 417: B. Ran, D. E. Boyce, Dynamic Urban Transpor- Vol. 444: S. J~rgensen, G. Zaccour (Eds.), Dynamic
tation Network Models. XV, 391 pages. 1994. Competitive Analysis in Marketing. X, 285 pages. 1996.
Vol. 418: P. Bogetoft, Non-Cooperative Planning Theory. Vol. 445: A. H. Christer, S. Osaki, L. C. Thomas (Eds.),
XI, 309 pages. 1994. Stochastic Modelling in Innovative Manufactoring. X, 361
Vol. 419: T. Maruyama, W. Takahashi (Eds.), Nonlinear and pages. 1997.
Convex Analysis in Economic Theory. VIII, 306 pages. 1995. Vol. 446: G. Dhaene, Encompassing. X, 160 pages. 1997.
Vol. 420: M. Peeters, Time-To-Build. Interrelated Invest- Vol. 447: A. Artale, Rings in Auctions. X, 172 pages. 1997.
ment and Labour Demand Modelling. With Applications to
Six OECD Countries. IX, 204 pages. 1995. Vol. 448: G. Fandel, T. Gal (Eds.), Multiple Criteria Decision
Making. XII, 678 pages. 1997.
Vol. 421: C. Dang, Triangulations and Simplicial Methods.
IX, 196 pages. 1995. Vol. 449: F. Fang, M. Sanglier (Eds.), Complexity and Self-
Organization in Social and Economic Systems. IX, 317
Vol. 422: D. S. Bridges, G. B. Mehta, Representations of pages, 1997.
Preference Orderings. X, 165 pages. 1995.
Vol. 450: P. M. Pardalos, D. W. Hearn, W. W. Hager, (Eds.),
Vol. 423: K. Marti, P. Kall (Eds.), Stochastic Programming. Network Optimization. VIII, 485 pages, 1997.
Numerical Techniques and Engineering Applications. VIII,
351 pages. 1995. Vol. 451: M. Salge, Rational Bubbles. Theoretical Basis,
Economic Relevance, and Empirical Evidence with a Special
Vol. 424: G. A. Heuer, U. Leopold-Wildburger, Silverman's Emphasis on the German Stock Market.JX, 265 pages. 1997.
Game. X, 283 pages. 1995.
Vol. 452: P. Gritzmann, R. Horst, E. Sachs, R. Tichatschke
Vol. 425: J. Kohlas, P.-A. Monney, A Mathematical Theory (Eds.), Recent Advances in Optimization. VIII, 379 pages.
of Hints. XIII, 419 pages, 1995. 1997.
Vol. 426: B. Finkenstadt, Nonlinear Dynamics in Eco- Vol. 453: A. S. Tangian, J. Gruber (Eds.), Constructing
nomics.IX, 156 pages. 1995. Scalar-Valued Objective Functions. VTII, 298 pages. 1997.
Vol. 427: F. W. van Tongeren, Microsimulation Modelling Vol. 454: H.-M. Krolzig, Markov-Switching Vector Auto-
of the Corporate Firm. XVII, 275 pages. 1995. regressions. XIV, 358 pages. 1997.
Vol. 428: A. A. Powell, Ch. W. Murphy, Inside a Modem Vol. 455: R. Caballero, F. Ruiz, R. E. Steuer (Eds.), Advances
Macroeconometric Model. XVIII, 424 pages. 1995. in Multiple Objective and Goal Programming. VIII, 391
Vol. 429: R. Durier, C. Michelot, Recent Developments in pages. 1997.
Optimization. VIII, 356 pages. 1995. Vol. 456: R. Conte, R. Hegselmann, P. Terna (Eds.), Simu-
Vol. 430: J. R. Daduna,l. Branco, J. M. Pinto Paixilo (Eds.), lating Social Phenomena. VllI, 536 pages. 1997.
Computer-Aided Transit Scheduling. XIV, 374 pages. 1995.
Vol. 457: C. Hsu, Volume and the Nonlinear Dynamics of Vol. 483: M. Konigstein, Equity, Efficiency and Evolutionary
Stock Returns. VIII, 133 pages. 1998. Stability in Bargaining Games with Joint Production. XII,
Vol. 458: K. Marti, P. Kall (Eds.), Stochastic Programming 197 pages. 2000.
Methods and Technical Applications. X, 437 pages. 1998. Vol. 484: D. D. Gatti, M. Gallegati, A. Kirman, Interaction
Vol. 459: H. K. Ryu, D. J. Slottje, Measuring Trends in U.S. and Market Structure. VI, 298 pages. 2000.
Income Inequality. XI, 195 pages. 1998. Vol. 485: A. Garnaev, Search Games and Other Applications
Vol. 460: B. Fleischmann, J. A. E. E. van Nunen, M. G. of Game Theory. VIII, 145 pages. 2000.
Speranza, P. Stahly, Advances in Distribution Logistic. XI, Vol. 486: M. Neugart, Nonlinear Labor Market Dynamics.
535 pages. 1998. X, 175 pages. 2000.
Vol. 461: U. Schmidt, Axiomatic Utility Theory under Risk. Vol. 487: Y. Y. Haimes, R. E. Steuer (Eds.), Research and
XV, 201 pages. 1998. Practice in Multiple Criteria Decision Making. XVII,
Vol. 462: L. von Auer, Dynamic Preferences, Choice 553 pages. 2000.
Mechanisms, and Welfare. XII, 226 pages. 1998. Vol. 488: B. Schmolck, Ommitted Variable Tests and
Vol. 463: G. Abraham-Frois (Ed.), Non-Linear Dynamics Dynamic Specification. X, 144 pages. 2000.
and Endogenous Cycles. VI, 204 pages. 1998. Vol. 489: T. Steger, Transitional Dynamics and Economic
Vol. 464: A. Aulin, The Impact of Science on Economic Growth in Developing Countries. VIII, 151 pages. 2000.
Growth and its Cycles. IX, 204 pages. 1998. Vol. 490: S. Minner, Strategic Safety Stocks in Supply
Vol. 465: T. J. Stewart, R. C. van den Honert (Eds.), Trends Chains. XI, 214 pages. 2000.
in Multicriteria Decision Making. X, 448 pages. 1998. Vol. 491: M. Ehrgott, Multicriteria Optimization. VIII, 242
Vol. 466: A. Sadrieh, The Alternating Double Auction pages. 2000.
Market. VII, 350 pages. 1998. Vol. 492: T. Phan Huy, Constraint Propagation in Flexible
Vol. 467: H. Hennig-Schmidt, Bargaining in a Video Ex- MlIIUIfac(uring. IX, 258. pages. 2000.
periment. Determinants of Boundedly Rational Behavillr. Vo.1A93: .J. Zhu, Modular Pricing of Options. X, 170 pages.
XII, 221 pages. 1999. 2000.
Vol. 468: A. Ziegler, A Game Theory Analysis of Options. Vol. 494: D. "Franzen, Design of Master Agreements for OTC
XIV, 145 pages. 1999. Derivatives. VIII, 175 pages. 2001.
Vol. 469: M. P. Vogel, Environmental Kuznets Curves. xm, VOl. 495: I Konnov, Combined Relaxation Methods for
197 pages. 1999. Variational Inequalities. Xl, 181 pages. 200 I.
Vol. 470: M. Ammann, Pricing Derivative Credit Risk. XII, Vol. 496: P. WeiS, Unemployment in Open Economies. XII,
228 pages. 1999. 226 pages. 200 I.
Vol. 471: N. H. M. Wilson (Ed.), Computer-Aided Transit Vol. 497: J. Inkmann, Conditional Moment Estimation of
Scheduling. Xl, 444 pages. 1999. Nonlinear Equation Systems. VIII, 214 pages. 2001.
Vol. 472: J.-R. Tyran, Money Illusion and Strategic Vol. 498: M. Reutter, A Macroeconomic Model of West
Complementarity as Causes of Monetary Non-Neutrality. German Unemployment. X, 125 pages. 2001.
X, 228 pages. 1999. Vol. 499: A. Casajus, Focal Points in Framed Games. Xl,
Vol. 473: S. Helber, Performance Analysis of Flow Lines 131 pages. 200 I.
with Non-Linear Flow of Material. IX, 280 pages. 1999. Vol. 500: F. Nardini, Technical Progress and Economic
Vol. 474: U. Schwalbe, The Core of Economies with Growth. XVII, 191 pages. 2001.
Asymmetric Information. IX, 141 pages. 1999. Vol. 501: M. Fleischmann, Quantitative Models for Reverse
Vol. 475: L. Kaas, Dynamic Macroeconomics with Imperfect Logistics. XI, 181 pages. 2001.
Competition. XI, 155 pages. 1999. Vol. 502: N. Hadjisavvas, J. E. Martinez-Legaz, J.-P. Penot
Vol. 476: R. Demel, Fiscal Policy, Public Debt and the (Eds.), Generalized Convexity and Generalized Mono-
Term Structure of Interest Rates. X, 279 pages. 1999. tonicity. IX, 410 pages. 2001.
Vol. 477: M. Th~ra, R. Tichatschke (Eds.), III-posed Vol. 503: A. Kirman, J.-B. Zimmermann (Eds.), Economics
Variational Problems and Regularization Techniques. VIII, with Heterogenous,lnteracting Agents. VII, 343 pages. 200 I.
274 pages. 1999. Vol. 504: P.-Y. Moix (Ed.),The Measurement of Market Risk.
Vol. 478: S. Hartmann, Project Scheduling under Limited Xl, 272 pages. 200 I.
Resources. XII, 221 pages. 1999. Vol. 505: S. VoS, J. R. Daduna (Eds.), Computer-Aided
Vol. 479: L. v. Thadden, Money, Inflation, and Capital Scheduling of Public Transport. XI, 466 pages. 2001.
Formation. IX, 192 pages. 1999. Vol. 506: B. P. Kellerhals, Financial Pricing Models in
Vol. 480: M. Grazia Speranza, P. Stlihly (Eds.), New Trends Continuous Time and Kalman Filtering. XIV, 247 pages.
in Distribution Logistics. X, 336 pages. 1999. 2001.
Vol. 481: V. H. Nguyen, J. J. Strodiot, P. Tossings (Eds.).
Optimation. IX, 498 pages. 2000.
Vol. 482: W. B. Zhang, A Theory of International Trade.
XI, 192 pages. 2000.