For many years researchers agreed that the first three principal components (PC) of bond yields contain all the relevant information in bond prices. As a result, a large class of affine term structure models assumes that only those three factors determine the evolution of bond prices. Recent studies, however, suggest that other variables contain additional information that is useful for predicting bond returns. This evidence comes from predictive regressions of bond returns on variables other than (and possibly including) the first three PC.1
In an influential paper, Cochrane and Piazzesi (2005) (CP), argue that a linear combination of forward rates, labeled the “return forecasting factor” (RFF), predicts 1-year excess bond returns of each maturity with an R2 of about 0.35.2 One implication of this observation is that a single factor may capture all the economically relevant variations in expected bond returns. As a result, Cochrane and Piazzesi (2008) propose to augment the three-factor affine term structure model by adding a fourth (RFF) factor as the main driver of risk premia.
The RFF is the linear combination of forwards rates that maximizes the predictive ability in 1-year excess returns regressions. Yet, if the RFF is the only driver of bond risk premia, then it should also summarize all the relevant variations in excess returns considering holding periods different from 1 year. We find that it does not. The 1-year RFF outperforms the first three PC only for 1-year excess return regressions. In regressions of excess returns for holding periods longer than a year, the first three PC provide a better fit of excess returns. And for sufficiently long holding periods, the PC have more than three times the explanatory power of the RFF.
We also reject the null hypothesis that there is a single factor structure in expected excess returns over different holding horizons. Although CP also reject the single factor model, they claim that the RFF captures all the economically relevant variations in 1-year excess returns. When considering different holding horizons, however, we find that it does not.
Finally, if the RFF is a relevant state variable in a term structure model, it should also be a significant explanatory variable in unrestricted regressions of bond yields on the RFF and the first three PC. Using the robust t-test developed by Ibragimov and Müller (2010), we reject the significance of the RFF on such regressions, while the first three PC are highly significant.
Taken together, these findings contradict the single factor structure in bond returns and imply that including the RFF as the main driver of bond risk premia in a term structure model is not supported by the data.
2 The return forecasting factor and excess bond returns
We first define notation. Let
denote the annualized h-year log holding return from buying an n-period zero coupon bond at time t and selling it at time t + h as an (n − h)-year zero coupon bond. The excess return from holding an n-year zero coupon bond for h years is defined as
Finally, we denote the average h-year excess return across maturities by
We use monthly data on US bond prices for the period Jan-1964 through December-2003. We focus on the same sample period used by Cochrane and Piazzesi (2005) to make sure that our results do not change because the sample is different. We gather data from two sources. First, we use the Fama-Bliss data (Fama and Bliss, 1987) of 1 through 5-year zero coupon bond prices. From these data, we construct five 1-year forward rates
One limitation of the Fama-Bliss data is that it is not possible to construct excess returns for holding periods different from those. For that reason, we also use the Gürkaynak, Sack, and Wright (2007) dataset (GSW) which contains monthly records of zero coupon treasury yields with yearly maturities
2.2 The predictive 1-year excess return regressions
Cochrane and Piazzesi (2005) run predictive regressions of 1-year excess returns on government bonds on the five Fama-Bliss forward rates at the beginning of the period. Their main result is that a single linear combination of the five Fama-Bliss forward rates (the RFF) predicts 1-year excess returns of bonds of each maturity with an R2 of about 0.35, and up to 0.44 if one includes lagged forward rates as regressors. Furthermore, the RFF outperforms the first three PC extracted from the yield curve in terms of forecasting ability.
CP consider a single-factor predictive regression of excess returns of the form
where n = 2, 3, 4, 5 and
where the coefficients 𝜷n are estimated independently for each maturity. CP show that estimating the one-factor model (1) leads to almost identical R2 s than estimating the unrestricted regressions (2). Moreover, the RFF has more predictive power than the first three PC from bond yields.
Table 1 reports the R2s of regressions of the average 1-year excess returns (across maturities) on the return forecasting factor and the principal components. The left panel of the table, labeled Fama-Bliss data, uses the five Fama-Bliss discount prices to perform the regressions. The right panel, labeled GSW data, displays the same results but using ten 1-year forward rates
R2s of 1-year excess return forecast regressions.
|Fama-Bliss data||GSW data|
|Regressor||Individual||Joint ||F-test||Individual||Joint ||F-test|
The table report the R2s of regressions of the average 1-year excess returns
These results suggest that variations in 1-year excess returns may be driven by a single factor since the R2 of the regressions using only the RFF is at least as good as those using the principal components. Therefore, Cochrane and Piazzesi (2005) conclude that conventional three-factor models of the yield curve miss an important component of bond returns and that researchers should augment the model by adding a fourth (RFF) factor as the main driver of bond risk premia.
2.3 The predictive h-year excess return regressions
To be regarded as the main driver of risk premia, the (1-year) RFF should also summarize variations in expected bond returns considering holding periods different from 1 year. In this section we show that it does not. The 1-year RFF outperforms the first three PC only in 1-year excess return regressions such as those in Table 1. In regressions of excess returns for holding period longer than a year, the first three PC provide a better fit of excess returns. And for sufficiently long holding periods, the PC have three times the explanatory power of the RFF.
We first regress the average 1-year excess return
Next, we run predictive regressions of the expected h-year excess return averaged across available maturities
Figure 1 displays the R2 of the regressions of the excess holding returns on the return forecasting factors using the Fama-Bliss dataset. For example, the top left panel of the figure reports the R2 of the regressing the 1- through 4-year ahead excess returns on the 1-year-ahead return forecasting factor (dotted-circled line) and on the first three principal components (dotted-squared line) as a function of the holding period
In a similar fashion, the upper right panel of Figure 1 shows the R2 of regressions of the average h-period excess holding returns (
This pattern is reinforced when we consider the GSW database. Figure 2 is identical to Figure 1 but using a larger set of holding periods h. For example, the upper left panel of the figure shows the R2 of regressions of the average excess holding returns
The observation that the 1-year-ahead return forecasting factor performs poorly when attempting to explain excess holding returns for other horizons can be rationalized in several ways. The 1-year expected excess return obtained from the regressions can be thought as a model-free version of the expected excess return obtained from a structural model or an affine model of bond prices. For example, Bansal and Shaliastovich (2013) propose an equilibrium long-run risk model with Epstein-Zin preferences in which bond risk premia depend on expected economic growth and inflation. This model is able to reproduce 1-year excess returns remarkably well. Yet, since in the model investors prefer early resolution of uncertainty, risks that affect returns in the short run can be very different from those that affect long-run returns.
Then we can understand the previous results in two alternative ways. The first is that the determinants of short and long run risks are different. Therefore, trying to fit short and long period excess returns using a single 1-year return forecasting factor is most likely deemed to fail. Alternatively, since the Bansal and Shaliastovich (2013) model predicts a positive relation between bond risk premia and the preference for early resolution of uncertainty, estimating the model using data for different holding horizons will most likely generate different measures of risk premia. The empirical counterpart of this observation is that the h-year return forecasting factor will be a good predictor of the h-year excess return but not necessarily for other holding horizons. Trying to fit a 4-year excess return using the 1-year return forecasting factor is akin to wrongly calibrating the model.
A possible exception to the previous argument is when there exists a single factor structure in all excess returns, a property that is not supported by the data as we show in Section 2.4. But before doing that we evaluate the out-of-sample forecast properties of the empirical models discussed above.
We now analyze the relative merits of the 1-year-ahead RFF, the first three PCs, and a constant in terms of their ability to forecast excess returns.7 We compute the out-of-sample forecasts based on a series of recursive forecasts beginning in 2003M12 and extending through 2014M12 for a forecast horizon of 1 month and for holding periods of 12, 24, 36, and 48 months.8
Table 2 reports the results of this exercise. The forecasting ability of the two empirical models (using the RFF and the PCs) is poor, and both are outperformed by a constant for all holding periods. Which of the two models (RFF or PCs) forecast returns better depend on which holding horizon is forecast. Perhaps surprisingly, the principal components outperform the 1-year-ahead RFF only for a holding period of 1 year. Figure 3 displays the out-of-sample forecasts of the excess returns for the four holding periods. Overall, these results show the poor forecasting performance of both models.
Out-of-sample forecast accuracy.
|Horizon||1-year RFF||Principal components||Constant|
|h = 12||3.305||3.065||1.163|
|h = 24||1.944||2.018||1.600|
|h = 36||1.136||1.315||1.007|
|h = 48||0.658||0.814||0.593|
The table shows the root-mean-squared errors of 1-month-ahead forecasts of excess bond returns for different holding periods
2.4 A single factor structure in excess returns?
Cochrane and Piazzesi (2008) propose to augment the standard affine three-factor model by including an additional (return forecasting) factor that drives fluctuations in all excess returns. The principal components would then only serve to match the cross-section of yields but not to influence bond returns. The relevance of this result is that, to the extent that this dimensionality reduction of the data holds for all the possible manifestations of bond risk premia, models of the term structure of interest rates could be greatly simplified: while the usual first three principal components may accurately explain the cross-section of bond yields, compensation for risk would depend only on a single (return forecasting) factor.
Although CP actually statistically reject the single factor model in 1-year excess return, they claim that, nevertheless, the RFF captures all the economically relevant movements in 1-year excess returns.9 Not surprisingly, we also reject a single factor model in excess returns. But we also find that the single factor model does not capture all the economically relevant variations in excess returns once we consider different holding periods.
Using average excess returns constructed from the Fama-Bliss data, we estimate by maximum likelihood an unrestricted system of equations
and a restricted, one-factor, system of equations given by
assuming normally distributed errors. A likelihood ratio test strongly rejects the single factor model with a likelihood ratio statistic of 133.2 (distribututed as chi-square with 12 degrees of freedom).10
Moreover, the statistical rejection is economically relevant. Performing an eigenvalue decomposition of the expected returns obtained from equation (5), we find that the first principal component explains less than 90 percent of the variance of expected returns, much smaller than the 99.5 percent reported in CP when considering only 1-year excess returns. In other words, a single factor model fails to capture all the economically relevant variations in expected returns.
3 The return forecasting factor and bond yields
In this section we run regressions of bond yields on the RFF and the first three PC. If the RFF is a relevant state variable in a term structure model, we would also expect it to be a significant predictor of bond yields.
We estimate the regressions
To deal with this potential bias, we use the t-tests developed by Ibragimov and Müller (2010) and recommended by Bauer and Hamilton (2018) in the context of return regressions. We partition the sample into q subsamples and estimate regression (7) in each of them. If the estimates across subsamples are approximately independent, we can test significance of individual coefficients using a t-test with q − 1 degrees of freedom. In our setting, however, this independence assumption is unlikely to hold because the PC and RFF are both constructed using the entire sample. Hence, we divide the sample into q = 8 equally sized subsamples but construct the PC and RFF separately for each subsample. We then estimate regression (7) for each subsample, and perform the t-test of the null hypothesis that the coefficient on the RFF is zero.
The upper left panel of Figure 4 plots the t-statistics of the null hypothesis that
4 Concluding remarks
Cochrane and Piazzesi (2005) claim that variations in 1-year excess bond returns of all maturities can be summarized by a single common factor. We reevaluate this dimensionality reduction in bond risk premia. What differentiates our analysis from other papers is that we consider different holding horizons. That is, we analyze different manifestations of bond risk premia besides 1-year excess bond returns.
We find that a given return forecasting factor only outperforms the first principal components for the horizon for which it has been constructed, and that the first three principal components predict bond returns better than the return forecasting factor, particularly for long holding horizons. In addition, our results suggest that there is no single factor model that captures all the economically relevant variations in bond returns. Moreover, the coefficient of the return forecasting factor in unrestricted regressions of bond yields is insignificant for all bond maturities after controlling for the principal components. When considering out-of-sample forecasts of excess returns, however, we find that the return forecasting factor and the principal components predict poorly and are dominated by a constant for all holding horizons.
Taken together, our results suggest that the restrictions necessary to summarize bond risk premia with a single factor are not supported by the data.
This appendix describes the algorithm that we use to compute the recursive forecasts of excess returns. For the forecasting exercise we consider only bonds with yearly maturities
We consider two models of average excess returns (across bond maturities). The first uses on-year ahead return forecasting factor of Cochrane and Piazzesi (2005); the second, the first three principal components of bond yields,
Let ft to denote the
We now describe how to produce recursive out-of-sample forecasts for the empirical models (9) and (10). We focus exclusively on the Fama-Bliss data set and consider four holding horizons,
We first fit in-sample a vector autoregression of the principal components with two lags (selected using the Schwarz information criterion),
To produce out-of-sample forecasts of the principal components, we reestimate recursively the vector autoregression adding one observation at a time. Begin by setting
- Construct the forecast of the excess return
- –Given data for
, construct the forecast of the principal components,
- –Using the forecast of the principal components, compute the forecast of the RFF,
- –Compute the forecast of the excess returns using the RFF and the first three PCs,
- –Given data for
- Update sample and parameters:Add one observation to the sample, so that the sample period is now
. Using the new sample, reestimate the principal components, the factor loadings , and update the coefficients , , , , , , , , , and for estimating the regression equations
- Set k → k + 1 and return to step 1. The algorithm stops when
Additional figures and tables
Figure 5 displays the estimated coefficients
A single factor structure in excess returns?
|Parameters of the unrestricted excess return regression|
|Log-likelihood = 2208.896 Akaike = 10.337 Schwarz = 10.564|
|Parameters of the restricted single-factor excess return regression|
|Log-likelihood = 2275.493 Akaike = 10.590 Schwarz = 10.703|
|LR statistic = 133.193 p-value = 0.0000|
Bansal, R., and I. Shaliastovich. 2013. “A Long-Run Risks Explanation of Predictability Puzzles in Bond and Currency Markets.” The Review of Financial Studies 26: 1–33.
Bauer, M. D., and J. D. Hamilton. 2018. “Robust Bond Risk Premia.” The Review of Financial Studies 31: 399–448.
Cochrane, J. H., and M. Piazzesi. 2005. “Bond Risk Premia.” American Economic Review 95: 138–160.
Cochrane, J. H., and M. Piazzesi. 2008. “Decomposing the Yield Curve.” Technical Report 18, Society for Economic Dynamics.
Fama, E., and R. R. Bliss. 1987. “The Information in Long-Maturity Forward Rates.” American Economic Review 77: 680–692.
Gürkaynak, R. S., B. Sack, and J. H. Wright. 2007. “The U.S. Treasury yield Curve: 1961 to the Present”. Journal of Monetary Economics 54: 2291–2304.
Ibragimov, R., and U. K. Müller. 2010. “t-Statistic Based Correlation and Heterogeneity Robust Inference.” Journal of Business & Economic Statistics 28: 453–468.
Müller, U. K. 2014. “Hac Corrections for Strongly Autocorrelated Time Series.” Journal of Business & Economic Statistics 32: 311–322.
Svensson, L. E. 1994. “Estimating and Interpreting Forward Interest Rates: Sweden 1992–1994.” Working Paper 4871, National Bureau of Economic Research.
The online version of this article offers supplementary material (DOI:
Bauer and Hamilton (2018) define the “spanning hypothesis” as the claim that the first three PC contain all the relevant information for predicting future interest rates and bond returns. That paper also contains an extensive list of papers that claim to find violations of the spanning hypothesis. However, they have questioned the value of variables other than the first three PC as predictors of bond returns. The conventional tests that support including those variables in predictive regressions are subject to small-sample distortions that lead to oversized tests.
The traditional Fama and Bliss (1987) regressions of excess returns on forwards rates attain R2s of about 0.17, while regressions on the first three principal components deliver R2s of about 0.27.
Cochrane and Piazzesi (2005) construct the return forecasting factor
In particular, consider the regressions
There are two differences between Figure 1 and Figure 2. First, while the Fama-Bliss dataset only allows us to construct excess returns for holding periods of 1, 2, 3 and 4 years, the GSW data allows us to consider a richer set of holding periods, ranging from 1 month (1/12 years) through 108 months (9 years). Hence, relative to Figure 1, the average excess holding returns used in Figure 2 are constructed using a larger set of bond maturities and the figure displays excess returns for a finer and longer set of holding periods. The results using the GSW data restricted to holding periods of 1 through 4 years are very similar to those in Figure 1 (results available upon request).
We include a constant as a benchmark to capture no predictability of excess returns.
From the unrestricted regressions (2), CP compute an eigenvalue decomposition of 1-year expected excess returns
In effect, using GSW data we find that the RFF is correlated with the missing PCs.