The Federal Reserve Board eagle logo links to home page

Skip to: [Printable Version (PDF)] [Bibliography] [Footnotes]
Finance and Economics Discussion Series: 2010-14 Screen Reader version

Variance Risk Premia, Asset Predictability Puzzles, and Macroeconomic Uncertainty1

Hao Zhou2
Federal Reserve Board
First Draft: May 2009
This Version: January 2010

Keywords: Short-run predictability, variance premium dynamics, equity premium puzzle, bond risk premia, credit spread puzzle, macroeconomic uncertainty, recursive preference.

Abstract:

This paper presents predictability evidence from the difference between implied and expected variances or variance risk premium that: (1) the variance difference measure predicts a significant positive risk premium across equity, bond, and credit markets; (2) the predictability is short-run, in that it peaks around one to four months and dies out as the horizon increases; and (3) such a short-run predictability is complementary to that of the standard predictor variables--P/E ratio, forward spread, and short rate. These findings are potentially justifiable by a general equilibrium model with recursive preference that incorporates stochastic economic uncertainty. Calibration evidence suggests that such a framework is capable of reproducing the variance premium dynamics, especially its high skewness and kurtosis, without introducing jumps. The calibrated model can also qualitatively explain the equity premium puzzle and the bond risk premia in short horizons.

JEL classification: G12, G13, G14.



1 Introduction

Option implied volatility, such as the Chicago Board Option Exchange's VIX index, is widely viewed by investors as the market gauge of fear (Whaley, 2000).3 In recent research, the difference between the implied and expected volatilities has been interpreted as an indicator of the representative agent's risk aversion (Rosenberg and Engle, 2002; Bollerslev, Gibson, and Zhou, 2008a; Bakshi and Madan, 2006). An alternative interpretation is that the implied-expected variance difference, as a proxy for variance risk premium, is due to the macroeconomic uncertainty risk (Bollerslev, Tauchen, and Zhou, 2009; Drechsler and Yaron, 2008). Such an approach relies on the non-standard recursive utility framework of Epstein and Zin (1991) and Weil (1989), such that the consumption uncertainty risk commands a time-varying risk premium. This method follows the spirit of the long-run risks models as pioneered by Bansal and Yaron (2004), but focuses on the independent consumption volatility as a primary source of the time variation in market risk premia, while shutting down the long-run risk channel completely.

This paper demonstrates that the difference between implied and expected variances, as a measure for variance risk premium, provides a significant predictability for short-run equity returns, bond returns, and credit spreads. The documented return predictability peaks around one-to-four month horizons across these markets, and then dies out as the forecasting horizon increases. More importantly, such a short-term forecastability of market risk premia is complementary to the standard established predictors--P/E ratio, forward spread, and short rate, which are indicated by the asset pricing theory or empirical evidence; in that when combined together, the statistical significance of the variance risk premium proxy is rather increased or at least stable, instead of being crowded out by the standard predictor variables. This constitutes an important evidence that risk premia across major financial markets comove in the short-run, and such a common component seems to be intimately related to the variance risk premium, which is constructed from the high quality derivatives market data and high-frequency underlying market prices.

This type of common short-run risk factor may be a proxy for stochastic economic uncertainty or consumption volatility risk that varies independently with the consumption growth rate--the latter being the main focus of the long-run risk models (Bansal and Yaron, 2004).4 These empirical results may be consistent with a self-contained general equilibrium model incorporating the effects of such a time-varying economic uncertainty, where the uncertainty risk is priced only under the recursive utility function. Calibration evidence shows that such a framework can replicate quantitatively the observed skewness and kurtosis in variance premium dynamics, without introducing jumps into the endowment process and/or the volatility process (as in, e.g., Eraker and Shaliastovich, 2008; Drechsler and Yaron, 2008). More importantly, within such a calibration parameter setting, the model can qualitatively explain the equity premium puzzle and bond risk premia in short-horizons.

There is a fundamental link between the notion of option-implied volatility risk premium and the notion of variance risk premium embedded in underlying assets. Within the arbitrage pricing framework, stochastic volatility of equity market can only be priced, if its innovation is correlated with the market return innovation (Bates, 1996; Bakshi and Kapadia, 2003; Heston, 1993). There is a great deal of empirical evidence that equity volatility is negatively correlated with the equity returns, such that the negative volatility risk premium embedded in equity options provides a vital hedging service for the average investors. However, within a consumption-based asset pricing framework, without assuming any arbitrary statistical correlation between the volatility and consumption innovations, one need to endow the economic agents with a preference for an earlier resolution of uncertainty and a stochastic volatility-of-volatility. Under such a setup the variance risk premium embedded in equity, bond, and credit markets must be positive, as more risk requires more return compensation. And the positive variance risk premium embedded in underlying assets is entirely consistent with the negative volatility risk premium implied from option prices.

Economic uncertainty and its effect on asset prices can be examined with other techniques under the recursive preference structure. Bansal and Shaliastovich (2008a) and Shaliastovich (2008) introduce information learning into the long-run risk model, such that endogenously asset prices requires a compensation for jump risks. Chen and Pakos (2008) model the endowment growth rate as a Markov switching process with a constant volatility, where learning brings about an endogenous uncertainty premium. Drechsler (2008) applies the Knightian uncertainty about model misspecification with realistic asset dynamics to explain the observed option pricing puzzles. Lettau, Ludvigson, and Wachter (2008) use a Markov switching model to describe the long-swing changes in the constant consumption volatility--Great Moderation--and to draw implications for the declining equity risk premiums.5

In contrast, the approach taken here shuts down the the long-run component in consumption growth, and attributes the higher order time-variation in financial market risk premia to the stochastic volatility-of-volatility in consumption growth. Conceptually, the root of macroeconomic or consumption uncertainty can be traced back to the volatility shocks in total factor productivity and firm-specific technology in a real business cycle model (Bloom, Floetotto, and Jaimovich, 2009; Bloom, 2009). Such an approach treats the short-run economic uncertainty as a fundamental priced risk factor, without relying on informational learning, heterogeneous beliefs, or behavioral assumptions. The emphasis here is also complementary to the long-run risk perspective (Bansal and Yaron, 2004) or rare disaster angle (Gabaix, 2009; Barro, 2006), and is orthogonal to the channel of time-varying risk aversion (Campbell and Cochrane, 1999). Lastly, but not the least, this framework has sharp empirical testable hypotheses, since the economic uncertainty factor or stochastic volatility-of-volatility is uniquely identified by the variance risk premium.6

The rest of the paper will be organized as the follows, the next section defines the variance risk premium and describes its empirical measurement; Section 3 presents the main empirical evidence of the short-run predictability from implied-expected variance difference for risk premia in asset markets; the following section discusses a general equilibrium model with stochastic economic uncertainty and provides some calibration implications for explaining the short-run risk premia dynamics; and Section 5 concludes.


2 Variance Risk Premium and Empirical Measurement

The central empirical finding of this paper is that market risk premia have a common short-run component--variance risk premium--that is not directly observable. However, an empirical proxy can be constructed from the difference between model-free option-implied variance and the conditional expectation of model-free realized variance.

2.1 Variance Risk Premium: Definition and Measurement

To define the procedure in quantifying the model-free implied variance, let  C_t(T,K) denote the price of a European call option maturing at time  T with strike price  K, and  B(t,T) denote the price of a time  t zero-coupon bond maturing at time  T. As shown by Carr and Madan (1998); Britten-Jones and Neuberger (2000) among others, the market's risk-neutral expectation of the return variance between time  t and  t+1 conditional on time  t information, or the implied variance  IV_{t,t+1}, can be expressed in a "model-free" fashion as the a portfolio of European calls,

\displaystyle IV_{t,t+1} \ \equiv \ \operatorname{E}_t^Q \left( \mbox{Var} _{t,t+1}\right) \ = \ 2 \int_{0}^{\infty}{C_t\left(t+1,\frac{K}{B(t,t+1)}\right) - C_t\left(t,{K}\right) \over K^2}\,dK, (1)

which relies on an ever increasing number of calls with strikes spanning from zero to infinity.7

This equation follows directly from the classical result in Breeden and Litzenberger (1978), that the second derivative of the option call price with respect to strike equals the risk-neutral density, such that all risk neutral moments payoff can be replicated by the basic option prices (Bakshi and Madan, 2000). In practice,  IV_{t,t+1} must be constructed on the basis of a finite number of strikes; which turns out to be a fairly accurate approximation to the true (unobserved) risk-neutral expectation of the future market variance, under reasonable assumptions about the underlying asset dynamics (Carr and Wu, 2008; Bollerslev, Gibson, and Zhou, 2008a; Jiang and Tian, 2005).

In order to define the measure in quantifying the actual return variation, let  p_t denote the logarithmic price of the asset. The realized variance over the discrete  t to  t+1 time interval can be measured in a "model-free" fashion by

\displaystyle RV_{t,t+1} \equiv \sum_{j=1}^{n}\left[p_{t+\frac{j}{n}} - p_{t+\frac{j-1}{n}}\right]^2 \longrightarrow   Var\displaystyle _{t,t+1}, (2)

where the convergence relies on  n\rightarrow\infty; i.e., an increasing number of within period price observations.

As demonstrated in the literature (see, e.g., Barndorff-Nielsen and Shephard, 2002; Andersen, Bollerslev, Diebold, and Ebens, 2001), this "model-free" realized variance measure based on high-frequency intraday data offers a much more accurate ex-post observation of the true (unobserved) return variance than the traditional ones based on daily or coarser frequency returns. In practice, various market microstructure frictions invariably limit the highest sampling frequency that may be used in reliably estimating  RV_{t,t+1}.

The variance risk premium is defined as the difference between the ex-ante risk neutral expectation of the future return variance and the objective or statistical expectation of the return variance over the  [t,t+1] time interval,

\displaystyle VRP_{t} ~ \equiv ~ \operatorname{E}_t^Q \left( \mbox{Var} _{t,t+1} \right) ~ - ~ \operatorname{E}_t^P\left( \mbox{Var} _{t,t+1}\right), (3)

which is not directly observable in practice. To construct an empirical proxy for such a variance risk premium concept (3), one need to estimate various reduced-form counterparts of the risk neutral and physical expectations, i.e.,  \widehat{VRP}_{t} ~ \equiv ~ \widehat{\operatorname{E}}_t^Q \left( \mbox{Var} _{t,t+1} \right) ~ - ~ \widehat{\operatorname{E}}_t^P\left( \mbox{Var} _{t,t+1}\right). In practice, the risk-neutral expectation  \widehat{\operatorname{E}}_t^Q \left( \mbox{Var} _{t,t+1} \right) is typically replaced by the CBOE implied variance or VIX ^2 and the true variance    Var _{t,t+1} is replaced by its discretized realization  RV_{t,t+1}. However, methods for constructing the objective expectation  \widehat{\operatorname{E}}_t^P \left(\cdot\right) vary in practice.

One approach is to estimate a reduced-form multi-frequency auto-regression with multiple lags for  \widehat{\operatorname{E}}_t^P\left(RV_{t,t+1}\right) (Bollerslev, Tauchen, and Zhou, 2009). For more specific structural jump-diffusion processes, one could use the model-implied objective expectation (Todorov, 2009). Based on the evidence of implied-realized variance forecast efficiency and a potential two-factor structure in variance dynamics, Drechsler and Yaron (2008) use lagged implied and realized variances to forecast the expected variance. For forecasting purposes only, one could use time- t realized variance  RV_{t-1,t} (Bollerslev, Tauchen, and Zhou, 2009), which ensures that the variance risk premium proxy is in the time  t information set. Following the common practice in the variance swap market, Carr and Wu (2008) use the ex-post forward realized variance to substitute for the expected variance in characterizing the variance risk premium. When presenting the empirical findings, I will focus on the method that uses the twelve lag auto-regressive estimate.

2.2 Data Description and Summary Statistics

For the option-implied variance of the S&P500 market return, I use the end-of-month data of the Chicago Board of Options Exchange (CBOE) volatility index VIX ^2 on a monthly basis, as a risk-neutral expectation of return variance for the next 30 days. Following the literature, the monthly realized variance measure for the S&P500 index is the summation of the 78 within day five-minute squared returns covering the normal trading hours from 9:30am to 4:00pm plus the close-to-open return.

Here I consider three market risk premium measures with their traditional predictor variables. Specifically, monthly P/E ratios and index returns for the S&P500 are obtained from Standard & Poor's, bond returns and forward rates are from the monthly CRSP Fama t-bill data set with 1 to 6 month maturities, and AAA and BAA corporate bond spreads are from Moody's with Fama-Bliss risk-free interest rates from CRSP. The empirical analysis here is for the sample period from January 1990 to December 2008, when the new VIX index based on S&P500 index becomes available.

To give a visual illustration, Figure 1 plots the monthly time series of variance risk premium, implied variance, and realized variance.The variance risk premium proxy is moderately high during the 1990 and 2001 recessions, and much higher around the 1997-1998 Asia-Russia-LTCM crisis and the 2002-2003 corporate accounting scandals. The variance spikes during October 2008 already surpasses the initial shock of the Great Depression in October 1929. The huge run-up of variance risk premium in the fourth quarter of 2008 actually leads an equity market bottom in the first quarter of 2009.

Table 1 Panel A compares the summary statistics of different variance risk premium proxies based on alternative ways to estimate the conditional expectation of realized variances. The mean level of variance risk premium is around 17 to 22 (in percentage-squared, monthly basis) across five different estimates, with a standard deviation around 22-28. Not surprisingly, the variance risk premiums based on current and lagged realized variance have the highest kurtosis of 44 to 46, while the full sample AR(12) estimate has the lowest of 17. Also noteworthy is that the variance risk premium estimates based on raw current and lag realized variance has a skewness of -3, while others are all positive skewed. The negative skewness is entirely driven by one monthly observation of a negative spike in October 2008 (Figure 2 upper two panels) and by not using the expected variance in constructing the variance risk premium. Finally, the auto-regressive coefficient of order one is generally low between 0.26 and 0.76, with the full sample AR (12) achieves the lowest value. Figure 2 also shows the variance risk premia based on other estimates of the expected variance, where the recursive AR(12) and MA(12) approaches both suggest that the variance risk premium had achieved the unprecedented historical level around October-November 2008.

Basic summary statistics for the monthly returns and predictor variables are given in Table 1 Panels B to D. The mean excess return on the S&P500 over the sample period equals 3.58 percent annually, reflecting the significantly lowered market returns during the 2007-2008 financial crisis and economic downturn. The one month holding period returns for 2-6 month t-bills are ranging from 0.44 to 0.86 percent annually, and the credit spread for Moody's AAA rating is 1.25 percent and BAA 2.14 percent. The sample means for the variance risk premium is about 18.30 (in percentages squared). P/E ratio and short rate are very persistent with first order autocorrelations 0.97 and 0.99. While forward spread and variance risk premium are rather stationary, with a serial correlation between -0.10 to 0.39. Of particular interests, is that the variance risk premium variable generally has very small correlations with standard long-run predictor variables--0.07 with P/E ratio, around 0.04 to 0.06 with forward spreads, and -0.09 with short rate--which may partially explain why the short-run predictability of the variance premium variable is complementary to those of the established standard predictors.


3 Short-Run Predictability Puzzles of Financial Assets

This section presents new predictability evidence of the variance risk premium proxy for equity returns, bond return, and credit spreads, with and without the standard predictor variables--P/E ratio, forward spread, and short rate. Data are monthly observations with horizons up to one year. All of the reported  t-statistics are based on heteroskedasticity and serial correlation consistent standard errors (Newey and West, 1987). The discussion focuses on the estimated slope coefficients and their statistical significance as determined by the robust  t-statistics. The forecasts accuracy of the regressions are also measured by the corresponding adjusted  R^2's.8

3.1 Equity Returns

For equity returns, I focus on the regression of S&P500 returns on a long-run predictor--P/E ratio and a short-run predictor--variance risk premium,

\displaystyle xr_{t+h} ~ = ~ b_0 (h) ~ + ~ b_1 (h) \ {VRP_t} ~ + ~ b_2 (h) \ \log \left(P_t/E_t\right) ~ + ~ u_{t+h,t}, (4)

where  xr_{t+h} is the horizon-scaled market excess return and the horizon  h goes out to 12 months. The presentation will be brief here, as detailed results on equity return predictability in comparison with most established predictor variables are discussed in Bollerslev and Zhou (2007).

Table 2 top row shows that the predictability or  R^2 of the variance risk premium  VRP_t, starts out at -0.43 percent at monthly, peaks around 8 percent at four month, and then gradually decreases toward zero with longer forecasting horizons. The robust  t-statistic is the highest at four month at 3.56 and remains marginally significant at the ten month horizon. On the other hand, as shown in the middle row in Table 2, the usual long-run predictor,  \log P_t/E_t ratio, starts out barely significant at 10 percent level from one to three month, and then becomes insignificant. However, the adjusted  R^2 of P/E ratio monotonically increases from 0.92 percent at one month to 6.34 percent at twelve month.

Turning to the joint regressions reported in the bottom row of Table 2, it is clear that combining the variance premium with the P/E ratio results in an even higher  R^2 of 12.59 percent at four month horizon, which is higher than the sum of two  R^2's in the respective univariate regressions. The  t-statistics for  VRP_t and  \log P_t/E_t are also somewhat higher at four month, 3.61 and -1.60, than their univariate counterparts respectively.

Figure 3 visualizes such a short-run predictability pattern in  R^2 and  t-statistics. The predictability of variance premium has a tent shape pattern maximizes at four month horizon. While P/E ratio has no statistical significance at one-to-twelve month horizons.

3.2 Bond Returns

The failure of the Expectations Hypothesis (EH) of interest rates can be best characterized as that bond excess return, estimated from forward rates, is largely predictable and time-varying countercyclically (Fama, 1984,1986; Stambaugh, 1988).9Here I adopt the conventional forward rate setup as in Fama (1984) and augment it with the variance risk premium variable,

\displaystyle xhpr_{t+h}^n \ = \ b_0^n(h) \ + \ b_1^n(h) \ {VRP_t} \ + \ b_2^n(h)\ \left[f_{t-1}(n-h,h)-y_{t-1}(h)\right] \ + \ u_{t+h,t}^n, (5)

where  xhpr_{t+h}^n is the excess holding period return of zero coupon bonds with hold period  h=1,2,3,4,5 month and maturity  n=2,3,4,5,6 month (in excess of the yield on a  h-month zero coupon bond);  f_{t-1}(n-h,h) is the forward rate for a contract  h-month ahead with  n-h-month length; and  y_{t-1}(h) is the  h-month zero coupon bond yield.10

As shown in Table 3, the variance risk premium can significantly forecast the one month holding period excess returns of the two-six month t-bills, with a positive slope coefficient around 0.006 to 0.013. Considering the average level of variance risk premium of 18.30, this magnitude translates to an average bond risk premium induced by variance risk around 11 to 24 basis points. More importantly, the Newey-West  t-statistics are all well above the 1 percent significance level, with an  R^2 around 2.77 to 4.57 percent. Moving to the two month holding period, the  t-ratios reduce to a marginal significance of 1.22 to 2.05, and the  R^2 decreases to 0.86 to 3.54 percent. The predictability of variance premium converges to zero as the holding period increase to three-five months.

As Table 4 indicates, the forward rate (spread) is indeed a powerful predictor for excess bond returns for two-to-six month bonds with one-to-five month holding periods-- t-statistics all above 1 percent level and  R^2 between 2.06 and 30.88 percent. Another pattern is that the magnitudes of  t-statistics and  R^2 are generally higher at the one-month horizon and lower toward the five-month horizon. Overall, the predictability of forward spreads for short-term bills reported here are similar to those reported by Fama (1984).

More importantly, when variance risk premium is combined with forward rates, as shown in Table 5, the predictability of the variance risk premium remains intact. For example, at one month horizon, the  t-statistics are much higher for 3, 5, 6 month t-bills; and slightly lower for 2 and 4 month t-bills. Note that the adjusted  R^2's for the one month horizon with both variance premium and forward spread are all higher than the ones with forward spreads alone (Table 4 top row), suggesting that variance premium variable indeed contributes to the short-run bond return predictability, independent of that provided by the forward spread.

These results can be visualized in Figure 4, where for six month t-bill returns, the predictability of the variance premium variable is significant at one month but is monotonically decreasing with the holding period. While for the forward spread variable, although being significant at all one-to-five month horizons, its predictability seems to have a tent shape pattern that peaks at three month. This result suggests that the variance difference and forward spread are proxies for different components in bond risk premia.

3.3 Credit Spreads

The relatively large and time-varying credit spread on corporate bond has long been viewed as an anomaly in the literature (see, e.g., Huang and Huang, 2003). Here I provide some new evidence that, in additional to the standard predictor, namely the interest rate level (Longstaff and Schwartz, 1995), the variance risk premium proxy also helps to explain the short-run movement in credit spreads, with the following standard forecasting regression,

CS\displaystyle _{t+h} \ = \ b_0(h) \ + \ b_1(h) \ {VRP_t} \ + \ b_2(h) \ r_{f,t} \ + \ u_{t+h,t}, (6)

where the credit spread CS _{t+h} of  h month ahead is being forecasted by the short rate  r_{f,t} and the variance risk premium  VRP_t.

As shown in Table 6, short term interest rate is indeed a predominant predictor of the future credit spread levels, with  t-statistics of -3.94 for investment grade (Moody's AAA rating) and -2.94 for speculative grade (Moody's BAA rating). The adjusted  R^2 is around 32 percent, and the negative sign of the slope coefficient is consistent with the risk-neutral drift interpretation in Longstaff and Schwartz (1995).11 Although the significance of the short rate level extends to the six-month horizon, it is a very persistent variable with an AR(1) coefficient of 0.99 (Panel E in Table 1). As shown in the lower two panels of Figure 5, the predictability of short rate for credit spread shows a monotonic pattern of decreasing, even though the slope coefficient value seems to be increasing.

Note that if we include the variance risk premium alone in the forecasting regressions, its statistical significance is above the 1 percent level at one month horizon, with  t-statistic being 2.49 for AAA grade and 2.35 for BAA grade. Given the average level of variance risk premium of 18.30, that translates into an average effect on credit spread in the order of 9 to 15 basis points. Once the forecasting horizons increase to 3, 6, 9, and 12 months, the  t-statistic for the variance difference variable become insignificant or marginal. As shown in the upper two panels of Figure 5, the predictability of variance premium for credit spread shows a hump-shape pattern, which peaks at two month horizon and then becomes generally insignificant.

When the variance risk premium is combined with the short rate, both become more significant at the short horizons. For example, at one month horizon, the  t-statistics for short rate is -4.59 for AAA and -3.53 for BAA; while for variance premium it is 4.26 and 3.77 respectively. In fact, the variance risk premium variable maintains at least a marginal significance in the joint regressions even at the 12 month horizon, even though the short rate drops out as insignificant beyond the 6 month horizon. Judging from the  R^2 perspective, e.g., at the one month horizon, the univariate  R^2 for variance risk premium is about 5-7 percent, but its contribution to the joint  R^2 is about 9-10 percent, on top of what the short rate level has already achieved--32 percent.

This is an important finding, in that the variance risk premium or implied-expected variance difference captures an important component in credit risk premium that is independent with and complementary to the systematic risk risk factor that has already been captured by the short-term interest rate. Section 4 tries to provide some economic interpretation for such an effect.

3.4 Correlations between Different Markets

There is a valid concern that if the different financial markets exhibit a high degree of (positive) correlation among the residuals, then as a joint statistical inference issue, the conventional  t-statistics of the regression slope coefficient may have to be discounted, even if they are highly significant in the univariate regressions (Boudoukh, Richardson, and Whitelaw, 2008). This concern can be alleviated by the fact that the variance risk premium variable is not persistent at all, with an AR(1) coefficient being 0.26; and that neither the  t-statistics nor the  R^2's for variance premium are monotonically increasing with the forecasting horizons, although the P/E ratio does show such a pattern in its adjusted  R^2's.

One still need to address the remaining concern on whether the residuals of different markets are heavily correlated with each other. As shown in Panel A of Table 7, the raw return correlations between these three markets are generally low--ranging from -0.18 to 0.09--and are mostly close to zero. However, the concern about join statistical inference does bite for the securities in the same market--t-bills returns have correlations between 0.61 and 0.95 and two credit spreads have a correlation of 0.91. One should note that this paper focuses on the short-run predictability pattern of the variance premium variable across different markets, but not for the instruments in the same market.

Furthermore, once regressed on the variance risk premium variable, as shown in Panel B of Table 7, the residual correlations between different markets are either closer to zero or become slightly more negative in all cases, except one where the correlation between equity return and 2 month t-bill return increases from 0.04 to 0.08, albeit still a very small number indeed. Therefore, it seems that no strong positive comovement among these regressions residuals can be detected. Nonetheless, the high positive correlations among t-bill returns and between credit spreads remain largely unchanged.

Even if there is a high degree of positive correlation among the same instruments of a particular market, it does not necessarily mean that the univariate significant  t-statistics need to be discounted. An alternative interpretation could be that there is another common risk factor missing, in addition to the short-run variance premium variable; which should be incorporated in explaining the time-variations in these risk premiums.12


4 A Model of Macroeconomic Uncertainty

It is very challenging to provide a conceptual framework to jointly explain the risk premium dynamics in equity, bond, and credit markets. Here I draw from a self-contained general equilibrium model with stochastic consumption volatility-of-volatility (Bollerslev, Tauchen, and Zhou, 2009), and try to give a unified qualitative interpretation of these short-term risk premia puzzles. One should admit upfront that such a stylized model cannot provide satisfactory quantitative explanation for various aspects of asset pricing puzzles, within the same parameter setting and being constrained by matching the consumption dynamics.13

To be more specific, I will try to calibrate a model with macroeconomic uncertainty, by matching the equity risk premium as in Bansal and Yaron (2004) and short-run equity risk premium as in Bollerslev, Tauchen, and Zhou (2009), constrained by reasonable preference structure and consumption dynamics. And then I explore how far such a unified parameter setting can go for jointly explaining some salient features in variance risk premium and bond risk premium. One can learn by knowing why such a model is successful for explaining some risk premium dynamics across markets and why it cannot simultaneously explain the short-run predictability pattern in these markets.

4.1 Model Assumptions and Price-Dividend Ratio

The representative agent in the economy is equipped with Epstein-Zin-Weil recursive preferences, and has the value function  V_t of her life-time utility as

\displaystyle V_t ~ = ~ \left[\left(1-\delta\right)C_t^{\frac{1-\gamma}{\theta}} +\delta\left(E_t\left[V_{t+1}^{1-\gamma}\right]\right)^{\frac{1}{\theta}}\right]^{\frac{\theta}{1-\gamma}}, (7)

where  C_t is consumption at time  t,  \delta denotes the subjective discount factor,  \gamma refers to the coefficient of risk aversion,  \theta = \frac{1-\gamma}{1 - \frac{1}{\psi}}, and  \psi equals the intertemporal elasticity of substitution (IES). The agent maximizes her utility subject to the budget constraint,  W_{t+1}=(W_t-C_t) \times r_{t+1}, where  W_t is the wealth of the agent and  r_{t} is the return on the consumption asset. The key assumptions are that  \gamma>1, implying that the agents are more risk averse than the log utility investor; and  \psi>1 hence  \theta <0, implying that agents prefer an earlier resolution of economic uncertainty. These restrictions ensure that the uncertainty or volatility risk in asset markets carries a positive risk premia.

Suppose that log consumption growth and its volatility follow the joint dynamics

\displaystyle g_{t+1} \displaystyle = \displaystyle \mu_{g} + \sigma_{g,t} z_{g,t+1} , (8)
\displaystyle \sigma^2_{g,t+1} \displaystyle = \displaystyle a_{\sigma} + \rho_{\sigma} \sigma^2_{g,t} + \sqrt{q_t} z_{\sigma, t+1} , (9)
\displaystyle q_{t+1} \displaystyle = \displaystyle a_q + \rho_q q_t + \varphi_q \sqrt{q_t}z_{q,t+1} , (10)

where  \mu_{g}>0 denotes the constant mean growth rate,  \sigma^2_{g,t+1} represents time-varying volatility in consumption growth, and  q_t introduces the volatility uncertainty process in the consumption growth process. The parameters satisfy  a_{\sigma} > 0, a_q > 0,  \vert\rho_{\sigma}\vert < 1,  \vert\rho_{q}\vert < 1,  \varphi_q>0; and  \{z_{g,t}\},  \{z_{\sigma,t}\} and  \{z_{q,t}\} are iid Normal (0,1) processes jointly independent with each other. The time-variation in  \sigma^2_{g,t+1} is one of the two components that drives the equity risk premium, or the "consumption risk"; while the time-variation in  q_t is not only responsible for the "uncertainty risk" component in equity risk premium, but also constitutes the main driver of bond and variance risk premia as explained bellow.

Let  w_t denote the logarithm of the price-dividend or wealth-consumption ratio, of the asset that pays the consumption endowment,  \{C_{t+i}\}^\infty_{i=1}; and conjecture a solution for  w_t as an affine function of the state variables,  \sigma^2_{g,t} and  q_t,

\displaystyle w_t = A_0 + A_\sigma \sigma^2_{g,t} + A_q q_t. (11)

One can solve for the coefficients  A_0,\, A_\sigma and  A_q (as in Bollerslev, Tauchen, and Zhou, 2009), using the standard Campbell and Shiller (1988) approximation  r_{t+1} = \kappa_0 + \kappa_1 w_{t+1} - w_{t} + g_{t+1}. The aforementioned restrictions that  \gamma>1 and  \psi>1, hence  \theta <0, readily imply that the impact coefficient associated with both consumption and volatility state variables are negative; i.e.,  A_{\sigma}<0 and  A_q<0. So if consumption and uncertainty risks are high, the price-dividend ratio is low, hence risk premia are high. In response to high economic uncertainty risks, agents sell risky assets, and consequently the wealth-consumption ratio falls; so that risk premiums rise.

4.2 Model-Implied Equity, Variance, Bond Risk Premia

Given the solution of the price-dividend ratio, one can easily solve for the variables of interest, like equity return and risk-free rate (see, Bollerslev, Tauchen, and Zhou, 2009, for details), the model-implied equity risk premium can be shown as

\displaystyle \operatorname{E}_t\left(r_{t+1}\right)-r_{f,t} ~ =~ \gamma \sigma_{g,t}^2 \, + \, (1-\theta)\kappa_1^2(A_q^2 \varphi_q^2 + A_{\sigma}^2)q_t \, > \, 0. (12)

The premium is composed of two separate terms. The first term,  \gamma \sigma_{g,t}^2, is compensating for the classic consumption risk term as in a standard consumption based CAPM model. The second term,  (1-\theta)\kappa_1^2(A_q^2 \varphi_q^2 + A_{\sigma}^2)q_t, represents a true premium for variance risk. The existence of the variance or uncertainty risk premium depends crucially on the dual assumptions of recursive utility, or  \theta \neq 1, as uncertainty would otherwise not be a priced factor; and time varying volatility-of-volatility, in the form of the  q_t process. The restrictions that  \gamma>1 and  \psi>1 implies that the variance risk premium embedded in the equity risk premium is always positive by construction. And since the variance risk premium embedded in equity returns loads on the same uncertainty risk factor,  q_t, as in the variance risk premium (demonstrated bellow), the latter becomes a perfect predictor for equity premium variation that is induced by the stochastic economic uncertainty.

The conditional variance of the time  t to  t+1 return,  \sigma^2_{r,t} \equiv   Var _t( r_{t+1} ), can be shown as  \sigma^2_{r,t} = \sigma_{g,t}^2 + \kappa_1^2 \left( A_{\sigma}^2 + A_q^2 \varphi_q^2 \right) q_t. The variance risk premium can be defined as the difference between risk-neutral and objective expectations of the return variance,

\displaystyle VRP_{t} \displaystyle \equiv \displaystyle \operatorname{E}_t^Q\left(\sigma_{r,t+1}^2\right)-\operatorname{E}_t^P\left(\sigma_{r,t+1}^2\right)  
  \displaystyle \approx \displaystyle (\theta-1)\kappa_1\left[ A_{\sigma} + A_q \kappa_1^2 \left( A_{\sigma}^2 + A_q^2 \varphi_q^2 \right) \varphi_q^2\right]q_t ~ >0~, (13)

where the approximation comes from the fact that the model-implied risk-neutral conditional expectation cannot be computed in closed form, and a log-linear approximation is applied. One key observation here is that any temporal variation in the endogenously generated variance risk premium, is due solely to the volatility-of-volatility or economic uncertainty risk,  q_t, but not the consumption growth risk,  \sigma^2_{g,t+1}. Moreover, provided that  \theta <0,  A_{\sigma}<0, and  A_q<0, as would be implied by the agents' preference of an earlier resolution of economic uncertainty (intertemporal elasticity of substitution--IES--bigger than one), this difference between the risk-neutral and objective expectations of return variances is guaranteed to be positive. If  \varphi_q=0, and therefore  q_t = q is constant--without stochastic volatility-of-volatility, the variance premium reduces to a constant  (\theta-1)\kappa_1A_{\sigma}q, and one cannot replicate the large skewness and kurtosis in the observed variance risk premium series.

The bond yield in this economy can be shown as an affine function of the state variables,  y_t(n)= - \frac{1}{n} \left[\begin{array}{ccc} A(n) & B(n) & C(n) \end{array} \right] \left[\begin{array}{ccc} 1 & \sigma^2_{g,t} & q_t \end{array} \right]', where the coefficients  A(n),  B(n), and  C(n) are given in Zhou (2008). Let  rx_{t+1}(n-1) be the bond excess return from  t to  t+1 for an  n-period bond holding for one period, then its expected value  rp_t^n or bond risk premium is given by,

\displaystyle rp_t^n = \left[B(n-1)(\theta-1)\kappa_1 A_{\sigma} +C(n-1)(\theta-1)\kappa_1 A_q \varphi_q^2 \right]q_t \, > \, 0, (14)

where the risk premia genuinely has two time-varying components--consumption risk and uncertainty risk, but they are co-linear in the same state variable  q_t. This is driven by the fact that the variances of both the volatility process and the volatility-of-volatility process are loading on the same state variable,  q_t. Recall that in the standard consumption-based asset pricing model, bond risk premium is zero by construction, because there is no time-varying volatility. The more interesting observation here is that, even if the consumption volatility is time-varying, as long as the volatility-of-volatility is constant, then the bond risk premia must be constant. Therefore the current modeling framework can qualitatively explains the existence of real bond risk premium variations, without introducing exogenous inflation or monetary distortion.

4.3 Calibrating Equity, Variance, and Bond Risk Premia

To more directly gauge how the model adopted here can explain the documented risk premium dynamics, I perform a limited calibration exercise. The basic strategy is to find the preference and distribution parameters that are constrained by reasonable consumption dynamics and can simultaneously match the observed equity, bond, and variance premia as much as possible.

4.3.1 Calibration Design and Parameter Setting

As shown in Table 8, the benchmark parameter settings are adapted from Bansal and Yaron (2004) such that the consumption growth rate (  \mu_g = 0.0015) is 2.4 percent annually and the consumption volatility is highly persistent (  \rho_\sigma = 0.978 at monthly), same across the two calibration settings in both this paper and Bollerslev, Tauchen, and Zhou (2009). I also use the same intertemporal elasticity of substitution (IES) parameter (  \psi = 1.5) such that agents have the same preference for an earlier resolution of economic uncertainty. The time preference is also the same as  \delta = 0.997 across two scenarios. The Campbell-Shiller approximation constants are chosen as  \kappa_1 = 0.9 hence  \kappa_0 = 0.3251, similar as those in Bansal and Yaron (2004).

However, these two papers differ dramatically in terms of both risk aversion and volatility risk. In Bollerslev, Tauchen, and Zhou (2009), the risk aversion coefficient is  \gamma= 10, same as in Bansal and Yaron (2004); but in this paper, I choose  \gamma= 2. On the other hand, in Bollerslev, Tauchen, and Zhou (2009), the consumption volatility is  \operatorname{E}( \sigma_g ) = 0.0078^2, similar as that in Bansal and Yaron (2004), while here it increases to  E( \sigma_g ) = (0.0078\times4)^2. A material implication of these parameter values is that the consumption volatility is increased to 4 times of the 2.7 percent annually, which may be justifiable if one "leverages" up the dividend shocks several times larger than the consumption shocks (Abel, 1999).14

Furthermore, Bollerslev, Tauchen, and Zhou (2009) choose the stochastic economic uncertainty process,  q_t, to have monthly persistence level of  \rho_q = 0.80, long-run mean of  \operatorname{E}(q) = 1 \times 10.0^{-6}, and volatility-of-volatility parameter as  \varphi_q = 0.001; but here I choose  \rho_q = 0.95,  \operatorname{E}(q) = 3.65 \times 10.0^{-4}, and  \varphi_q = 0.008, implying that the  q_t process is not only magnified but also more persistent and volatile. As demonstrated bellow, such a modification is critical in balancing the needs of jointly matching the equity, variance, and bond risk premia in the short-run.

4.3.2 Calibrated Equity, Variance, and Bond Risk Premia

The calibration results in Table 9 indicate that one can achieve a reasonable compromise to simultaneously match the equity, variance, bond risk premium puzzles, without relying on introducing jumps into the endowment and/or its volatility process (as in Eraker and Shaliastovich, 2008; Drechsler and Yaron, 2008) or introducing exogenous inflation process for the nominal bonds (Bansal and Shaliastovich, 2008b; Gallmeyer, Hollifield, Palomino, and Zin, 2008; Wachter, 2006). It can be achieved if one would be willing to "levered" up the volatility-of-volatility process and to lower risk aversion to a moderate level. Admittedly such a "short-run" risk model lacks the long-run component to match the consumption growth predictability, which is subject to some debate due to the consumption data measurement problem.

The resulting equity risk premium is 5.11 percent and real interest rate 1.86 percent, which are different than Bollerslev, Tauchen, and Zhou (2009) (7.79 and 0.69 percent) but similar to the observed premiums of 3.58 and 1.13 percent for the sample period 1990-2008 (as opposed to 7.84 and 0.87 percent for 1930-2008). The parameter choice here has a advantage of getting a better equity volatility of 13.24 percent, as opposed to 4 percent in Bollerslev, Tauchen, and Zhou (2009), comparing to 14.60 percent observed in the recent period. On the other hand, this paper overfits the volatility of risk-free rate at 8.83 percent, versus 3.37 percent in recent data and 2.95 percent in Bollerslev, Tauchen, and Zhou (2009). This sacrifice in matching the equity risk premium has gained in explaining the variance and bond risk premia. However, the underfitting of equity risk premium (5.11 percent) and the overfitting of short rate volatility (7.83 percent) disappear if one considers the earlier sample period of 1891-1949 (4.97 and 7.83 percent respectively).

The model-implied variance risk premium has a mean of 18.30 and a standard deviation of 25.12, which are very close to observed values of 18.30 and 22.69. More importantly, the model produces quite reasonable values in skewness (2.48) and kurtosis (13.18), which are slightly less than the observed values of 2.79 and 16.62. On the other hand, the parameter specification in Bollerslev, Tauchen, and Zhou (2009) would completely miss the level and standard deviation of variance premium and underfit the the skewness and kurtosis. These results are non-trivial in that a stochastic volatility-of-volatility model can generate realistic skewness and kurtosis in variance risk premium, similar to the result provided in Drechsler and Yaron (2008), where the common jumps have been introduced into the consumption growth and its volatility processes. Finally, the model-implied persistence level--equal to  \rho_q--is as high as 0.95, as opposed to the low value of 0.26 reflected in the data and 0.80 in Bollerslev, Tauchen, and Zhou (2009). It turns out that such a high level of persistence is necessary to generate high enough variance risk premium, without sacrificing too much in matching the equity risk premium and risk-free rate, but with overfitting of the interest rate volatility and bond risk premia.

The observed bond risk premium for 2-6 month t-bills holding one month are positive but only mildly upward sloping around 0.44 to 0.86 percent. If one chooses the specification in Bollerslev, Tauchen, and Zhou (2009), the bond risk premia would be extremely hump-shaped with about 3 percent at 2 month, 7 percent at 4 month, and  -1.06\times10^6 percent at 6 month.15 In contrast, the specification chosen in this paper would produce a term structure of bond risk premia from 2.92 percent to 13.90 percent for 2-6 month t-bills holding one month, which is certainty more steep than the observed data, but is quite closer to the observed term structure or bond risk premiums than Bollerslev, Tauchen, and Zhou (2009). Again, this is a non-trivial result because one has not introduced the inflation dynamics or monetary distortion yet into the model.

4.3.3 Challenge for Explaining the Term Structure of Predictability

The current modeling framework has been shown (Bollerslev, Tauchen, and Zhou, 2009, Figure 1) to be able to replicate the short-run predictability pattern in equity returns (top two panels in Figure 3). However, such a model implies that the bond risk premium is entirely driven by the uncertainty risk factor,  q_t, (equation 14), which is also the only state variable in the variance risk premium (equation 13). Therefore, as shown in Figure 6, the prediction  R^2's should be always equal to one, as opposed to the low single digit; and the (deterministic) slope coefficients, although upward sloping, are many times larger than their empirical counterparts. Therefore it remains a challenge to simultaneously reproduce the short-run predictability patterns of the variance risk premium for both equity and bond returns, and for credit spreads the current model still lacks a defaultable sector.

4.4 Extension to Credit Spread Puzzles

It is also challenging to incorporate the default risk of a representative firm into the current modeling framework. The strategy could follow Chen (2008) and Bhamra, Kuehn, and Strebulaev (2009), where the recursive preference plus macroeconomic uncertainty generate richer dynamics in the credit spread dynamics. To fix the idea, assume that in a Merton (1974) type model as in Chen, Collin-Dufresne, and Goldstein (2008), the credit spread of a discount bond for a defaultable firm with  T maturity can be shown as

CS\displaystyle _t(T) = - \frac{1}{T} \log\left\{1-\mbox{LGD} \times \mbox{Normal} \left[\mbox{Normal}^{-1} \left(\mbox{PD}\right)+\lambda\sigma \sqrt{T}\right] \right\}, (15)

where LGD is the loss given default, PD is the real default probability,  \lambda is the market price of asset risk, and  \sigma is the asset return volatility. All these important variables are constants or deterministic in the original Merton model. It is well known that such a simplified model cannot explain the high credit spread level and its time variation (Huang and Huang, 2003).

The equilibrium structural approach (as in, e.g., Chen, 2008; Bhamra, Kuehn, and Strebulaev, 2009, among others) can be viewed as letting the real default probability PD _t to be time-varying and countercyclical, with possible business cycle fluctuations of the firm's refinancing decision or default barrier. It is also possible to model the recovery rate  1-LGD _t as stochastic to help explain the cyclical behavior in credit spreads, but the quantitative improvement could be marginal in this angle. Chen, Collin-Dufresne, and Goldstein (2008) take a novel approach to allow for the market price of risk  \lambda_t to be driven by a countercyclical risk aversion motivated by the Campbell and Cochrane (1999) external habit model. Finally, one can allow the asset volatility  \sigma_t to be time-varying and countercyclical. Such an extension may be consistent with the empirical and calibration evidence in Zhang, Zhou, and Zhu (2009), where the stochastic asset volatility can help structural models to explain the credit spread puzzles.16


5 Conclusion

The implied-expected variance difference can be viewed as a measure for the variance risk premium. This paper provides consistent empirical evidence that the variance risk premium can significantly predict short-run equity returns, bond returns, and credit spreads. The documented return predictability peaks around one-to-four month and decline with the forecasting horizon. Importantly, such a short-term forecastability of risk premia is complementary to the established predictors--P/E ratio, forward spread, and short rate. This constitutes an important evidence that risk premia across major financial markets co-vary in short-term, and such a comovement seems to be driven by a common risk factor, measured by the implied-expected variance difference.

Such a common short-run risk factor may be a proxy for the macroeconomic uncertainty or consumption volatility risk that varies independently with the consumption growth risk--the latter being the main focus of the long-run risks model (Bansal and Yaron, 2004). The empirical results may be consistent with a general equilibrium model incorporating the effects of such a time-varying economic uncertainty component, where the uncertainty risk is priced only under the recursive preference. This paper provides calibration evidence that the equity premium puzzle, variance premium dynamics, and bond risk premia in short-horizon may be qualitatively explained by the proposed model with the same calibration parameter setting. Extension to credit spread puzzle requires a defaultable sector in the modeling framework.

Although the stylized model examined here can provide qualitative justification for the short-run predictability of asset market returns from the variance risk premium, it is not rich enough to simultaneously explain such effects within the same parameter setting. More importantly, to jointly interpret the long-run and short-run comovements in asset markets, a consumption growth factor may be needed to quantitatively replicate various predictability puzzles established in the literature. Finally, the short-run forecastability of variance risk premium documented here as in the time-series domain need to be reconciled with the cross-sectional evidence of asset market returns. I leave these challenging issues for future research.


Bibliography

Abel, Andrew B. (1999),
"Risk Premia and Term Premia in General Equilibrium," Journal of Monetary Economics, vol. 43, 3-33.
Andersen, Torben G., Tim Bollerslev, Francis X. Diebold, and Heiko Ebens (2001),
"The Distribution of Realized Stock Return Volatility," Journal of Financial Economics, vol. 61, 43-76.
Bakshi, Gurdip and Nikunj Kapadia (2003),
"Delta-Hedged Gains and the Negative Market Volatility Risk Premium," Review of Financial Studies, vol. 16, 527-566.
Bakshi, Gurdip and Dilip Madan (2000),
"Spanning and Derivative-Security Valuation," Journal of Financial Economics, vol. 55, 205-238.
Bakshi, Gurdip and Dilip Madan (2006),
"A Theory of Volatility Spread," Management Science, vol. 52, 1945-1956.
Bansal, Ravi, Dana Kiku, and Amir Yaron (2007),
"A Note on the Economics and Statistics of Predictability: A Long Run Risks Perspective," Working Paper, Duke University and University of Pennsylvania.
Bansal, Ravi, Dana Kiku, and Amir Yaron (2009),
"An Empirical Evaluation of the Long-Run Risks Model for Asset Prices," Duke University and University of Pennsylvania.
Bansal, Ravi and Ivan Shaliastovich (2008a),
"Learning and Asset-Price Jumps," Working Paper, Duke University.
Bansal, Ravi and Ivan Shaliastovich (2008b),
"A Long-Run Risks Explanation of Preditability Puzzles in Bond and Currency Markets," Working Paper, Duke University.
Bansal, Ravi and Amir Yaron (2004),
"Risks for the Long Run: A Potential Resolution of Asset Pricing Puzzles," Journal of Finance, vol. 59, 1481-1509.
Barndorff-Nielsen, Ole and Neil Shephard (2002),
"Econometric Analysis of Realised Volatility and Its Use in Estimating Stochastic Volatility Models," Journal of Royal Statistical Society, Series B, vol. 64, 253-280.
Barro, Robert J. (2006),
"Rare Disasters and Asset Markets in the Twentieth Century," Quarterly Journal of Economics, vol. 121, 823-866.
Bates, David S. (1996),
"Jumps and Stochastic Volatility: Exchange Rate Process Implicit in Deutsche Mark Options," The Review of Financial Studies, vol. 9, 69-107.
Beeler, Jason and John Y. Campbell (2009),
"The Long-Run Risks Model and Aggregate Asset Prices: An Empirical Assessment," Working Paper, Harvard University.
Bekaert, Geert, Marie Hoerova, and Martin Scheicher (2009),
"What Do Asset Prices Have to Say about Risk Appetite and Uncertainty?" Working Paper, European Central Bank.
Benzoni, Luca, Collin-Dufresne, and Robert S. Goldstein (2008),
"Explaining Pre- and Post-1987 Crash Asset Prices within a Unified General Equilibrium Framework," Working Paper, Federal Reserve Bank of Chicago.
Bhamra, Harjoat S., Lars-Alexander Kuehn, and Ilya A. Strebulaev (2009),
"The Levered Equity Risk Premium and Credit Spreads: A Unified Framework," Review of Financial Studies, forthcoming.
Bloom, Nicholas (2009),
"The Impact of Uncertainty Shocks," Econometrics, vol. 77, 623-685.
Bloom, Nicholas, Max Floetotto, and Nir Jaimovich (2009),
"Really Uncertain Business Cycles," Working Paper, Stanford University.
Bollerslev, Tim, Mike Gibson, and Hao Zhou (2008a),
"Dynamic Estimation of Volatility Risk Premia and Investor Risk Aversion from Option-Implied and Realized Volatilities," Journal of Econometrics, forthcoming.
Bollerslev, Tim, Natalia Sizova, and George Tauchen (2008b),
"Volatility in Equilibrium: Asymmetries and Dynamics Dependencies," Working Paper, Department of Economics, Duke University.
Bollerslev, Tim, George Tauchen, and Hao Zhou (2009),
"Expected Stock Returns and Variance Risk Premia," Review of Financial Studies, forthcoming.
Bollerslev, Tim and Hao Zhou (2007),
"Expected Stock Returns and Variance Risk Premia," Finance and Economics Disucussion Series 2007-11, Federal Reserve Board.
Boudoukh, Jacob, Matthew Richardson, and Robert F. Whitelaw (2008),
"The Myth of Long-Horizon Predictability," Review of Financial Studies, vol. 21, 1577-1605.
Breeden, Douglas and Robert Litzenberger (1978),
"Prices of State-Contingent Claims Implicit in Option Prices," Journal of Business, vol. 51, 621-651.
Britten-Jones, Mark and Anthony Neuberger (2000),
"Option Prices, Implied Price Processes, and Stochastic Volatility," Journal of Finance, vol. 55, 839-866.
Buraschi, Andrea, Fabio Trojani, and Andrea Vedolin (2009),
"When Uncertainty Blows in the Orchard: Comovement and Equilibrium Volatility Risk Premia," Working Paper, Imperial College London and University of Lugano.
Campbell, John Y. and John H. Cochrane (1999),
"By Force of Habit: A Consumption Based Explanation of Aggregate Stock Market Behavior," Journal of Political Economy, vol. 107, 205-251.
Campbell, John Y. and Robert J. Shiller (1988),
"The Dividend-Price Ratio and Expectations of Future Dividends and Discount Factors," Review of Financial Studies, vol. 1, 195-228.
Campbell, John Y. and Motohiro Yogo (2006),
"Efficient Tests of Stock Return Predictability," Journal of Financial Economics, vol. 81, 27-60.
Carr, Peter and Dilip Madan (1998),
"Towards a Theory of Volatility Trading," in Volatility: New Estimation Techniques for Pricing Derivatives, chap.29, 417-427, Robert Jarrow (ed.). London: Risk Books.
Carr, Peter and Liuren Wu (2008),
"Variance Risk Premiums," Review of Financial Studies, forthcoming.
Chen, Hui (2008),
"Macroeconomic Conditions and the Puzzles of Credit Spreads and Capital Structure," Journal of Finance, forthcoming.
Chen, Hui and Michal Pakos (2008),
"Asset Pricing with Uncertainty for the Long Run," Working Paper, MIT.
Chen, Long, Pierre Collin-Dufresne, and Robert S. Goldstein (2008),
"On the Relation Between the Credit Spread Puzzle and the Equity Premium Puzzle," Review of Financial Studies, forthcoming.
Cochrane, John H. and Monika Piazzesi (2005),
"Bond Risk Premia," American Economic Review, vol. 95, 138-160.
Cremers, Martijn, Joost Driessen, and Pascal Maenhout (2007),
"Explaining the Level of Credit Spreads: Option-Implied Jump Risk Premia in a Firm Value Model," Review of Financial Studies, forthcoming.
David, Alexander (2008),
"Inflation Uncertainty, Asset Valuations, and the Credit Spreads Puzzle," Reveiew of Financial Studies, vol. 21, 2487-2534.
Drechsler, Itamar (2008),
"Uncertainty, Time-Varying Fear, and Asset Prices," Working Paper, University of Pennsylvania.
Drechsler, Itamar and Amir Yaron (2008),
"What's Vol Got to Do With It," Working Paper, University of Pennsylvania.
Epstein, Larry G. and Stanley E. Zin (1991),
"Substitution, Risk Aversion, and the Temporal Behavior of Consumption and Asset Returns: An Empirical Analysis," Journal of Political Economy, vol. 99, 263-286.
Eraker, Bjørn and Ivan Shaliastovich (2008),
"An Equilibrium Guide to Designing Affine Pricing Models," Mathematical Finance, vol. 18, 519-543.
Fama, Eugene F. (1984),
"The Information in the Term Structure," Journal of Financial Economics, vol. 13, 509-528.
Fama, Eugene F. (1986),
"Term Premiums and Default Premiums in Money Markets," Journal of Financial Economics, vol. 17, 175-196.
Fama, Eugene F. and Robert T. Bliss (1987),
"The Information in Long-Maturity Forward Rates," American Economic Review, vol. 77, 680-692.
Gabaix, Xavier (2009),
"Variable Rare Disasters: An Exactly Solved Framework for Ten Puzzles in Macro-Finance," Working Paper, NYU Stern.
Gallmeyer, Michael F., Burton Hollifield, Francisco Palomino, and Stanley E. Zin (2008),
"Term Premium Dynamics and the Taylor Rule," The David A. Tepper School of Business, Carnegie Mellon University.
Goyal, Amit and Ivo Welch (2008),
"A Comprehensive Look at the Empirical Performance of Equity Premium Prediction," Review of Financial Studies, forthcoming.
Gvozdeva, Evgenia and Praveen Kumar (2009),
"Structural Uncertainty, Learning, and Asset Pricing," Working Paper, University of Houston.
Heston, Steven (1993),
"A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options," Review of Financial Studies, vol. 6, 327-343.
Huang, Jing-Zhi and Ming Huang (2003),
"How Much of the Corporate-Treasury Yield Spread Is Due to Credit Risk?" Working Paper, Penn State University.
Jiang, George and Yisong Tian (2005),
"Model-Free Implied Volatility and Its Information Content," Review of Financial Studies, vol. 18, 1305-1342.
Lettau, Martin, Sydney Ludvigson, and Jessica Wachter (2008),
"The Declining Equity Premium: What Role does Macroeconomic Risk PLay?" Review of Financial Studies, forthcoming.
Longstaff, Francis A. and Eduardo S. Schwartz (1995),
"A Simple Approach to Valuing Risky Fixed and Floating Rate Debt," Journal of Finance, vol. 50, 789-820.
Ludvigson, Sydney C. and Serena Ng (2008),
"Macro Factors in Bond Risk Premia," Review of Financial Studies, forthcoming.
Mehra, Rajnish and Edward C. Prescott (1985),
"The Equity Premium: A Puzzle," Journal of Monetary Economics, vol. 15, 145-161.
Merton, Robert C. (1974),
"On the Pricing of Corporate Debt: The Risk Structure of Interest Rates," Journal of Finance, vol. 29, 449-470.
Newey, Whitney K. and Kenneth D. West (1987),
"A Simple Positive Semi-Definite, Heteroskedasticity and Autocorrelation Consistent Covariance Matrix," Econometrica, vol. 55, 703-708.
Pástor, Lubos and Robert F. Stambaugh (2009),
"Are Stocks Really Less Volatile in the Long Run?" The University of Chicago Booth School of Business.
Rosenberg, Joshua V. and Robert F. Engle (2002),
"Empirical Pricing Kernels," Journal of Financial Economics, vol. 64, 341-372.
Shaliastovich, Ivan (2008),
"Learning, Confidence, and Option Price," Working Paper, Duke University.
Stambaugh, Robert F. (1988),
"The Information in Forward Rates: Implications for Models of the Term Structure," Journal of Financial Economics, vol. 21, 41-70.
Stambaugh, Robert F. (1999),
"Predictive Regressions," Journal of Financial Economics, vol. 54, 375-421.
Todorov, Viktor (2009),
"Variance Risk Premium Dynamics: The Role of Jumps," Review of Financial Studies, forthcoming.
Wachter, Jessica A. (2006),
"A Consumption-Based Model of the Term Structure of Interest Rates," Journal of Financial Economics, vol. 79, 365-399.
Weil, Philippe (1989),
"The Equity Premium Puzzle and the Risk Free Rate Puzzle," Journal of Monetary Economics, vol. 24, 401-421.
Whaley, Robert E. (2000),
"The Investor Fear Gauge," Journal of Portfolio Management, vol. 26, 12-17.
Wright, Jonathan and Hao Zhou (2009),
"Bond Risk Premia and Realized Jump Risk," Journal of Banking and Finance, forthcoming.
Zhang, Benjamin Yibin, Hao Zhou, and Haibin Zhu (2009),
"Explaining Credit Default Swap Spreads with the Equity Volatility and Jump Risks of Individual Firms," Review of Financial Studies, forthcoming.
Zhou, Guofu and Yingzi Zhu (2009),
"A Long-Run Risks Model with Long- and Short-Run Volatilities: Explaining Predictability and Volatility Risk Premium," Working Paper, Olin School of Business, Washington University.
Zhou, Hao (2008),
"Stochastic Economic Uncertainty and Asset Pricing Puzzles--Analytical Solutions," Unpublished Techinical Note, Federal Reserve Board.


Table 1: Summary statistics

The sample period extends from January 1990 to December 2008. All variables are reported in annualized percentage form whenever appropriate. Panel A reports the variance risk premiums using different methods for estimating the physical expectation of the realized variance: (1) ex post realized variance, (2) lagged realized variance, (3) 12 month moving average, (4) recursive AR(12) forecast, and (5) full sample AR(12) forecast. Panel B reports S&P500 market returns, log price-earning ratio. Panel C reports the 1 month excess holding period return and forward minus spot rate spreads for Treasury zero coupon bond with 1-6 month maturities. Panel D reports the Moody's AAA and BAA credit spread indices with the US short rate level.


Table 1, Panel A: Comparison of  VRP_t with different  RV_t forecasts
 RV Lag  RV  MA(12) Recursive  AR(12) Full Sample  AR(12)
Summary Statistics: Mean 17.22 17.07 20.49 21.81 18.30
Summary Statistics: Std Dev 20.27 19.99 25.03 27.86 22.69
Summary Statistics: Skewness -3.12 -3.32 4.19 4.39 2.79
Summary Statistics: Kurtosis 44.20 46.41 27.51 34.15 16.62
Summary Statistics: AR(1) 0.28 0.76 0.66 0.76 0.26
Correlation Matrix:  RV 1.00 0.28 0.11 0.23 0.04
Correlation Matrix: Lag  RV   1.00 0.06 0.12 0.18
Correlation Matrix:  MA(12)     1.00 0.63 0.75
Correlation Matrix: Recursive  AR(12)       1.00 0.21
Correlation Matrix: Full Sample  AR(12)         1.00


Table 1, Panel B: Equity
 r_{t}-r_{f,t}  VRP  \log(P_t/E_t)
Summary Statistics: Mean 3.58 18.30 3.12
Summary Statistics: Std Dev 50.58 22.69 0.25
Summary Statistics: Skewness -0.90 2.79 0.48
Summary Statistics: Kurtosis 50.08 16.62 2.53
Summary Statistics: AR(1) 0.07 0.26 0.97
Correlation Matrix:  r_{t}-r_{f,t} 1.00 0.01 -0.11
Correlation Matrix:  VRP   1.00 0.07
Correlation Matrix:  \log(P_t/E_t)     1.00


Table 1, Panel C: Treasury Bill
               f_t(1,1)  f_t(2,1)  f_t(3,1)  f_t(4,1)  f_t(5,1)
   xhpr_{t+1}^2  xhpr_{t+1}^3  xhpr_{t+1}^4  xhpr_{t+1}^5  xhpr_{t+1}^6  VRP  -y_t(1)  -y_t(1)  -y_t(1)  -y_t(1)  -y_t(1)
Summary Statistics: Mean 0.44 0.57 0.56 0.78 0.86 18.30 0.41 0.46 0.46 0.50 0.54
Summary Statistics: Std Dev 0.63 0.70 0.86 1.03 1.33 22.69 0.96 0.76 0.69 0.69 0.69
Summary Statistics: Skewness 0.91 0.74 1.16 1.18 0.65 2.79 1.15 0.51 0.97 0.88 0.74
Summary Statistics: Kurtosis 4.85 4.59 6.05 5.88 8.87 16.62 5.31 4.05 6.26 5.29 5.19
Summary Statistics: AR(1) 0.47 0.36 0.29 0.32 0.30 0.26 -0.10 0.26 0.35 0.37 0.39
Correlation Matrix:  xhpr_{t+1}^2 1.00 0.91 0.80 0.73 0.61 0.22 0.71 0.76 0.75 0.74 0.71
Correlation Matrix:  xhpr_{t+1}^3   1.00 0.92 0.88 0.77 0.18 0.57 0.66 0.64 0.62 0.60
Correlation Matrix:  xhpr_{t+1}^4     1.00 0.95 0.86 0.18 0.46 0.52 0.53 0.51 0.50
Correlation Matrix:  xhpr_{t+1}^5       1.00 0.92 0.21 0.37 0.44 0.43 0.44 0.43
Correlation Matrix:  xhpr_{t+1}^6         1.00 0.22 0.27 0.33 0.33 0.33 0.38
Correlation Matrix:  VRP           1.00 0.06 0.04 0.05 0.05 0.06
Correlation Matrix:  f_t(1,1)-y_t(1)             1.00 0.88 0.88 0.88 0.84
Correlation Matrix:  f_t(2,1)-y_t(1)               1.00 0.94 0.94 0.93
Correlation Matrix:  f_t(3,1)-y_t(1)                 1.00 0.97 0.95
Correlation Matrix:  f_t(4,1)-y_t(1)                   1.00 0.97
Correlation Matrix:  f_t(5,1)-y_t(1)                     1.00


Table 1, Panel D: Credit Spread
AAA BAA  VRP Short Rate
Summary Statistics: Mean 1.25 2.14 18.30 4.24
Summary Statistics: Std Dev 0.45 0.68 22.69 1.83
Summary Statistics: Skewness 0.89 2.15 2.79 -0.20
Summary Statistics: Kurtosis 3.18 10.50 16.62 2.37
Summary Statistics: AR(1) 0.98 1.05 0.26 0.99
Correlation Matrix: AAA 1.00 0.91 0.16 -0.61
Correlation Matrix: BAA   1.00 0.17 -0.62
Correlation Matrix:  VRP     1.00 0.10
Correlation Matrix: Short Rate       1.00



Table 2: Equity Returns, Variance Risk Premia, and P/E ratios

The sample period extends from January 1990 to December 2008. All of the regressions are based on monthly observations. Robust  t-statistics following Newey and West (1987) with 24 lags are reported in parentheses. All variable definitions are identical to Table 1.
Monthly Horizon 1 2 3 4 5 6 7 8 9 10 11 12
(1) Constant 3.13 1.74 -0.69 -2.90 -0.43 0.44 0.82 1.72 2.60 2.90 3.56 3.33
(1) Constant (t-statistic) (0.55) (0.34) (-0.14) (-0.56) (-0.09) (0.10) (0.19) (0.41) (0.64) (0.72) (0.93) (0.86)
(1)  VRP_t 0.02 0.20 0.24 0.40 0.27 0.23 0.21 0.17 0.13 0.13 0.10 0.11
(1)  VRP_t (t-statistic) (0.07) (0.85) (1.47) (3.56) (2.54) (2.36) (2.50) (2.13) (1.80) (2.04) (1.51) (2.38)
(1) Adj.  R^2 (%) -0.43 1.06 2.29 8.11 4.50 3.85 3.97 2.97 1.83 1.73 0.83 1.51
(2) Constant 75.94 80.65 72.44 67.22 61.90 58.25 54.26 52.42 52.36 52.70 53.83 55.11
(2) Constant (t-statistic) (1.84) (1.95) (1.73) (1.56) (1.48) (1.42) (1.34) (1.31) (1.34) (1.39) (1.48) (1.57)
(2) log(P/E) -23.18 -24.06 -22.00 -20.22 -18.44 -17.23 -15.91 -15.24 -15.15 -15.21 -15.53 -15.90
(2) log(P/E) (t-statistic) (-1.75) (-1.82) (-1.64) (-1.47) (-1.37) (-1.30) (-1.22) (-1.18) (-1.20) (-1.24) (-1.32) (-1.40)
(2) Adj.  R^2 (%) 0.92 2.28 2.92 3.40 3.74 3.96 3.95 4.21 4.70 5.13 5.69 6.34
(3) Constant 75.94 79.84 72.98 67.96 62.20 58.53 54.60 52.65 52.44 52.90 53.87 55.38
(3) Constant (t-statistic) (1.83) (1.90) (1.70) (1.50) (1.43) (1.38) (1.31) (1.29) (1.32) (1.38) (1.47) (1.56)
(3) log(P/E) -23.42 -25.10 -23.72 -22.83 -20.17 -18.71 (-17.31) -16.40 -16.04 -16.10 -16.18 -16.74
(3) log(P/E) (t-statistic) (-1.65) (-1.80) (-1.70) (-1.60) (-1.44) (-1.37) (-1.28) (-1.24) (-1.24) (-1.29) (-1.36) (-1.46)
(3)  VRP_t 0.04 3.59 0.27 0.42 0.29 0.24 0.23 0.19 0.15 0.14 0.11 0.13
(3)  VRP_t (t-statistic) (0.12) (1.50) (1.51) (3.61) (2.67) (2.52) (2.61) (2.21) (1.86) (2.17) (1.69) (2.77)
(3) Adj.  R^2 (%) 0.51 5.59 5.76 12.59 9.08 8.60 8.75 7.91 7.15 7.52 7.04 8.57


Table 3: Bond Returns and Variance Risk Premia

The sample period extends from January 1990 to December 2008. All of the regressions are based on monthly observations. The regression takes the form
\displaystyle xhpr_{t+h}^n \ = \ b_0^n(h) \ + \ b_1^n(h) \ VRP_t \ + \ u_{t+h,t}^n,    

where  h=1,2,3,4,5 month and  n=2,3,4,5,6 month. Robust  t-statistics following Newey and West (1987) with 12 lags are reported in parentheses. All variable definitions are identical to Table 1 Panel C.
Holding Period 2 Month Bill 3 Month Bill 4 Month Bill 5 Month Bill 6 Month Bill
1 Month/Const. 0.33 0.47 0.43 0.60 0.62
1 Month/Const. (t-statistic) (4.00) (5.66) (4.67) (5.35) (4.34)
1 Month/Const.  VRP_t 5.96e-3 5.55e-3 6.88e-3 9.85e-3 13.12e-3
1 Month/Const.  VRP_t (t-statistic) (2.38) (2.33) (2.86) (2.76) (3.23)
1 Month/Const. Adj.  R^2 4.20 2.77 2.86 4.24 4.57
2 Months/Const.   0.25 0.22 0.22 0.29
2 Months/Const. (t-statistic)   (5.69) (5.69) (5.23) (4.53)
2 Months/Const.  VRP_t   1.76e-3 2.58e-3 4.26e-3 7.17e-3
2 Months/Const.  VRP_t (t-statistic)   (1.22) (2.05) (1.95) (2.03)
2 Months/Const. Adj.  R^2   0.86 1.34 2.23 3.54
3 Months/Const.     0.16 0.26 0.32
3 Months/Const. (t-statistic)     (4.97) (5.39) (4.23)
3 Months/Const.  VRP_t     1.28e-3 1.98e-3 3.70e-3
3 Months/Const.  VRP_t (t-statistic)     (1.19) (1.24) (1.42)
3 Months/Const. Adj.  R^2     0.70 0.93 1.59
4 Months/Const.       0.18 0.28
4 Months/Const. (t-statistic)       (5.74) (4.74)
4 Months/Const.  VRP_t       0.57e-3 0.84e-3
4 Months/Const.  VRP_t (t-statistic)       (0.50) (0.50)
4 Months/Const. Adj.  R^2       -0.17 -0.22
5 Months/Const.         0.16
5 Months/Const. (t-statistic)         (4.31)
5 Months/Const.  VRP_t         0.44e-3
5 Months/Const.  VRP_t (t-statistic)         (0.38)
5 Months/Const. Adj.  R^2         -0.28


Table 4: Bond Returns and Lagged Forward Spreads

The sample period extends from January 1990 to December 2008. All of the regressions are based on monthly observations. Regression takes the form
\displaystyle xhpr_{t+h}^n \ = \ b_0^n(h) \ + \ b_2^n(h)\ \left[f_{t-1}(n-h,h)-y_{t-1}(h)\right] \ + \ u_{t+h,t}^n,    

where  h=1,2,3,4,5 month and  n=2,3,4,5,6 month. Robust  t-statistics following Newey and West (1987) with 12 lags are reported in parentheses. All variable definitions are identical to Table 1 Panel C.
Holding Period 2 Month Bill 3 Month Bill 4 Month Bill 5 Month Bill 6 Month Bill
1 Month/Const. 0.23 0.36 0.34 0.51 0.59
1 Month/Const. (t-statistic) (4.90) (5.91) (4.35) (4.28) (4.54)
1 Month/Const. Forward Spread 0.47 0.39 0.44 0.41 0.37
1 Month/Const. Forward Spread (t-statistic) (7.22) (7.93) (5.59) (4.83) (4.23)
1 Month/Const. Adj.  R^2 21.64 30.88 17.99 11.75 8.05
2 Months/Const.   0.22 0.26 0.27 0.39
2 Months/Const. (t-statistic)   (5.86) (5.74) (4.02) (4.11)
2 Months/Const. Forward Spread   0.24 0.29 0.47 0.43
2 Months/Const. Forward Spread (t-statistic)   (3.45) (4.92) (5.57) (3.16)
2 Months/Const. Adj.  R^2   5.59 9.13 18.07 8.67
3 Months/Const.     0.14 0.20 0.25
3 Months/Const. (t-statistic)     (4.54) (4.60) (3.85)
3 Months/Const. Forward Spread     0.23 0.36 0.44
3 Months/Const. Forward Spread (t-statistic)     (3.25) (5.96) (3.87)
3 Months/Const. Adj.  R^2     6.70 12.01 9.76
4 Months/Const.       0.15 0.23
4 Months/Const. (t-statistic)       (5.51) (5.03)
4 Months/Const. Forward Spread       0.18 0.24
4 Months/Const. Forward Spread (t-statistic)       (3.18) (5.60)
4 Months/Const. Adj.  R^2       4.53 5.59
5 Months/Const.         0.14
5 Months/Const. (t-statistic)         (5.01)
5 Months/Const. Forward Spread         0.12
5 Months/Const. Forward Spread (t-statistic)         (2.26)
5 Months/Const. Adj.  R^2         2.06


Table 5: Bond Returns, Variance Risk Premia, and Lagged Forward Spreads

The sample period extends from January 1990 to December 2008. All of the regressions are based on monthly observations. Regression takes the form
\displaystyle xhpr_{t+h}^n \ = \ b_0^n(h) \ + \ b_1^n(h) \ VRP_t \ + \ b_2^n(h)\ \left[f_{t-1}(n-h,h)-y_{t-1}(h)\right] \ + \ u_{t+h,t}^n,    

where  h=1,2,3,4,5 month and  n=2,3,4,5,6 month. Robust  t-statistics following Newey and West (1987) with 12 lags are reported in parentheses. All variable definitions are identical to Table 1 Panel C.
Holding Period 2 Month Bill 3 Month Bill 4 Month Bill 5 Month Bill 6 Month Bill
1 Month/Const. 0.18 0.27 0.26 0.34 0.41
1 Month/Const. (t-statistic) (3.62) (4.51) (3.17) (3.37) (3.13)
1 Month/Const. Forward Spread 0.44 0.39 0.42 0.40 0.34
1 Month/Const. Forward Spread (t-statistic) (7.08) (7.99) (5.29) (4.79) (4.06)
1 Month/Const.  VRP_t 3.65e-3 5.04e-3 4.85e-3 9.23e-3 11.10e-3
1 Month/Const.  VRP_t (t-statistic) (2.07) (3.02) (2.75) (2.99) (3.29)
1 Month/Const. Adj.  R^2 22.98 33.23 19.25 15.47 11.18
2 Months/Const.   0.20 0.22 0.22 0.29
2 Months/Const. (t-statistic)   (4.86) (4.52) (3.35) (3.08)
2 Months/Const. Forward Spread   0.23 0.28 0.46 0.40
2 Months/Const. Forward Spread (t-statistic)   (3.46) (4.84) (5.67) (3.34)
2 Months/Const.  VRP_t   1.25e-3 2.00e-3 3.11e-3 5.82e-3
2 Months/Const.  VRP_t (t-statistic)   (0.99) (1.92) (1.89) (2.11)
2 Months/Const. Adj.  R^2   5.82 9.80 19.11 10.86
3 Months/Const.     0.13 0.18 0.21
3 Months/Const. (t-statistic)     (4.31) (4.07) (2.86)
3 Months/Const. Forward Spread     0.22 0.35 0.42
3 Months/Const. Forward Spread (t-statistic)     (3.24) (5.98) (4.03)
3 Months/Const.  VRP_t     0.58e-3 1.37e-3 2.85e-3
3 Months/Const.  VRP_t (t-statistic)     (0.67) (1.08) (1.28)
3 Months/Const. Adj.  R^2     6.51 12.27 10.55
4 Months/Const.       0.15 0.22
4 Months/Const. (t-statistic)       (4.85) (3.91)
4 Months/Const. Forward Spread       0.18 0.24
4 Months/Const. Forward Spread (t-statistic)       (3.38) (5.72)
4 Months/Const.  VRP_t       0.18e-3 0.52e-3
4 Months/Const.  VRP_t (t-statistic)       (0.18) (0.33)
4 Months/Const. Adj.  R^2       4.13 5.25
5 Months/Const.         0.14
5 Months/Const. (t-statistic)         (3.18)
5 Months/Const. Forward Spread         0.12
5 Months/Const. Forward Spread (t-statistic)         (2.45)
5 Months/Const.  VRP_t         0.15e-3
5 Months/Const.  VRP_t (t-statistic)         (0.15)
5 Months/Const. Adj.  R^2         1.64


Table 6, Part 1 of 5: Credit Spreads, Variance Risk Premia, and Interest Rates - 1 Month Horizon

The sample period extends from January 1999 to December 2008. All of the regressions are based on monthly observations. Regression takes the form
CS\displaystyle _{t+h} \ = \ b_0(h) \ + \ b_1(h) \ VRP_t \ + \ b_2(h) \ r_{f,t} \ + \ u_{t+h,t},    

where the credit spread of  h month ahead is being forecasted. Robust  t-statistics following Newey and West (1987) with 24 lags are reported in parentheses. All variable definitions are identical to Table 1 Panel D.

Moody's AAA Bond Yield Spread: Constant

Moody's AAA Bond Yield Spread:  r_{f,t} Moody's AAA Bond Yield Spread:  VRP_t Moody's AAA Bond Yield Spread: Adj.  R^2 Moody's BAA Bond Yield Spread: Constant Moody's BAA Bond Yield Spread:  r_{f,t} Moody's BAA Bond Yield Spread:  VRP_t Moody's BAA Bond Yield Spread: Adj.  R^2
1.86 -0.14   32.68 3.06 -0.22   32.15
(10.56) (-3.94)     (8.18) (-2.94)    
1.17   4.68e-3 5.19 2.00   7.99e-3 6.64
(11.19)   (2.49)   (14.21)   (2.35)  
1.78 -0.15 5.85e-3 41.13 2.93 -0.23 9.76e-3 42.37
(11.74) (-4.59) (4.26)   (9.10) (-3.53) (3.77)  


Table 6, Part 2 of 5: Credit Spreads, Variance Risk Premia, and Interest Rates - 3 Month Horizon
Moody's AAA Bond Yield Spread: Constant Moody's AAA Bond Yield Spread:  r_{f,t} Moody's AAA Bond Yield Spread:  VRP_t Moody's AAA Bond Yield Spread: Adj.  R^2 Moody's BAA Bond Yield Spread: Constant Moody's BAA Bond Yield Spread:  r_{f,t} Moody's BAA Bond Yield Spread:  VRP_t Moody's BAA Bond Yield Spread: Adj.  R^2
1.77 -0.12   22.33 2.92 -0.18   21.83
(9.33) (-3.11)     (7.66) (-2.52)    
1.19   3.79e-3 2.61 2.01   7.37e-3 4.52
(10.63)   (1.84)   (13.56)   (2.08)  
1.71 -0.13 5.24e-3 27.75 2.81 -0.19 9.58e-3 29.77
(10.00) (-3.49) (3.09)   (8.39) (-3.00) (2.87)  


Table 6, Part 3 of 5: Credit Spreads, Variance Risk Premia, and Interest Rates - 6 Month Horzon
Moody's AAA Bond Yield Spread: Constant Moody's AAA Bond Yield Spread:  r_{f,t} Moody's AAA Bond Yield Spread:  VRP_t Moody's AAA Bond Yield Spread: Adj.  R^2 Moody's BAA Bond Yield Spread: Constant Moody's BAA Bond Yield Spread:  r_{f,t} Moody's BAA Bond Yield Spread:  VRP_t Moody's BAA Bond Yield Spread: Adj.  R^2
1.66 -0.09   12.78 2.75 -0.14   12.32
(7.80) (-2.17)     (6.92) (-1.93)    
1.21   2.72e-3 0.93 2.08   3.93e-3 0.77
(10.37)   (1.19)   (12.31)   (1.35)  
1.61 -0.10 4.26e-3 15.71 2.69 -0.15 6.25e-3 14.95
(7.94) (-2.43) (2.17)   (7.08) (-2.21) (2.28)  


Table 6, Part 4 of 5: Credit Spreads, Variance Risk Premia, and Interest Rates - 9 Month Horizon
Moody's AAA Bond Yield Spread: Constant Moody's AAA Bond Yield Spread:  r_{f,t} Moody's AAA Bond Yield Spread:  VRP_t Moody's AAA Bond Yield Spread: Adj.  R^2 Moody's BAA Bond Yield Spread: Constant Moody's BAA Bond Yield Spread:  r_{f,t} Moody's BAA Bond Yield Spread:  VRP_t Moody's BAA Bond Yield Spread: Adj.  R^2
1.53 -0.06   4.97 2.55 -0.09   5.00
(6.65) (-1.30)     (6.59) (-1.35)    
1.22   2.94e-3 1.17 2.07   5.15 1.64
(9.89)   (1.29)   (12.36)   (1.48)  
1.49 -0.07 3.90e-3 7.35 2.49 -0.10 6.67e-3 8.00
(6.71) (-1.49) (1.86)   (6.92) (-1.60) (1.85)  


Table 6, Part 5 of 5: Credit Spreads, Variance Risk Premia, and Interest Rates
Moody's AAA Bond Yield Spread: Constant Moody's AAA Bond Yield Spread:  r_{f,t} Moody's AAA Bond Yield Spread:  VRP_t Moody's AAA Bond Yield Spread: Adj.  R^2 Moody's BAA Bond Yield Spread: Constant Moody's BAA Bond Yield Spread:  r_{f,t} Moody's BAA Bond Yield Spread:  VRP_t Moody's BAA Bond Yield Spread: Adj.  R^2
1.37 -0.02   0.25 2.30 -0.03   0.23
(5.82) (-0.45)     (6.92) (-0.57)    
1.20   4.20e-3 2.85 2.06   5.49e-3 1.88
(9.66)   (2.05)   11.49   (1.92)  
1.33 -0.03 4.67e-3 3.78 2.24 -0.04 6.19e-3 2.67
(5.88) (-0.66) (2.50)   (7.22) (-0.80) (2.10)  


Table 7: Correlations between Different Markets

This table reports the correlations between equity, bond, and credit markets for the raw excess returns and the residual excess returns after regressing on the variance risk premium (VRP) variable. The correlations are for the monthly sample period 1990-2008.
Table 7, Panel A: Raw Correlations
  Equity 2 M. Bill 3 M. Bill 4 M. Bill 5 M. Bill 6 M. Bill AAA BAA
Equity 1.00 0.04 0.03 -0.06 -0.06 -0.12 -0.15 -0.18
2 M. Bill   1.00 0.91 0.80 0.73 0.61 -0.17 -0.16
3 M. Bill     1.00 0.92 0.88 0.77 -0.11 -0.09
4 M. Bill       1.00 0.95 0.86 -0.07 -0.01
5 M. Bill         1.00 0.92 -0.09 0.06
6 M. Bill           1.00 0.01 0.09
AAA             1.00 0.91
BAA               1.00

Table 7, Panel B: Residual Correlations from Regression on  VRP
  Equity 2 M. Bill 3 M. Bill 4 M. Bill 5 M. Bill 6 M. Bill AAA BAA
Equity 1.00 0.08 0.03 -0.04 -0.09 -0.08 -0.21 -0.29
2 M. Bill   1.00 0.91 0.79 0.71 0.59 -0.24 -0.25
3 M. Bill     1.00 0.92 0.87 0.76 -0.18 -0.16
4 M. Bill       1.00 0.95 0.86 -0.12 -0.07
5 M. Bill         1.00 0.92 -0.07 0.01
6 M. Bill           1.00 -0.05 0.02
AAA             1.00 0.90
BAA               1.00


Table 8: Model Calibration Parameter Setting

This table reports the calibration parameter values for the stochastic volatility-of-volatility model used in this paper. BTZ2009 refers to the calibration setting of Bollerslev, Tauchen, and Zhou (2009), with an emphasis on equity risk premium and its short-run predictability. The Campbell-Shiller linearization constants are  \kappa_1 = 0.9 and  \kappa_0 = 0.3251.
Calibration Parameters BTZ2009 This Paper
Preference Parameters:  \delta = 0.997

 \gamma= 10

 \psi = 1.5

 \delta = 0.997

 \gamma= 2

 \psi = 1.5

Endowment Parameters:  \mu_g = 0.0015

 a_{\sigma}=1.34\times10^{-6}

 \rho_{\sigma}=0.978

 a_q=2\times10^{-7}

 \rho_q=0.8

 \varphi_q = 0.001

 \mu_g = 0.0015

 a_{\sigma}=2.14\times10^{-5}

 \rho_{\sigma}=0.978

 a_q=1.825\times10^{-5}

 \rho_q = 0.95

 \varphi_q = 0.008


Table 9: Calibrated Equity, Variance, and Bond Risk Premia

This table reports the calibration output values for the stochastic volatility-of-volatility model used in this paper. BTZ2009 refers to the calibration setting of Bollerslev, Tauchen, and Zhou (2009). For equity risk premiums, the long historical sample is based on the annual CRSP data from 1930 to 2008 (Bansal, Kiku, and Yaron, 2009), and the short historical sample is based on the annual S&P data from 1891 to 1949 (Gvozdeva and Kumar, 2009). While the short data sample is based on the monthly data from 1990 to 2008 used in this paper. The calibrated equity premiums and bond risk premia are analytical results, while the variance premium moments are based a simulated sample of 1,000,000 months.
Calibration Target Variables 1930-2008 1891-1949 1990-2008 BTZ2009 This Paper
Equity Risk Premium (%): Equity Risk Premium 7.84 4.97 3.58 7.79 5.11
Equity Risk Premium (%): Equity Premium Volatility 20.16 20.36 14.60 4.02 13.24
Equity Risk Premium (%): Risk-Free Rate 0.86 1.96 1.13 0.69 1.86
Equity Risk Premium (%): Risk-Free Rate Volatility 1.74 7.83 3.37 2.95 8.83
Variance Risk Premium(% ^2): Mean --- --- 18.30 3.70 18.30
Variance Risk Premium(% ^2): Std Dev --- --- 22.69 7.29 25.12
Variance Risk Premium(% ^2): Skewness --- --- 2.79 1.70 2.48
Variance Risk Premium(% ^2): Kurtosis --- --- 16.62 11.42 13.18
Variance Risk Premium(% ^2): AR(1) --- --- 0.26 0.80 0.95
Bond Risk Premium (%): 2 Month Bill --- --- 0.44 2.99 2.92
Bond Risk Premium (%): 3 Month Bill --- --- 0.57 5.81 5.78
Bond Risk Premium (%): 4 Month Bill --- --- 0.56 7.08 8.56
Bond Risk Premium (%): 5 Month Bill --- --- 0.78 -1.41  \times10^2 11.27
Bond Risk Premium (%): 6 Month Bill --- --- 0.86 -1.06  \times10^6 13.90



Figure 1: Variance Risk Premium, Implied and Realized Variances
Figure 1: Variance Risk Premium, Implied and Realized Variances. This figure plots the variance risk premium or implied-expected variance difference (top panel), the implied variance (middle panel), and the realized variance (bottom panel) for the S&P500 market index from January 1990 to December 2008. The variance risk premium is based on the realized variance forecast from a full sample AR(12). The shaded areas represent NBER recessions.
This figure plots the variance risk premium or implied-expected variance difference (top panel), the implied variance (middle panel), and the realized variance (bottom panel) for the S&P500 market index from January 1990 to December 2008. The variance risk premium is based on the realized variance forecast from a full sample AR(12). The shaded areas represent NBER recessions.


Figure 2: Variance Risk Premiums with Alternative  RV_t Forecasts
Figure 2: Variance Risk Premiums with Alternative RV_t Forecasts. The figure plots the variance risk premium series constructed using alternative ways to forecast the realized variance: (1) ex post $RV_t$, (2) lagged $RV_{t-1}$, (3) MA(12) estimate, and (4) recursive AR(12) estimate. The shaded areas represent NBER recessions.
The figure plots the variance risk premium series constructed using alternative ways to forecast the realized variance: (1) ex post  RV_t, (2) lagged  RV_{t-1}, (3) MA(12) estimate, and (4) recursive AR(12) estimate. The shaded areas represent NBER recessions.


Figure 3: Estimated Slopes and  R^2's of Equity Returns
Figure 3: Estimated Slopes and R^2's of Equity Returns. The figure shows the estimated slope coefficients and pointwise 95 percent confidence intervals, along with the corresponding adjusted $R^2$'s from the regressions of the scaled $h$-month S&P500 excess returns on the variance risk premium and P/E ratio. All of the regressions are based on monthly observations from January 1990 to December 2008.
The figure shows the estimated slope coefficients and pointwise 95 percent confidence intervals, along with the corresponding adjusted  R^2's from the regressions of the scaled  h-month S&P500 excess returns on the variance risk premium and P/E ratio. All of the regressions are based on monthly observations from January 1990 to December 2008.


Figure 4: Estimated Slopes and  R^2's of 6 Month T-Bill Returns
Figure 4: Estimated Slopes and R^2's of 6 Month T-Bill Returns. The figure shows the estimated slope coefficients and pointwise 95 percent confidence intervals, along with the corresponding adjusted $R^2$'s from the regressions of the 6 month t-bill excess returns with one-to-five month holding periods on the variance risk premium and lagged forward spread. All of the regressions are based on monthly observations from January 1990 to December 2008.
The figure shows the estimated slope coefficients and pointwise 95 percent confidence intervals, along with the corresponding adjusted  R^2's from the regressions of the 6 month t-bill excess returns with one-to-five month holding periods on the variance risk premium and lagged forward spread. All of the regressions are based on monthly observations from January 1990 to December 2008.


Figure 5: Estimated Slopes and  R^2's of BAA Credit Spreads
Figure 5: Estimated Slopes and R^2's of BAA Credit Spreads. The figure shows the estimated slope coefficients and pointwise 95 percent confidence intervals, along with the corresponding adjusted $R^2$'s from the regressions of the $h$-month ahead credit spread on the variance risk premium and short rate. All of the regressions are based on monthly observations from January 1990 to December 2008.
The figure shows the estimated slope coefficients and pointwise 95 percent confidence intervals, along with the corresponding adjusted  R^2's from the regressions of the  h-month ahead credit spread on the variance risk premium and short rate. All of the regressions are based on monthly observations from January 1990 to December 2008.


Figure 6: Model-Implied and Estimated Slopes and  R^2's for 2-6 Month T-Bills
Figure 6: Model-Implied and Estimated Slopes and R^2's for 2-6 Month T-Bills. The figure shows the calibrated model-implied slope coefficients and adjusted $R^2$'s for regressing the 2-6 months t-bill excess returns on variance risk premium, along with their estimated empirical counterparts from Table 3.
The figure shows the calibrated model-implied slope coefficients and adjusted  R^2's for regressing the 2-6 months t-bill excess returns on variance risk premium, along with their estimated empirical counterparts from Table 3.




Footnotes

1. I benefited from helpful discussions with Gurdip Bakshi, Ravi Bansal, Tim Bollerslev, John Campbell, Peter Christoffersen, Pete Kyle, George Tauchen, Jonathan Wright, Hong Yan, Moto Yogo, Harold Zhang, and Guofu Zhou. I am also grateful for the comments received from seminar participants at BIS--Hong Kong, Tsinghua University, FRB--San Francisco, Renmin University, University of Texas at Dallas, Duke University, and Purdue University. The author acknowledges the recognitions from CQA Academic Competition (3rd place) and Crowell Prize Finalist by PanAgora Asset Managment. James Marrone and Paul Reverdy provided excellent research assistance. The views presented here are solely those of the author and do not necessarily represent those of the Federal Reserve Board or its staff. Return to Text
2. Risk Analysis Section, Federal Reserve Board, Mail Stop 91, Washington DC 20551, USA; E-mail [email protected], Phone 202-452-3360, Fax 202-728-5887. Return to Text
3. For example, in the final quarter of 2008, the VIX index has closed above 50 percent for almost twelve weeks and peaked around 90 percent. As reported by the Wall Street Journal on November 12, 2008, if market volatility continues to remain above 50 percent for just over five weeks, it would have surpassed the Great Depression in the 1930s; and "such a high volatility signifies all those unknowns that are a greater cloud of what we call Uncertainty". Return to Text
4. Beeler and Campbell (2009) show that the calibration setting in Bansal, Kiku, and Yaron (2007) puts more emphasis on the persistent volatility channel as opposed to the long-run risk channel as in Bansal and Yaron (2004). Return to Text
5. Pástor and Stambaugh (2009) use the uncertainties of expected future and current returns to argue that long-run stock returns are indeed more volatile. Bekaert, Hoerova, and Scheicher (2009) use VIX and VDAX to assess the relative importance of time-varying risk aversion versus economic uncertainty for both US and Germany. Buraschi, Trojani, and Vedolin (2009) examine uncertainty's effect on volatility risk premium in a heterogeneous beliefs setting with independent Lucas trees. Return to Text
6. There are quite a few extensions following the similar direction: Zhou and Zhu (2009) incorporate two volatility factors into the long-run risks model to explain equity return predictability and volatility risk premium; Benzoni, Collin-Dufresne, and Goldstein (2008) use a similar stochastic volatility model to explain the option pricing puzzles pre- and post-1987 stock market crash; and Bollerslev, Sizova, and Tauchen (2008b) apply the same framework to interpret the volatility asymmetry and dynamic dependency puzzles. Return to Text
7. Such a characterization abstracts from the realistic economic environment that allows for (1) lumpy dividend payment, (2) stochastic interest rate, (3) underlying asset jumps, and (4) limited number and range of option strikes--discretization and truncation errors. See Jiang and Tian (2005) for detailed discussions. Return to Text
8. For the highly persistent predictors like P/E ratio and short rate, the conventional  t-statistics and  R^2's for the overlapping multi-period return regressions need to be interpreted with caution (Boudoukh, Richardson, and Whitelaw, 2008; Goyal and Welch, 2008; Stambaugh, 1999; Campbell and Yogo, 2006). Return to Text
9. The forward rate regression is recently extended by Cochrane and Piazzesi (2005) to multiple forward rates, by Ludvigson and Ng (2008) to incorporate extracted macroeconomic factors, and by Wright and Zhou (2009) to augment with a realized jump risk measure. However, these studies use 2-5 year zero coupon bonds with a one year holding period, where the variance risk premium variable has virtually zero forecasting power of the bond risk premia. Return to Text
10. Note that Fama and Bliss (1987) suggested using the lagged forward spread to break the potential first order serial correlation in the market microstructure error, which may be more relevant to our short-run t-bill predictability regressions. The results based on the current forward spread are similar with higher  R^2's and are available upon request. Return to Text
11. If one includes the term spread alone, it is marginally significant with  R^2 of 8-9 percent and  t-statistics of 1.73 and 2.09. However, when short rate and term spread are combined together, term spread is driven out as insignificant with  t-statistics being -0.94 and -0.93. This result in tabular form is available upon request. Return to Text
12. For example, this is indeed the case with the t-bill returns in this paper--once regressed on both variance premium and forward spread, the residual correlations are significantly decreased in all cases from a few percentage points to as much as 29 percentage points. This result in tabular form is available upon request. Return to Text
13. In fact, traditional consumption-based asset pricing models cannot even qualitatively explain the existence of risk premia in bond and credit markets; and can only explain a negligible portion of the observed equity premium (Mehra and Prescott, 1985). Return to Text
14. For example, the dividend volatility is levered up to 5.96 times of the consumption volatility in Bansal, Kiku, and Yaron (2007). Return to Text
15. In fact, both Bollerslev, Tauchen, and Zhou (2009) and Bansal and Yaron (2004) would produce a five-year real rate near negative infinity; which, as argued by Bansal, Kiku, and Yaron (2009), would not hamper the resolution of the Expectations Hypothesis puzzle in the nominal interest rates, if an appropriate inflation process is attached to the long-run risks model. Return to Text
16. This idea may also be observationally equivalent to the approaches based on option-implied jump risk premia (Cremers, Driessen, and Maenhout, 2007) and macroeconomic risk induced by the inflation uncertainty with heterogeneous beliefs (David, 2008). Return to Text

This version is optimized for use by screen readers. Descriptions for all mathematical expressions are provided in LaTex format. A printable pdf version is available. Return to Text