An Econometric Model of the Term Structure of InterestRate Swap Yields
by
Kenneth J. Singleton, Darrell Duffie
Citation
Title:
An Econometric Model of the Term Structure of InterestRate Swap Yields
Author:
Kenneth J. Singleton, Darrell Duffie
Year:
1997
Publication:
The Journal of Finance
Volume:
52
Issue:
4
Start Page:
1287
End Page:
1321
Publisher:
Language:
English
URL:
Select license:
Select License
DOI:
PMID:
ISSN:
Updated: February 4th, 2013
Abstract:
An Econometric Model of the Term Structure of
InterestRate Swap Yields
DARRELL DUFFIE and KENNETH J. SINGLETON*
ABSTRACT
This article develops a multifactor econometric model of the term structure of interestrate swap yields. The model accommodates the possibility of counterparty default, and any differences in the liquidities of the Treasury and Swap markets. By parameterizing a model of swap rates directly, we are able to compute modelbased estimates of the defaultable zerocoupon bond rates implicit in the swap market without having to specify a priori the dependence of these rates on default hazard or recovery rates. The time series analysis of spreads between zerocoupon swap and treasury yields reveals that both credit and liquidity factors were important sources of variation in swap spreads over the past decade.
ALTHOUGHPLAIN VANILLA interestrate swaps comprise a major
FIXEDforfloatingsegment of the fixedincome derivative market, notably few econometric models for pricing swaps have been developed in the literature. Perhaps the primary reasons for this are: (i) swap contracts embody default risk and hence equilibrium or arbitragefree term structure models developed for defaultfree government bond markets are not directly applicable to the swap market; (ii) empirical modeling of the default event underlying credit spreads on defaultable bonds and swaps has met with limited success at explaining the timeseries properties of spreads; and (iii) swap spreads are likely to depend on other factors such as liquidity that are not directly related to default events. Also, until recently, data have not been widely available. In this article we develop a multifactor econometric model of the term structure of U.S. fixedforfloating interestrate swap yields that accommodates many of the institutional features of swap markets. Specifically, using results in Duffie and Singleton (1996), we show that the fixed payment rate of a swap, assuming that the floating rate is London Interbank Offering Rate (LIBOR), can be expressed in terms of present values of net cash flows of the swap contract
" Graduate School of Business, Stanford University. We are grateful for the research assistance of Qiang Dai, Stephen Gray, and Raj Tewari, and comments from Mark Fisher, Dilip Madan, Ming Huang, Rene Stulz, an anonymous referee, and seminar participants at the National Bureau of Economic Research, the University of Wisconsin Finance Symposium, the University of California, Berkeley, the University of Chicago, Duke University, the University of Arizona, and the University of California at San Diego, and from an anonymous referee. Data were kindly provided by Goldman Sachs and Co. Financial support was provided by the Stanford GSB Financial Research Initiative. These results appeared in preliminary form under the title "Econometric Modeling of Term Structures of Defaultable Bonds." An extended version of the valuation models from that article now appears in Duffie and Singleton (1996).
discounted by a default and liquidityadjusted instantaneous short rate. In other words, there is an adjusted short rate process that allows us to develop a term structure model for the swap market in much the same way that models have been developed for government yield curves. Default and liquidity risks are "collapsed" into a riskadjusted short rate for computing the present values of future risky cash flows.
Our formulation provides a modelbased alternative to swap valuation mod els used by many investment and commercial banks. Many financial institu tions use forward interest rates obtained by interpolating between fixed ma turity points on the swap yield curve to derive discount rates for future net cash payments between the swap counterparties. Instead, we develop an arbitragefree multifactor term structure model of swap yields that leads to an implied, modelbased zero curve that can be used in valuing offtherun swaps, caps, floors, and swaptions (options to enter into a swap contract at some future date). Moreover, unlike interpolation schemes, our model provides a framework for dynamically hedging swaps or derivatives, such as caps and swaptions.
The literature on swap rates has typically focused on the spreads between swap rates and the corresponding point on the U.S. Treasury (defaultfree) term structure. Sun, Sundaresan, and Wang (1993) examine the average swap spreads to Treasuries and the bidlask spreads by the credit class of the swap counterparties. Brown, Harlow, and Smith (1994) regress swap spreads on various contemporaneous measures of credit risk and the hedging costs of market makers. They find that both sets of variables are correlated with swap spreads. Neither study develops a dynamic swappricing model.
Grinblatt (1995) argues that liquidity risk is a more plausible explanation for swap spreads than credit risk. Liquidity enters his model as a convenience yield to Treasuries associated with their relative liquidity and potential to go "on special" in markets for repurchase agreements (repo). Holders of Treasury bonds that go on special can effectively receive a special dividend by borrowing at below market rates, using Treasuries as collateral. (See also Duffie (1996) and Jordan and Jordan (1997).) Using the onemonth LIBORTreasury bill spread as a proxy for the convenience yield of holding government securities, Grinblatt calibrates a Vasicek (1977) representation of riskless and conve nienceyield processes. He finds that this model explains about 35 to 40 percent of the variation in swap spreads for maturities of two through ten years.
Instead of focusing on swap spreads to the defaultfree term structure, we focus on swap yields directly. Our goal is to develop a model of the swap market without needing to be precise a priori about the economic mechanisms that generate swap spreads. Subsequent to fitting the model, we can study the properties of the defaultable zerocoupon yields implied by the swap market with the goal of a better understanding of the economic factors that determine their spreads to Treasury zerocoupon yields. In Section I we show that, under the assumption that the counterparties have symmetric probabilities of de fault, a swap is "priced" by the present value of its cash flows discounted by a risk and liquidityadjusted shortrate process. Once one adopts a parameter ization of this shortrate process, the parameters of the model can be estimated without having to specify a priori the functional forms for default probabilities, liquidity premiums, and so forth. Virtually any of the models examined pre viously for government yield curves can be used to model the defaultadjusted shortrate process, including affine processes for riskless rates and Heath JarrowMorton (1992) type models of forward rates.1
In Section 11, we discuss the econometric models of defaultable swap yields. The econometric model of swaps studied in Section I11 presumes that the liquidity and defaultadjusted shortrate process is the sum of two indepen dent squareroot diffusions (a 2factor model). Assuming that the model ex actly prices swaps at two points along the swap yield curve, and using the fact that the conditional distribution of discretely sampled data from a squareroot diffusion has a noncentral chisquare density (e.g., Cox, Ingersoll, and Ross (1985) (CIR)), we estimate our model using the joint likelihood of swap yields at several maturities.2 The distributions of the fitted swap yields, evaluated at the maximum likelihood (ML) estimates, have sample moments that are similar to those of the corresponding actual yields. Moreover, deviations be tween the actual and fitted swap yields are on average zero, with standard deviations between four and seven basis points for seven years of weekly data. (By comparison, the bidlask spreads in the swap market averaged about four basis points during our sample period.)
In light of the close fit of the model to swap yields, we proceed to compute the implied risky zerocoupon bond yields, evaluated at the ML estimates, and compute their spreads to the corresponding U.S. Treasury zerocoupon yields. In Section IV the properties of these swap zerocoupon yield spreads are studied in the context of a multivariate vector autoregression WAR) in an attempt to shed some light on the relative importance of liquidity and credit factors in the determination of swap spreads. Included in the VAR, along with swap zero spreads, are variables proxying for corporate credit risk and liquid ity in the U.S. Treasury market. The results suggest that both liquidity and credit factors affect the temporal behavior of swap zero spreads, but that the responses of swap spreads to changes in these factors follow very different time paths. Liquidity effects are shortlived, whereas responses to credit shocks are weak initially and then increase in importance over a horizon of several months. Concluding remarks are presented in Section V.
I. Valuation of Swaps
Consider a set of M plain vanilla fixedforfloating swaps. The mth swap has rm years to maturity. The floating side is reset semiannually to the sixmonth LIBOR rate from six months prior. The fixed side pays a coupon cmat the reset
See Duffie and Singleton (1996) for a discussion of alternative formulations of risky discount rate processes for valuing defaultable bonds. This likelihood function has the same form as that studied by Chen and Scott (1993) and Pearson and Sun (1994) in their studies of the U.S. Treasury market.
dates. Let rf denote the LIBOR rate set at date t for loans maturing six months in the future and PV(t, t + 7,) denote the present value of the promised net cash flows between the parties to the swap agreement. We assume that, at the inception date of the swap,
where Eg denotes expectation under an equivalent martingale measure for the information sets (9,: t 10) commonly available to investors, and where R is an instantaneous discountrate process defined below. In other words, the fixedside coupon c," of the swap is set at date t so that the present value of the net cash flows exchanged by the counterparties at the reset dates is zero at the inception of the swap, as in equation (1).
We let
be the discount factor at time t for maturity 7 associated with the short term discountrate process R. We further assume that risky zerocoupon bonds are priced at the appropriate LIBOR rate in the interbank lending market, or
Using equations (1)to (3), the present value of the floating rate payments is 1B:m. Combining these observations, the coupon rate c," on the fixed side of the swap can be expressed as
In deriving the valuation model (1)(4), we have implicitly made several important assumptions. Since the assumptions needed for equation (4) to hold are significantly weaker than those typically made (see, for example, Litzen berger (1992) for a review of the literature), we briefly examine in more depth the key underlying assumptions.
The discounting in equations (1)and (2) takes the same form as the dis counting of defaultfree cash flows by the riskless, instantaneous interest rate in riskneutral representations of standard termstructure models. However, models for pricing defaultfree cash flows are not directly applicable to swaps and LIBOR contracts, because they are defaultable instruments. Nevertheless, using results in Duffie and Singleton (1996), we can interpret R as a default adjusted discount rate and, under this interpretation plus additional assump tions outlined below, equation (1)correctly prices a defaultable swap.
More precisely, we define a defaultable claim to be a pair ((X, T),(X', T')) of contingent claims. The underlying claim (X, T)is the obligation of the issuer. The secondary claim (X',T') defines the stopping time T' at which the issuer defaults, and the payment X' to be received at default. This means that the actual claim (2, TI generated by a defaultable claim ((X, T), (X',T')) is defined
by X ifT<Tf, 7 = min(T, T'); Z = {X' ifT?Tf,
where 7 is the stopping time at which Z is actually paid. The exdividend price process V of any given contingent claim (2, r) is defined by V, = 0 for t 2 T and
Evaluation of the pricing formula (6) is complicated in practice by the possi bility of default reflected in the payoff2 in equation (5) and that the probability of default will in general be correlated with the shortrate process, r. To circumvent these difficulties, we follow Duffie and Singleton (1996), and as sume that the secondary claim (X', T') is defined by a hazard rate process h and a fractional default loss process A.The hazard rate h can be thought of as the arrival intensity of a Poisson process whose first jump occurs at default. The state dependent process h defines the fraction of market value of the claim that is lost upon default.
Under mild technical regularity conditions given in Duffie and Singleton (1996), the valuation of this defaultable claim can proceed as if the promised payoff X is defaultfree; however, with discounting at a defaultadjusted dis count rate R instead of the riskless rate. Consequently, under those conditions,
where
We interpret each promised cash flow of a swap contract as one of the promised payments X described above. By repeatedly applying the same logic used to derive equation (7) to the sequence of promised cash flows of a swap, we get equation (1). This approach to swap valuation avoids a more complicated structural model based on knowledge of the assetliability structure of the swap counterparties, in the style of Merton (1974), and as applied to the swaps
market by Rendleman (1992). The development of a structural model seems impractical given our data on generic market swap rates. ~e~ond
default, a factor that may affect swap spreads to Treasuries is the relative liquidities of the two markets (Grinblatt (1995), Bansal and Coleman (1996)). Therefore, we also include a convenience yield I, that allows for the effect of differences in liquidity and rep0 specialness between the Treasury and swap markets. With this modification, the swap rate cm at inception solves equation (1)with the adjusted discount rate,
The processes r, I, h, and h are assumed to be adapted to the investors' information sets (9,: t 2 0) and thus, may be state dependent and mutually correlated.3
For the case of swaps, modeling the default time as an inaccessible stopping time, such as a Poisson arrival, seems reasonable because default events, when they do occur, are rarely fully anticipated even a short time before the default. Changing expectations concerning the likelihood of default are captured by the stochastic properties of the hazard rate process h. Indeed, under the risk neutral measure &, the conditional probability at time t of default over the next "instant" of time of length At is approximately h,At.
The cash payment upon default of a swap contract is based on the market value of the remaining obligations under the terms of the swap and negotia tions between the counterparties. If PV(td, td + 7) denotes the value of a 7year swap just prior to the default time td, then (1htd)PV(td, td + 7) is the present value at td of the cash flows ultimately generated by the negotiated settlement. The liability, if any, according to standard International Swap Dealers Asso ciation (ISDA) swap contracts, is based on midmarket quotes for similar swaps.
One implication of these observations is that the assumption of no default made by Smith, Smithson, and Wakeman (1988) and Sun, Sundaresan, and Wang (1993), among others, in justifying expressions like equation (4) is unnecessarily strong. The market yields on Arated LIBOR issues are occa sionally at large spreads to U.S. Treasuries, and these spreads fluctuate substantially over time, as we show in Section 111. The pricing relation (4) implicitly allows this spread to arise from a spread between R and the riskless rate r.
There are still, however, some important implicit assumptions remaining in the valuation model (1)(4).
A. Exogenous Default Risk
The valuation model (l), for a given "credit spread" process Ah, implicitly presumes that revaluation of the swap is not in itself a significant determinant
These processes are also assumed to be jointly measurable in state and time and to satisfy mild integrability conditions described in Duffie and Singleton (1996).
of default likelihood or losses on default. (This endogeneity would exist, for example, with a swap that constitutes the major part of the liabilities of one of the counterparties.) If we were to compensate by allowing Ah to depend endogenously on the market value of the swap, the linear valuation model (1) would be replaced with a nonlinear model described by Duffie and Singleton (1996). This would significantly complicate the econometric model. The fact that we are using generic market quotes for swap rates, and do not depend on the valuation of a swap between two particular counterparties, presumably mitigates the impact of any endogeneity of the default spread A,h,. Moreover, based on the numerical results noted below, the impact of moderate depen dence of Ah on the market value of the swap is minimal.
B. "Refreshed" AQuality Counterparties
Discounting the net cash flows at all maturities with the same risky discount rate R in equation (1)presumes that the counterparties maintain the credit rating underlying generic (say, Arated) swaps for the life of the swap contract. More precisely, at the inception of the swap, the counterparties are presumed to have a credit rating of A. Subsequently, there is the possibility that either counterparty could be upgraded or downgraded (or suffer a change in liquidity) during the life of the swap. This itself is not a problem for pricing in terms of equation (I), since the potential of a change in quality can be captured in the stochastic process for the discount rate R. However, we do not have data on the yields for specific swaps as they mature, but rather on new Arated swaps with constant maturity. Thus, we are valuing a hypothetical swap priced by dealers who presume that the counterparties will maintain the quality of newly issued Arated debt over the life of the swap. This problem is not unique to our model of swaps, but rather is inherent in any model of newissue rates on defaultable debt from a fixed credit class. Based on the numerical calculations of Duffie and Huang (1996) and Li (19951, described below, there is a minimal impact on swap rates of reasonably anticipated variation in the credit quality of the counterparties over the life of the swap.4
C. Symmetric CreditQuality
Related to the point above is the possibility that the two counterparties to a given swap may have different credit quality, either at the inception of the swap, or subsequently. Extending our model to allow for asymmetric credit qualities of the counterparties would add substantial complexity to the pricing model. For example, even if the creditspreads AAhA of counterparty A and ABhBof counterparty B are each given exogenously, and do not depend on the market value of the swap, in keeping with the exogenous default risk assump tion above, there is nevertheless the effect of an endogenous dependence of the
For alternative models of the impact of default risk on swap rates, see Abken (19931, Cooper and Mello (1991), Hull and White (1992), Jarrow and Turnbull (19951, Li (1995), Rendleman (1992), Solnik (1990), Sorensen and Bollier (19941, and Sundaresan (1991).
credit spread on market value if hAhA f h,h,. This arises from the fact that the current market value of the swap determines the firm whose default risk currently "matters," based on netting provisions and the implications of one way or twoway ("no fault") payment schemes that are usually adopted in standard swap agreements. (If the market value of the swap to counterparty A is positive at a given point in time, then it is the default risk of counterparty B that "matters" at that moment.) Therefore, the relevant credit spread would become dependent on the market value of the swap, and the nonlinear valu ation model discussed in Duffie and Singleton (1996) would apply. Duffie and Huang (1996) and Li (1995) explore theoretically and numerically this aspect of swap contracts, and find that the degree of asymmetry in credit quality is a relatively minor determinant of swap rates for typical interestrate swaps. (For example, with less than a 10year maturity and under typical parameters, there is a correction of roughly 1basis point or less in swap yields for a credit risk asymmetry generating a 100 basis point difference in bond yields.) Thus, there is likely to be negligible misspecification error from proceeding under the assumption that swaps are priced as if there is symmetric counterparty risk.
D. Homogeneous LIBORSwapMarket Credit Quality
There is no reason in theory that LIBOR corporate bond rates, including the floatingrate swap payments rk, should be determined by discounting at the same liquidity and credit adjusted shortrate R used for discounting net swap payments in equation (1). The default scenarios may be different in the two markets, recovery rates may differ, and the liquidities of the two markets are also typically different. That is, it is a nontrivial assumption that equation (3) holds.5
Although these remarks have focused on the default component of swap spreads, there is no presumption in equation (9) that default is the only or even the primary determinant of swap spreads. Empirically, I,, which represents other factors that determine effective carrying costs such as liquidity, may be as (or more) important a source of variation in R, r,.
11. Econometric Models of Defaultable Swap Yields
In this section we discuss two alternative formulations of the adjusted shortrate R that may lead to econometrically identified models of the term structure of swap rates. One approach is to focus directly on R and assume that R, = p(Y,), where Y is a Markov state vector determining the defaultadjusted shortrate. Under this approach, no attempt is made to distinguish between
There are many alternative variablerate payments streams which, in place of LIBOR floating rate payments, would justify equation (4). It is the fact that the LIBOR floating rate payments are made in arrears that forces us to make equation (3)as an assumption. Specifically, if equations (11, (21, and (4) are to hold at all maturities, and the floating rate payment at time t is known at time t 0.5, then one can show by induction, using the law of iterated expectations, that the floating rate payment at time t + 0.5 must be (BP5)11.
the contributions of the riskless rate r and the premium z, = 1, + htht. The second approach parameterizes both r and the mean loss rate z.~
A. Formulations of the Adjusted Shortrate R
Consider again a generic, defaultable, contingent claim with a promised payoff of X at maturity date T. Suppose that there is a Markov (under the equivalent martingale measure Q) statevariable process Y, such that the promised contingent claim is of the formX =g(YT), for some functiong, and the default and liquidity adjusted shortrate process R is of the form p(Y,), for some function p(.). Then equation (7) implies that the claim to payment of g(YT) at time T has a price at time t, assuming that the claim has not defaulted by time t, of
We emphasize that the riskless shortrate r, and mean loss rate z, = hth, 1, do not enter directly into this pricing model (lo), but rather enter implicitly through the defaultadjusted shortrate R, = p(Y,). For a given defaultable bond market, parameterizing R directly, as opposed to separate parameteriza tion of r and z, provides less information concerning the meanlossrate process
z. On the other hand, this formulation permits empirical characterizations of the adjusted shortrate process R without a need to commit to a formulation of the credit spread. One reason that this may be attractive is that there may be nondefault factors that determine spreads between bonds of various classes. Many of these nondefault factors can be accommodated in our model by appropriate reinterpretation of the adjusted shortrate process. This robust ness is an attractive feature of this econometric modeling strategy when one is most concerned with characterizing the distribution of R, say for the purpose of computing a zerocoupon yield curve implied by a defaultable yield curve. On the other hand, the robustness highlights the fact that nothing can be learned directly about the default processes h and h from this approach.
Pursuing this approach further, suppose Y, = (Y,,, . . . , Y,,)', for some n, solves a stochastic differential equation of the form
where B is a standard Brownian motion in Rn under Q, and where p and a are well behaved functions on Rn into Rn and RnXn;respectively. Then we know from the "FeynmanKac formula" that, under technical conditions (Friedman
Even if liquidity effects are absent, it is in general not possible to separately identify the components h, and A, ofz, using swap data (or corporate bond data) alone. See Duffie and Singleton (1996) for a more indepth discussion of this point and a demonstration that the hazard and recovery rates can be identified from price information on credit derivatives.
(1975), Krylov (1980)), equation (11)implies that V solves the backward Kolmogorov partial differential equation
Qi","V(y, t) + p(y)V(y, t) = 0, (y, t) E RnX [O, TI, (12)
with the boundary condition
where
The identification problem for this modeling strategy is identical to that of standard termstructure models for defaultfree yield curves. All of the models for the shortrate process r that have been successfully studied for defaultfree term structures are also identified, at least in principle, in the case of default able bonds, simply by replacing r with R. For instance, consider the special case Y, = (Y,,, . . . , Y,,)', where Y,, Y,, . . . , Yn are independent (under riskneutral probabilities) squareroot processes. That is, equation (11)applies with pi(y) = ~~(0~yi), aii(y) = cia,for positive constants K~,
Oi, and ifi, and aJy) = 0, for i # j. By taking p to be affine, one can then apply the CIR solution to equations (13) and (14) withg(y) = 1, allowing defaultable zerocoupon bond prices to be computed in closed form. This is the approach taken in Section 111.
The same valuation mode1 (12)(13) applies, under mild technical conditions given in Duffie and Singleton (1996), when the underlying statevariable process Y is a jumpdiffusion, or more general continuoustime Markov pro cess. One merely replaces %@,"with the infinitesimal generator of Y. Allowing for "jumps" in credit quality may be useful if one has in mind the potential for unusually large creditquality events.
B. Modeling the Mean Loss Rate Process
A second modeling strategy is to parameterize the behavior of the joint process (r,2)' for the defaultfree shortrate and the riskneutral meanloss rate processes, respectively.7 Perhaps the simplest example of this strategy would be a case in which one studies the joint distribution of the returns on a defaultable bond and the associated reference Treasury bond used in pricing the defaultable bond. To be concrete, suppose that w = (r, 2)' follows an "affine" diffusion
dw, = (a,+ cu,w,)dt + p(w,)dB,, (15)
In this section we abstract from the nondefault factors discussed in Section V. If one has observable proxies for liquidity, for example, then the following discussion is extendable to accommodate a liquidity premium in R.
where
B is a standard Brownian motion in LQ2 under an equivalent martingale
measure Q;
a, E LQ2;a1 is a 2 X 2 matrix;
For each i andj, there is a fixed constant yoij and fixed vector y,,. in LQ2 such
that [P(w)P(w)'li,.= yoij + ylij w; and,
a,, al, and p satisfy joint restrictions8 for the existence of solutions to
equation (15).
Let G: and C: denote the prices of defaultfree and defaultable zerocoupon bonds of maturity n, respectively. The identification of z as the instantaneous default premium comes from the simultaneous estimation of the implied pricing equations for G: and C:, and the imposition of identifying restrictions on a, and p(.). The latter are necessary because, as shown in Dai and Singleton (1996),there will in general be a family of observationally equivalent affine models in the absence of normalizations in a, and P.
As an illustrative example of how information about z can be inferred from data on G: and C:, consider the special case (15) in which
The defaultfree term structure is driven by CIRstyle factor r in equation (16). Thus, the parameters a,,, a,,, yo,,, and y,,, are identifiable from information on one or more government bond prices GF. The parameters of r and the default process z are therefore identifiable from information on defaultfree and defaultable bond prices (see Dai and Singleton (1996) and Duffie and Singleton (1996)). Note that the Brownian motions associated with r, and z, may be correlated, though a,, and v,,, must satisfy certain existence conditions.
Two special cases of equation (16)and (17) are: (i)R, = Br, + st,with E, and r, being independent diffusions, and (ii)either h, or A, is specified a priori at a fixed value with the other assumed to follow a given affine diffusion. In the first case, if the issuer is exposed to interest rate risk, then this could be captured by the assumption that B > 1. In the second case, the strategy of modeling h and h simplifies to this case of modeling z. One example of the latter simplification is the model of Nielsen and Ronn (1995),who assume that A, is fixed at either 0.0 or 0.5, the drift of z is independent of r, and the diffusion coefficient for z is constant over time. The affine example (17) illustrates how these special cases can easily be extended to allow the meanlossrate process z to depend on the riskless interest rate. This seems like a potentially important extension in light of the documented cyclicality in credit spreads.
These restrictions are basically that the state vector does not enter a region where the square root of a negative number would have to be taken, as for example with a negative drift at zero in the CoxIngersollRoss (1985) model. See Duffle and Kan (1996) for details.
More generally, credit spreads are related to businesscycle variables, as shown by Bernanke (1990), Friedman and Kuttner (1993), Jaffee (1975), and Stock and Watson (1989). Among the more prominent variables are general stock market returns, measures of output growth, panel data on consumer sentiment, and levels of capital investment. There may be value in pursuing empirical models in which the riskneutral expected loss rate process z is linked with such macroeconomic variables. The state vector in equation (15) could be expanded to introduce macro information, although the observation frequency of many macro series might limit the applications of such models.
These examples lead to closed or nearly closedform expressions for the defaultable zerocoupon bond prices in terms of r, and z,. However, there is no need to restrict attention to analytic solutions for zerocoupon prices. As long as it is computationally feasible to compute the zerocoupon bond prices nu merically at the same time that the objective function defining the estimator of the unknown parameters governing the process (r, z) is being optimized, our preceding comments continue to apply. Of course, the correlation structure of r, and z, must be such that the parameters are identified.
111. An Econometric Model of Swap Yields
In this section we implement the first estimation strategy under the assump tion that the default and liquidity adjusted shortrate process R is a linear combination of independent squareroot diffusion models. To be concrete, let Y, be a Markov state vector that determines the current riskadjusted shortrate
R. Let %(Y,, Po) denote a vector of M1 prices of defaultable zerocoupon bonds and %(Y,, p,) denote a vector of M, yields on newly issued swaps (M = MI + M,) implied by the termstructure model. The parameter vector governing the probability model of the state process Y is Po.Also, let E, denote an Mvector of measurement errors contaminating the observed counterparts of the zero prices and yields, B, and C,, that are independent of the state vector Y. Then the econometric model takes the form
where, for ease of notation, we combine the parameter vectors governing the probability laws of Y and E into Po. The special case of equation (18) that we will study is the multifactor square root model, with
where (W1, W2) is a standard Brownian motion in 9i2.The positive scalars K~, Oi,and aihave interpretations in terms of meanreversion, steadystate mean, and volatility, respectively, that have been developed by Feller (1951) and, in the context of termstructure models, by CIR. In our empirical analysis we take R, = Y: + Y:, in which case equation (2) and the independence of {Y1) and {Y2) under an equivalent martingale measure Q leave us with
where pi(y, T) = U~(T)~~Z(~)Y,
for
and
a(i) = 2~~0~lg,yL = ((K~ + hi)' + 24)'12, and hi denotes a riskpremium coefficient explained by CIR. That is, eachpi is the form of a zerocoupon bond price in a univariate squareroot diffusion model, and the price of a defaultable zerocoupon bond in our setting is the product of the univariate bond price formulas.
The state vector Y is unobservable. Therefore we proceed, as do Chen and Scott (1993), under the assumption that two elements of the yield vector (93,%) (corresponding to the number of state variables) are measured without error. The pricing model can then be inverted to express the state variables as functions of these swap and LIBOR yields. Alternatively, we could allow all of the swap and LIBOR yields to be measured with error and then estimate the model using the simulated method of moments (Duffie and Singleton (1993)). While this strategy would permit us to obtain consistent estimates of the parameters, the state variables would be unknown. In principle, the latent state variables could be estimated using filtering methods. However, the approach taken here has several potential advantages for pricing, including having direct observations of the state variables and forcing the model to fit a subset of the swap yields exactly. These are important considerations for valuing derivative claims based on swap or LIBOR yields. The model we examine assumes that the two and tenyear swap rates satisfy
for T = 2, 10, with no measurement error^.^ Since equation (23) is a nonlinear function of the state variables, inferring (Y:, Y:) from (c:, c:') must be done numerically for each observation.
We also examine a model in which the price of a sixmonth LIBOR instrument satisfies
lo year yield
6

2 I
z
2 year110 year slope I'TL
,I""'\/"" \/w,v
2: 'A <\
1,.. cr*
'*c,,
, \ .
,">*&
1
.\\'.\* /*\,*'J,*~
i d
'% 1O \I
2 , I 1 uu.p 8 8 8 9 90 9 1 92 9 3 94
Date
Figure 1. Level and Slope of Swap Yield Curve, January 8,1988 to October 28,1994.
We make one additional modification to the model. The squareroot param eterization precludes negative state variables. As can be seen from Figure 1, the swap data exhibit substantial changes in the level of rates (as measured by the 10year yield) and in the magnitude and sign of the slope of the swap curve (as measured by the 2to10year yield spread). We are unable to find admis sible parameters for equations (19) and (20) for which the implied state variables Y1and Y2 are positive for the entire sample period. Therefore, we modify the discount process to be
where is a positive constant.10 With this modification, the zerocoupon bond price defined by equation (20) is replaced with
and in which equation (23) holds for T = 10. We will comment briefly in the next section on the comparative fit of this alternative model.
While the inclusion of? solves the problem of negative state variables, the adjusted discount rate R may be negative. In this respect, the modified model differs from the original squareroot diffusion specification. Pearson and Sun (1994) study the same specification of a twofactor model of the riskless rate r underlying U.S. Treasury yields (their "extended model").
The conditional densities of the state variables {Y1, y2) are well known to be noncentral chisquare (e.g., CIR). Through a change of variables, these con ditional densities can be rewritten as functions of observed swap yields mul tiplied by the Jacobian of the transformations (23). The Jacobian is nonlinear and time dependent.
Measurement errors are included in virtually all econometric models of the term structure because the state vector Y usually has low dimension (say, three or less) relative to the number of yields M.ll Thus, without additional sources of uncertainty, the models would imply deterministic relations among prices or yields that are clearly violated in the data. In constructing the likelihood function, we assume that the (nonzero) measurement errors {E,) for the swap yields follow univariate AR(1) processes, with innovations that are normally distributed and that may have nonzero correlation. The loglikeli hood function of the swap yield data is the sum of the logdensity of the noncentral chi squares of the squareroot processes, adjusted for j, and the logdensity of the multivariate normal associated with the measurement er rors.12 We estimate the swappricing model using weekly data from January 4, 1988 through October 28, 1994. The swap yields are constructed as follows. Weekly data on constantmaturity U.S. Treasury bond yields are constructed by concatenating yields on the current, ontherun Treasury bonds for matu rities two, three, five, seven, and ten years. The average of the quoted bidlask swap spreads are then added to these Treasury series to obtain the swap yield data.l3 The shortterm, sixmonth rate is taken to be the dollar LIBOR.
Initially, several onefactor models are fitted with models indexed by the point on the swap curve used to extract the single state variable Y'. In all of these models, the yield curves evaluated at the ML estimates tended to be too flat on average compared to the actual yield curves. Additionally, the one factor models are unable to fit simultaneously the volatilities of changes in yields at the long and short ends of the swap curve. Finally, the deviations between actual and modelimplied slopes of the swap curve exceeds 50 basis points for several extended periods. These poor results motivate our focus on a twofactor model. The state variables for the twofactor model are extracted from the two and tenyear swap rates. The rates assumed to be measured with error are those on the three, five, and sevenyear swaps. LIBOR is excluded from the econometric model, but is used subsequently in assessing the fit of the
l1 Pearson and Sun (1994) proceed instead by forming portfolios of bond yields and using as many portfolios as there are unobserved state variables. Thus, their model is presumed to fit the portfolio yields exactly. We have chosen to work directly with a crosssection of maturities of swap yields in order to assess the fit at specific points along,the swap yield curve, and because a crosssection of swap yields embodies more information about the yield curve than portfolios do.
l" See Chen and Scott (1993) and Pearson and Sun (1994) for further discussion of the likelihood function for squareroot diffusions. Our likelihood function differs from that used by Chen and Scott because of the nonlinear Jacobian for swap yields. We are grateful to Qiang Dai for developing the approximation to a noncentral chi square distribution used in our numerical routines.
l3The swap data are taken from the Telerate brokers screens and represent average bid and ask rates quoted by several large dealers.
Table I
Estimates of the 2Factor Model Weekly Data, January 8, 1988 to
October 28,1994
The parameters K~,Bi,and uzgovern the diffusion for the ith state variable Yi,
and hiis the associated risk premium. The ps are the autocorrelations of the measurement errors. Two sets of standard errors are given in parentheses below each estimate: the first is based on the outer product of the score and the second is based on the usual Hessian of the likelihood function.
Panel A: Parameters of Diffusion Y1
0.544  0.374  0.023  0.036 
(0.052)  (0.444)  (0.014)  (0.046) 
(0.029)  (0.017)  (0.0005)  (0.012) 
Panel B: Parameters of Diffusion Y
Panel C: Autocorrelations of Measurement Errors
twofactor model.14 Thus, the model has twelve parameters: four each for the two state variables (including the risk premia), the adjustment parameter 7, and the autocorrelations of the measurement errors for the three, five, and sevenyear swap rates. The parameter estimates are displayed in Table I, with two sets of standard errors. The first is computed from the outer product of the score vector, and the second is based on the usual Hessian matrix of the loglikelihood function. As is commonly the case in small samples, the two estimates differ.
The estimate of the adjustment factor 7 is 0.58, with estimated standard errors of 0.58 and 0.003. Together with the estimates of the ai,this implies an
l4 Initially, we keep LIBOR in the econometric model as a defaultable zero rate measured with error, but find that the model does not fit well relative to the version with LIBOR omitted altogether. Therefore, we report estimates for the version excluding LIBOR, and then examine the fit of the implied sixmonth swap rate to LIBOR. The previous draft of this article also presents estimates for a twofactor model in which the state variables are extracted from LIBOR and the tenyear swap rate. The fitting errors for the various swap rates are large and often over 60 basis points.
estimated longrun mean R = 0, + 0, of R, of 5.2 percent.15 The risk premium coefficients A, and A, are both negative, which implies that the term premiums are, on average, positive as maturity increases. In the context of pricing defaultable bonds, excess returns over the instantaneous rate R may reflect term premiums in the underlying riskless term structure andlor an increasing term structure of average credit spreads. The term structure of credit spreads for the zerocoupon bond yields implied by our swappricing model are subsequently examined in more depth.
The estimates of the meanreversion parameters K, and K, + A, are much larger than the corresponding estimates for the second factor. Indeed, K, + A, is very close to zero, suggesting that there is at most weak mean reversion in the second state variable. To interpret these findings, it is informative to examine the relations between the state variables extracted from the swap yields (Y' and d)and the swap yield curve.16 The sample correlation of aP1 with changes in the slope of the swap curve (10 year2 year) is approximately 0.99, whereas the correlation of aY1with changes in the tenyear swap yield is 0.02. The correlations between APand changes in the two, five, and tenyear swap yields are 0.60, 0.78, and 0.93, respectively. Thus, the first factor behaves like the negative of the slope of the swap curve, and the second factor behaves like the level of the tenyear swap yield.
Further insights come from examination of the coefficients of the zerocoupon bond yields, In BYT (see Fig. 2). From equation (20), these coefficients are In ai(r)/r and ~,(T)/T for the two state variables and j = 0.580 for the scale factor in equation (25). Across the maturity T, the coefficients ~,(T)/T
are nearly constant at unity and the coefficients 1n ~,(T)/T
are essentially zero. It follows that the second factor represents a parallel shift in the entire zerocoupon yield curve induced by changes in Y;.In light of the correlations of swap yields and state variables, this parallelshift factor is well proxied empirically by the tenyear swap yield. In particular, the likelihood function is not maximized by selecting a shortterm swap rate as the "level" risk factor. Figure 2 also shows that ~,(T)/T such that positive shifts in Y:induce
declines with maturity T flattenings in the slope of the swap curve. The findings that AY: and slope of the swap curve were nearly perfectly negatively correlated, and that ~,(T)/T
is large for small values of T, suggest that slope changes were associated with greater variation in shortmaturity yields during the sample period. The rapid decline in ~,(T)/T
as T increases is induced in the squareroot model largely by a fast rate of mean reversion of the first factor (large K,). The terms In ~,(T)/T contribute to an inherent upward slope
and In ~,(T)/T in the zerocoupon yield curve, as the former is upward sloping, while the latter
l5 We are grateful to Mark Fisher for pointing out that the longrun mean of the R was implausible in an earlier version of this article. Upon exploring the likelihood frontier further, we found a higher value of the likelihood function with the parameter estimates reported in Table 1 and, in particular, a plausible value for R. The timeseries properties of the fitted swap rates are unchanged compared to the previous results.
l6 More precisely, time series on the two state variables are computed by inverting the pricing model evaluated at the ML estimates for each date of the sample.
Maturity (Years)
Figure 2. Coefficients of Zero Coupon Bond Yields. The price of a zerocoupon bond with maturity T is
where Y1(t)and Y2(t)are the two state variables. Therefore, the yield on a Tyear zero is
This figure displays these coefficients of zero yields as functions of maturity 7.
is essentially zero for all maturities. Thus, y2represents nearly a pure parallel shift. On the other hand, at (Y1 = 0, y2= O), the yield curve is upward sloping due to the contribution of In ~,(T)/T.
It follows that the slope of the yield curve is matched by adjusting Y', and that matching an inverted swap curve may require a large value of Y1 to offset the upward slope induced by In u,(T)/T. In such cases, the resultant large value of Y;~,(T)/T (the mean of Y: is 0.37) would clearly overstate the level of zerocoupon yields, especially for short maturities. Thus, in order to simultaneously explain the temporal behavior of the slope of the swap curve and fit the average level of the curve, an adjustment is necessary. This explains the economically significant value of the constant adjustment of 3.
Descriptive statistics of the swap yields implied by the models are presented in Table 11. Panel A displays the sample means of the levels and first differ ences of the historical swap yields ("Sample"), as well as the swap yields implied by the twofactor model. The model matches the average three, five,
Table I1 Descriptive Statistics from the 2Factor Model
Panel A displays the sample means of the historical (Sample = c7) and modelimplied (Model = t') rates on the London Interbank Offering Rate (LIBOR) and nyear swap contracts for n = 3, 5,and
7. Results for the level of rates (c7) and changes in rates !Ac7i are presented. Panel B presents the analogous results for the sample standard deviations of swap yields. Panel C restricts attention to changes in swap rates (dc') and displays sample historical and modelimplied skewness and kurtosis statistics. Panel D displays the standard deviations of the differences between the actual and modelimplied swap rates (c7 t7),and of (Ac' At7).
Panel A: Means of Swap Yields
T =  LIBOR  3 Yr  5 Yr  7 Yr  

Sample  c'  
AcT  
Model  t  
AtT 
Panel B: Standard Deviations of Swap Yields
Sample c '
ac7
Model tAtT
Panel C: Skewness/Kurtosis for Changes in Swap Yields
Sample Skewness 0.316 0.156 0.218 0.216
Kurtosis 4.357 3.001 3.071 3.148 Model Skewness 0.052 0.141 0.152 0.145
Kurtosis 3.096 2.875 2.878 2.931
Panel D: Standard Deviations of Fitting Errors (Basis Points)
and sevenyear swap yields, as well as the standard deviations of yield changes (Panel B), to the third decimal place. The skewness and kurtosis statistics of the changes in implied swap yields displayed in Panel C are slightly smaller than the corresponding statistics for their sample counterparts.
The two and tenyear swap yields are fit exactly by construction. As such, perhaps the most challenging rate to fit is the fiveyear rate, because changes in the curvature of the swap curve may result in movements of the fiveyear swap rate that are independent of the longer and shorter ends of the swap curve. Figure 3 shows that the model fits the fiveyear swap yield to within 20 basis points (bp) over the sevenyear sample period. The mean error is only 0.3 basis points. Moreover, the sample standard deviation of the measurement error for the level of the fiveyear swap rate is 7.2 bp (see Table 11,Panel D), which was approximately the bidlask spread during this sample period.
15!1
88 89 90 91 92 93 94
Date
Figure 3. Deviation from Fit of Swap Yields. The dashed line is the deviation, in basis points, between the actual and fitted fiveyear swap rates. The solid line is the corresponding deviation for the sevenyear minus threeyear swap yield spread.
Figure 3 shows that there is also a close correspondence between movements in the slope of the swap yield curve (7 year3 year) and the implied slope from the model. The maximal deviation is less than 16 basis points over the sample period, with an average fitting error of 1.1basis points. The most challenging period to fit seems to have been the trough in swap rates during late 1993 and early 1994.
As additional evidence on goodnessoffit, we report in Table 111,Panel A, the results from regressing the changes in the actual swap yields against their fitted counterparts. If the twofactor model describes the data well, then we would expect that the estimated intercept (h) and slope (p)to satisfy p = 1and h = 0, and that the R~Swould be close to one. Regressions are run with differenced yields, because of the high degree of persistence in levels and our desire to evaluate the model's explanatory power for changes. Only the results for LIBOR differ markedly from those predicted by the model. Indeed, it is striking that the model explains approximately 95 percent of the variation of the changes in individual swap rates over the seven year sample period. Table I1 expresses this finding in terms of basis points: The standard deviations of the fitting errors for the swap yields are between 4.48 bp and 7.16 bp.
Table I11
Correlations of Actual and Fitted Yields
Panel A displays the results from regressing the changes in actual swap rates (LC;) on the changes
in fitted rates from the model (At;). Standard errors of the estimated intercept (iu) and slope (P)are
given in parentheses. R2is the coefficient of determination and s.e.e. is the standard error of the
residual. Results are reported for the London Interbank Offering Rate (LIBOR) rate and three
year, five year, and seven year swap rates. Panel B displays the corresponding results for a model
fit with Treasury bond yields instead of swap rates (i.e., c; is the yield on a Tyear Treasury bond).
Maturity  a  fi  R2 s.e.e. 
Panel A: Actual Against Fitted Swap Yiel  d Regressions  
LIBOR 
3 Year
5 Year
7 Year
Panel B: Actual Against Fitted U.S. Treasury Yield Regressions
3 Year
5 Year
7 Year
For comparison, we reestimate our twofactor squareroot model using con stantmaturity U.S. Treasury yields for the same sample period.17 As with swaps, the state variables are extracted from the two and tenyear maturity instruments. The results from regressing actual on fitted changes in Treasury yields are displayed in Table 111, Panel B. The Treasury model also fits reasonably well during this sample period, although the coefficients are significantly different from one at conventional significance 1evels.la Also, the standard error of the estimate (s.e.e.) for the fiveyear rate is larger for the Treasury model than for the swap model.
Although LIBOR is not used in estimation, the sixmonth rate implied by the twofactor model can be computed at the ML estimates and compared to LIBOR. The results in Table I1 show that the fitted LIBOR is about 18 basis
The constant maturity Treasury data are not exactly what one would want in order to fit a model for the U.S. Treasury curve, since the maturities of, say, the tenyear note change slightly and the ontherun issue is not always at par between auctions. However, for our purposes of providing a benchmark for assessing fit, the data seem adequate.
l8 Standard errors are computed using Hansen's (1982) correction for serial correlation and heteroskedasticity with five noncontemporaneous terms in the asymptotic covariance matrix.
Date
Figure 4. Deviation from Fit to LIBOR. This figure displays the deviation in percent between the actual and fitted sixmonth LIBOR rates.
points too small on average and is more volatile than the actual LIBOR. Moreover, the standard deviation of the fitting error is over 32 bp. The sample skewness of changes in LIBOR is negative and its kurtosis is larger than three, whereas the sample estimates of skewness for the changes in swap rates are positive and the sample estimates of their kurtoses are close to three. Not surprisingly, the model does not explain the negative skewness and excess kurtosis of LIBOR. Together, these findings suggest that the sixmonth LIBOR series has distinctive characteristics that are not shared by the longer end of the swap yield curve. One interpretation of these results is that LIBOR loans are of a somewhat lower average perceived quality than multiyear swaps of the same credit rating. This is consistent with the evidence presented in Sun, Sunderson, and Wang (1993). Alternatively, it may be that additional non credit factors are needed to describe simultaneously the distribution~ of the long and short ends of the swap curve.
The poor fit for LIBOR is displayed in Figure 4. The deviations between actual and implied LIBOR fluctuate substantially and often exceed 50 basis points in absolute value. Moreover, there is evidently a seasonal pattern to the spreads between actual and fitted LIBOR. In every year, there are large changes in the deviations around the calendar year end. (The vertical lines in Figure 4 occur on the last Friday of each calendar year that markets were open.) For endofyear 1989, 1990, and 1991, the actualfitted LIBOR differ ences fell substantially. On the other hand, for endofyear 1992 and 1993,
there is evidence of an increase, although in the latter case the increase was quickly followed by a steady decline in the actualfitted LIBOR deviations. These calendar effects may well be related to balancesheet adjustments by financial institutions concerned about capital requirements or the risk profile of their securities positions at year end. Whatever the source, these results suggest that using a shortterm rate as the state variable describing the entire yield curve may lead to misleading conclusions about the shapes of the distri butions of longterm swap yields.
How can the poor fit for LIBOR and relatively good fit for the swap rates be reconciled with the assumption of homogeneous LIBORswap market credit quality? One reconciliation comes from noting that equations (4) and (18) can be used to fit the modelbased discount factors B,' implicit in the swap market without reference to the LIBOR market. If, as seems common in practice, swap traders also extract estimates B,' for pricing swaps from swap rates (and Eurodollar futures) without reference to LIBOR markets, then our economet ric analysis provides a modelbased construction of the traders' implied dis count factors, B,', and pricing rules for the longer end of the swap curve. At the short end of the swap curve, additional factors are evidently necessary for our model to generate consistent pricing of Eurodollar futures and LIBOR con tracts, that is, to reproduce the discount factor at the short end of the swap zerocoupon curve.
Regressing the deviation between actual and fitted LIBOR (DEVLIB) on the first and second lagged values of itself giveslg
DEVLIB, = 0.002 + 0.787 DEVLIB,I + 0.173 DEVLIB,,,
(0.006) (0.05) (0.05)
with an R2of 0.90 and with a standard error of the residual of 11.5 basis points. The finding that a loworder autoregression explains most of the within sample variation in the misfitting of LIBOR suggests that the twofactor model could be modified to include a third factor to accommodate the dynamics of the entire swapLIBOR yield curve. Evidence consistent with this conjecture is presented in Dai and Singleton (1996).
IV. Analysis of Implied Swap Zerocoupon Yields
We turn next to an examination of the properties of the zerocoupon bond yields implied by our swappricing model. Specifically, we compute the yields (log B,'lr) x 100 (i.e., continuously compounded yields) with Bi, given by equation (20), evaluated at the ML estimates. Our focus is on the spreads of these zerocoupon bond yields to their counterparts in the U.S. Treasury market. In light of the small residuals from fitting swap yields, the defaultable zerocoupon yields implied by our model should serve as reasonable proxies for
''While supported by a different model and different data, our characterization is largely consistent in this regard with independent work by Brown and Schaefer (1993), who study U.S. Treasury yields.
Mean Zero Spreads
0'  I  
0  1  2  3  4  5  6  7  8  9  10 
Maturity (Years) 
Figure 5. Means and Standard Deviations of Swap Zero Spreads, Sample Period: Jan uary 15, 1988 to October 2, 1994. The squares and diamonds are the standard deviations and means, respectively, of the spreads between zero rates implied by the swap and Treasury yield curves. The mean and standard deviation of the actual LIBOR spread to the sixmonth Treasury bill rate are displayed as a triangle and x, respectively.
the defaultable zero yields implicit in the swap curve. For the U.S. Treasury zerocoupon (continuously compounded) yields, we use the discount function implied by a statistical spline model.20 Whereas the swap curve is sparsely reported and a modelbased interpolation scheme seems desirable, the sub stantial data on coupon U.S. Treasury yields and strip prices are used to compute Treasury zero rates.
Figure 5 displays the term structure of sample means and standard devia tions of spreads between the swap and U.S. Treasury zerocoupon bond yields. The term structure of zerocoupon default spreads is, on average, upward sloping during this sample period from about 20 basis points for the sixmonth spread, up to about 40 basis points for the fiveyear spread. Beyond five years to maturity, the spreads are nearly constant at about 38 to 40 basis points. The large difference between the means of the fitted and actual spreads of LIBOR to the sixmonth Treasury rate is consistent with the presence of money market effects noted earlier that are not captured by the twofactor swap model.
20 The U.S. Treasury zerocoupon yields, provided to us by Goldman Sachs, are computed from the coupon yields and strip prices.
Spread volatilities tend to be increasing beyond two years to maturity, although volatilities are roughly flat between four and seven years. (The volatilities of the zerocoupon bond yields decline with maturity for both the
U.S. Treasury and swap markets.) As maturity declines below two years, the spread volatilities increase. The actual sixmonth spread volatility in the money market is much lower than the implied spread volatility from the swap model.
To assess the relative importance of the liquidity (1,) and credit (h,h,) components of R, in determining swap zero spreads, vector autoregressions (VAR) are estimated for zero spreads, and for proxies for credit and liquidity. For squareroot diffusions, the zero yields are linear in the state variables Y1 and y2.If the U.S. Treasury curve can similarly be described by a discount function that is the sum of squareroot diffusions (Table I11 suggests that this is approximately so during our sample), then these VARs can be interpreted as linear projections of differences in linear combinations of the state variables driving the riskadjusted and Treasury interest rate processes onto the vari ables in the VAR.21 The variables included in the VARs are the sixmonth Treasury zero yield (TB6), the spread between the generic threemonth rep0 rate for the tenyear Treasury note and the rep0 rate of the current ontherun Treasury note (REPOSP), the spread between rates on BAA and AAArated commercial paper (CPS), and the spread between the tenyear zero rates implied by the swap and Treasury markets (ZEROSP10).22 All VARs are fit with eight lags.z3
The inclusion of TB6 is intended to capture the effect of the level of riskless interest rates on spreads for defaultable zerocoupon bonds. CPS captures default risk. Liquidity differences may also contribute to a spread between BAA and AAArated corporate paper, but we expect these differences to be small compared to the effects of the different credit rating.
REPOSP captures the specialness of the Treasury coupon note with the same maturity as the zero being studied. Of course, there is not a simple linear mapping between the specialness of a coupon bond and the spread in the underlying zero markets. However, the specialness in the coupon market is the only set of rep0 data available, and it should shed some light on the importance of rep0 effects on swap spreads. An increase in REPOSP (the rep0 rate on the
Although, under the assumptions of our model, the innovations in the VARs are not normal, the model does imply linear expectations conditional on past Ys. Therefore, we expect that the VAR analysis will provide suggestive descriptive evidence about the contributions of liquidity and credit factors to swap spreads.
22 We also fit VARs including the spread between Moody's corporate bond yield indices for BAA and AAArated credits as an additional conditioning variable to proxy for the longterm corporate spreads. However, this spread is not a significant explanatory variable for zero swap spreads, and the proportion of variance of ZEROSPlO explained by the Moody's series is essentially zero.
23 A Bayesian prior on the lag distribution, which has the own first lag of each variable equal to unity and all of the other coefficients equal to zero, is also imposed. Specifically, in the notation of the manual of the RATS statistical computer package, we use a symmetric prior with TIGHT = 0.15, and 0.5 as the weight on the other variables in the equation. The results are not substantially different with and without the prior.
Table IV FTests of Exclusion Restrictions in VARs
This table displays the values of the Fstatistics for testing the hypothesis that all of the coefficients on lagged values of the variable indicated under "Variable Excluded are zero in the equation with the dependent variable given in the lefthand column. Marginal significance levels of the statistics are given in parentheses. Two VARs are estimated: VAR4 includes the 6month Treasury bill rate (TB6), the spread between generic and ontherun repo rates for tenyear Treasury bonds (REPOSP), the spread between BAA and AAArated commercial paper (CPS), and the spread between the tenyear zero rates implied by the swap and Treasury markets (ZEROSP10); VAR3 includes TB6, REPOSP, and ZEROSP. The * indicates exclusion of CPS from the VAR3 model.
Variable Excluded
TB6  REPOSP  CPS  ZEROSPlO  

Dependent  
Variable  VAR4 VAR3 VAR4  VAR3  VAR4 VAR3 VAR4  VAR3  
TB6  476.5  
(0.000)  
REPOSP  0.54  
(0.827)  
CPS  0.80  
(0.602)  
ZEROSPlO  1.42  
(0.187) 
current bond falls relative to the generic rep0 rate) implies that holders of the currentissue bond receive an extra "dividend." Therefore, as argued by Grin blatt (19951, when the ontherun Treasury notes go on special, Treasury prices tend to rise and swap spreads tend to increase.24
Table IV presents the Ftests and their marginal significance levels for the exclusion of all eight lags for each of the explanatory variables in the VAR. Two sets of results are displayed: those for a fourvariable VAR (VAR4) and those for a threevariable VAR (VAR3) that excludes the commercial paper spread CPS. The test statistics for the ZEROSPlO equation, presented in the last row of Table IV, suggest that the histories of ZEROSPlO and REPOSP have significant predictive power for ZEROSP10. The Fstatistic for TB6 is smaller, but significant at the 10 percent level in VAR3. The corporate spread CPS has little explanatory power for ZEROSPlO in VAR4, which is why we also exam ine VAR3. Interestingly, REPOSP is significant or borderline significant at conventional significance levels in all equations of both VARs.
24 Alternatively, Evans and ParenteBales (1991) and Brown, Harlow, and Smith (1994) argue that the costs to dealers of hedging net swap exposure is a key factor underlying changes in swap spreads. Brown, Harlow, and Smith (1994) use the repo rate as a proxy for a dealer's costs of hedging a swap exposure with Treasury bonds. A net exposure to the pay fixed side of a swap would be hedged by taking a long position in government bonds, which are financed in the repo market. Thus, they argue that a decrease in the repo rate, by lowering hedging costs, reduces swap spreads. The sign of this effect depends, however, on whether in the aggregate dealers are net pay or receivefixed counterparties in the market, so this hedging effect may change in sign over time.
Within a VAR system, the effects of, say, CPS on the future course of ZEROSPlO depends not only on the direct effects tested by the Fstatistics in Table IV, but also on the indirect effects of CPS on the other variables in the VAR. CPS has significant predictive power, for example in the equation for TB6. In light of these feedback effects within the multivariate system, it is instructive to examine the impulseresponse functions for the VARs. Each function traces out the effects of an innovation in one of the variables on the future course of ZEROSP10, taking into account the dynamic interactions among the variables. The shock is positive with magnitude equal to the estimated standard deviation of the variable's own residual in the VAR. The results for VAR4 are displayed in Figure 6.25
The impulseresponse patterns are traced out over one year (52 weeks) subsequent to the impulse. In order to assess the significance of the responses, standarderror bands are computed using Monte Carlo methods. Specifically, assuming that residuals in the VARs are i.i.d. normal random variables, draws are made from the posterior distribution of the VAR coefficients and the impulseresponse functions are computed based on these coefficients. These calculations are repeated 1000 times to derive upper and lower twostandard error bands for the responses over 52 weeks.26
The cumulative effects of these shocks are summarized in Table V, which presents the decomposition of variance. For a given VAR, each row must sum to one as the sum represents the total variation in ZEROSP10. For example, after 4 weeks, about 6.29 percent of the variation of the change in ZEROSPlO due to shocks in the four variables in VAR4 is due to TB6, about 2.74 percent is due to REPOSP, and so on.
The largest response of ZEROSPlO over the first weeks following the shocks is due to its own shock (Figure 6d). However, the effects of the own shocks die out relatively quickly. An increase in the level of interest rates (TB6) implies an increase in the tenyear zero spread (Figure 6a). This is consistent with the view that credit spreads widen during market selloffs and narrow during rallies. Given the strong positive correlation between generic rep0 and TB6, this pattern is also consistent with an explanation based on dealer hedging costs under the presumption that dealers have more payfixed than receive fixed swaps on their books. Interest rate "level" effects peak after about six
25 The interpretation of the patterns is subject to the usual caveat that the VAR residuals are transformed to orthogonal shocks before computing the impulse responses. The variables are ordered in VAR4 as presented in Figure 6: TB6, REPOSP, CPS, and ZEROSP10. So contempora neous correlation between the residual for TB6 and the residuals for the other three variables is attributed to variation in TB6, etc.
26 In order to make determination of the posterior distribution tractable, we work directly with the unconstrained VAR and assume normally distributed innovations. So, in particular, the Bayesian priors imposed in calculating the results in Tables IV and V are not imposed here. As noted previously, the results are similar with and without the priors imposed. Of course, these are estimated standarderror bands based on normal innovations. Their use is therefore subiect to the usual caveats that the sample may be small or the distribution for the Monte Carlo analysis may be misspecified.
a: Shocks to TB6
3.2
2.4
1.6 .......................................
0.8
0.0 ..................................
0.8
1.6 7
b: Shocks to REPOSP
3.2
2.4
1.6 ______________
0.8
0.0 .
1
0.8 I
c: Shocks to CPS
3.2
2.4
1.6
d: Shocks to ZEROSPlO
3.2 I
(
Figure 6. Impulse Responses for ZEROSPlO in VAR4. A fourvariable VAR is estimated with the sixmonth Treasury bill rate (TB6), the spread between generic and ontherun repo rates for tenyear Treasury bonds (REPOSP), the spread between BAAand AAArated commercial paper (CPS), the spread between the tenyear zero rates implied by the swap and Treasury markets (ZEROSP10). Then each variable is perturbed by a onestandard deviation shock in its (orthogo nalized) innovation. This figure displays the responses of ZEROSPlO to these shocks in basis points, over a period of 52 weeks following the impulses. The dashed lines represent plus and minus twostandard error bands around the impulse responses, estimated by Monte Carlo.
months (Table V, column TB6),and over a twoyear horizon explain about 11 percent of the variation in ZEROSP10.
A shock to REPOSP initially has a small negative effect, but after the first couple of weeks the effect turns positive and statistically significant. This positive effect is predicted by Grinblatt's (1995) model of liquidity premiums, with REPOSP representing the convenience yield associated with treasury bonds. Over an eightweek horizon, REPOSP explains about 11.9 percent of the variation in ZEROSP10. Subsequently, variation in ZEROSPlO due to REPOSP levels off at about 20 percent.
The immediate impact of an increase in CPS on the zerocoupon spreads is essentially zero. After a few weeks, the zerocoupon spreads tend to decrease (Figure 6c);a widening of the credit spread between BAA and AAA commercial
Table V Decompositions of Variance of TenYear Swap Zero Spreads
Each column displays the proportion of the variance of the error from forecasting ZEROSPlO k weeks ahead, due to the impulse in the variable indicated in the column heading. k ranges from 1to 104 weeks (2 years). Two VARs are estimated: VAR4 includes the 6month Treasury bill rate (TB61, the spread between generic and ontherun repo rates for tenyear Treasury bonds (RE POSP), the spread between BAA and AAArated commercial paper (CPS), and the spread between the tenyear zero rates implied by the swap and Treasury markets (ZEROSP10); VAR3 includes TB6, REPOSP, and ZEROSPlO. The * indicates exclusion of CPS from the VAR3 model.
Percentage Due to Impulse From
TB6  REPOSP  CPS  ZEROSPlO  

Weeks  
Ahead  VAR4  VAR3  VAR4  VAR3  VAR4  VAR3  VAR4  VAR3 
paper leads to a narrowing of the spread between the swap zero and Treasury zero yields. A potential explanation for these patterns is that relative credit spreads (BAA AAA) tend to widen during recessions when default probabil ities increase, even though commercial paper spreads to Treasury bill rates tend to narrow. Thus, a widening of CPS might be expected to eventually lead to a narrowing of ZEROSPlO if the former reflects a weakening of the economy. That commercial paper spreads to Treasuries are informative leading indica tors of the U.S. business cycle has been documented by Friedman and Kuttner (19931, among others. CPS might reasonably be expected to have a similar property.
Interpreting the effects of REPOSP and CPS on ZEROSPlO as liquidity and credit shocks, respectively, the patterns suggest that liquidity shocks are relatively important over short horizons of a few months. After about 9 weeks, the response to REPOSP shocks is insignificant out to about six months. Credit shocks have little impact on ZEROSPlO over the first couple of months. Over longer horizons the effect of CPS grows in importance to the point that after two years it accounts for the second largest percentage (20.3 percent) in the variance decomposition for ZEROSP10. The estimated standard error bands for CPS include zero for all of the horizons out to one year, however, so credit effects on swap spreads may be weaker than the numbers in Table V suggest.
For comparison, we also estimate a threevariable VAR including TB6, REPOSP, and ZEROSP10, in that order. The impulse responses are displayed in Figure 7. The effect of TB6 is larger than in VAR4 (compare Figures 6a and 7a). This is reflected in the variance decompositions (Table V), where 42 percent of the variation in ZEROSPlO over a oneyear horizon is attributable
a: Shocks to TB6
3.5
3.0
2.5
2.0 .______
/."I
1.5
1.0
0.5 ...........................................
0.0 0.5 1.0
6"3"6"9"12 15'18 21 24 27 30 33 36 39 42 45 48 51
b: Shocks to REPOSP
3.5
3.0
2.5
c: Shocks to ZEROSPlO
3.5 , I
Figure 7. Impulse Responses for ZEROSPlO in VAR3. A threevariable VAR is estimated with the sixmonth Treasury bill rate (TB6), the spread between generic and ontherun repo rates for tenyear Treasury bonds (REPOSP), the spread between the tenyear zero rates implied by the swap and Treasury markets (ZEROSP10). Then each variable is perturbed by a onestandard deviation shock in its (orthogonalized) innovation. This figure displays the responses of ZEROSPlO to these shocks in basis points, over a period of 52 weeks following the impulses. The dashed lines represent plus and minus twostandard error bands around the impulse responses, estimated by Monte Carlo.
to shocks in TB6 in VAR3, compared with 12.8 percent in VAR4. Evidently, omission of CPS leads to the variation in ZEROSPlO attributed to CPS in VAR4 being attributed to TB6 in VAR3.27 An explanation for these results can be gleaned from the results for CPS in VAR4. The Fstatistics in Table IV suggest that there is statistically significant feedback from CPS to TB6. Furthermore, the impulse responses for TB6 show a negative effect of a shock to CPS on TB6: A widening of the commercial paper credit spread is associated
27 Notice also that the variation in ZEROSPlO due to its own shock over a twoyear horizon declines in VAR3, and this reduction is attributed largely to TB6. However, in light of the wide standard error bands around the impulse response function of ZEROSPlO over horizons beyond six months (Figure 71, this change may not be of much significance.
0 ++tt++i++.'~ 9 0 91 92 93 94 Date
Figure 8. Repo Rates for OntheRun Ten and FiveYear U.S. Treasuries.
with a subsequent decline in TB6. Associating widening relative credit spreads within the commercial paper market with a weakening economy, it appears that, as the economy weakens, Treasury rates tend to decline, relative com mercial paper spreads tend to widen, and the tenyear swap zero spread tends to narrow. Thus, both changes in Treasury bill rates and relative commercial paper credit spreads reflect the weakness of the economy and any induced changes in the likelihood of default by swap counterparties. When CPS is omitted from the VAR, the predictive power of changing relative commercial paper credit spreads for ZEROSPlO is largely captured by TB6.
The effects of shocks in tenyear rep0 specialness, as proxied by REPOSP, are striking and much larger than when the corresponding rep0 spread for fiveyear Treasury bonds is used instead. Upon further exploration, we find a strong seasonal pattern in the ontherun tenyear rep0 rate (actual rate level) that is not present in the fiveyear rate (see Figure 8).Especially since 1992, there is a strong quarterly seasonal effect, and the tenyear rep0 rate (REPO10) generally is below the fiveyear rep0 rate. This suggests that the liquidity premium (convenience yield) is typically larger in the tenyear sector and more seasonally variable.
Several factors contribute to the quarterly seasonality in REPO10. First, the peaks and troughs in REP010 correspond approximately to the calendar endofquarter dates. Efforts by Treasury bond dealers to cover their short positions at the ends of quarters puts pressure on bond borrowing and, hence, downward pressure on rep0 rates (Wall Street Journal, September 29, 1992, page C20). Short covering is encouraged by the treatment of shortsales as liabilities on dealers' quarterly balance sheet reports. Additionally, the Trea sury note futures contracts expire on the last date of the end of each calendar quarter. Borrowing bonds for delivery could contribute to this pressure in the rep0 market.28
However, the squeezes in the bondborrowing market do not occur every quarter. Indeed, during 1992 and 1993, the largest dips in REPOlO were during March and September. The explanation for this pattern appears to be the reopening of tenyear note issues by the U.S. Treasury. Specifically, in May and November of 1992 and 1993 the U.S. Treasury reopened the bond from the previous auction instead of auctioning a new note. The presence of an out standing stock of the reopened note on the reopening date implied that these issues did not go on special as much as the typical newly auctioned notes. Consequently, the rep0 rate for the ontherun tenyear Treasury did not fall as much as for regular auction dates (February and August). The induced sea sonal patterns in REPOlO appear to be associated with significant widenings of tenyear swap spreads after regular auctions. These Treasury "experiments" show a strong liquidity effect on swap spreads.
Within the fourvariable VAR examined, nearly half of the variations of swap spreads over twoyear horizons are explained by their own shocks. This suggests that an understanding of the timeseries properties of swap zero spreads may require a deeper investigation of swapspecific market activity. For example, segmentation of supply or demand pressures by maturity sector due to institutional or accounting considerations may be important.
The behavior of swap spreads is also complicated by the state tax abatement on Treasuries, which would widen spreads, but in a manner whose dynamics are not clear or easily measured.
V. Conclusion
We have shown that, under the assumption of symmetric counterparty credit risk, swaps can be priced using standard termstructure models based on a risk and liquidityadjusted discount rate. For the purpose of our econo metric analysis, we assume that the adjusted discount rate can be expressed as the sum of two independent squareroot diffusions. Upon maximizing the likelihood function for about seven years of weekly data, we find that this twofactor model fits many aspects of the swap term structure well. The primary exception is the very short end of the swap curve represented by sixmonth LIBOR. The deviations between the actual and the modelimplied LIBORs evidence strong seasonality and tend to be particularly large near the calendar yearends. The autocorrelation structure of the deviations suggests that a third factor might allow one to fit the shortterm LIBOR curve and the swap curve simultaneously. This possibility, as well as the implications of using two versus threefactor models for valuing LIBORbased derivatives, are left for future research.
28 Other factors include the heavy use of tenyear treasury notes as hedges against positions in mortgagerelated securities and the more frequent auction cycle of five versus tenyear notes.
Using the estimated model, we then examine the dynamic properties of the spreads between the zerocoupon bond yields implicit in the swap curve and their Treasury counterparts. Three notable findings are as follows. First, the specialness of the ontherun tenyear Treasury note has a positive effect on zero spreads that tends to peak in the first few weeks following the impulse, and over long horizons explains about 20 percent of the variation in tenyear zero spreads. These patterns suggest that liquidity advantages to trading in Treasury markets, and Treasury rep0 specials, have a substantial effect on swap spreads, consistent with Grinblatt's model, but that this liquidityrep0 effect is not the only, or indeed the primary, determinant of swapTreasury zero spreads. The same conclusion applies to the actual swap spreads as well.
This is not surprising given that the convenience yield 1, captures relative liquidity and rep0 effects in swap and treasury markets. The demands for the pay and receivefixed sides of swaps change over time and along the maturity spectrum. These demands may well have important effects on swap spreads. Indeed, the second conclusion that we draw from the VAR analysis is that, for the explanatory variables considered, between 35 percent and 48 percent of the variation in the tenyear zero spread over a twoyear horizon is explained by its own shocks. That is, after accounting for the proxies of hedging costs, liquidity, and credit effects, a substantial fraction of the variation in swap spreads is left unexplained. Brown, Harlow, and Smith (1994) find weak evidence of an effect on the supply of corporate debt in the long end of the swap curve. Further exploration of these and other supply effects in the swap market, as well as the impact of asymmetric tax treatments, seems worthwhile.
Finally, interpreting the spreads between AAA and BAA commercial paper rates as a credit spread, we find that changes in the market's perception of credit risk have large effects on zero spreads at the tenyear maturity. The response of swap spreads to credit shocks is very different from the response to liquidity shocks. The credit effects are weak initially, then peak about six to seven months after the impulse. In terms of contribution to variation in zero spreads, the relative importance of credit shocks increases for over two years, reaching 20 percent.
However, the confidence intervals for the impulse responses for credit shocks are wide. Moreover, the channel by which CP credit spreads are correlated with swap spreads appears to be through the shortterm riskless rate. The patterns are consistent with economic downturns leading to both a deterioration in credit quality (widening of relative credit spreads) and lower interest rates. Although an increasing hazard rate of default during economic downturns is one interpretation of these patterns, further analysis of the links between default hazard rates, the levels of riskless interest rates, and swap spreads seems necessary to conclude more definitively that businesscycle impacts on credit quality represent an economically significant credit compo nent to variation in swap spreads.
These findings will prove useful in guiding a more complete model of swap spreads. Within our theoretical framework, the next step would be to param eterize the process generating the riskless rate (r,),the liquidity convenience yield (I,), and the credit adjustment (h,h,). The VAR analysis suggests that the parameterizations adopted should be sufficiently flexible to allow very differ ent dynamic structures for I, and h,h,, and allow for correlation between default hazard rates and the riskless interest rate. The VAR analysis may also provide potentially useful insights into the dynamic relations among our liquidity and credit proxies that may prove useful in parameterizing the correlations between the processes for r,, I,, and h,h,. As elaborated by Duffie and Singleton (1996), one is unable to disentangle the separate influences of A, and h, using only yield data, without the benefit of default incidence or loss data, or of data on the prices of instruments that respond nonlinearly with the prices of the underlying bonds at default, such as creditspread options.
REFERENCES
Abken, Peter, 1993, Valuation of defaultrisky interestrate swaps, Advances In Futures And Options Research 6, 93116. Bansal, Ravi, and William Coleman 11,1996, A monetary explanation of the equity premium, term premium, and riskfree rate puzzles, Journal of Political Economy 104, 11351171. Bernanke, Ben S., 1990, Clearing and settlement during the crash, Review ofFinancial Studies 3, 133151. Brown, Kevin, W. V. Harlow, and Donald J. Smith, 1994, An empirical analysis of interestrate swap spreads, Journal of Fixed Income 3, 6178.
Brown, Roger, and Steven Schaefer, 1993, Interest rate volatility and the term structure of interest rates, Philosophical Transactions of the Royal Society: Physics, Science, and Engi neering 347, 563576.
Chen, RenRaw, and Louis Scott, 1993, Maximum likelihood estimation for a multifactor equilib rium model of the term structure of interest rates, Journal of Fixed Income 3, 1431. Cooper, Ian A., and Antonio S. Mello, 1991, The default risk of swaps, Journal of Finance 46, 597620. Cox, John C., Jonathan E. Ingersoll, Jr., and Stephen A. Ross, 1985, A theory of the term structure of interest rates, Econometrica 53, 385407. Dai, Qiang, and Kenneth J. Singleton, 1996, Specification analysis of affine term structure models,
Working paper, Graduate School of Business, Stanford University. DufEe, Darrell, 1996, Special repo rates, Journal of Finance 51, 493526. DufEe, Darrell, and Ming Huang, 1996, Swap rates and credit quality, Journal of Finance 51,
921949. Duffie, Darrell, and Rui Kan, 1996, A yieldfactor model of interest rates, Mathematical Finance 6, 379406. Duffie, Darrell, and Kenneth J. Singleton, 1993, Simulated moments estimation of Markov models of asset prices, Econometrica 61, 929952. Duffie, Darrell, and Kenneth J. Singleton, 1996, Modeling term structures of defaultable bonds, Working paper, Graduate School of Business, Stanford University. Evans, Ellen, and Gioia ParenteBales, 1991, What drives interest rate swap spreads, in Carl R.
Beidleman, Ed.: InterestRate Swaps (Business One Irwin Pubs., Homewood, Illinois). Feller, William, 1951, Two singular diffusion problems, Annals of Mathematics 54, 173182. Friedman, Benjamin, and Kenneth Kuttner, 1993, Why does the paperbill spread predict real
economic activity? in NBER Macro Annual, 213253. Grinblatt, Mark, 1995, An analytic solution for interestrate swap spreads, Working paper, UCLA Anderson Graduate School of Management. Hansen, Lars, 1982, Large sample properties of generalized method of moments estimators, Econometrica 61, 929335. Heath, David, Robert Jarrow, and Andrew Morton, 1994, Bond pricing and the term structure of interest rates: A new methodology for contingent claims valuation, Econometrica 60, 77106.
A Model of Term Structure of Interest Rate Swap Yields 1321
Hull, John, and Alan White, 1992, The impact of default risk on the prices of options and other derivative securities, Journal of Banking and Finance 19, 299322. Jaffee, Dwight, 1975, Cyclical varations in the risk structure of interest rates, Journal of Monetary Economics 1,309325. Jarrow, Robert, and Steven Turnbull, 1995, Pricing options on financial securities subject to default risk, Journal of Finance 50, 5386. Jordan, Brad, and Susan Jordan, 1997, Special repo ratesan empirical analysis, Journal of Finance, Forthcoming.
Krylov, N., 1980, Controlled Diffusion Processes (Springer Verlag, New York).
Li, Haitao, 1995, Pricing of swaps with default risk, Working paper, Yale University, School of Organization and Management. Litzenberger, Robert H., 1992, Swaps: Plain and fanciful, Journal of Finance 47, 831850. Merton, Robert, 1974, On the pricing of corporate debt: The risk structure of interest rates,
Journal of Finance 29, 449470. Nielsen, S., and Ehud Ronn, 1995, The valuation of default risk, in Corporate Bonds and Interest Rate Swaps (University of Texas at Austin, Austin).
Pearson, Neil, and TongSheng Sun, 1994, An empirical examination of Cox, Ingersoll, and Ross model of the term structure of interest rates using the method of maximum likelihood, Journal of Finance 49, 929359.
Rendleman, Robert J., Jr., 1992,How risks are shared in interest rateswaps, Journal of Financial Services Research 7, 534. Smith, Clifford W., Jr., Charles W. Smithson, and Lee MacDonald Wakeman, 1988,The market for interest rate swaps, Financial Management 17, 3444. Solnik, Bruno, 1990, Swappricing and default risk: A note, Journal of International Financial Management and Accounting 2, 7991. Sorensen, Eric H., and Thierry F. Bollier, 1994, Pricing swap default risk, Financial Analysts Journal 50, 2333. Stock, James, and Mark Watson, 1989,New indexes of coincident and leading economic indicators, in NBER Macro Anual, 351395. Sun, TongSheng, Suresh Sundaresan, and Ching Wang, 1993, Interest rate swapsan empirical investigation, Journal of Financial Economics 34, 7799. Sundaresan, Suresh, 1991, Valuation of swaps, in Recent Development in International Banking (Elsevier, North Holland). Vasicek, Oldrich, 1977, An equilibrium characterization of the term structure, Journal of Finan cial Economics 5. 177188.
Comments