Dynamic estimation of volatility risk premia and investor risk aversion from option-implied and realized volatilities

This paper proposes a method for constructing a volatility risk premium, or investor risk aversion, index. The method is intuitive and simple to implement, relying on the sample moments of the recently popularized model-free realized and option-implied volatility measures. A small-scale Monte Carlo experiment confirms that the procedure works well in practice. Implementing the procedure with actual S&P500 option-implied volatilities and high-frequency five-minute-based realized volatilities indicates significant temporal dependencies in the estimated stochastic volatility risk premium, which we in turn relate to a set of macro-finance state variables. We also find that the extracted volatility risk premium helps predict future stock market returns.     

Jump tails,extreme dependencies,and the distribution of stock returns

We provide a new framework for estimating the systematic and idiosyncratic jump tail risks in financial asset prices. Our estimates are based on in-fill asymptotics for directly identifying the jumps, together with Extreme Value Theory (EVT) approximations and methods-of-moments for assessing the tail decay parameters and tail dependencies. On implementing the procedures with a panel of intraday prices for a large cross-section of individual stocks and the S&P 500 market portfolio, we find that the distributions of the systematic and idiosyncratic jumps are both generally heavy-tailed and close to symmetric, and show how the jump tail dependencies deduced from the high-frequency data together with the day-to-day variation in the diffusive volatility account for the “extreme” joint dependencies observed at the daily level.     

Variance trading and market price of variance risk

This paper develops a new approach for variance trading. We show that the discretely-sampled realized variance can be robustly replicated under very general conditions, including when the price can jump. The replication strategy specifies the exact timing for rebalancing in the underlying. The deviations from the optimal schedule can lead to surprisingly large hedging errors. In the empirical application, we synthesize the prices of the variance contract on S&P 500 index over the period from 01/1990 to 12/2009. We find that the market variance risk is priced, its risk premium is negative and economically very large. The variance risk premium cannot be explained by the known risk factors and option returns.     

Hermite polynomial based expansion of European option prices

We seek a closed-form series approximation of European option prices under a variety of diffusion models. The proposed convergent series are derived using the Hermite polynomial approach. Departing from the usual option pricing routine in the literature, our model assumptions have no requirements for affine dynamics or explicit characteristic functions. Moreover, convergent expansions provide a distinct insight into how and on which order the model parameters affect option prices, in contrast with small-time asymptotic expansions in the literature. With closed-form expansions, we explicitly translate model features into option prices, such as mean-reverting drift and self-exciting or skewed jumps. Numerical examples illustrate the accuracy of this approach and its advantage over alternative expansion methods.     

Chi-squared tests for evaluation and comparison of asset pricing models

This paper presents a general statistical framework for estimation, testing and comparison of asset pricing models using the unconstrained distance measure of Hansen and Jagannathan (1997). The limiting results cover both linear and nonlinear models that could be correctly specified or misspecified. We propose modified versions of the existing model selection tests and new pivotal specification and model comparison tests with improved finite-sample properties. In addition, we provide formal tests of multiple model comparison. The excellent size and power properties of the proposed tests are demonstrated using simulated data from linear and nonlinear asset pricing models.     

Jumps and betas: A new framework for disentangling and estimating systematic risks

We provide a new theoretical framework for disentangling and estimating the sensitivity towards systematic diffusive and jump risks in the context of factor models. Our estimates of the sensitivities towards systematic risks, or betas, are based on the notion of increasingly finer sampled returns over fixed time intervals. We show consistency and derive the asymptotic distributions of our estimators. In an empirical application of the new procedures involving high-frequency data for forty individual stocks, we find that the estimated monthly diffusive and jump betas with respect to an aggregate market portfolio differ substantially for some of the stocks in the sample.     

Mutual excitation in Eurozone sovereign CDS

We study self- and cross-excitation of shocks in the Eurozone sovereign CDS market. We adopt a multivariate setting with credit default intensities driven by mutually exciting jump processes, to capture the salient features observed in the data, in particular, the clustering of high default probabilities both in time (over days) and in space (across countries). The feedback between jump events and the intensity of these jumps is the key element of the model. We derive closed-form formulae for CDS prices, and estimate the model by matching theoretical prices to their empirical counterparts. We find evidence of self-excitation and asymmetric cross-excitation. Using impulse-response analysis, we assess the impact of shocks and a potential policy intervention not just on a single country under scrutiny but also, through the effect on cross-excitation risk which generates systemic sovereign risk, on other interconnected countries.     

Spanning tests in return and stochastic discount factor mean–variance frontiers: A unifying approach

We propose new spanning tests that assess if the initial and additional assets share the economically meaningful cost and mean representing portfolios. We prove their asymptotic equivalence to existing tests under local alternatives. We also show that unlike two-step or iterated procedures, single-step methods such as continuously updated GMM yield numerically identical overidentifying restrictions test, so there is arguably a single spanning test. To prove these results, we extend optimal GMM inference to deal with singularities in the long run second moment matrix of the influence functions. Finally, we test for spanning using size and book-to-market sorted US stock portfolios.     

Realized Laplace transforms for estimation of jump diffusive volatility models

We develop an efficient and analytically tractable method for estimation of parametric volatility models that is robust to price-level jumps. The method entails first integrating intra-day data into the Realized Laplace Transform of volatility, which is a model-free estimate of the daily integrated empirical Laplace transform of the unobservable volatility. The estimation is then done by matching moments of the integrated joint Laplace transform with those implied by the parametric volatility model. In the empirical application, the best fitting volatility model is a non-diffusive two-factor model where low activity jumps drive its persistent component and more active jumps drive the transient one.     

Econometric analysis of jump-driven stochastic volatility models

This paper introduces and studies the econometric properties of a general new class of models, which I refer to as jump-driven stochastic volatility models, in which the volatility is a moving average of past jumps. I focus attention on two particular semiparametric classes of jump-driven stochastic volatility models. In the first, the price has a continuous component with time-varying volatility and time-homogeneous jumps. The second jump-driven stochastic volatility model analyzed here has only jumps in the price, which have time-varying size. In the empirical application I model the memory of the stochastic variance with a CARMA(2,1) kernel and set the jumps in the variance to be proportional to the squared price jumps. The estimation, which is based on matching moments of certain realized power variation statistics calculated from high-frequency foreign exchange data, shows that the jump-driven stochastic volatility model containing continuous component in the price performs best. It outperforms a standard two-factor affine jump–diffusion model, but also the pure-jump jump-driven stochastic volatility model for the particular jump specification.     

Equilibrium price dynamics in an overlapping-generations exchange economy

An extension to the Yaari (1965)–Blanchard (1985) continuous time overlapping generations model for an endowment Arrow–Debreu economy with an age-structured population is presented. It is proved that Arrow–Debreu equilibrium prices are represented by a double linear integral equation, and depend on the age-distributions of population and endowments. For an economy with a balanced growth, and logarithmic utility, we prove that bubbles may exist if endowments are distributed earlier than some critical age.     

The common and specific components of dynamic volatility

This paper develops a dynamic approximate factor model in which returns are time-series heteroskedastic. The heteroskedasticity has three components: a factor-related component, a common asset-specific component, and a purely asset-specific component. We develop a new multivariate GARCH model for the factor-related component. We develop a univariate stochastic volatility model linked to a cross-sectional series of individual GARCH models for the common asset-specific component and the purely asset-specific component. We apply the analysis to monthly US equity returns for the period January 1926 to December 2000. We find that all three components contribute to the heteroskedasticity of individual equity returns. Factor volatility and the common component in asset-specific volatility have long-term secular trends as well as short-term autocorrelation. Factor volatility has correlation with interest rates and the business cycle.     

Linear and nonlinear regression with stable errors

In this paper we describe methods and evaluate programs for linear regression by maximum likelihood when the errors have a heavy tailed stable distribution. The asymptotic Fisher information matrix for both the regression coefficients and the error distribution parameters are derived, giving large sample confidence intervals for all parameters. Simulated examples are shown where the errors are stably distributed and also where the errors are heavy tailed but are not stable, as well as a real example using financial data. The results are then extended to nonlinear models and to non-homogeneous error terms.     

Modelling volatility by variance decomposition

In this paper, we propose two parametric alternatives to the standard GJR-GARCH model of Glosten et al. (1993), based on additive and multiplicative decompositions of the variance. They allow the variance of the model to have a smooth time-varying structure. The suggested parameterizations describe structural change in the conditional and unconditional variances where the transition between regimes over time is smooth. The main focus is on the multiplicative decomposition of the variance into an unconditional and conditional components. Estimation of the multiplicative model is discussed in detail. An empirical application to daily stock returns illustrates the functioning of the model. The results show that the ‘long memory type behaviour’ of the sample autocorrelation functions of the absolute returns can also be explained by deterministic changes in the unconditional variance.     

On the properties of the coefficient of determination in regression models with infinite variance variables

We examine the asymptotic properties of the coefficient of determination, R²$R^2$, in models with α-stable$\alpha \text{-stable}$   random variables. If the regressor and error term share the same index of stability α<2$\alpha <2$, we show that the R²$R^2$  statistic does not converge to a constant but has a nondegenerate distribution on the entire [0,1]$[0,1]$ interval. We provide closed-form expressions for the cumulative distribution function and probability density function of this limit random variable, and we show that the density function is unbounded at 0 and 1. If the indices of stability of the regressor and error term are unequal, we show that the coefficient of determination converges in probability to either 0 or 1, depending on which variable has the smaller index of stability, irrespective of the value of the slope coefficient. In an empirical application, we revisit the Fama and MacBeth (1973) two-stage regression and demonstrate that in the infinite-variance case the R²$R^2$  statistic of the second-stage regression converges to 0 in probability even if the slope coefficient is nonzero. We deduce that a small value of the R²$R^2$  statistic should not, in itself, be used to reject the usefulness of a regression model.     

Term structure models and the zero bound: An empirical investigation of Japanese yields

When Japanese short-term bond yields were near their zero bound, yields on long-term bonds showed substantial fluctuation, and there was a strong positive relationship between the level of interest rates and yield volatilities/risk premiums. We explore whether several families of dynamic term structure models that enforce a zero lower bound on short rates imply conditional distributions of Japanese bond yields consistent with these patterns. Multi-factor “shadow-rate” and quadratic-Gaussian models, evaluated at their maximum likelihood estimates, capture many features of the data. Furthermore, model-implied risk premiums track realized excess returns during extended periods of near-zero short rates. In contrast, the conditional distributions implied by non-negative affine models do not match their sample counterparts, and standard Gaussian affine models generate implausibly large negative risk premiums.     

No-arbitrage macroeconomic determinants of the yield curve

No-arbitrage macro-finance models use variance decompositions to gauge the extent of association between the macro variables and yields. We show that results generated by this approach are sensitive to the order of variables in the recursive identification scheme. In a four-factor model, one may obtain 18 different sets of answers out of 24 possible. We propose an alternative measure that is based on levels of macro variables as opposed to shocks. We account for the correlation between the macro and latent factors via projection of the latter onto the former. As a result, the association between macro variables and yields can be computed uniquely via an R²$R^2$. Macro variables explain 80% of the variation in the short rate and 50% of the slope, and 54% to 68% of the term premia.     

Density approximations for multivariate affine jump-diffusion processes

We introduce closed-form transition density expansions for multivariate affine jump-diffusion processes. The expansions rely on a general approximation theory which we develop in weighted Hilbert spaces for random variables which possess all polynomial moments. We establish parametric conditions which guarantee existence and differentiability of transition densities of affine models and show how they naturally fit into the approximation framework. Empirical applications in option pricing, credit risk, and likelihood inference highlight the usefulness of our expansions. The approximations are extremely fast to evaluate, and they perform very accurately and numerically stable.