Yield Curve Construction: Draft 2
|
|
- Ira Taylor
- 5 years ago
- Views:
Transcription
1 Yield Curve Construction: Draft 2 Stefanus Lie Risk Management Institute National University of Singapore Abstract In this work, several combination of models will be explored. Analysis, comparison, and stress testing will firstly be done for set of government bonds, before possible future implementation for policy banks and corporate bonds. Construction of our own yield curves could be useful for two reasons: Firstly, satisfy the desire to remove black-box part from Chinabond yield curve, for example in bonds valuation using Hull-White model. Secondly, to be implemented in future pricing project, or even to make it a commercial product. This Draft 2 contribution to the literature is comprehensive analysis, comparison, and stress testing among different combinations of filtering, models, knot points, and weight (these terms will be explained later) in Chinese bond market. Keywords: Spot Rate, Forward Rate, Discounting Rate, Spline, Filtering, Smoothness, Stress Testing, Knot Points, Global Optimization Preprint submitted to Risk Management Institute December 16, 2016
2 1. Preliminaries 1.1. Introduction and Notation Define s(t) to be t-year spot rate, f(t) to be t-year forward rate, and d(t) to be t-year discounting factor (or price of risk-free zero-coupon bond paying 1 in t year), for all t 0. Recall the following relations: 1. d(t) = exp( ts(t)). 2. d(t) = exp( t f(u)du) f(t) = ts (t) + s(t). If we know either one of s(t), f(t), d(t), we know all. Term Structure of Interest Rates and Yield Curve of a class of bonds could be used interchangeably, and they may mean: 1. Yield to Maturity with respect to time to maturity t Spot rate with respect to time to maturity t Forward rate with respect to time to maturity t Discounting factor with respect to time to maturity t 0. If we know either one of (2), (3), or (4), we know all Government Bond Market in China Most of government bond in China are zero-coupon bond or fixed-coupon bond with no option. The term to maturity could be 0.25 years, 0.5 years, 0.75 years, 1 year, 2 years, 3 years, 5 years, 7 years, 8 years (only once in 1999), 9 years (only twice in 1983 and 1984), 10 years, 15 years, 20 years, 30 years, and 50 years. According to coupon payment structure, a bond could be classified as zero-coupon bond (if the issuance price is 100 and there is final payment of 100 plus coupon at maturity), discounted bond (if the issuance price is less than 100 and there is final payment of 100 at maturity), and coupon bond (if there are some coupon payments between issuance and maturity, besides payment of coupon at maturity). Here are the summary about the bond market: 2
3 Figure 1: Summary of government bonds in Chinese market Note that bonds are unique to its short name, not symbol (because symbols are different for different market: Interbank, Shanghai Exchange, Shenzhen Exchange, Bank Counter, and Others). Note also that for some bonds, after the first issuance, the government could issue the bonds with same property in primary market, and thus increase the amount outstanding. The second issuance (and so on) could have different symbol than the first bond, although it s on the same market. 3
4 1.3. Description of Data Currently, there are two sets of data in the database as basis of our project: 1. Static: short name, face value, symbol, market, issue end date, issue amount, carry date, maturity date, coupon frequency, coupon at issue, coupon description, interest reference, tax exemption, tax rate, coupon date description, subordinated, benchmark interest name, special provisions, option-embedded, option-embedded description, spread at issue, redemption date, min coupon adjustment, max coupon adjustment, benchmark method. 2. Dynamic: date, short name, market, open price, previous close price, close price, high price, low price, volume, amount, average, quotes, best bid clean price, best bid yield, best ask clean price, best ask yield, last trading day, settlement data (transaction, average, high, close, open). The blue variables are those that are currently used in the model. The static data is obtained by combining data from Wind Data Explorer (active, delisted, and matured bonds) to excel, then upload the corresponding data to database. This is only updated to 31 October The dynamic data is obtained from static data and Wind code generator in Python, directly uploaded to database. Hence, static and dynamic data are only updated to 31 October. The problem to our daily update is reduced to the problem of updating the static data automatically, i.e., how to know recently issued bonds? 1.4. Introduction to Models Two steps are crucial in building yield curve for a given day: 1. Filtering: Selection of data from market, i.e., which data of prices/yields should be included? 2. Modeling: Run the model that will generate yield curve from filtered data. The empirical model for yield curve construction, from my perspective, is always compromise between fitting current market data and imposing theoretical/ parametric assumptions. By strictly assuming that the constructed yield curve follows the current market data, we may have to relax some theoretical/ parametric assumptions (no-arbitrage, economic meaning, shape, 4
5 and smoothness), while by assuming some theoretical/parametric assumptions, we have to accept some noise in market data in the constructed yield curve and have less parameters. Call the extremes to be left and right extremes. It should be noted that, for pricing and hedging purposes, we should not get too much to the right, while for policy making, most countries uses right models which has less number of parametes, as in Bank of International Settlements (2005). Most common examples of right models are Nelson-Siegel (1987) - having only 4 parameters or less to be estimated - and Svennson (1994, 1995) - having only 6 paramaters or less to be estimated. 2. Literature Review In this sections, I will provide selected historical perspective of yield curve construction. By selected, I mean: Only contributions that are relevant to our modelling purpose will be discussed in more details Pre-1970s era Report from Durand (1942) and US Treasury Department (1966), among others, consist only of first type yield curve (yield to maturity vs time). The fitting methods are by hand-fitting: not quantitative and not objective. An attempt was made by Cohen, Kramer, and Waugh (1966), by using regression with one dependent variables and two independent variables. The dependent variable is yield to maturity or log of yield to maturity, the first independent variable is time to maturity of squared time to maturity, and the second independent variable is squared log of time to maturity. Yet, approximating YTM vs time to maturity directly is not rigorous, because of coupon effect: bonds with equal maturity but with different coupon could not be seen having the same YTM. The spot/forward/discounting curve should be more appropiate Introduction to Spline There are several definitions for spline, and for consistency, throughout this project, we use the following definition (if some paper by authors is quoted, we adjust their definition to this definition): Definition 1: For given n > 1 and increasing sequence t 1 < t 2 < < t k called knot points, a function f : [t 1, t k ] R is called n-spline if f (n) is constant on each intervals (t 1, t 2 ), (t 2, t 3 ),... (t k 1, t k ), and f (n 1) is continuous on [t 1, t k ]. 5
6 The following propositions could be shown by induction: Proposition 1: For n-spline defined above, set of n-spline is also a vector space over R with dimension n + k 1. Define this vector space to be V (n, t 1, t 2,..., t k ). Proposition 2: Let S be subset of {1, 2,... k}. Define W (S, n, t 1, t 2,..., t k ) = {f V (n, t 1, t 2,..., t k ) f(t s ) = 0 s S}. Then W is subspace of V with dim(w ) = n + k S 1. In particular, the space of 3-spline (cubic spline) f s such that f(t s ) = 0 s = 1, 2,... k has dimension 2. Proposition 3 (consequence of well-known result about vector space): Let S be subset of {1, 2,... k}, t 1 < t 2 < < t k be knot points, and f j : [t 1, t k ] R(1 j n + k S 1) be linearly independent n-spline corresponding to the knot points such that f j (t s ) = 0 s S. Then the set {f j 1 j k + 2 S } becomes the basis of vector space W (S, n, t 1, t 2,..., t k ) = {f V (n, t 1, t 2,..., t k ) f(t s ) = 0 s S}. There are several ways to define the basis of the spline. A popular one is B-spline, discussed comprehensively in DeBoor (1978). We could also use McCulloch definition of basis cubic spline. Steeley (1991) suggested that using B-spline may be preferrable if we want to impose some constraints of the curve: Say, limit the first/second derivative, impose monotonicity, etc Breakthrough by McCulloch(1971, 1975) Recall the discounting function d(t). In McCulloch (1971), d(t) is assumed to be 2-spline based on knot points 0 = t 1 < t 2 < < t k such that d(0) = 1. Therefore the family of function d 1 (d minus constant function 1) is a vector space with dimension 2 + k 1 1 = k, according to Proposition 2. This means, all possible d satisfying condition above is set of { k } 1 + a i f i (t) a i R i=1 where f i (t) : [0, t k ] R are k + 1 functions as a choice of basis. One of the disadvantage of assuming the discounting function to be 2-spline is that the forward curves will have some knuckles at knot points, since first derivative of d is not differentiable at knot points. In McCulloch (1975), d(t) is assumed to be cubic spline based on knot points 0 = t 1 < t 2 < < t k such that d(0) = 1. Therefore the family of 6
7 function d 1 (d minus constant function 1) is a vector space with dimension 3 + k 1 1 = k + 1, according to Proposition 2. This means, all possible d satisfying condition above is set of { k+1 } 1 + a i f i (t) a i R i=1 where f i (t) : [0, t k ] R are k + 1 functions as a choice of basis. McCulloch (1975) proposed a particular choice of (k + 1) basis functions. Now, it is assumed that n P = A j d(s j ) + (w(p )ɛ) for all bonds in the bond sample set. This is equivalent to P = P = j=1 j=1 j=1 n k+1 ) A j (1 + a i f i (s j ) + (w(p )ɛ) i=1 i=1 n k+1 ( n ) A j + a i A j f i (s j ) + (w(p )ɛ) j=1 The parameters need to be estimated are a 1, a 2, a k+1, and this could be estimated by weighted linear regression. McCulloch (1971) proposed that the weight is broker fee plus half of bid-ask price difference, using the quotes data. Now, the natural questions to be asked: How many knot points there should be? How to choose the knot points? Suppose that the ordered bond maturities in the filtered set are m 1 < m 2 < < m n. McCulloch (1975) suggested that knot points are chosen such that the number of knot points is around n 1. Let the number of knot points be k. The choice of knot points 0 = t 1 < t 2 <... t k = m n depends on the sequence m 1 < m 2 < < m n : Suppose we wish to compute t i for 1 i k, then define f(i) = n(i 1)/8. Denote g(i) and h(i) to be integer part and non-integer part of f(i), respectively. Then define t i = m g(i) + h(i)(m g(i)+1 m g(i) ) for 2 i k 1. This will ensure almost equal number of bonds with maturity in each intervals (t i, t i+1 ). We could see that by having more knot points, we will have more parameters to be estimated, and the model goes to left. Moreover, ensuring the 7
8 equal number of bonds in each intervals is to make the curve between adjacent knot points stable: not sensitive to price change/error in one particular bond for which the maturity is in the interval Contribution by Carleton & Cooper (1976): Direct discrete discounting estimation Carleton & Cooper (1976) assumes that price of coupon bonds could be seen as the sum of several zero-coupon bonds, and hence, linear combination of discounting factors. During 1970s, most of US Government Bonds had coupon payments only in Feb 15, May 15, Aug 15, and Nov 15. Thus, by selecting about 40 bonds maturing no later than 4 years, price of all bonds is linear combination of 16 discounting factors. Standard linear regression with bid, ask, and mean prices for each bond, are used to estimate the discounting factors. Afterwards, spot yield curve and forward yield curve are constructed by linear interpolation of forward rates between two discrete dates. This model is not implementable in our project, since the coupon payment date is irregular in Chinese market, so we have more variables than equations Pioneer of a good Parametric Model: Cooper (1977) Recall that Cohen, Cramer, Waugh (1966) use parametric model, and then regression to construct YTM yield curve. The seed of parametric model could be traced in Echols and Eliott (1976), who assumed that simple forward rate follows f(t) = A exp(bt). From this, it is implied that simple spot rate is linear combination of t, 1/t, and 1. By considering the coupon effect on spot rate, Cohen, Cramer, Waugh (1966) use regression of simple spot rate with dependent variables: 1, t, 1/t, c (coupon). However, this is still a raw estimation. A good contribution was made by Cooper (1977), who assumes that f(t) = A + (B A) exp( Ct), implying that s(t) = A + ((B A)/Ct)(1 exp( Ct)). The parameters A, B, C are estimated that give smallest residual standard errors in price Attempt to improve McCulloch (1971, 1975): Vasicek-Fong (1982) Vasicek-Fong (1982) criticized McCulloch s assumption that discounting factor is cubic spline, while actually it should be exponential-like, then suggested transformation of discounting function. For t 0, define x = 1 exp( αt), and then define G(x) = d(t), for 0 x < 1. This will make α to be long-term forward rate. Since d(0) = 1, d is decreasing, and d(t) 8
9 as t, G(0) = 1, G(1) = 0 (where G(1) is defined as left limit), and G is decreasing. Now, G is assumed to be cubic spline, with converted knot points 0 = x 1 < x 2 < < x k < x k+1 < 1, where x k = 1 exp( αt k ). Therefore the family of function G 1 + x (G minus constant function 1 plus function y = x) is a vector space with dimension 3 + (k + 1) 2 1 = k + 1. Now, all possible G satisfying conditions above is set of { k+1 } 1 x + a i f i (t) a i R i=1 where f i (t) : [0, 1] R are k + 1 functions as a choice of basis. The parameters need to be estimated are a i s. Define y j = 1 exp( αs j ). Now, it is assumed that n P = A j g(y j ) + (w(p )ɛ) for all bonds in the bond sample set. This is equivalent to P = P = j=1 n A j j=1 j=1 n k+1 ) A j (1 y j + a i f i (y j ) + (w(p )ɛ) j=1 i=1 i=1 n k+1 ( n ) y j A j + a i A j f i (y j ) + (w(p )ɛ) j=1 The parameters need to be estimated are a 1, a 2, a k+1, and this could be estimated by weighted linear regression. Vasicek & Fong (1982) suggested using elasticity (first derivative of price with respect to yield) as weight in linear regression. This paper doesn t provide any empirical test Extension to McCulloch (1971, 1975) and Vasicek & Fong (1982): Smoot (1983) Smoot (1983) in his Ph.D thesis proposes a new method of applying spline directly on spot rate. There are several layers of consideration in this paper: Knot placement, tax models, weight for price error, and the model itself (McCulloch/Vasicek/Smoot). After deciding to choose a particular knot placement and tax model, data from is tested using all combinations of weight and model. For each choice of weight, the lowest MSE is produced by using the new model: spline on spot rate. However, I doubt 9
10 that this result could be generalized to our project directly: Firstly, the data could be very different, in terms of time and place. Secondly, Smoot (1983) had already zoom in to a specific knot points and tax models. Therefore, we still consider McCulloch and Vasicek-Fong models for consideration. In the rest of this section, more explanation about this model will be given: Spot rate s(t) is assumed to be cubic spline based on knot points 0 = t 1 < t 2 < < t k. Therefore the family of spot function s(t) is a vector space with dimension 3 + k 1 = k + 1, according to Proposition 2. This means, all possible s satisfying condition above is set of { k+2 } a i f i (t) a i R Now, note that the discounting function is ( k+2 ) d(t) = exp( ts(t)) = exp t a i f i (t) j=1 j=1 i=1 Let there are m bonds in the sample set, denote P l (1 l m) to be its dirty prices, and denote n 1, n 2,..., n m to be number of future payment times of each bond. For l-th bond, let s l,1 < s l,2 < < s l,nl to be time to next coupon/final payments, and A l,1, A l,2,..., A l,nl to be corresponding payment amounts. Note that if we assume spot rate as above, the price of l-th bond should be: n l n l k+2 ) n l k+2 A l,j d(s l,j ) = A l,j exp ( s l,j a i f i (s l,j ) = A l,j exp( s l,j f i (s l,j )) a i Therefore, our aim to minimize sum of weighted square of error: ( m ( nl k+2 ) A l,j exp( s l,j f i (s l,j )) a i P l w(p l ) l=1 j=1 i=1 i=1 where w(p l ) is weight associated with l-th bond Chambers, Carleton, Waldman (1984) The spot rate is assumed to be polynomial: s(t) = J x j t j 1 j=1 10 i=1 j=1 i=1 ) 2
11 This paper do empirical testing on using 1, 2, 3, 4, 5 degree polynomial, minimizing MSE by nonlinear regression with equal error weight. Futhermore, by arguing that MSE-minimizing do not produce good results (heterodascity), maximum likelihood methods are proposed considering power of time to maturity as weight of variance. In my opinion, the contribution of this paper is not much. Assuming spot rate to be polynomials everywhere may not be good since it may not capture the bend in longer term, contrast to spline methods Nelson-Siegel (1987) and Svennson(1994, 1995) These models are more extreme to the right. NS model assumes that the forward rate follows: f(t) = β 0 + β 1 exp( t/τ) + β 2 (t/τ) exp( t/τ) The first term is relevent for long-term, the second term for short term, and the third term for the medium (hump) term. NSS model is generalization of NS model: f(t) = β 0 + β 1 exp( t/τ 1 ) + β 2 (t/τ 1 ) exp( t/τ 1 ) + β 3 (t/τ 2 ) exp( t/τ 2 ) The last term would be for the second hump Yield Curve Smoothness Concept by Adams and van Deventer (1994) This paper is the pioneer in applying smoothness concept to the forward curve, borrowing the concept from numerical analysis, i.e. a curve f is smooth if T (f (s)) 2 ds is relatively small. Given discounting at discrete 0 points d(t 1 ), d(t 2 ),..., d(t k ), a closed formula is for forward curve is proposed such that: (1) It match with the discounting at discrete points; (2) Among all that satisfy (1), the curve is the smoothest one. The computation in this paper is further corrected by Lim & Xiao (2002). Although this paper only suggested solution given discounting at discrete points, which could only be relevant for Carleton & Cooper (1976), the way to quantify the smoothness would be useful and used in subsequent papers Fisher, Nychka, and Zervos (1994): Estimating and Smoothing simultaneously While Adams and van Deventer (1994) proposes smoothing after estimating the discounting at several discrete points, Fisher, Nychka, and Zervos (1994) is the first to combine the estimation and smoothness. This paper 11
12 listed three possible estimation: Spline with smoothing on discounting, spline with smoothing on log discounting, and spline with smoothing on forward rates. This entity is to be minimized: Sum of squared difference between model and observed price (with equal weight), in addition to λ T 0 (g (s)) 2 ds, where g could be discounting, log discounting, or forward curves. λ is chosen using Generalized Cross Validation formula. When trying to apply the three methods to data from , this paper says that fitting forward curve provides the most accurate fit. However, this result is very limited by the choice of data set and, more importantly, using only equal weight for price error. Therefore, we could still consider those three methods (the first one is originally from McCulloch (1971, 1975)). On the rest of the section would be the implementation details for cubic spline of forward rates and on log discounting. Spline forward rate: Forward rate f(t) is assumed to be cubic spline based on knot points 0 = t 1 < t 2 < < t k. Therefore the family of forward function f(t) is a vector space with dimension 3 + k 1 = k + 1, according to Proposition 2. This means, all possible f satisfying condition above is set of { k+2 } a i f i (t) a i R i=1 Define g i (t) = t 0 f i(u)du. Now, note that the discounting function is ( d(t) = exp t 0 ) ( k+2 ) f(u)du = exp a i g i (t) Let there are m bonds in the sample set, denote P l (1 l m) to be its dirty prices, and denote n 1, n 2,..., n m to be number of future payment times of each bond. For l-th bond, let s l,1 < s l,2 < < s l,nl to be time to next coupon/final payments, and A l,1, A l,2,..., A l,nl to be corresponding payment amounts. Note that if we assume forward rate as above, the price of l-th bond should be: n l j=1 n l A l,j d(s l,j ) = j=1 ( k+2 A l,j exp i=1 ) a i g i (s l,j ) = i=1 n l k+2 A l,j j=1 i=1 exp( g i (s l,j )) a i 12
13 Therefore, our aim to minimize sum of weighted square of error: ( m ( nl l=1 j=1 A l,j k+2 exp( g i (s l,j )) a i P l )/w(p l ) i=1 where w(p l ) is weight associated with l-th bond. Spline of log discounting: c(t) = log d(t) is assumed to be cubic spline based on knot points 0 = t 1 < t 2 < < t k such that c(0) = 0. Therefore the family of function c is a vector space with dimension 3+k 1 1 = k +1, according to Proposition 2. This means, all possible d satisfying condition above is set of { k+1 } a i c i (t) a i R i=1 where f i (t) : [0, t k ] R are k + 1 functions as a choice of basis. Now, it is assumed that n P = A j exp(c(s j )) + (w(p )ɛ) j=1 for all bonds in the bond sample set. This is equivalent to P = n ( k+1 ) A j exp(c i (s j )) a i + (w(p )ɛ) j=1 i=1 Therefore, our aim to minimize sum of weighted square of error: ( m ( n ( k+1 ) ) A j exp(c i (s j )) a i + (w(p )ɛ) P l w(p l ) l=1 j=1 i=1 where w(p l ) is weight associated with l-th bond Waggoner(1997) Waggoner generalized the smoothness penalty to T 0 λ(s)(g (s)) 2 ds. In the paper, it is suggested that λ(s) = 0.1 for 0 s 1, λ(s) = 100 for 1 s 100, and λ(s) = for 10 s, forcing more smoothness on longer term and less smoothness on shorter term. For details on the reasoning, see the paper. 13 ) 2 ) 2
14 2.13. Lapshin & Wang (2013) This paper proposes a method to fit forward curve such that it is always positive (no-arbitrage), so in this case f(t) = g(t) 2, and g(t) is the curve to be estimated, minimizing weighted sum of square error, in addition to smoothness penalty λ T 0 (g (s)) 2 ds. The result is tested on Chinese market and compared with Chinabond yield curve A Short Reflection Filtering data is a very market-contextual and subjective task, while modelling filtered data to yield curve is a more discussable task: Many authors propose model, others accept/reject by theoretical/empirical results, others improve or combine models in the past, etc. Currently (up to 16 December 2016), we will zoom in to implementation for cubic spline methods without smoothness penalty. It could be summarized in one sentence with four blanks: After filtering the data with (criteria), to construct yield curve, we assume cubic spline on (model), with choice of knot points to be (knot method), minimizing sum of squared difference of observed and model price weighted by (weight method) Regarding error weight, and knot points The choice for error weight also depends on the sample data. In our sample, there are data from volume and quotes. Therefore, using weight in McCulloch (bid-ask price difference) is not available. Using equal weight is standard. Using price as weight is suggested in Litzenberger (1982). Using elasticity or duration is suggested in Vasicek & Fong (1982), and several countries central bank in BIS Report Implementation and Analysis: Part 1 Following previous discussions, the choices are: 1. Criteria: 1, 2, 3 2. Model: McCulloch (spline on discounting), Vasicek-Fong (exponential spline on discounting) 3. Knot method: McCulloch (in short, mc), fixed following Chinabond, custom 14
15 4. Weight method: equal, price, duration, elasticity Let m 1 < m 2 < < m n be list of time to maturity of all bonds in the bond sample, for a given day. There are about 100 bonds in bond sample, for a given day. Regarding the choice of knot points, currently, there are three alternatives: 1. McCulloch (1975) suggestion: Knot points are chosen such that the number of knot points is around n 1. Therefore, we assume that the number of knot points is 9. The choice of knot points 0 = t 1 < t 2 <... t 9 = m n depends on the sequence m 1 < m 2 < < m n : Suppose we wish to compute t i for 1 1 9, then define f(i) = n(i 1)/8. Denote g(i) and h(i) to be integer part and non-integer part of f(i), respectively. Then define m g(i) + h(i)(m g(i)+1 m g(i) ). This will ensure almost equal number of bonds with maturity in each intervals (t i, t i+1 ). 2. Follow Chinabond key terms, so that the knot points are fixed at 0, 0.167, 0.25, 0.5, 0.75, 1, 2, 3, 5, 7, 10, 15, 20, 30 years. 3. Use custom fixed knot points. For criteria 1, use 0, 0.2, 0.5, 2, 4, 5, 7, 9, 50, while for criteria 3, use 0, 0.167, 0.5, 2, 3.5, 4.5, 6, 8.5, 50 years. By using equal weight, we minimize (P P ) 2 By using price as weight weight, we minimize (P P ) 2 P By using duration as weight, we minimize (P P ) 2 D(P ) By using elasticity as weight, we minimize Note that (P P ) 2 E(P ) E(P ) = P/ y D(P ) = P ( P/ y) 15
16 3.1. Metric Function In this section, metric functions are introduced to have quantitative measure of distance between function. Definition 2: For two functions f, g : [a, b] R and [c, d] [a, b], define metric of f and g on [c, d] to be m [c,d] (f, g) = d c f(x) g(x) dx Definition 3: For two functions f, g : [a, b] R and [c, d] [a, b], define adjusted metric of f and g on [c, d] to be am [c,d] (f, g) = m [c,d] (f, g)/(d c) Definition 4: For two functions f, g : [a, b] R and [c, d] [a, b], define max-metric of f and g on [c, d] to be 3.2. Regarding the Criteria Criteria 1: mm [ c, d](f, g) = max{f(x) g(x) x [c, d]} 1. Suppose that we want to build yield curve for a given transaction day. Collect all zero-coupon/fixed-coupon bonds with no option from set of active bonds in the category (government/adbc/cdb/export- Import) in that day. 2. For each bond in the collection, if the day is issuance end date of the bond, define the value of the bond to be the issue price in Interbank or Exchange Market; 3. Otherwise, if there is volume in Interbank Market, define the value of the bond to be the average transaction dirty price; 4. Otherwise, if there are quotes in Interbank Market, define the value of the bond to be the average of best bid and best ask dirty price, where the computation of accrued interest assumes settlement on the day of quotes; 5. Otherwise, remove the bond from our collection. 6. In the end, we have collection of bonds with corresponding value. This will serves as daily input in our model. 16
17 Now, we try to do robustness test on current criteria, i.e., how stable is the daily set of bonds with respect to the criteria? The testing period is May 2016 to October There are 125 trading days in this period, with one anomaly on 12 June 2016, where there are only about 20 bonds with trading/quotes recorded on Wind, possibly because quotes data are missing on that day. Figure 2: Criteria 1 samples Now, define B 1, B 2,... B 125 be be set of bonds included in 1st, 2nd,..., 125th transaction day in the testing period. Define T NY i = B i \ B i 1 for 2 i 125 (number of bonds in i-th day but not in (i 1)-th day). Furthermore, define T NT i = B i \ B i+1 for 1 i 124 (number of bonds in i-th day but not in (i + 1)-th day). Because of one anomaly, there will be one extreme outliers on each sequences T NY i and T NT i. We take average of sequence T NY i, resulting in 3.179, and for T NT i is This is far above average number of new issuance per day which is This means, by comparing set of bonds on a given day to previous day, about three bonds have been excluded and three new bonds have been included. By looking at Criteria 1, four disadvantages could be seen: 17
18 It may include an illiquid bond if there is volume on a given day, although there had been no trading for long time before. It may exclude a liquid bond if there is no volume/quotes on a given day, although there had been trading for long time before. It may include bonds with price anomaly: Large bid-ask price difference, or suddenly-jump transaction price. From observation, there is usually no volume/quotes record on the transaction day after the issuance end time, although afterwards, the bond may be very liquid. Define Criteria 2 as follows: 1. Suppose that we want to build yield curve for a given transaction day. Collect all zero-coupon/fixed-coupon bonds with no option from set of active bonds in the category (government/adbc/cdb/export- Import) in that day. 2. For each bond in the collection, if the day or previous trading day is issuance end date of the bond, define the value of the bond to be the corresponding dirty price of issue price in Interbank or Exchange Market; 3. Otherwise, if there is volume in Interbank Market on that day and previous day and (one of previous-previous or previous-previous-previous day), define the value of the bond to be the average transaction dirty price; 4. Otherwise, if there are quotes in Interbank Market, define the value of the bond to be the average of best bid and best ask dirty price, where the computation of accrued interest assumes settlement on the day of quotes; 5. Otherwise, remove the bond from our collection. 6. In the end, we have collection of bonds with corresponding value. This will serves as daily input in our model. 18
19 Figure 3: Criteria 2 samples We now compute the average (removing extreme outlier) of sequence T NY i, resulting in 1.236, and for T NT i is The average number of new issuance per day which is We could see that the robustness improved compared to Criteria 1. Define Criteria 3 as follows: 1. Suppose that we want to build yield curve for a given transaction day. Collect all zero-coupon/fixed-coupon bonds with no option from set of active bonds in the category (government/adbc/cdb/export- Import) in that day. 2. For each bond in the collection, if the day or previous trading day is issuance end date of the bond, define the value of the bond to be the corresponding dirty price of issue price in Interbank or Exchange Market; 3. Otherwise, if there is volume in Interbank Market on that day and previous day and (one of previous-previous or previous-previous-previous day), define the value of the bond to be the average transaction dirty price; 4. Otherwise, if there are quotes in Interbank Market for which the bidask yield spread not exceeding 50 bps, define the value of the bond to be 19
20 the average of best bid and best ask dirty price, where the computation of accrued interest assumes settlement on the day of quotes; 5. Otherwise, remove the bond from our collection. 6. In the end, we have collection of bonds with corresponding value. This will serves as daily input in our model. Figure 4: Criteria 3 samples We now compute the average (removing extreme outlier) of sequence T NY i, resulting in 1.667, and for T NT i is The average number of new issuance per day which is We could see that the robustness improved compared to Criteria 1. Furthermore, stress test is done on the yield curve, using a hypothetical, but realistic condition. The question to be asked is: Suppose three bonds (one with maturity about 0.01, one with maturity about 1.2, and one with maturity about 6) are removed from the set of bonds using Criteria 1 on a given day. The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate derived by removing 3 bonds (all figures are in bps): 20
21 Figure 5: Criteria 1 remove 3 bonds maturity about 0.01, 1.2, 6 From observations: 1. When removing 3 bonds from the sample list, the shift will be relatively local and not that much (less than 1 bps). Theoretical explanation: Sample size is large 2. The shift effects become even smaller when using duration or elasticity as weight (less than 0.5 bps). Theoretical explanation: More weight put into shorter term bonds, and there are many of them, so removing 3 will not affect much. 3. The shift effects become relatively smaller when using custom or mc as knot methods. Theoretical explanation: using fixed knot methods, sparse data of maturities in one interval make it sensitive. Furthermore, stress test is done on the yield curve, using a hypothetical, but realistic condition. The question to be asked is: Suppose one bond with maturity around 0.5 year are removed from the set of bonds using Criteria 3 on a given day. The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate derived by removing 1 bond (all figures are in bps): 21
22 Figure 6: Criteria 3 remove 1 bond maturity about 0.5 From observations: 1. When removing 1 bond from the sample list, the shift will be relatively local and not that much (less than 0.5 bps). Theoretical explanation: Sample size is large. 2. The shift effects become even smaller when using duration or elasticity as weight. Theoretical explanation: More weight put into shorter term bonds, and there are many of them, so removing 3 will not affect much. 3. Theoretically, since the input is more stable for Criteria 3, the yield curve will be more stable. Comparison of Criteria 1 vs Criteria 3: The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate when changing from Criteria 1 to Criteria 3: Figure 7: Criteria 1 vs Criteria 3 From observations: 22
23 1. Changing the Criteria 1 to Criteria 3 shifts the curve much especially in for maturity larger than 10 years. Theoretical explanation: Little trading/quotes data for long-term bonds, some may be outlier addressed in Criteria Regarding the knot methods Firstly, let us see how many bond s maturities in each intervals of Chinabond knot points, for each Criteria. The table is taken by taking average from May to Oct 2016: Figure 8: Average number of maturities in Chinabond knot points We could see that some intervals are very sparse: 0-2 data available, making the curve on that intervals very sensitive to the sample data. This is one of disadvantange of adopting Chinabond fixed knot points, and could be solved by McCulloch knot points: choosing knot points such that there are almost equal number of bond s maturities on each intervals. Let me start with giving example on how Mcculloch knot method work. On 31 Oct 2016, this is the distribution of maturities on bonds in the sample: Figure 9: Distribution 0-20 years Figure 10: Distribution years In this case, the knot points would be: 0, 0.23, 0.863, 2.479, 3.693, 4.978, 6.559, 9.51, If the distribution of maturities change, the knot points will also change. So the knot points would change daily. Possibly this makes 23
24 the curve using Mcculloch knot points unstable from day to day (although this needs further research). I am then motivated to define custom knot points for each criteria, such that it is fixed and there are almost equal number of bond s maturity on each intervals. For criteria 1, use 0, 0.2, 0.5, 2, 4, 5, 7, 9, 50, while for criteria 3, use 0, 0.167, 0.5, 2, 3.5, 4.5, 6, 8.5, 50 years. Here are the average for May-Oct 2016: Figure 11: Custom knot points We use 9 knot points in McCulloch choice of knot points currently. The question now: What if we change to 10 points? The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate when changing from 9 knot points to 10 knot points: Figure 12: Stress Testing McCulloch knot points from 9 to 10 24
25 From observations: 1. The shift is quite large for 0-1 year, and moderate for other intervals. 2. The shift becomes very small using duration or elasticity as weight. How does three knot methods compare to each other? The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate, for pairwise comparison among Chinabond fixed, Mcculloch, and custom knot methods: 25
26 Figure 13: Comparison among 3 knot methods From observations: 1. Moderate to large shift when changing from fixed to mc or fixed to custom. 2. Low to moderate shift when changing from mc to custom. 3. The shift become small when using duration or elasticity as weight. 26
27 3.4. Regarding the weight methods By giving more weight for a bond, we allow more uncertainty/noise of the bond prices. Let us look at bid-ask price difference for quotes recorded in Criteria 1 on one day. The trend is: As maturity increases, price spread increases, allowing for more uncertainty. Figure 14: Bid-ask price difference vs time Let us look at bid-ask yield difference for quotes recorded in Criteria 1 on one day. There is no trend that could be seen here. Figure 15: Bid-ask yield difference vs time Therefore, theoretically, using duration or elasticity as weight is better 27
28 than using equal or price as weight. Indeed, by using elasticity, the error, neglecting convexity effect, is approximately the error term in the regression will be approximately yield error. The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate, for pairwise comparison equal vs duration, equal vs price, and duration vs elasticity: Figure 16: Weight comparisons From observations: 1. Not much shift when changing from using equal to price as weight, or duration to elasticity as weight, except on some intervals in Criteria Much shift when changing from equal to duration, especially for shorter term. 28
29 3.5. Stress Testing alpha in Vasicek-Fong What if we change long-term forward rate α in Vasicek-Fong from to or 0.040? The testing period is October 2016Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate, for pairwise comparison between α vs and vs 0.040: From observations: Figure 17: Stress Testing alpha 1. In all cases, choosing different alpha will not change the curve much Regarding the model How does Mcculloch compare to Vasicek-Fong? The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate, for McCulloch vs Vasicek-Fong: 29
30 Figure 18: Models comparison From observations: 1. Changing mcculloch to vasicek-fong method will not shift the curve much, but the shift become larger as maturities become larger. Theoretical explanation: alpha (long-term forward rate) will pull the curve from infinity. 2. Advantage of vasicek-fong: We could determine long-term forwards rates, and this choice of alpha will not be sensitive, until 30 years Smooothness The testing period is October Here are the summary of results of smoothness calculated by Adams & van Deventer (1994) measure: 30
31 From observations: Figure 19: Smoothness comparison 1. More smooth when using criteria More smooth when using duration or elasticity as weight. 3. More smooth when using mc or custom as knot methods Conclusions From observations, using equal or price as weight could be rejected (with theoretical and empirical reasons). Furthermore, using fixed Chinabond knot points will also be rejected (with theoretical and empirical reasons). Although it could be seen that Criteria 3 improves from Criteria 1, Criteria 3 also has disadvantages: It may exclude data from volume but still including the quotes. We will still consider using Criteria 1 and 3 in further testing. 4. Implementation and Analysis: Part 2 Following previous discussions, the choices are: 31
32 1. Criteria: 1, 3 2. Model: McCulloch (spline on discounting), Vasicek-Fong (exponential spline on discounting), Smoot (spline of spot), spline on log discounting, spline on forward. 3. Knot method: McCulloch (in short, mc), custom according to the criteria 4. Weight method: duration, elasticity 4.1. Regarding the Criteria Stress test is done on the yield curve, using a hypothetical, but realistic condition. The question to be asked is: Suppose three bonds (one with maturity about 0.01, one with maturity about 1.2, and one with maturity about 6) are removed from the set of bonds using Criteria 1 on a given day. The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate derived by removing 3 bonds (all figures are in bps): Figure 20: Criteria 1 remove 3 bonds maturity about 0.01, 1.2, 6 From observations: 1. When removing 3 bonds from the sample list, the shift will be relatively local and not that much (less than 1 bps). Theoretical explanation: Sample size is large. Suppose one bond with maturity around 0.5 year are removed from the set of bonds using Criteria 3 on a given day. The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate derived by removing 1 bond (all figures are in bps): 32
33 Figure 21: Criteria 3 remove 1 bond maturity about 0.5 From observations: 1. When removing 1 bond from the sample list, the shift will be relatively local and not that much (less than 0.3 bps). Theoretical explanation: Sample size is large. 2. Theoretically, since the input is more stable for Criteria 3, the yield curve will be more stable. Comparison of Criteria 1 vs Criteria 3: The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate when changing from Criteria 1 to Criteria 3: Figure 22: Criteria 1 vs Criteria 3 From observations: 1. Changing the Criteria 1 to Criteria 3 shifts the curve much especially in for maturity 10 years. Theoretical explanation: Little trading/quotes data for long-term bonds, some may be outlier addressed in Criteria 3. 33
34 4.2. Regarding the knot methods We use 9 knot points in McCulloch choice of knot points currently. The question now: What if we change to 10 points? The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate when changing from 9 knot points to 10 knot points: Figure 23: Stress Testing McCulloch knot points from 9 to 10 From observations: 1. The shift is moderate to small in all cases How does three knot methods compare to each other? The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate, for pairwise comparison among Chinabond fixed, Mcculloch, and custom knot methods: 34
35 Figure 24: Comparison among 3 knot methods From observations: 1. For Criteria 1, moderate shift for short term and small shift for longer term when changing knot methods from mc to custom. 2. For Criteria 3, small shift for all term when changing knot methods from mc to custom. 3. Using custom could appeal because of knot points stability on daily basis, but this need further research Regarding the weight methods The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate, using weight duration vs elasticity: 35
36 Figure 25: Weight comparisons From observations: 1. Not much shift when changing from elasticity to duration as weight, except for longer maturity in Criteria Stress Testing alpha in Vasicek-Fong What if we change long-term forward rate α in Vasicek-Fong from to or 0.040? The testing period is October 2016Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate, for pairwise comparison between α vs and vs 0.040: 36
37 Figure 26: Stress Testing alpha From observations: 1. In all cases, choosing different alpha will not change the curve much Regarding the model How does Mcculloch compare to Vasicek-Fong? The testing period is October Here are the summary of results of adjusted metric and maximum metric for different intervals of spot rate, for McCulloch vs Vasicek-Fong: 37
38 Figure 27: Models comparison 38
39 Figure 28: Models comparison From observations: 1. Not much shift when using either mc, vf, logdisc. 2. Larger shift for longer maturity when using forward model. 39
40 4.6. Stress Testing by chaning YTM of one bond ±1% The testing period is Oct Here are the average shift when changing YTM of one bond with maturity about 0.2 years, by one percent: Figure 29: YTM one bond (maturity 0.2 years) plus one percent Here are the average shift when changing YTM of one bond with maturity about 0.2 years, by minus one percent: Figure 30: YTM one bond (maturity 0.2 years) minus one percent 40
41 Here are the average shift when changing YTM of one bond with maturity about 4 years, by one percent: Figure 31: YTM one bond (maturity 4 years) plus one percent Here are the average shift when changing YTM of one bond with maturity about 4 years, by minus one percent: Figure 32: YTM one bond (maturity 4 years) minus one percent Here are the average shift when changing YTM of one bond with largest maturity in sample, by one percent: 41
42 Figure 33: YTM one bond (maturity largest in sample) plus one percent Here are the average shift when changing YTM of one bond with maturity about 0.2 years, by minus one percent: Figure 34: YTM one bond (maturity largest in sample) minus one percent From observations: 1. In all cases, stressing YTM of a sample only affects the curve locally. 42
43 4.7. Smooothness The testing period is October Here are the summary of results of smoothness calculated by Adams & van Deventer (1994) measure: Figure 35: Smoothness comparison From observations: 1. Less smoothness when using smoot or forward may be inconclusive, could be affected by shorter term bend. Use Waggoner (1997) for further research. 2. More smooth when using criteria Conclusions There are tradeoffs of using McCulloch knot points and custom knot points. Using Criteria 3 will exclude outlier in volume, but still include quotes data. Further improvements could be made on the Criteria. Not much difference when using duration vs elasticity. 43
44 5. Discussions and Future Work In the future, the criteria could be improved, using ranking process as in Chinamoney/SSI yield curves. Furthermore, smoothness penalty could be implemented for some curves, but the optimization would be non-trivial anymore and would take longer time. Imposing some conditions could also be done: Assuming spot/forward rate to be fixed at several points. Or, we could try to do the same test for policy banks/corporate bonds first? Exploration for other models could also be done. By comparing price using the model curve and Chinabond yield curve in the future, we could learn why there is consistent difference between both of them Optimization Problem In this section, the following optimization problem is stated, and will be used repeatedly afterwards.given positive integers n 1, n 2,..., n m and real numbers P 1, P 2,... P m. Besides these, given positive real numbers a i,j,l for 1 i n l, 1 l m and 0 j k, and we try to find real numbers (β 1, β 2,..., β k ) that minimize m n l ( a i,0,l a β 1 i,1,l aβ 2 i,2,l aβ k i,k,l P l) 2 l=1 i=1 Currently, this optimization problem is solved using scipy.optimize.minimize, in particular, Newton-Conjugate-Gradient algorithm. Currently, this method could not give confidence that the optimized point would be global optimum (using the above method with zero vector as starting point), although when some cases are tried using different starting points, it indeed converge to the same point. Further research need to be done on this issue How Chinabond construct the spot yield curve? Lapshin and Wang (2013) and other sources noted that there are 4 steps: 1. Filtering of data on several bonds with maturity near key terms: Removing outliers from quotes significantly different from history, etc. 2. Augmenting trading data with expert estimates and sometimes historical data. 44
45 3. Key yields are chosen for interpolation, then estimate YTM at this key terms. Then there will be approximation of spot rate from YTM at those key terms. 4. Use monotone cubic spline interpolation (this will be unique) from key points before to obtain the whole curve. The way they estimate YTM first rather that estimating the spot rates directly is a non-rigorous step. Their advantage is having more data (for example, details of market maker and transaction/quotes time). 6. Bibliography [1] Adams, K.J., van Deventer, D.R., Fitting Yield Curves and Forward Rate Curves with Maximum Smoothness. Journal of Fixed Income, June [2] Bank of International Setllements, Zero-Coupon Yield Curves: Technical Documentation. BIS Papers No 25, October [3] Carleton, W.T., Cooper, I.A., Estimation and Uses of the Term Structure of Interest Rates. The Journal of Finance, Vol XXXI, No 4, September [4] Chambers, D.R., Carleton, W.T., Waldman, D.W., A New Approach to the Estimatiomn of the Term Structure of Interest Rates. Jpurnal of Financial and Quantitative Analysis, Vol 19, No 3, September [5] China Central Depository and Clearing Co., Ltd., Chinabond Pricing System [6] Cohen, K.J., Kramer, R.L., Waugh, W.H., Regressin Yield Curves for U.S. Government Securities. Management Science, December [7] Cooper, I.A., Asset Values, Interest-Rate Changes, and Duration. The Journal of Financial and Quantitative Analysis, Vol 12, No 5, December [8] de Boor, C., A Practical Guide to Splines. Applied Mathematical Sciences Volume 27, first edition 1978, second edition [9] Durand, D., Basic Yields of Corporate Bonds, National Bureau of Economic Research, June
Estimating A Smooth Term Structure of Interest Rates
E STIMATING A SMOOTH LTA 2/98 TERM STRUCTURE P. 159 177 OF INTEREST RATES JARI KÄPPI 1 Estimating A Smooth Term Structure of Interest Rates ABSTRACT This paper extends the literature of the term structure
More informationChinese Bond Market: A Need for Sound Estimation of Term Structure Interest Rates
World Applied Sciences Journal 4 (3): 358-363, 3 ISSN 88-495 IDOSI Publications, 3 DOI:.589/idosi.wasj.3.4.3.35 Chinese Bond Market: A Need for Sound Estimation of Term Structure Interest Rates Victor
More informationProblems and Solutions
1 CHAPTER 1 Problems 1.1 Problems on Bonds Exercise 1.1 On 12/04/01, consider a fixed-coupon bond whose features are the following: face value: $1,000 coupon rate: 8% coupon frequency: semiannual maturity:
More informationSpline Methods for Extracting Interest Rate Curves from Coupon Bond Prices
Spline Methods for Extracting Interest Rate Curves from Coupon Bond Prices Daniel F. Waggoner Federal Reserve Bank of Atlanta Working Paper 97-0 November 997 Abstract: Cubic splines have long been used
More informationEstimating Maximum Smoothness and Maximum. Flatness Forward Rate Curve
Estimating Maximum Smoothness and Maximum Flatness Forward Rate Curve Lim Kian Guan & Qin Xiao 1 January 21, 22 1 Both authors are from the National University of Singapore, Centre for Financial Engineering.
More informationIndian Sovereign Yield Curve using Nelson-Siegel-Svensson Model
Indian Sovereign Yield Curve using Nelson-Siegel-Svensson Model Of the three methods of valuing a Fixed Income Security Current Yield, YTM and the Coupon, the most common method followed is the Yield To
More informationWorking paper. An approach to setting inflation and discount rates
Working paper An approach to setting inflation and discount rates Hugh Miller & Tim Yip 1 Introduction Setting inflation and discount assumptions is a core part of many actuarial tasks. AASB 1023 requires
More informationImproving Nelson-Siegel term structure model under zero / super-low interest rate policy
Improving Nelson-Siegel term structure model under zero / super-low interest rate policy July 14th, 2015 Koji Inui School of Interdisciplinary Mathematical Sciences, Meiji University 4-21-1 Nakano Nakano-ku,
More informationWe consider three zero-coupon bonds (strips) with the following features: Bond Maturity (years) Price Bond Bond Bond
15 3 CHAPTER 3 Problems Exercise 3.1 We consider three zero-coupon bonds (strips) with the following features: Each strip delivers $100 at maturity. Bond Maturity (years) Price Bond 1 1 96.43 Bond 2 2
More informationInstantaneous Error Term and Yield Curve Estimation
Instantaneous Error Term and Yield Curve Estimation 1 Ubukata, M. and 2 M. Fukushige 1,2 Graduate School of Economics, Osaka University 2 56-43, Machikaneyama, Toyonaka, Osaka, Japan. E-Mail: mfuku@econ.osaka-u.ac.jp
More informationINTERPOLATING YIELD CURVE DATA IN A MANNER THAT ENSURES POSITIVE AND CONTINUOUS FORWARD CURVES
SAJEMS NS 16 (2013) No 4:395-406 395 INTERPOLATING YIELD CURVE DATA IN A MANNER THAT ENSURES POSITIVE AND CONTINUOUS FORWARD CURVES Paul F du Preez Johannesburg Stock Exchange Eben Maré Department of Mathematics
More informationLiquidity Premium in Emerging Debt Markets
Liquidity Premium in Emerging Debt Markets Abstract Developed markets are currently beset with credit risk though there is not much of a liquidity risk in these markets. However, it is the other way round
More informationThe Term Structure of Expected Inflation Rates
The Term Structure of Expected Inflation Rates by HANS-JüRG BüTTLER Swiss National Bank and University of Zurich Switzerland 0 Introduction 1 Preliminaries 2 Term Structure of Nominal Interest Rates 3
More informationSmooth estimation of yield curves by Laguerre functions
Smooth estimation of yield curves by Laguerre functions A.S. Hurn 1, K.A. Lindsay 2 and V. Pavlov 1 1 School of Economics and Finance, Queensland University of Technology 2 Department of Mathematics, University
More informationThe Duration Derby: A Comparison of Duration Based Strategies in Asset Liability Management
The Duration Derby: A Comparison of Duration Based Strategies in Asset Liability Management H. Zheng Department of Mathematics, Imperial College London SW7 2BZ, UK h.zheng@ic.ac.uk L. C. Thomas School
More informationForward Rate Curve Smoothing
Forward Rate Curve Smoothing Robert A Jarrow June 4, 2014 Abstract This paper reviews the forward rate curve smoothing literature The key contribution of this review is to link the static curve fitting
More informationEstimating Term Structure of U.S. Treasury Securities: An Interpolation Approach
Estimating Term Structure of U.S. Treasury Securities: An Interpolation Approach Feng Guo J. Huston McCulloch Our Task Empirical TS are unobservable. Without a continuous spectrum of zero-coupon securities;
More informationAn investigation into popular methods for constructing yield curves
An investigation into popular methods for constructing yield curves by Paul Fourie du Preez Submitted in partial fulfillment of the requirements for the degree Magister Scientiae in the Department of Mathematics
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationIntroduction to Bonds The Bond Instrument p. 3 The Time Value of Money p. 4 Basic Features and Definitions p. 5 Present Value and Discounting p.
Foreword p. xv Preface p. xvii Introduction to Bonds The Bond Instrument p. 3 The Time Value of Money p. 4 Basic Features and Definitions p. 5 Present Value and Discounting p. 6 Discount Factors p. 12
More informationThe Fixed Income Valuation Course. Sanjay K. Nawalkha Gloria M. Soto Natalia A. Beliaeva
Interest Rate Risk Modeling The Fixed Income Valuation Course Sanjay K. Nawalkha Gloria M. Soto Natalia A. Beliaeva Interest t Rate Risk Modeling : The Fixed Income Valuation Course. Sanjay K. Nawalkha,
More informationFX Smile Modelling. 9 September September 9, 2008
FX Smile Modelling 9 September 008 September 9, 008 Contents 1 FX Implied Volatility 1 Interpolation.1 Parametrisation............................. Pure Interpolation.......................... Abstract
More informationThe Fixed Income Valuation Course. Sanjay K. Nawalkha Gloria M. Soto Natalia A. Beliaeva
Interest Rate Risk Modeling The Fixed Income Valuation Course Sanjay K. Nawalkha Gloria M. Soto Natalia A. Beliaeva Interest t Rate Risk Modeling : The Fixed Income Valuation Course. Sanjay K. Nawalkha,
More informationThe duration derby : a comparison of duration based strategies in asset liability management
Edith Cowan University Research Online ECU Publications Pre. 2011 2001 The duration derby : a comparison of duration based strategies in asset liability management Harry Zheng David E. Allen Lyn C. Thomas
More informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationTerm Par Swap Rate Term Par Swap Rate 2Y 2.70% 15Y 4.80% 5Y 3.60% 20Y 4.80% 10Y 4.60% 25Y 4.75%
Revisiting The Art and Science of Curve Building FINCAD has added curve building features (enhanced linear forward rates and quadratic forward rates) in Version 9 that further enable you to fine tune the
More informationThe Information Content of the Yield Curve
The Information Content of the Yield Curve by HANS-JüRG BüTTLER Swiss National Bank and University of Zurich Switzerland 0 Introduction 1 Basic Relationships 2 The CIR Model 3 Estimation: Pooled Time-series
More informationApplying the Principles of Quantitative Finance to the Construction of Model-Free Volatility Indices
Applying the Principles of Quantitative Finance to the Construction of Model-Free Volatility Indices Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg
More informationAppendix A Financial Calculations
Derivatives Demystified: A Step-by-Step Guide to Forwards, Futures, Swaps and Options, Second Edition By Andrew M. Chisholm 010 John Wiley & Sons, Ltd. Appendix A Financial Calculations TIME VALUE OF MONEY
More informationLecture Quantitative Finance Spring Term 2015
implied Lecture Quantitative Finance Spring Term 2015 : May 7, 2015 1 / 28 implied 1 implied 2 / 28 Motivation and setup implied the goal of this chapter is to treat the implied which requires an algorithm
More informationThe Nelson-Siegel-Svensson Model for U.S. Treasury Securities and Its Interpretation
1 The Nelson-Siegel-Svensson Model for U.S. Treasury Securities and Its Interpretation By Lisa Patrick 1 Introduction Whether you are an investor in equities, bonds, real estate, or other financial securities,
More informationCzech Government Bond yields under FX pressure
Czech Government Bond yields under FX pressure Mojmír Simerský May 2018 Abstract This paper presents some results of the yield curve (YC) estimation method proposed in [5]. We focus on the Czech Government
More informationMULTISTAGE PORTFOLIO OPTIMIZATION AS A STOCHASTIC OPTIMAL CONTROL PROBLEM
K Y B E R N E T I K A M A N U S C R I P T P R E V I E W MULTISTAGE PORTFOLIO OPTIMIZATION AS A STOCHASTIC OPTIMAL CONTROL PROBLEM Martin Lauko Each portfolio optimization problem is a trade off between
More informationEXAMINATION II: Fixed Income Valuation and Analysis. Derivatives Valuation and Analysis. Portfolio Management
EXAMINATION II: Fixed Income Valuation and Analysis Derivatives Valuation and Analysis Portfolio Management Questions Final Examination March 2011 Question 1: Fixed Income Valuation and Analysis (43 points)
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationCredit Risk. June 2014
Credit Risk Dr. Sudheer Chava Professor of Finance Director, Quantitative and Computational Finance Georgia Tech, Ernest Scheller Jr. College of Business June 2014 The views expressed in the following
More informationOne-Factor Models { 1 Key features of one-factor (equilibrium) models: { All bond prices are a function of a single state variable, the short rate. {
Fixed Income Analysis Term-Structure Models in Continuous Time Multi-factor equilibrium models (general theory) The Brennan and Schwartz model Exponential-ane models Jesper Lund April 14, 1998 1 Outline
More informationA nonparametric zero coupon yield curve
A nonparametric zero coupon yield curve Mojmír Simerský January 016 Keywords: yield curve estimation, nonparametric regression, optimal smoothness Abstract The paper adds to a volume of publications on
More informationarxiv: v2 [q-fin.pr] 23 Nov 2017
VALUATION OF EQUITY WARRANTS FOR UNCERTAIN FINANCIAL MARKET FOAD SHOKROLLAHI arxiv:17118356v2 [q-finpr] 23 Nov 217 Department of Mathematics and Statistics, University of Vaasa, PO Box 7, FIN-6511 Vaasa,
More informationMartingale Pricing Theory in Discrete-Time and Discrete-Space Models
IEOR E4707: Foundations of Financial Engineering c 206 by Martin Haugh Martingale Pricing Theory in Discrete-Time and Discrete-Space Models These notes develop the theory of martingale pricing in a discrete-time,
More informationF A S C I C U L I M A T H E M A T I C I
F A S C I C U L I M A T H E M A T I C I Nr 38 27 Piotr P luciennik A MODIFIED CORRADO-MILLER IMPLIED VOLATILITY ESTIMATOR Abstract. The implied volatility, i.e. volatility calculated on the basis of option
More informationTHE NEW EURO AREA YIELD CURVES
THE NEW EURO AREA YIELD CURVES Yield describe the relationship between the residual maturity of fi nancial instruments and their associated interest rates. This article describes the various ways of presenting
More informationIntroduction to Bond Markets
1 Introduction to Bond Markets 1.1 Bonds A bond is a securitized form of loan. The buyer of a bond lends the issuer an initial price P in return for a predetermined sequence of payments. These payments
More informationStatistics and Finance
David Ruppert Statistics and Finance An Introduction Springer Notation... xxi 1 Introduction... 1 1.1 References... 5 2 Probability and Statistical Models... 7 2.1 Introduction... 7 2.2 Axioms of Probability...
More informationMFE8812 Bond Portfolio Management
MFE8812 Bond Portfolio Management William C. H. Leon Nanyang Business School January 16, 2018 1 / 63 William C. H. Leon MFE8812 Bond Portfolio Management 1 Overview Value of Cash Flows Value of a Bond
More informationStatistical Models and Methods for Financial Markets
Tze Leung Lai/ Haipeng Xing Statistical Models and Methods for Financial Markets B 374756 4Q Springer Preface \ vii Part I Basic Statistical Methods and Financial Applications 1 Linear Regression Models
More informationA Non-Parametric Technique of Option Pricing
1 A Non-Parametric Technique of Option Pricing In our quest for a proper option-pricing model, we have so far relied on making assumptions regarding the dynamics of the underlying asset (more or less realistic)
More informationMarket interest-rate models
Market interest-rate models Marco Marchioro www.marchioro.org November 24 th, 2012 Market interest-rate models 1 Lecture Summary No-arbitrage models Detailed example: Hull-White Monte Carlo simulations
More informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
More informationCrashcourse Interest Rate Models
Crashcourse Interest Rate Models Stefan Gerhold August 30, 2006 Interest Rate Models Model the evolution of the yield curve Can be used for forecasting the future yield curve or for pricing interest rate
More informationModelling the Zero Coupon Yield Curve:
Modelling the Zero Coupon Yield Curve: A regression based approach February,2010 12 th Global Conference of Actuaries Srijan Sengupta Section 1: Introduction What is the zero coupon yield curve? Its importance
More informationTutorial 4 - Pigouvian Taxes and Pollution Permits II. Corrections
Johannes Emmerling Natural resources and environmental economics, TSE Tutorial 4 - Pigouvian Taxes and Pollution Permits II Corrections Q 1: Write the environmental agency problem as a constrained minimization
More informationMTH6154 Financial Mathematics I Interest Rates and Present Value Analysis
16 MTH6154 Financial Mathematics I Interest Rates and Present Value Analysis Contents 2 Interest Rates 16 2.1 Definitions.................................... 16 2.1.1 Rate of Return..............................
More informationNo arbitrage conditions in HJM multiple curve term structure models
No arbitrage conditions in HJM multiple curve term structure models Zorana Grbac LPMA, Université Paris Diderot Joint work with W. Runggaldier 7th General AMaMeF and Swissquote Conference Lausanne, 7-10
More informationMortality Rates Estimation Using Whittaker-Henderson Graduation Technique
MATIMYÁS MATEMATIKA Journal of the Mathematical Society of the Philippines ISSN 0115-6926 Vol. 39 Special Issue (2016) pp. 7-16 Mortality Rates Estimation Using Whittaker-Henderson Graduation Technique
More informationA Study on the Risk Regulation of Financial Investment Market Based on Quantitative
80 Journal of Advanced Statistics, Vol. 3, No. 4, December 2018 https://dx.doi.org/10.22606/jas.2018.34004 A Study on the Risk Regulation of Financial Investment Market Based on Quantitative Xinfeng Li
More informationInterest Rate Basis Curve Construction and Bootstrapping Guide
Interest Rate Basis Curve Construction and Bootstrapping Guide Michael Taylor FinPricing The term structure of an interest rate basis curve is defined as the relationship between the basis zero rate and
More informationPricing Implied Volatility
Pricing Implied Volatility Expected future volatility plays a central role in finance theory. Consequently, accurate estimation of this parameter is crucial to meaningful financial decision-making. Researchers
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationFixed-Income Options
Fixed-Income Options Consider a two-year 99 European call on the three-year, 5% Treasury. Assume the Treasury pays annual interest. From p. 852 the three-year Treasury s price minus the $5 interest could
More informationOULU BUSINESS SCHOOL. Ilkka Rahikainen DIRECT METHODOLOGY FOR ESTIMATING THE RISK NEUTRAL PROBABILITY DENSITY FUNCTION
OULU BUSINESS SCHOOL Ilkka Rahikainen DIRECT METHODOLOGY FOR ESTIMATING THE RISK NEUTRAL PROBABILITY DENSITY FUNCTION Master s Thesis Finance March 2014 UNIVERSITY OF OULU Oulu Business School ABSTRACT
More informationMorningstar Fixed-Income Style Box TM
? Morningstar Fixed-Income Style Box TM Morningstar Methodology Effective Apr. 30, 2019 Contents 1 Fixed-Income Style Box 4 Source of Data 5 Appendix A 10 Recent Changes Introduction The Morningstar Style
More informationMFE8812 Bond Portfolio Management
MFE8812 Bond Portfolio Management William C. H. Leon Nanyang Business School January 23, 2018 1 / 118 William C. H. Leon MFE8812 Bond Portfolio Management 1 Overview Shapes of the Term Structure Evolution
More informationA SUMMARY OF OUR APPROACHES TO THE SABR MODEL
Contents 1 The need for a stochastic volatility model 1 2 Building the model 2 3 Calibrating the model 2 4 SABR in the risk process 5 A SUMMARY OF OUR APPROACHES TO THE SABR MODEL Financial Modelling Agency
More informationMTH6154 Financial Mathematics I Interest Rates and Present Value Analysis
16 MTH6154 Financial Mathematics I Interest Rates and Present Value Analysis Contents 2 Interest Rates and Present Value Analysis 16 2.1 Definitions.................................... 16 2.1.1 Rate of
More informationOutline. 1 Introduction. 2 Algorithms. 3 Examples. Algorithm 1 General coordinate minimization framework. 1: Choose x 0 R n and set k 0.
Outline Coordinate Minimization Daniel P. Robinson Department of Applied Mathematics and Statistics Johns Hopkins University November 27, 208 Introduction 2 Algorithms Cyclic order with exact minimization
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Consider
More informationCharacterization of the Optimum
ECO 317 Economics of Uncertainty Fall Term 2009 Notes for lectures 5. Portfolio Allocation with One Riskless, One Risky Asset Characterization of the Optimum Consider a risk-averse, expected-utility-maximizing
More informationChapter 2: BASICS OF FIXED INCOME SECURITIES
Chapter 2: BASICS OF FIXED INCOME SECURITIES 2.1 DISCOUNT FACTORS 2.1.1 Discount Factors across Maturities 2.1.2 Discount Factors over Time 2.1 DISCOUNT FACTORS The discount factor between two dates, t
More informationMarkowitz portfolio theory
Markowitz portfolio theory Farhad Amu, Marcus Millegård February 9, 2009 1 Introduction Optimizing a portfolio is a major area in nance. The objective is to maximize the yield and simultaneously minimize
More informationJaime Frade Dr. Niu Interest rate modeling
Interest rate modeling Abstract In this paper, three models were used to forecast short term interest rates for the 3 month LIBOR. Each of the models, regression time series, GARCH, and Cox, Ingersoll,
More informationINTEREST RATES AND FX MODELS
INTEREST RATES AND FX MODELS 7. Risk Management Andrew Lesniewski Courant Institute of Mathematical Sciences New York University New York March 8, 2012 2 Interest Rates & FX Models Contents 1 Introduction
More informationUtility Indifference Pricing and Dynamic Programming Algorithm
Chapter 8 Utility Indifference ricing and Dynamic rogramming Algorithm In the Black-Scholes framework, we can perfectly replicate an option s payoff. However, it may not be true beyond the Black-Scholes
More informationCalculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the
VaR Pro and Contra Pro: Easy to calculate and to understand. It is a common language of communication within the organizations as well as outside (e.g. regulators, auditors, shareholders). It is not really
More informationForecast Horizons for Production Planning with Stochastic Demand
Forecast Horizons for Production Planning with Stochastic Demand Alfredo Garcia and Robert L. Smith Department of Industrial and Operations Engineering Universityof Michigan, Ann Arbor MI 48109 December
More information1 Unemployment Insurance
1 Unemployment Insurance 1.1 Introduction Unemployment Insurance (UI) is a federal program that is adminstered by the states in which taxes are used to pay for bene ts to workers laid o by rms. UI started
More informationValuation and Optimal Exercise of Dutch Mortgage Loans with Prepayment Restrictions
Bart Kuijpers Peter Schotman Valuation and Optimal Exercise of Dutch Mortgage Loans with Prepayment Restrictions Discussion Paper 03/2006-037 March 23, 2006 Valuation and Optimal Exercise of Dutch Mortgage
More informationarxiv: v1 [math.st] 6 Jun 2014
Strong noise estimation in cubic splines A. Dermoune a, A. El Kaabouchi b arxiv:1406.1629v1 [math.st] 6 Jun 2014 a Laboratoire Paul Painlevé, USTL-UMR-CNRS 8524. UFR de Mathématiques, Bât. M2, 59655 Villeneuve
More informationMarket risk measurement in practice
Lecture notes on risk management, public policy, and the financial system Allan M. Malz Columbia University 2018 Allan M. Malz Last updated: October 23, 2018 2/32 Outline Nonlinearity in market risk Market
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationThe Yield Envelope: Price Ranges for Fixed Income Products
The Yield Envelope: Price Ranges for Fixed Income Products by David Epstein (LINK:www.maths.ox.ac.uk/users/epstein) Mathematical Institute (LINK:www.maths.ox.ac.uk) Oxford Paul Wilmott (LINK:www.oxfordfinancial.co.uk/pw)
More informationA Study on the Relationship between Monetary Policy Variables and Stock Market
International Journal of Business and Management; Vol. 13, No. 1; 2018 ISSN 1833-3850 E-ISSN 1833-8119 Published by Canadian Center of Science and Education A Study on the Relationship between Monetary
More information1 The Solow Growth Model
1 The Solow Growth Model The Solow growth model is constructed around 3 building blocks: 1. The aggregate production function: = ( ()) which it is assumed to satisfy a series of technical conditions: (a)
More informationEstimating the Interest Rate Term Structures of Treasury and Corporate Debt with Bayesian Penalized Splines
Journal of Data Science 3(25), 223-24 Estimating the Interest Rate Term Structures of Treasury and Corporate Debt with Bayesian Penalized Splines Min Li 1 and Yan Yu 2 1 California State University, Sacramento
More informationWhere Has All the Value Gone? Portfolio risk optimization using CVaR
Where Has All the Value Gone? Portfolio risk optimization using CVaR Jonathan Sterbanz April 27, 2005 1 Introduction Corporate securities are widely used as a means to boost the value of asset portfolios;
More informationDiscrete time interest rate models
slides for the course Interest rate theory, University of Ljubljana, 2012-13/I, part II József Gáll University of Debrecen, Faculty of Economics Nov. 2012 Jan. 2013, Ljubljana Introduction to discrete
More informationSimultaneous optimization for wind derivatives based on prediction errors
2008 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 2008 WeA10.4 Simultaneous optimization for wind derivatives based on prediction errors Yuji Yamada Abstract Wind
More informationCash Balance Plans: Valuation and Risk Management Cash Balance Plans: Valuation and Risk Management
w w w. I C A 2 0 1 4. o r g Cash Balance Plans: Valuation and Risk Management Cash Balance Plans: Valuation and Risk Management Mary Hardy, with David Saunders, Mike X Zhu University Mary of Hardy Waterloo
More informationPricing in markets modeled by general processes with independent increments
Pricing in markets modeled by general processes with independent increments Tom Hurd Financial Mathematics at McMaster www.phimac.org Thanks to Tahir Choulli and Shui Feng Financial Mathematics Seminar
More informationMarket Risk Analysis Volume I
Market Risk Analysis Volume I Quantitative Methods in Finance Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume I xiii xvi xvii xix xxiii
More informationBayesian Finance. Christa Cuchiero, Irene Klein, Josef Teichmann. Obergurgl 2017
Bayesian Finance Christa Cuchiero, Irene Klein, Josef Teichmann Obergurgl 2017 C. Cuchiero, I. Klein, and J. Teichmann Bayesian Finance Obergurgl 2017 1 / 23 1 Calibrating a Bayesian model: a first trial
More informationRISKMETRICS. Dr Philip Symes
1 RISKMETRICS Dr Philip Symes 1. Introduction 2 RiskMetrics is JP Morgan's risk management methodology. It was released in 1994 This was to standardise risk analysis in the industry. Scenarios are generated
More informationBOND ANALYTICS. Aditya Vyas IDFC Ltd.
BOND ANALYTICS Aditya Vyas IDFC Ltd. Bond Valuation-Basics The basic components of valuing any asset are: An estimate of the future cash flow stream from owning the asset The required rate of return for
More informationUniversité de Montréal. Rapport de recherche. Empirical Analysis of Jumps Contribution to Volatility Forecasting Using High Frequency Data
Université de Montréal Rapport de recherche Empirical Analysis of Jumps Contribution to Volatility Forecasting Using High Frequency Data Rédigé par : Imhof, Adolfo Dirigé par : Kalnina, Ilze Département
More informationEffectiveness of CPPI Strategies under Discrete Time Trading
Effectiveness of CPPI Strategies under Discrete Time Trading S. Balder, M. Brandl 1, Antje Mahayni 2 1 Department of Banking and Finance, University of Bonn 2 Department of Accounting and Finance, Mercator
More informationInt. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p approach
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p.5901 What drives short rate dynamics? approach A functional gradient descent Audrino, Francesco University
More informationStructural Models of Credit Risk and Some Applications
Structural Models of Credit Risk and Some Applications Albert Cohen Actuarial Science Program Department of Mathematics Department of Statistics and Probability albert@math.msu.edu August 29, 2018 Outline
More informationP2.T5. Market Risk Measurement & Management. Bruce Tuckman, Fixed Income Securities, 3rd Edition
P2.T5. Market Risk Measurement & Management Bruce Tuckman, Fixed Income Securities, 3rd Edition Bionic Turtle FRM Study Notes By David Harper, CFA FRM CIPM www.bionicturtle.com Tuckman, Chapter 6: Empirical
More informationBloomberg. Portfolio Value-at-Risk. Sridhar Gollamudi & Bryan Weber. September 22, Version 1.0
Portfolio Value-at-Risk Sridhar Gollamudi & Bryan Weber September 22, 2011 Version 1.0 Table of Contents 1 Portfolio Value-at-Risk 2 2 Fundamental Factor Models 3 3 Valuation methodology 5 3.1 Linear factor
More information