Estimating affine multifactor term structure models using closed-form likelihood expansions

We develop and implement a technique for closed-form maximum likelihood estimation (MLE) of multifactor affine yield models. We derive closed-form approximations to likelihoods for nine Dai and Singleton (2000) affine models. Simulations show our technique very accurately approximates true (but infeasible) MLE. Using US Treasury data, we estimate nine affine yield models with different market price of risk specifications. MLE allows non-nested model comparison using likelihood ratio tests; the preferred model depends on the market price of risk. Estimation with simulated and real data suggests our technique is much closer to true MLE than Euler and quasi-maximum likelihood (QML) methods.     

Non-linear Gaussian sovereign CDS pricing models

Prior literature indicates that quadratic models and the Black–Karasinski model are very promising for CDS pricing. This paper extends these models and the Black [J. Finance 1995, 50, 1371–1376] model for pricing sovereign CDS’s. For all 10 sovereigns in the sample quadratic models best fit CDS spreads in-sample, and a four factor quadratic model can account for the joint effects on CDS spreads of default risk, default loss risk and liquidity risk with no restriction to factors correlation. Liquidity risk appears to affect sovereign CDS spreads. However, quadratic models tend to over-fit some CDS maturities at the expense of other maturities, while the BK model is particularly immune from this tendency. The Black model seems preferable because its out-of-sample performance in the time series dimension is the best.     

Approximate pricing of swaptions in affine and quadratic models

This paper proposes new bounds on the prices of European-style swaptions for affine and quadratic interest rate models. These bounds are computable whenever the joint characteristic function of the state variables is known. In particular, our lower bound involves the computation of a one-dimensional Fourier transform independently of the swap length. In addition, we control the error of our method by providing a new upper bound on swaption price that is applicable to all considered models. We test our bounds on different affine models and on a quadratic Gaussian model. We also apply our procedure to the multiple curve framework. The bounds are found to be accurate and computationally efficient.     

Do Bonds Span Volatility Risk in the U.S. Treasury Market? A Specification Test for Affine Term Structure Models

We propose using model‐free yield quadratic variation measures computed from intraday data as a tool for specification testing and selection of dynamic term structure models. We find that the yield curve fails to span realized yield volatility in the U.S. Treasury market, as the systematic volatility factors are largely unrelated to the cross‐section of yields. We conclude that a broad class of affine diffusive, quadratic Gaussian, and affine jump‐diffusive models cannot accommodate the observed yield volatility dynamics. Hence, the Treasury market per se is incomplete, as yield volatility risk cannot be hedged solely through Treasury securities.     

Discounting earnings with stochastic discount rates

This paper presents new equity valuation formulae in closed form that extend the abnormal earnings growth (AEG) valuation of Ohlson [2005. “On Accounting-Based Valuation Formulae.” Review of Accounting Studies 10: 323–347] to the cases of time-varying or stochastic cost of capital as in Ang and Liu [2004. “How to Discount Cash Flows with Time-Varying Expected Returns.” Journal of Finance 59 (6): 2745–2783] or to cases of stochastic interest rates as in Ang and Liu [2001. “A General Affine Earnings Valuation Model.” Review of Accounting Studies 6: 397–425]. Interest rates are modelled by quadratic term structure models, which are not hindered by restrictions to factors correlation or by other shortcomings of affine term structure models in discounting long-term earnings. This is crucial since valuation can be very sensitive to the correlation between the factors driving earnings and interest rates. Positive correlation reduces price-earnings ratios according to US data. Valuation is also sensitive to the ‘volatility’ of abnormal earnings growth. The residual earnings risk-neutral valuation of Ang and Liu (2001) is adapted to quadratic term structure models.     

On break-even correlation: the way to price structured credit derivatives by replication

We consider the pricing of European-style structured credit pay-off under the Gaussian Copula Model (GCM). When no sudden jump-to-default events occur, the perfect replication of these pay-offs under the GCM is obtained if and only if the underlying single-name credit spreads follow a particular family of dynamics and if the pricing parameters are given by so-called ‘break-even’ correlations. We exhibit a class of Merton-style models that are consistent with this result. We calculate break-even correlations explicitly to price nth-to-default baskets under the GCM. Finally, we illustrate the usefulness of this concept as a relative-value tool.     

Term structure dynamics with macro-factors using high frequency data

This paper empirically studies the role of macro-factors in explaining and predicting daily bond yields. In general, macro-finance models use low-frequency data to match with macroeconomic variables available only at low frequencies. To deal with this, we construct and estimate a tractable no-arbitrage affine model with both conventional latent factors and macro-factors by imposing cross-equation restrictions on the daily yields of bonds with different maturities, credit risks, and inflation indexation. The estimation results using both the US and the UK data show that the estimated macro-factors significantly predict actual inflation and the output gap. In addition, our daily macro-term structure model forecasts better than no-arbitrage models with only latent factors as well as other statistical models.     

A New Approach for Using Lévy Processes for Determining High-Frequency Value-at-Risk Predictions

A new approach for using Lévy processes to compute value-at-risk (VaR) using high-frequency data is presented in this paper. The approach is a parametric model using an ARMA(1,1)-GARCH(1,1) model where the tail events are modelled using fractional Lévy stable noise and Lévy stable distribution. Using high-frequency data for the German DAX Index, the VaR estimates from this approach are compared to those of a standard nonparametric estimation method that captures the empirical distribution function, and with models where tail events are modelled using Gaussian distribution and fractional Gaussian noise. The results suggest that the proposed parametric approach yields superior predictive performance.     

What does a term structure model imply about very long-term interest rates?

We address two empirical issues related to the long end of the yield curve based on euro swap rates. First, for maturities longer than 20 years we find evidence for an ‘excess’ downward slope that cannot be explained by convexity. Second, volatility at the very long end of the yield curve is larger than predicted by no-arbitrage models. We construct a model-based arbitrage-free extrapolation of the yield-curve and compare it to the regulatory discount curve. Because of near-zero mean reversion, there is no convergence towards an ‘ultimate forward rate’ and convexity effects cause the arbitrage-free extrapolations to have slightly downward sloping curves. The low level of mean-reversion also implies that the volatility of long-term rates does not decline relative to the 20-year volatility. Therefore, we conclude that the mean-reversion and resulting smoothing adopted by the regulatory curve is much too strong.     

The affine Heston model with correlated Gaussian interest rates for pricing hybrid derivatives

In this article we define a multi-factor equity–interest rate hybrid model with non-zero correlation between the stock and interest rate. The equity part is modeled by the Heston model and we use a Gaussian multi-factor short-rate process. By construction, the model fits in the framework of affine diffusion processes, allowing fast calibration to plain vanilla options. We also provide an efficient Monte Carlo simulation scheme.     

“Extended Black” term structure models

This paper examines “Extended Black” term structure models (EBTSM), which are multi-factor extensions of the one-factor Black model (Black, F., 1995. Interest rates as options. Journal of Finance 50, 1371-1376). EBTSM are not affected by the admissibility restrictions that plague canonical affine models. EBTSM encompass quadratic models, but unlike in quadratic models bond yields are sufficient statistics to infer the latent factors driving the short interest rate. EBTSM are amenable to econometric estimation despite the need to solve bond pricing equations through finite difference numerical methods. Estimation through the Iterated Extended Kalman filter reveals that a two-factor EBTSM fit well the observed cross section and time series of Japanese Government bond yields. A three-factor EBTSM is also proposed.     

Quadratic Hawkes processes for financial prices

We introduce and establish the main properties of QHawkes (‘Quadratic’ Hawkes) models. QHawkes models generalize the Hawkes price models introduced in Bacry and Muzy [Quant. Finance, 2014, 14(7), 1147–1166], by allowing feedback effects in the jump intensity that are linear and quadratic in past returns. Our model exhibits two main properties that we believe are crucial in the modelling and the understanding of the volatility process: first, the model is time-reversal asymmetric, similar to financial markets whose time evolution has a preferred direction. Second, it generates a multiplicative, fat-tailed volatility process, that we characterize in detail in the case of exponentially decaying kernels, and which is linked to Pearson diffusions in the continuous limit. Several other interesting properties of QHawkes processes are discussed, in particular the fact that they can generate long memory without necessarily being at the critical point. A non-parametric fit of the QHawkes model on NYSE stock data shows that the off-diagonal component of the quadratic kernel indeed has a structure that standard Hawkes models fail to reproduce. We provide numerical simulations of our calibrated QHawkes model which is indeed seen to reproduce, with only a small amount of quadratic non-linearity, the correct magnitude of fat-tails and time reversal asymmetry seen in empirical time series.     

A note on “Modelling exchange rate returns: which flexible distribution to use?”

Corlu and Corlu [Quant. Finance, 2014, doi: 10.1080/14697688.2014.942231] provided a novel modelling of exchange rate data for nine currencies using five flexible distributions. They stated that the generalized lambda, skew t and normal inverse Gaussian distributions ‘do a good job’. Here, we reanalyse the data and show that a distribution simpler than all of these fits at least as well as these distributions. We also find that the normal inverse Gaussian distribution provides good fits for only one of the data-sets.     

Can interest rate volatility be extracted from the cross section of bond yields?

Most affine models of the term structure with stochastic volatility predict that the variance of the short rate should play a ‘dual role’ in that it should also equal a linear combination of yields. However, we find that estimation of a standard affine three-factor model results in a variance state variable that, while instrumental in explaining the shape of the yield curve, is essentially unrelated to GARCH estimates of the quadratic variation of the spot rate process or to implied variances from options. We then investigate four-factor affine models. Of the models tested, only the model that exhibits ‘unspanned stochastic volatility’ (USV) generates both realistic short rate volatility estimates and a good cross-sectional fit. Our findings suggest that short rate volatility cannot be extracted from the cross-section of bond prices. In particular, short rate volatility and convexity are only weakly correlated.     

An empirical comparison of GARCH option pricing models

Recent empirical studies have shown that GARCH models can be successfully used to describe option prices. Pricing such contracts requires knowledge of the risk neutral cumulative return distribution. Since the analytical forms of these distributions are generally unknown, computationally intensive numerical schemes are required for pricing to proceed. Heston and Nandi (2000) consider a particular GARCH structure that permits analytical solutions for pricing European options and they provide empirical support for their model. The analytical tractability comes at a potential cost of realism in the underlying GARCH dynamics. In particular, their model falls in the affine family, whereas most GARCH models that have been examined fall in the non-affine family. This article takes a closer look at this model with the objective of establishing whether there is a cost to restricting focus to models in the affine family. We confirm Heston and Nandi's findings, namely that their model can explain a significant portion of the volatility smile. However, we show that a simple non affine NGARCH option model is superior in removing biases from pricing residuals for all moneyness and maturity categories especially for out-the-money contracts. The implications of this finding are examined. JEL Classification G13     

Unspanned stochastic volatility and fixed income derivatives pricing

We propose a parsimonious ‘unspanned stochastic volatility’ model of the term structure and study its implications for fixed-income option prices. The drift and quadratic variation of the short rate are affine in three state variables (the short rate, its long-term mean and variance) which follow a joint Markov (vector) process. Yet, bond prices are exponential affine functions of only two state variables, independent of the current interest rate volatility level. Because this result holds for an arbitrary volatility process, such a process can be calibrated to match fixed income derivative prices. Furthermore, this model can be ‘extended’ (by relaxing the time-homogeneity) to fit any arbitrary term structure. In its ‘HJM’ form, this model nests the analogous stochastic equity volatility model of Heston (1993) [Heston, S.L., 1993. A closed form solution for options with stochastic volatility. Review of Financial Studies 6, 327–343]. In particular, if the volatility process is specified to be affine, closed-form solutions for interest rate options obtain. We propose an efficient algorithm to compute these prices. An application using data on caps and floors shows that the model can capture very well the implied Black spot volatility surface, while simultaneously fitting the observed term structure.     

Diagnosing affine models of options pricing: Evidence from VIX

Affine jump-diffusion models have been the mainstream in options pricing because of their analytical tractability. Popular affine jump-diffusion models, however, are still unsatisfactory in describing the options data and the problem is often attributed to the diffusion term of the unobserved state variables. Using prices of variance-swaps (i.e., squared VIX) implied from options prices, we provide fresh evidence regarding the misspecification of affine jump-diffusion models, as variance-swap prices are affine functions of the state variables in a broader class of models that do not restrict the diffusion term of the state variables. We apply the nonparametric methodology used by Aït-Sahalia (1996b), supplemented with bootstrap tests and other parametric tests, to the S&P 500 index options data from January 1996 to September 2008. We find that, while the affine diffusion term of the state variables may contribute to the misspecification as the literature has suggested, the affine drift of the state variables, jump intensities, and risk premiums are also sources of misspecification.     

最优多因子仿射利率期限结构模型选择

在利用NS模型估计出市场即期利率的基础上,采用卡尔曼滤波方法对多因子Vasieck和CIR模型进行参数估计,最后运用蒙特卡罗模拟方法对交易所国债价格进行模拟,并与实际价格进行比较,进而确定了符合我们国债市场的最优多因子仿射利率期限结构模型。研究结果表明：多因子CIR模型对数据的拟合效果及对国债价格模拟效果要明显优于多因子Vasicek模型;对于多因子CIR模型而言,因子个数增加并没有提高模型的价格模拟效果;两因子CIR模型具有最优的国债价格模拟效果。     

Estimating and Testing Exponential-Affine Term Structure Models by Kalman Filter

This paper proposes a unified state-space formulation for parameter estimation of exponential-affine term structure models. The proposed method uses an approximate linear Kalman filter which only requires specifying the conditional mean and variance of the system in an approximate sense. The method allows for measurement errors in the observed yields to maturity, and can simultaneously deal with many yields on bonds with different maturities. An empirical analysis of two special cases of this general class of model is carried out: the Gaussian case (Vasicek 1977) and the non-Gaussian case (Cox Ingersoll and Ross 1985 and Chen and Scott 1992). Our test results indicate a strong rejection of these two cases. A Monte Carlo study indicates that the procedure is reliable for moderate sample sizes.     

Identification of Maximal Affine Term Structure Models

Building on Duffie and Kan (1996) , we propose a new representation of affine models in which the state vector comprises infinitesimal maturity yields and their quadratic covariations. Because these variables possess unambiguous economic interpretations, they generate a representation that is globally identifiable. Further, this representation has more identifiable parameters than the “maximal” model of Dai and Singleton (2000) . We implement this new representation for select three‐factor models and find that model‐independent estimates for the state vector can be estimated directly from yield curve data, which present advantages for the estimation and interpretation of multifactor models.