Bond risk premia and the return forecasting factor

Agustin Gutierrez; Constantino Hevia; Martin Sola

doi:10.1515/snde-2018-0009

Publicly Available Published by De Gruyter April 24, 2019

Bond risk premia and the return forecasting factor

Agustin Gutierrez , Constantino Hevia and Martin Sola

From the journal Studies in Nonlinear Dynamics & Econometrics

https://doi.org/10.1515/snde-2018-0009

Abstract

The return forecasting factor is a linear combination of forward rates that seems to predict 1-year excess bond returns of bond of all maturities better than traditional measures obtained from the yield curve. If this single factor actually captures all the relevant fluctuations in bond risk premia, then it should also summarize all the economically relevant variations in excess returns considering different holding periods. We find that it does not. We conclude that including the return forecasting factor as the main driver of risk premia in a term structure model, as has been suggested, is not supported by the data.

Keywords: affine term structure models; bond risk premia; excess returns; return forecasting factor

1 Introduction

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 R² 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 pt(n) denote the log price of an n-year zero coupon bond, where time and maturity are both expressed in years. The log yield is yt(n)=−pt(n)/n and the log forward rate for loans starting at t + n and maturing at t + n + h is ft(n,n+h)=pt(n)−pt(n+h). Let

r t → t + h ( n ) = ( p t + h ( n − h ) − p t ( n ) ) / h

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

r x t → t + h ( n ) = r t → t + h ( n ) − y t ( h ) .

Finally, we denote the average h-year excess return across maturities by rx¯t→t+h=1H∑n>hrxt→t+h(n), where H is the number of bonds with maturities longer than h periods ahead.

2.1 Data

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 (yt(1),ft(1,2),ft(2,3),ft(3,4),ft(4,5)) and 1- through 4-year excess holding returns.

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 (1/12,2/12,...,120/12). These yields are generated using the extended Nelson-Siegel parametrization of the yield curve proposed by Svensson (1994), whose coefficients are available on a daily basis. From the fitted curves we recover bond prices, yields, and forward rates of any maturity ranging from 1 month to 10 years.

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 R² 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

(1) r x t → t + 1 ( n ) = b n ( γ ′ f t ) + ϵ t + 1 ( n ) ,

where n = 2, 3, 4, 5 and ft=[1, yt(1), ft(1,2), ft(2,3), ft(3,4), ft(4,5)]′ is a vector with a constant and five Fama-Bliss forwards rates at time t. Equation (1) imposes a single factor structure on 1-year excess returns.^[3] The alternative, unrestricted, regressions are given by

(2) r x t → t + 1 ( n ) = β n ′ f t + ϵ t + 1 ( n ) ,

where the coefficients 𝜷_n are estimated independently for each maturity. CP show that estimating the one-factor model (1) leads to almost identical R² 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 R²s 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 ft=[1, yt(1), ft(1,2), ..., ft(9,10)]′ constructed using the Gürkaynak, Sack, and Wright (2007) data. In both panels, the column labeled “F-test” reports the p-values of the null hypothesis that the coefficient δ_j on the jth principal component is equal to zero.^[4]

Table 1:

R²s of 1-year excess return forecast regressions.

	Fama-Bliss data			GSW data
Regressor	Individual	Joint PC1:PCi	F-test	Individual	Joint PC1:PCi	F-test
RFF	0.351	−	−	0.358	−	−
PC₁	0.028	0.028	−	0.033	0.033	−
PC₂	0.213	0.242	0.0000	0.169	0.202	0.0000
PC₃	0.024	0.266	0.0001	0.059	0.261	0.0000
PC₄	0.085	0.351	0.0000	0.009	0.269	0.0195
PC₅	0.000	0.351	0.9140	0.021	0.291	0.0002
PC₆				0.026	0.318	0.0000
PC₇				0.013	0.331	0.0028
PC₈				0.014	0.344	0.0021
PC₉				0.010	0.355	0.0066
PC₁₀				0.003	0.358	0.1561

The table report the R²s of regressions of the average 1-year excess returns rx¯t→t+1 on the RFF and the PC of bond yields, and the p-values of F-tests of the null hypothesis that a model estimated with j − 1 principal components is as good as a model with j principal components for j=2,3,.... The columns labeled Fama-Bliss data are based on the five Fama-Bliss discount bond prices. The columns labeled GSW data uses the treasury yields constructed by Gürkaynak, Sack, and Wright (2007). The sample period is that used by Cochrane and Piazzesi (2005), Jan-1964 to Dec-2003.

These results suggest that variations in 1-year excess returns may be driven by a single factor since the R² 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 rx¯t→t+1 on the forward rates f_t using both the FB and GSW datasets. We call this the “1-year-ahead” RFF because it maximizes the forecasting performance in a regression of the average 1-year excess return. Likewise, we label the “m-year-ahead” RFF as the fitted value of a regression of the average m-year excess return rx¯t→t+m on the forward rates f_t. Let RFFt(m) denote the m-year-ahead return forecasting factor computed using this procedure (below we consider m=1,2,3,4).

Next, we run predictive regressions of the expected h-year excess return averaged across available maturities rx¯t→t+h, on the m-year return forecasting factor RFFt(m), and on the first three principal components PC1t, PC2t, and PC3t extracted from the same datasets,^[5]

(3) r x ¯ t → t + h = α 0 ( h ) + α 1 ( h ) R F F t ( m ) + ϵ t + h

(4) r x ¯ t → t + h = δ 0 ( h ) + δ 1 ( h ) P C 1 t + δ 2 ( h ) P C 2 t + δ 3 ( h ) P C 3 t + η t + h .

Figure 1 displays the R² 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 R² 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 h=1,2,3,4. While the RFF has more predictive ability than the PC at the 1-year holding horizon, the relation is reversed at longer holding periods: the RFF explains less than 10 percent of the variability of the 4-year excess return while the PC explain about 25 percent of the variance.

$Figure 1: R2 of the average h-year excess holding return regressions using Fama-Bliss yields. The figure shows the R2s of a regression of the average h-period excess holding returns (h = 1, 2, 3, 4) on the return forecasting factor and on the first three principal components $PC_{1t}$PC1t, $PC_{2t}$PC2t, and $PC_{3t}$PC3t together. In the upper left panel the return forecasting factor is computed using a regression of 1-year excess returns on the current forwards; in the upper right panel the return return forecasting factor is computed with a regression of 2-year excess returns on the current forwards, and so forth.$

Figure 1:

R² of the average h-year excess holding return regressions using Fama-Bliss yields. The figure shows the R²s of a regression of the average h-period excess holding returns (h = 1, 2, 3, 4) on the return forecasting factor and on the first three principal components PC1t, PC2t, and PC3t together.

In the upper left panel the return forecasting factor is computed using a regression of 1-year excess returns on the current forwards; in the upper right panel the return return forecasting factor is computed with a regression of 2-year excess returns on the current forwards, and so forth.

In a similar fashion, the upper right panel of Figure 1 shows the R² of regressions of the average h-period excess holding returns (h=1,2,3,4) on the 2-year ahead return forecasting factor, RFF(2), and the PC, and likewise for the lower two panels of the figure but using RFFt(3)andRFFt(4). Again, while the h-year-ahead RFF outperforms the PC at the h-year holding horizon, the PC tend to do better than the RFF at other holding horizons.

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 R² of regressions of the average excess holding returns rx¯t→t+h for h=1/12,2/12,...,108/12 on the 1-year-ahead return forecasting factor RFFt(1) and the first three principal components, as a function of the holding period h. Note that for longer holding horizons, the predictive power of the principal components is more than three times larger than that of the 1-year-ahead return forecasting factor and twice as large as that of the 2-year-ahead return forecasting factor. Finally, the other two return forecasting factors, RFFt(3) and RFFt(4), do not provide any predictive improvement over the principal components.^[6]

Figure 2:

R² of h-year excess holding return regressions using GSW data.

The figure shows R²s of regressions of the average h-year excess holding returns on the 1- though 4-year ahead RFF and on the first three PC. In the upper upper left panel the RFF is computed using a regression of 1-year excess returns on 1- through 10-year forwards; in the upper right panel the RFF is computed with a regression of 2-year excess returns on 1- through 10-year forwards, and so forth.

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.

Table 2:

Out-of-sample forecast accuracy.

	Root-mean-squared-error
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 h=12,24,36,48. The column 1-year-ahead RFF computes the forecasts using Cochrane and Piazzesi (2005) RFF, the column labeled Principal components does so using the first three principal components of bond yields. In the last column the forecast of the excess return is the unconditional mean.

Figure 3:

Out-of-sample performance of the empirical models.

One-month ahead out of sample forecasts of h-period excess returns of the 1-year-ahead RFF, the principal components, and a constant.

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

(5) [ r x ¯ t → t + 1 r x ¯ t → t + 2 r x ¯ t → t + 3 r x ¯ t → t + 4 ] = [ α 1 α 2 α 3 α 4 ] + [ β 1 , 1 β 1 , 2 β 1 , 3 β 1 , 4 β 1 , 5 β 2 , 1 β 2 , 2 β 2 , 3 β 2 , 4 β 2 , 5 β 3 , 1 β 3 , 2 β 3 , 3 β 3 , 4 β 3 , 5 β 4 , 1 β 4 , 2 β 4 , 3 β 4 , 4 β 4 , 5 ] [ y t ( 1 ) f t ( 1 , 2 ) f t ( 2 , 3 ) f t ( 3 , 4 ) f t ( 4 , 5 ) ] + [ ϵ t + 1 ( 1 ) ϵ t + 2 ( 2 ) ϵ t + 3 ( 3 ) ϵ t + 4 ( 4 ) ] ,

and a restricted, one-factor, system of equations given by

(6) [ r x ¯ t → t + 1 r x ¯ t → t + 2 r x ¯ t → t + 3 r x ¯ t → t + 4 ] = [ α 1 α 2 α 3 α 4 ] + [ γ 1 γ 2 γ 3 γ 4 γ 5 b 2 γ 1 b 2 γ 2 b 2 γ 3 b 2 γ 4 b 2 γ 5 b 3 γ 1 b 3 γ 2 b 3 γ 3 b 3 γ 4 b 3 γ 5 b 4 γ 1 b 4 γ 2 b 4 γ 3 b 4 γ 4 b 4 γ 5 ] [ y t ( 1 ) f t ( 1 , 2 ) f t ( 2 , 3 ) f t ( 3 , 4 ) f t ( 4 , 5 ) ] + [ ϵ t + 1 ( 1 ) ϵ t + 2 ( 2 ) ϵ t + 3 ( 3 ) ϵ t + 4 ( 4 ) ] ,

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

(7) y t ( n ) = ζ 0 ( n ) + λ ( n ) ( γ ^ ′ f t ) + ζ 1 ( n ) P C 1 t + ζ 2 ( n ) P C 2 t + ζ 3 ( n ) P C 3 t + η t ( n ) ,

for n=1/12,2/12...,120/12 and test the significance of the coefficient λ⁽ⁿ⁾. To test hypotheses, it is usual to perform a Newey-West correction due to the serial correlation of the regression residuals. This type of correction, however, may be inappropriate for two reasons. First, the RFF is likely to be correlated with the residuals: equation (7) omits PC that are relevant to explain yields, perhaps with a small but still significant coefficient. Those missing PC are in the residual. This omission is not a problem in a regression of yields on the first three PC since all PC are orthogonal. But if we add the RFF to the regression, which is constructed using the data from which we extract the PC, it is likely that it will be correlated with the missing PC and, hence, with the residuals.^[11] And second, the right hand variables of the regression are persistent and conventional HAC corrections work poorly with strongly dependent data (Müller 2014).

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 λ(n)=0 for each bond maturity. The tests do not reject the null hypothesis for any maturity. That is, the RFF is not a significant predictor of bond yields. In contrast, the first three PC are highly significant predictors of bond yields.

Figure 4:

Testing if factors do not explain yields using Ibrahimov-Müller t-test.

Robust t-statistics and 0.025, 0.975 critical values (dashed lines) of the null hypothesis that the coefficients of the RFF and of the first three PC are zero in regression (7).

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.

Appendix A

Out-of-sample forecasts

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 n=1,2,...,N. Let T be the total number of observations and T_s < T the number of observations used for estimation Thus, T − T_s observations are used for forecasting.

Let yt=[yt(1),yt(2),…,yt(N)] be a row vector of yields at time t, and Λ_s an N × N matrix with columns containing the eigenvectors covariance matrix of y. Here, the subindex s is used to denote that the loadings are constructed using in-sample data. In addition, let PC_t be a 1×N row vector with the principal components of yields at time t, that is, PCt=ytΛs.

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,

(8) r x ¯ t → t + h = α 0 ( h ) + α 1 ( h ( γ 0 + γ 1 y t ( 1 ) + γ 2 f t ( 1 , 2 ) + γ 3 f t ( 2 , 3 ) + ⋯ + γ n f t ( n − 1 , n ) ) + ϵ ¯ t + h

(9) r x ¯ t → t + h = δ 0 ( h ) + δ 1 ( h ) P C 1 t + δ 2 ( h ) P C 2 t + δ 3 ( h ) P C 3 t + η ¯ t + h .

Let f_t to denote the 1×N row vector of 1 year forwards rates. Since the forwards constitute N different linear combinations of the yields y_t, the set of forwards and the set of yields span the same information. Moreover, since the N principal component span the same information of the set of yields, we can express model (8) in terms of the vector of principal components,

(10) r x ¯ t → t + h = α 0 ( h ) + α 1 ( h ) ( ω 0 + ω 1 P C 1 t + ω 2 P C 2 t + ⋯ + ω N P C N t ) + ϵ ¯ t + h

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, h={12,24,36,48}. For the in-sample regressions, use observations from 1964M1 through 2003M12; to perform out-of-sample forecasts, we use data from 2004M1 through 2014M12.

We first fit in-sample a vector autoregression of the principal components with two lags (selected using the Schwarz information criterion),

P C t = Φ 0 + Φ 1 P C t − 1 + Φ 2 P C t − 2 + ξ t .

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 Φ^0,0=Φ0, Φ^1,0=Φ1, and Φ^2,0=Φ2, for the regression of the principal components; ω^j,0=ω0 for j=1,2,...,N, for the equation relating the RFF with the principal components; α^0,0(h)=α0(h) and α^1,0(h)=α1(h), for the equation relating the excess returns with the RFF; and δ^0,0(h)=δ0(h), δ^1,0(h)=δ1(h), and δ^2,0(h)=δ2(h), for the equation relating the excess returns with the first three principal components. Then, set k = 1 and for holding period h={12,24,36,48} do the following,

Construct the forecast of the excess return rx¯Ts+k→Ts+k+h at Ts+k:
1. Given data for t=1,...,Ts+k−1, construct the forecast of the principal components,
  
  P C ^ T s + k = Φ ^ 0 , k − 1 + Φ ^ 1 , k − 1 P C T s + k − 1 + Φ ^ 2 , k − 1 P C T s + k − 2 ,
2. Using the forecast of the principal components, compute the forecast of the RFF,
  
  R F F ^ T s + k = ω ^ 0 , k − 1 + ω ^ 1 , k − 1 P C ^ 1 , T s + k + ω ^ 2 , k − 1 P C ^ 2 , T s + k + ⋯ + ω ^ n , k − 1 P C ^ n , T s + k
3. Compute the forecast of the excess returns using the RFF and the first three PCs,
  
  r x ¯ ^ T s + k → T s + k + h = α ^ 0 , k − 1 ( h ) + α ^ 1 , k − 1 ( h ) R F F ^ T s + k r x ¯ ^ T s + k → T s + k + h = δ ^ 0 , k − 1 ( h ) + δ ^ 1 , k − 1 ( h ) P C ^ 1 , T s + k + δ ^ 2 , k − 1 ( h ) P C ^ 2 , T s + k + δ ^ 3 , k − 1 ( h ) P C ^ 3 , T s + k .
Update sample and parameters:

Add one observation to the sample, so that the sample period is now t=1,2,...Ts+k. Using the new sample, reestimate the principal components, the factor loadings Λs,k, and update the coefficients α^0,k(h), α^1,k(h), δ^0,k(h), δ^1,k(h), δ^2,k(h), δ^3,k(h), Φ^0,k, Φ^1,k, Φ^2,k, and ω^j,k(h) for j=0,1,...,N estimating the regression equations

r x ¯ t → t + 12 = ω ^ 0 , k ( h ) + ω ^ 1 , k ( h ) P C 1 t + . . . + ω ^ N , k ( h ) P C N t + ϵ ¯ t → t + 12 R F F t = ω ^ 0 , k ( h ) + ω ^ 1 , k ( h ) P C 1 t + . . . + ω ^ N , k ( h ) P C N t r x ¯ t → t + h = α ^ 0 , k ( h ) + α ^ 1 , k ( h ) R F F t + ϵ ¯ t + h r x ¯ t → t + h = δ ^ 0 , k ( h ) + δ ^ 1 , k ( h ) P C 1 t + δ ^ 2 , k ( h ) P C 2 t + δ ^ 3 , k ( h ) P C 3 t + η ¯ t + h P C t = Φ ^ 0 , k + Φ ^ 1 , k P C t − 1 + Φ ^ 2 , k P C t − 2 + ξ t .
Set k → k + 1 and return to step 1. The algorithm stops when Ts+k+12=T.

Appendix B

Additional figures and tables

Figure 5 displays the estimated coefficients α^1(h),δ^j(h),δ^2(h), and δ^3(h) of the regressions (3) and (4) using the Fama-Bliss dataset. Each line in the left panel represent the estimated coefficient as a function of the holding period h for a given return forecasting factor RFFt(m) for m=1,2,3,4. For example, the blue circled-dotted line represent the estimated coefficients α^1(h) as a function of the holding period h when using the 1-year-ahead return forecasting factor RFFt(1) as the regressor in equation (3). The other lines display the same coefficients but using the return forecasting factor RFFt(2),RFFt(3), and RFFt(4). Likewise, the right panel of the figure displays the estimated coefficients δ^1(h),δ^2(h), and δ^3(h) of the regression (3) on the principal components, for h=1,2,3,4.

Figure 5:

Estimated coefficients using the Fama-Bliss data set.

Figure 6:

Estimated coefficients using the GSW data set.

Figure 6 is identical to Figure 5 but using the GSW dataset and holding periods from 1 through 9 years ahead.

Table 3 reports the estimates of the constrained and unconstrained system of equations (5) and (6) that we use to test the null hypothesis of a single factor structure on excess returns.

Table 3:

A single factor structure in excess returns?

Parameters of the unrestricted excess return regression
α	β
[ 4.499 ( 0.747 ) − 3.336 ( 0.517 ) − 3.046 ( 0.459 ) − 0.977 ( 0.183 ) ]	[ 0.531 − 3.597 − 7.306 − 1.416 10.415 ( 0.079 ) ( 0.815 ) ( 1.078 ) ( 1.286 ) ( 1.361 ) 0.431 − 0.714 − 2.744 − 2.829 1.583 ( 0.062 ) ( 0.518 ) ( 0.752 ) ( 1.026 ) ( 1.048 ) 0.397 − 0.564 − 2.164 − 2.747 0.897 ( 0.056 ) ( 0.479 ) ( 0.694 ) ( 0.930 ) ( 0.943 ) 0.155 − 0.178 − 0.633 0.058 − 0.316 ( 0.025 ) ( 0.232 ) ( 0.386 ) ( 0.418 ) ( 0.417 ) ]
Log-likelihood = 2208.896 Akaike = 10.337 Schwarz = 10.564
Parameters of the restricted single-factor excess return regression
α	γ	b
[ − 5.712 ( 0.747 ) − 2.736 ( 0.429 ) − 2.382 ( 0.376 ) − 0.486 ( 0.134 ) ]	[0.698(0.082)−2.794(0.790)−6.336(0.989)−2.685(1.310)6.871(1.395)]	[ 0.509 ( 0.036 ) 0.447 ( 0.033 ) 0.116 ( 0.017 ) ]
Log-likelihood = 2275.493 Akaike = 10.590 Schwarz = 10.703
LR statistic = 133.193 p-value = 0.0000

References

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.10.3386/w18357Search in Google Scholar

Bauer, M. D., and J. D. Hamilton. 2018. “Robust Bond Risk Premia.” The Review of Financial Studies 31: 399–448.10.3386/w23480Search in Google Scholar

Cochrane, J. H., and M. Piazzesi. 2005. “Bond Risk Premia.” American Economic Review 95: 138–160.10.3386/w9178Search in Google Scholar

Cochrane, J. H., and M. Piazzesi. 2008. “Decomposing the Yield Curve.” Technical Report 18, Society for Economic Dynamics.10.2139/ssrn.1333274Search in Google Scholar

Fama, E., and R. R. Bliss. 1987. “The Information in Long-Maturity Forward Rates.” American Economic Review 77: 680–692.Search in Google Scholar

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.10.1016/j.jmoneco.2007.06.029Search in Google Scholar

Ibragimov, R., and U. K. Müller. 2010. “t-Statistic Based Correlation and Heterogeneity Robust Inference.” Journal of Business & Economic Statistics 28: 453–468.10.1198/jbes.2009.08046Search in Google Scholar

Müller, U. K. 2014. “Hac Corrections for Strongly Autocorrelated Time Series.” Journal of Business & Economic Statistics 32: 311–322.10.1080/07350015.2014.931238Search in Google Scholar

Svensson, L. E. 1994. “Estimating and Interpreting Forward Interest Rates: Sweden 1992–1994.” Working Paper 4871, National Bureau of Economic Research.10.3386/w4871Search in Google Scholar

Supplementary Material

The online version of this article offers supplementary material (DOI: https://doi.org/10.1515/snde-2018-0009).

Published Online: 2019-04-24

Bond risk premia and the return forecasting factor

Abstract

1 Introduction

2 The return forecasting factor and excess bond returns

2.1 Data

2.2 The predictive 1-year excess return regressions

2.3 The predictive h-year excess return regressions

2.4 A single factor structure in excess returns?

3 The return forecasting factor and bond yields

4 Concluding remarks

Appendix A

Out-of-sample forecasts

Appendix B

Additional figures and tables

References

Supplementary Material

Journal and Issue

Articles in the same Issue