ORIGINAL SCIENTIFIC PAPER RECEIVED: DECEMBER 2015 REVISED: MARCH 2016 ACCEPTED: MARCH 2016 DOI: 10.1515/ngoe-2016-0011 UDK: 303.725.35 JEL: E43, G12, C52 Citation: Lorencic, E. (2016). Testing the Performance of Cubic Splines and Nelson-Siegel Model for Estimating the Zero-coupon Yield Curve. Nase gospodarstvo/Our Economy, 62(2), 42-50. DOI: 10.1515/ngoe-2016-0011 Testing the Performance of Cubic Splines and Nelson-Siegel Model for Estimating the Zero-coupon Yield Curve Eva Lorenčič Ranca 35, 2211 Pesnica pri Mariboru, Slovenia eva.loren@gmail.com Abstract Understanding the relationship between interest rates and term to maturity of securities is a prerequisite for developing financial theory and evaluating whether it holds up in the real world; therefore, such an understanding lies at the heart of monetary and financial economics. Accurately fitting the term structure of interest rates is the backbone of a smoothly functioning financial market, which is why the testing of various models for estimating and predicting the term structure of interest rates is an important topic in finance that has received considerable attention for many decades. In this paper, we empirically contrast the performance of cubic splines and the Nelson-Siegel model by estimating the zero-coupon yields of Austrian government bonds. The main conclusion that can be drawn from the results of the calculations is that the Nelson-Siegel model outperforms cubic splines at the short end of the yield curve (up to 2 years), whereas for medium-term maturities (2 to 10 years) the fitting performance of both models is comparable. Keywords: Cubic splines; Nelson-Siegel; yield curve; zero-coupon bonds; term structure of interest rates 1 Introduction NG NAŠE GOSPODARSTVO OUR ECONOMY Vol. . 62 No. 2 2016 pp. . 42-50 The purpose of this paper is to construct a zero-coupon yield curve based on the data for coupon yields, price, and maturity of Austrian government bonds. The models used are the Nelson-Siegel and cubic splines. Forecasting the term structure of interest rates is a prerequisite for managing investment portfolios, pricing financial assets and their derivatives, calculating risk measures, valuing capital goods, managing pension funds, formulating economic policy, making household finance decisions, and managing fixed-income wealth. Prices of fixed-income securities such as swaps, bonds, and mortgage-backed securities depend on the yield curve. The yields on default-free government bonds that have different maturities, when examined together, reveal information about forward rates, which can predict real economic activity and are, hence, of interest for policymakers, market participants, and economists. For instance, forward rates are often inputs to pricing models and may indicate market expectations of the movement of inflation rates and currency appreciation/depreciation rates in the future. Understanding the relationship between interest rates and the term to maturity of securities is a prerequisite for developing and testing financial theory and, therefore, lies at the heart 24 Eva Lorencic: Testing the Performance of Cubic Splines and Nelson-Siegel Model for Estimating the Zero-coupon Yield Curve of monetary and financial economics. Accurately fitting the term structure of interest rates is the backbone of a smoothly functioning financial market, which is why refining yield curve modeling and forecasting methods is an important topic in finance that has received considerable attention for many decades (De Rezende & Ferreira, 2013; Diebold, Li, & Yue, 2008; Exterkate, van Dijk, Heij, & Groenen, 2013; Ioannides, 2003; Jordan & Mansi, 2003; Linton, Mammen, Nielsen, & Tanggaard, 2001; Rugengamanzi, 2013). The yield curve is a graphical representation of the term structure of interest rates (i.e., a one-to-one relationship between yields and corresponding maturities of default-free zero-coupon securities issued by sovereign lenders). The term structure of interest rates contains information about the yields of zero-coupon bonds1 of various maturities at a certain date. Constructing the term structure of interest rates is not a straightforward task due to the scarcity of zero-coupon bonds on the market, which represent the essential part of the term structure of interest rates. The majority of bonds traded in the market bear coupons. The yields to maturity on coupon-bearing bonds, whose maturities or coupons differ, are not immediately comparable. As a result, a uniform way of measuring the term structure of interest rates is needed: The spot interest rates2 (i.e., the yields earned on bonds that pay no coupon) must be estimated from coupon bond prices of bonds with different maturities by using interpolation methods, such as polynomial splines (e.g., cubic splines) and parsimonious functions (e.g., Nelson-Siegel). This is how the yield curve of zero-coupon bonds is constructed (Christensen, Diebold, & Rudebusch, 2011; Christofi, 1998; Gauthier & Simonato, 2012; Luo, Han, & Zhang, 2012; Teichmann & Wüthrich, 2013; Yu & Zivot, 2011). The most widely used models for estimating the zero-coupon yield curve are Nelson-Siegel and cubic splines. For instance, the central banks of Belgium, Finland, France, Germany, Italy, Norway, Spain, and Switzerland use the Nel-son-Sigel model or some type of its enhanced extension to fit and forecast yield curves (BIS, 2005). The European Central Bank uses the Sonderlind-Svensson model, an extension of the Nelson-Siegel model, to estimate yield curves in the Eurozone (Coroneo, Nyholm & Vidova-Koleva, 2011). The remainder of the paper is organized as follows. Section 2 gives an overview of the existing literature and relevant research studies. Section 3 presents the data; Section 4 lays out the methodology and the results. Section 5 concludes the paper. 1 Another name for "zero-coupon bond" is discount bond. 2 The terms "yield to maturity on a zero-coupon bond," "zero-cou- pon interest rate," "spot interest rate" and "zero-coupon yield" are synonyms; they all describe the same aspect of reality. 2 Literature Review The yield curve estimation methods originated in McCull-och's (1975) cubic splines and in Nelson and Siegel's (1987) parsimonious function. Most of the research studies on the term structure of interest rates build on these two methods and propose improvements and extensions. McCulloch (1975) modeled the discount curve with a spline. The fitted discount curve gives a poor fit of the yield curves, most notably at longer maturities. Nelson-Siegel's parsimonious function allows for various shapes of the yield curve. The forward rates are a solution to a second-order differential equation. The forward rate curve under McCulloch's method is not smooth, whereas the forward rate curve under Nelson-Siegel's method is smooth, but still unable to accurately price instruments at the longest end of the yield curve (Rugengamanzi, 2013). The Nelson-Siegel type of models are relatively efficient in capturing the general shapes of the yield curve, which is why they are extensively used by central banks and market practitioners. Nevertheless, the Nelson-Siegel type of models are still inferior to the dynamic term structure models, like the quadratic or affine term structure models or the forward rate-based arbitrage-free model introduced by Heath, Jarrow, and Morton (1992). Jordan and Mansi (2003) used five distinct yield curve-smoothing methods to derive spot rates from on-the-run treasuries. All methods used the bootstrapping3 technique either in discrete time or in continuous time. Yield curve-smoothing methods based on continuous-time bootstrapping deliver a superior approximation of the term structure of interest rates to those methods that employ the discrete-time bootstrapping technique. Of the five yield curve-smoothing methods, the linear spline interpolation produces the worst results; the second-worst method is cubic splines; and the two best methods are Mansi-Phillips and Nelson-Siegel. Diebold and Li (2006) developed a dynamic version of the Nelson-Siegel model. They showed that the three factors present in the Nelson-Siegel model can be interpreted as the level, slope, and curvature of the yield curve. They also corroborated that the dynamic model improves the forecasting accuracy. Diebold, Rudebusch, and Aruoba (2006) combined the Diebold-Li model with macroeconomic variables to analyze the relationship between the yield curve and the economy. Diebold et al. (2008) extended the dynamic Diebold-Li model to a global context. A large set of country yield curves was modeled in a setting that allows for both country-specific and global factors. The researchers found that the global yield level and slope factors exist and can 3 "Bootstrapping" refers to calculating spot rates from bond yields in an iterative manner (Jordan et al., 2003). 49 NAŠE GOSPODARSTVO / OUR ECONOMY Vol. 62 No. 2 / June 2016 explain a substantial proportion of variation in country bond yields. Christensen, Diebold, and Rudebusch (2009) proposed a generalized arbitrage-free Nelson-Siegel model using five factors. Christensen et al. (2011) added the arbitrage-free restriction to the Diebold-Li model, thereby creating an affine arbitrage-free Nelson-Siegel model (AFNS). The results show an improvement in the model's predictive performance. Yu and Zivot (2011) empirically tested Diebold and Li's dynamic Nelson-Siegel three factor model and found that the dynamic Diebold-Li factor AR(1) model is superior to other models on the out-of-sample forecast accuracy. Luo et al. (2012) compared the forecasting ability of the Diebold-Li, dynamic Svennson, and dynamic Bjork and Christensen models for the term structure of Chinese treasury yields. The results showed that all three models fit the data very well and that more flexible models produced superior in-sample fitting performance. Gauthier and Simonato (2012) developed linearized algorithms for estimating spot interest rate term structures. These algorithms converge much faster while retaining the important characteristics of the original approaches. These algorithms are superior in that they enable the inclusion of prior information about some of the parameters, thereby enhancing the precision of the estimated spot rate curves. Using the Brazilian yield curve data, De Rezende and Ferreira (2013) compared the in-sample adjustment and the out-of-sample forecasting performance of four different Nelson-Siegel type models: Nelson-Siegel, Bliss, Svensson, and a new five-factor model that is an extension of the Svensson model and could improve the fitting flexibility. The introduction of the fifth factor into the model produced the best in-sample fitting, but poor out-of-sample forecasting. Exterkate et al. (2013) investigated whether the inclusion of additional macroeconomic information into the Nelson-Siegel model results in improved yield curve forecasts. In general, the forecasts could not be improved in stable times (e.g., 1994-1998): When the yields are not volatile, the dynamic Nelson-Siegel model produces good yield curve forecasts. The inclusion of additional macroeconomic variables can substantially improve the yield curve forecast accuracy when the yields are volatile (e.g., 2008-2009). 3 Data We empirically tested the fitting performance of the Nelson-Siegel model and cubic splines on the data for Austrian government bonds. The data were retrieved from Bloomb-erg4 on October 8, 2013(see Table 1). All prices are in euros. Table 1. Data on Austrian Government Bonds (retrieved from Bloomberg on October 8, 2013). AT -Benchmark Maturity Coupon (%) Bid (EUR) Ask (EUR) 1YBUND 20.10.2014 3,40 103,32 103,35 2YBUND 15.07.2015 3,50 105,72 105,77 3YBUND 15.09.2016 4,00 110,14 110,18 4YBUND 15.09.2017 4,30 113,39 113,46 5YBUND 15.01.2018 4,65 115,53 115,58 6YBUND 15.03.2019 4,35 116,25 116,29 7YBUND 15.07.2020 3,90 114,85 114,91 8YBUND 15.09.2021 3,50 112,44 112,50 9YBUND 20.04.2022 3,65 113,35 113,42 10YBUND 20.10.2023 1,75 95,98 96,02 13YBUND 15.03.2026 4,85 125,63 125,81 4 We do not have access to more recent data, as Bloomberg terminal subscription costs around $20,000 per user per annum, which is why the license is mostly bought by and available in major institutions only. Source: Bloomberg. The raw data given in Table 1 was used to calculate the accrued interest and the dirty price (i.e., the market value, the present value). The dirty price is a sum of the clean price (i.e., the average of the "bid quote" and the "ask quote") and the accrued interest. The accrued interest is the interest that a bond holder would have obtained in theory between the last coupon date of each bond and the current date (which we assumed to be October 8, 2013—the day the data were extracted). The Nelson-Siegel model and cubic splines were applied to the data in Table 1 to estimate the zero-coupon bond yield curve. 4 Methodology and Results The two main categories of methods for estimating a yield curve are the spline methods and the parsimonious methods. The spline methods use a piecewise polynomial function (usually a cubic one) to approximate the yield curve. Cubic splines were first introduced by McCulloch (1975) and subsequently improved by Fisher, Nychka, and Zervos (1995), Waggoner (1997), and Anderson and Sleath (1999). The parsimonious methods (such as Nelson-Siegel model [1987]; Nelson-Siegel-Svensson model [Sevensson, 1995], and models described by Wiseman [1994] and Bjork and Christensen [1997]) approximate the yield curve by estimating the parameters in a single parametric function. With spline methods, it is possible to capture almost any shape of the yield curve, whereas parsimonious methods can capture only yield curves obeying certain financial constraints. Both groups of methods have specific advantages and disadvantages, but none of them seems to be able to consistently ~26 Eva Lorencic: Testing the Performance of Cubic Splines and Nelson-Siegel Model for Estimating the Zero-coupon Yield Curve outperform the other (Manousopoulos & Michalopoulos, 2009). In this paper, we employed the most representative type of each method (cubic splines and Nelson-Siegel) with the goal of estimating the zero-coupon yield curve. 4.1 Cubic splines model The cubic splines method, developed by McCulloch (1971) and McCulloch (1975), divides the zero-coupon yield curve into distinct intervals. In each of these intervals, a cubic spline acts as vertebra in the vertebrate spinal column. If the yield curve is divided into k-1 knots, then we need k parameters to describe the entire zero-coupon yield curve. The optimal parameters are obtained by constructing the matrix A, defined in the continuation of the paper. The optimal parameters and the optimal zero-coupon yield curve minimize the discrepancy (the error) between the model price and the market price of government bonds (McCulloch, 1975; Jankowitsch & Pichler, 2003; Rugengamanzi, 2013). To fit the observed market data for government bond yield curve, McCulloch (1971) used the discount function presented in (1): D(t) = 1 + Z/=i/}(0 1 aj, (1) a. When j = 1: d = d= 0 i i-i b. When j = k: f (t) = t (6) (7) In f. (t), t stands for time and j is a knot number (j = 1, ..., k-1). Because the number of knots is equal to k-1, we need k parameters. The knots are denoted as d.. We used four parameters (k = 4) and three knots (k-1 = 3). We set the first knot equal to zero (dl = 0), the second knot equal to five (d2 = 5), and the third knot equal to 13 (d3 = 13), which corresponds to the maturity of the last bond (McCulloch, 1975; Jankowitsch & Pichler, 2003; Rugengamanzi, 2013). The zero-coupon yields are calculated from the discount function defined by (1) as follows: (8) In order to calculate the parameters a. (which are needed in (1) and (8)), we first need to construct matrix A: A = (X-XTy1-XT-Y, (9) where a. are the parameters that need to be estimated. For j < k, f. (t) is a cubic polynomial defined as follows: I. When t < d i-i f (t) = 0 (2) II. When d., < t < d: i-i i m = (t - dj-i) 6 ■ (dj - dj^y (3) III. When d < t < d : i i+i /#CO = — + — + - 2 2 6-(dj+x~ (4) where c = d- d , and e = t - d.. i i-i i IV. When d < t: i+i fj(t) = (dj+1 - dj Two additional conditions have to be met: 2 ■ d]+1 - dj - d}_! ^ t - dj+1 where X is a matrix and Y is a matrix. X1 and Y are the transposed matrices of matrices X and Y. The matrix A has j rows and one column. The element a of matrix A (where j = 1, ..., k) is an optimal parameter a., needed to calculate the discount function, as defined by (1). The matrix X is given as follows: xij = YjZith)-fi{wi(h)) (10) h=1 for every i=1,...,n and for every}=1,...,k, where R. = the number of future cash flows of bond i, i ' Z. (h) = a future cash flow h of bond i, w. (h) = the time to maturity of a next cash flow h of bond i in years, and f (t) = a function of time t, defined by equations (j2) to (5). . (5) The matrix Y is given as follows: Ri yi = Pi - ^ Zi(h)for everyi = h.-.n (11) /1=1 49 NAŠE GOSPODARSTVO / OUR ECONOMY Vol. 62 No. 2 / June 2016 where P. = the dirty price (the present value, the market value) of bond i. The estimated zero-coupon yield curve is displayed in Figure 1. The dirty price (the present value, the market value) Pi of bond is defined as follows: P=clean price from Bloomberg+accrued interest. (12) The accrued interest is defined as follows: accrued _ the next cash flow interest (coupon amount) time that has passed between the last coupon payment and today time between two consecutive coupon payments (13) Table 2. Estimated Discount Factors (D(t)) and Zero-coupon Yields for t = 1, 2.....10. t D(t) r(t) r(t) in % 1 1,00056 -0,00056 -0,05606 2 0,99443 0,00279 0,27943 3 0,98238 0,00594 0,59424 4 0,96517 0,00890 0,89031 5 0,94355 0,01169 1,16896 6 0,91830 0,01431 1,43070 7 0,89024 0,01675 1,67481 8 0,86020 0,01900 1,90021 9 0,82901 0,02105 2,10539 10 0,79751 0,02288 2,28837 The final results obtained by empirically testing the cubic splines model (equations (1) to (13)) on data given in Table 1 (Austrian government bonds) are summarized in Table 2. Figure 1. The zero-coupon yield curve, estimated with the cubic splines model. 4.2 Nelson-Siegel model Nelson and Siegel (1987) developed a parsimonious function to model forward rates. The zero-coupon yield (the spot rate) can be derived as follows: r(t) =ß1+ß2 1 — e t_ t_ Tl ' + ß3 ' 1 — e t_ = ßi + (ß2+ß3) l-e Ti t ti ~ß3 t_ Tl ■ e Ti, (14) where t = time to maturity of a bond (in years), = parameter beta 1 (the level factor), P2 = parameter beta 2 (the slope factor), P3 = parameter beta 3 (the curvature factor), t = parameter tau 1 (the rate of exponential decay), and e = exponential function. The parameters P1, P2, P3, and tcan be calculated with the Excel add-in "Solver" by minimizing the sum of squared residuals between the dirty price (market value, present value) of the bonds and the model price of the bonds. (The dirty price is a sum of the clean price, retrieved from Bloomberg, and accrued interest.) The price of zero-coupon securities for time t is calculated as follows: -r(t)t p(t) = e «o . (15) The market value (MV) of a bond i according to the Nelson-Siegel model is calculated as follows: MV Nelson Siegel (bondt) = ZfcLiZt(h) ■ P(wj(7i)), (16) 2,5 2 1,5 g 1 0,5 0 -0,5 12 345 6 7 8 9 10 t where R. = the number of future cash flows of bond i, i ' Z. (h) = a future cash flow h of bond i, w. (h) = the time to maturity of a future cash flow h of bond i in years, and P(t) = the present value of a zero-coupon security with nominal value 1 and maturity in time t. The error of the Nelson-Siegel model for bond i is defined as follows: The error of the model for bond i = = (dirty price of bond. - MV Nelson Siegel of bond. )2 (17) ~28 X Eva Lorencic: Testing the Performance of Cubic Splines and Nelson-Siegel Model for Estimating the Zero-coupon Yield Curve The final results obtained by empirically testing the cubic splines model (equations (14) to (17)) on data given in Table 1 (Austrian government bonds) are summarized in Table 3. The estimated zero-coupon yield curve is displayed in Figure 2. 4.3 Comparison of the two models Our coupon bond price estimation results (Tables 4 and 5) are comparable with the in-sample coupon bond price estimation results of Jordan and Mansi (2003; see Table 2). Jordan and Mansi's (2003) results show that, at the short end of the yield curve (0 to 5 years), the Nelson-Siegel model outperforms cubic splines; the same is true for the intermediate range (5 to 10 years). Our findings are similar to Jordan and Mansi's (2003) findings in that the Nelson-Siegel model performed better than cubic splines at the short end of the yield curve (up to 5 years) and in the intermediate range (5-10 years); however, the supremacy of the Nelson-Siegel model over cubic splines was more pronounced for short-term maturities than for medium-term maturities. The cubic splines and Nelson-Siegel estimates of zero-coupon yields are summarized and compared in Table 6. The Nelson-Siegel model outperformed cubic splines at the short end of the yield curve (up to 2 years), whereas for medium-term maturities (2 to 10 years), the fitting performance of both models was comparable. t Table 4. Comparison of Cubic Splines and Nelson-Siegel Estimates of Bond Prices. Maturity of the bond Actual bond price Cubic splines bond price estimation Nelson-Siegel bond price estimation Cubic splines squared price error* Nelson-Siegel squared price error** 1 year 106.62 106.85 106.68 0.05061 0.00308 2 years 106.56 106.63 106.55 0.00532 0.00006 3 years 110.41 110.24 110.31 0.02882 0.01081 4 years 113.70 113.60 113.69 0.00949 0.00001 5 years 118.94 118.89 118.97 0.00301 0.00082 6 years 118.74 118.76 118.77 0.00058 0.00134 7 years 115.79 115.90 115.84 0.01271 0.00242 8 years 112.69 112.67 112.60 0.00038 0.00797 9 years 115.10 115.26 115.21 0.02767 0.01350 10 years 97.69 97.47 97.52 0.04878 0.02825 13 years 128.47 128.52 128.54 0.00213 0.00422 Sum: 0.18951 Sum: 0.07246 Notes. *Cubic splines squared price error is equal to the squared difference between the actual bond price and the cubic splines bond price estimation. **Nelson-Siegel squared price error is equal to the squared difference between the actual bond price and the Nelson-Siegel bond price estimation. Table 3. Estimated Zero-coupon Yields for t = 1, 2.....10. t r(t) in % 1 0.11 2 0.30 3 0.57 4 0.86 5 1.15 6 1.42 7 1.67 8 1.89 9 2.09 10 2.26 Figure 2. The zero-coupon yield curve, estimated with the Nelson-Siegel model. NAŠE GOSPODARSTVO / OUR ECONOMY Vol. 62 No. 2 / June 2016 Table 5. Comparison of Cubic Splines and Nelson-Siegel Squared Price Errors for Short-term and Medium-term Maturities. Sum of squared Sum of squared Maturity price errors, cubic price errors, splines Nelson-Siegel Up to 5 years 0.0972" 0.01477" > 5-10 years 0.09012"" 0.05348"" Notes. *Sum of squared price errors for maturities from 1 to 5 years. **Sum of squared price errors for maturities from 6 to 10 years. Table 6. Comparison of Cubic Splines and Nelson-Siegel Estimates of Zero-coupon Yields. Absolute Cubic splines Nelson-Siegel difference estimates of estimates of between the zero-coupon zero-coupon two model yields in % yield in % estimates in basis points l -0.06 0.11 17 2 0.28 0.30 2 3 0.59 0.57 2 4 0.89 0.86 3 5 1.17 1.15 2 6 1.43 1.42 1 7 1.67 1.67 0 s 1.90 1.89 1 9 2.11 2.09 2 lO 2.29 2.26 3 Sum of absolute differences in basis points: 33 The term structure of interest rates, as estimated by the two models, is displayed in Figure 3. Figure 3. The term structure of interest rates, as estimated with cubic splines (grey line) and Nelson-Siegel (black line). 2,5 2 1,5 g l r 0,5 O -0,5 l234 56 t 7 s 9 lO conclusion drawn from the results of the calculations was that the Nelson-Siegel model outperformed cubic splines at the short end of the yield curve (up to 2 years), whereas for medium-term maturities (2 to 10 years) the fitting performance of both models was comparable. In estimating the term structure of interest rates, we employed the simplest versions of the two models without any further extensions. In reality, Nelson-Siegel is one of the most widely used models for deriving the zero-coupon yield curve; however, the central banks and more sophisticated commercial banks use its enhanced versions. Our study is limited in that we test the fitting performance of the models on the data of government bonds from only one country (Austria). In order to further substantiate our findings, the fitting performance of the models could be tested on a wider set of data, such as government bond data from various countries. 5 Conclusion In this paper, we empirically contrasted the performance of cubic splines and the Nelson-Siegel model by estimating the zero-coupon yields of Austrian government bonds. The main As shown in the theoretical part of this paper, the researchers recently derived the dynamic and no-arbitrage improved versions of the Nelson-Siegel model, which are characterized by a greater estimate precision and forecast accuracy. Further research studies are warranted to evaluate and analyze the performance of the more recent versions of the model. References 1. Anderson, N., & Sleath, J. (1999). New estimates of the UK real and nominal yield curves. Bank of England Quarterly Bulletin, November, 384-392. 2. BIS. (2005). Zero-coupon yield curves: Technical documentation. Basle: Bank for International Settlements. 3. Bjork, T., & Christensen, B. (1997). Interest rate dynamics and consistent forward rate curves. University of Aarhus Working Paper, 1-38. 4. Christensen, J. H., Diebold, F. X., & Rudebusch, G. (2009). An arbitrage-free generalized Nelson-Siegel term structure model. The Econometrics Journal, 12(3), 33-64. http://dx.doi.org/10.1111/j.1368-423X.2008.00267.x 5. Christensen, J. H. E., Diebold, F. X., & Rudebusch, G. D. (2011). The affine arbitrage-free class of Nelson-Siegel term structure models. Journal of Econometrics, 164(1), 4-20. http://dx.doi.org/10.1016/jjeconom.2011.02.011 ~30~ Eva Lorencic: Testing the Performance of Cubic Splines and Nelson-Siegel Model for Estimating the Zero-coupon Yield Curve 6. Christofi,A. C. (1998). Estimation of the Nelson-Siegel parsimonious modeling of yield curves using an exponential GARCH process. Managerial Finance, 24(9/10), 1-19. http://dx.doi.org/10.1108/03074359810765813 7. Coroneo, L., Nyholm, K., & Vidova-Koleva, R. (2011). How arbitrage-free is the Nelson-Siegel model? Journal of Empirical Finance, 18(3), 393-407. http://dx.doi.org/10.1016/jjempfin.2011.03.002 8. De Rezende, R. B., & Ferreira, M. S. (2013). Modeling and forecasting the yield curve by an extended Nelson-Siegel class of models: A quantile autoregression approach. Journal of Forecasting, 32(2), 111-123. http://dx.doi.org/10.1002/for.1256 9. Diebold, F. X., & Li, C. (2006). Forecasting the term structure of government bond yields. Journal of Econometrics, 130(2), 337-364. http://dx.doi.org/10.1016/j.jeconom.2005.03.005 10. Diebold, F. X., Li, C., & Yue, V. Z. (2008). Global yield curve dynamics and interactions: A dynamic Nelson-Siegel approach. Journal of Econometrics, 146(2), 351-363. http://dx.doi.org/10.1016/jjeconom.2008.08.017 11. Diebold, F. X., Rudebusch, G. D., & Aruoba, B. (2006). The macroeconomy and the yield curve: a dynamic latent factor approach. Journal of Econometrics, 131(1-2), 309-338. http://dx.doi.org/10.1016/jjeconom.2005.01.011 12. Exterkate, P., van Dijk, D., Heij, C., & Groenen, P. J. F. (2013). Forecasting the yield curve in a data-rich environment using the factor-augmented Nelson-Siegel model. Journal of Forecasting, 32(3), 193-214. http://dx.doi.org/10.1002/for.1258 13. Fisher, M., Nychka, D., & Zervos, D. (1995). Fitting the term structure of interest rates with smoothing splines. Working Paper 95-1, Finance and Economics Discussion Series, Federal Reserve Board. 14. Gauthier, G., & Simonato, J-G. (2012). Linearized Nelson-Siegel and Svensson models for the estimation of spot interest rates. European Journal of Operational Research, 219(2), 442-451. http://dx.doi.org/10.10167j.ejor.2012.01.004 15. Heath, D., Jarrow, R., & Morton, A. (1992). Bond pricing and the term structure of interest rates: a new methodology for contingent claims valuation. Econometrica, 60(1), 77-105. http://dx.doi.org/10.2307/2951677 16. loannides, M. (2003). A comparison of yield curve estimation techniques using UK data. Journal of Banking & Finance, 27(1), 1-26. http://dx.doi.org/10.1016/S0378-4266(01)00217-5 17. Jankowitsch, R., & Pichler, S. (2003). Estimating zero-coupon yield curves for EMU government bonds. Work based on a research project funded by the Austrian National Bank (OeNB-Jubilaumsfondprojekt Nr. 8808). Vienna University of Economics and Business Administration, December 2003. 18. Jordan, J. V., & Mansi, S. A. (2003). Term structure estimation from on-the-run Treasuries. Journal of Banking & Finance, 27(8), 1487-1509. http://dx.doi.org/10.1016/S0378-4266(02)00273-X 19. Linton, O., Mammen, E., Nielsen, J. P., & Tanggaard, C. (2001). Yield curve estimation by kernel smoothing methods. Journal of Econometrics, 105(1), 185-223. http://dx.doi.org/10.1016/S0304-4076(01)00075-6 20. Luo, X., Han, H., & Zhang, J. E. (2012). Forecasting the term structure of Chinese Treasury yields. Pacific-Basin Finance Journal, 20(5), 639-659. http://dx.doi.org/10.1016/j~.pacfin.2012.02.002 21. Manousopoulos, P., & Michalopoulos, M. (2009). Comparison of non-linear optimization algorithms for yield curve estimation. European Journal of Operational Research, 192(2), 594-602. http://dx.doi.org/10.1016Zj.ejor.2007.09.017 22. McCulloch, J. H. (1971). Measuring the term structure of interest rates. The Journal of Business, 44, 19-31. 23. McCulloch, J. (1975). The tax-adjusted yield curve. Journal of Finance, 30, 811-830. http://dx.doi.org/10.1111/jl540-6261.1975. tb01852.x 24. Nelson, C., & Siegel, A. (1987). Parsimonious modeling of the yield curve. Journal of Business, 60(4), 473-489. 25. Rugengamanzi, M. N. (2013). Term structure estimation based on a generalized optimization framework (Doctoral dissertation, No. 1539). Sweden, Linkoping University, Department of Mathematics. 26. Svensson, L. (1995). Estimating forward interest rates with the extended Nelson and Siegel method. Sveriges Riksbank Quarterly Review, 3. 27. Teichmann, J., & Wüthrich, M. V. (2013). Consistent long-term yield curve prediction. (Under review). Retrieved from https://people. math.ethz.ch/~wueth/Papers/2013_Yield_Curve_Prediction.pdf 28. Waggoner, D. F. (1997). Spline methods for extracting interest rate curves from coupon bond prices. Federal Reserve Bank of Atlanta, Working Paper 97-10. 29. Wiseman, J. (1994). The exponential yield curve model. JPMorgan European Fixed Income Research. 30. Yu, W.-C., & Zivot, E. (2011). Forecasting the term structures of Treasury and corporate yields using dynamic Nelson-Siegel models. International Journal of Forecasting, 27(2), 579-591. http://dx.doi.org/10.1016/j.ijforecast.2010.04.002 Author Eva Lorencic holds a Master's Degree in Economics, awarded by the University of Maribor, Faculty of Economics and Business, and is currently a trainee at the European Central Bank in Frankfurt. She has authored and co-authored several original scientific papers, review papers, and professional papers. She interned with Erste Group Bank AG (in Vienna), Sberbank Europe AG (in Vienna), and Nova KBM (in Maribor). Her professional interests include corporate finance, banking regulation and supervision, monetary and fiscal policy, macroeconomics and international economics. 49 NAŠE GOSPODARSTVO / OUR ECONOMY Vol. 62 No. 2 / June 2016 Testiranje učinkovitosti modela kubičnih zlepkov in Nelson-Sieglovega modela pri ocenjevanju krivulje donosa brezkuponskih vrednostnih papirjev Izvleček Razumevanje razmerja med obrestnimi merami in časom do dospelosti vrednostnih papirjev je osnovni pogoj za razvoj in ovrednotenje pravilnosti finančne teorije. Ta tematika je zato v osrčju monetarne in finančne ekonomije. Natančno prilagajanje terminske strukture obrestnih mer je hrbtenica tekoče delujočega finančnega trga. To je razlog, da je testiranje različnih modelov, ki ocenjujejo in napovedujejo terminsko strukturo obrestnih mer, na področju financ pomembna vsebina, ki je že nekaj desetletij deležna precejšnje pozornosti. V tem članku empirično primerjamo učinkovitost modela kubičnih zlepkov in Nelson-Sieglovega modela, tako da ocenimo donosnost brezkuponskih avstrijskih državnih vrednostnih papirjev. Ključni sklep, ki ga lahko izpostavimo na podlagi dobljenih rezultatov, je, da Nelson-Sieglov model bolje aproksimira brezkuponsko krivuljo donosa na kratkem koncu (do dveh let), na srednjem delu krivulje donosa (od dveh do desetih let) pa med rezultati obeh modelov ni bistvenih razlik. Ključne besede: kubični zlepki, Nelson-Sieglov model, krivulja donosa, brezkuponske obveznice, terminska struktura obrestnih mer ~32