On the duality principle in option pricing: semimartingale setting

The purpose of this paper is to describe the appropriate mathematical framework for the study of the duality principle in option pricing. We consider models where prices evolve as general exponential semimartingales and provide a complete characterization of the dual process under the dual measure. Particular cases of these models are the ones driven by Brownian motions and by Lévy processes, which have been considered in several papers. Generally speaking, the duality principle states that the calculation of the price of a call option for a model with price process S=e ^H (with respect to the measure P) is equivalent to the calculation of the price of a put option for a suitable dual model S′=e ^H′ (with respect to the dual measure P′). More sophisticated duality results are derived for a broad spectrum of exotic options. The second named author acknowledges the financial support from the Deutsche Forschungsgemeinschaft (DFG, Eb 66/9-2). This research was carried out while the third named author was supported by the Alexander von Humboldt foundation.     

On the valuation of reverse mortgage insurance

This article presents a closed-form formula for calculating the loan-to-value (LTV) ratio in an adjusted-rate reverse mortgage (RM) with a lump sum payment. Previous literatures consider the pricing of RM in a constant interest rate assumption and price it on fixed-rate loans. This paper successfully considers the dynamic of interest rate and the adjustable-rate RM simultaneously. This paper also considers the housing price shock into the valuation model. Assuming that house prices follow a jump diffusion process with a stochastic interest rate and that the loan interest rate is adjusted instantaneously according to the short rate, we demonstrate that the LTV ratio is independent of the term structure of interest rates. This argument holds even when housing prices follow a general process: an exponential Lévy process. In addition, the HECM (Home Equity Conversion Mortgage) program may be not sustainable, especially for a higher level of housing price volatility. Finally, when the loan interest rate is adjusted periodically according to the LIBOR rate, our finding reveals that the LTV ratio is insensitive to the parameters characterizing the CIR model.     

Hedging derivative securities with genetic programming

One of the most recent applications of GP to finance is to use genetic programming to derive option pricing formulas. Earlier studies take the Black–Scholes model as the true model and use the artificial data generated by it to train and to test GP. The aim of this paper is to provide some initial evidence of the empirical relevance of GP to option pricing. By using the real data from S&P 500 index options, we train and test our GP by distinguishing the case in-the-money from the case out-of-the-money. Unlike most empirical studies, we do not evaluate the performance of GP in terms of its pricing accuracy. Instead, the derived GP tree is compared with the Black–Scholes model in its capability to hedge. To do so, a notion of tracking error is taken as the performance measure. Based on the post-sample performance, it is found that in approximately 20% of the 97 test paths GP has a lower tracking error than the Black–Scholes formula. We further compare our result with the ones obtained by radial basis functions and multilayer perceptrons and one-stage GP. Copyright © 1999 John Wiley & Sons, Ltd.     

Monte Carlo Option Pricing for Tempered Stable (CGMY) Processes

Lévy processes are popular models for stock price behavior since they allow to take into account jump risk and reproduce the implied volatility smile. In this paper, we focus on the tempered stable (also known as CGMY) processes, which form a flexible 6-parameter family of Lévy processes with infinite jump intensity. It is shown that under an appropriate equivalent probability measure a tempered stable process becomes a stable process whose increments can be simulated exactly. This provides a fast Monte Carlo algorithm for computing the expectation of any functional of tempered stable process. We use our method to price European options and compare the results to a recent approximate simulation method for tempered stable process by Madan and Yor (CGMY and Meixner Subordinators are absolutely continuous with respect to one sided stable subordinators, 2005).     

Multivariate asset models using Lévy processes and applications

In this paper, we propose a multivariate asset model based on Lévy processes for pricing of products written on more than one underlying asset. Our construction is based on a two-factor representation of the dynamics of the asset log-returns. We investigate the properties of the model and introduce a multivariate generalization of some processes which are quite common in financial applications, such as subordinated Brownian motions, jump-diffusion processes and time-changed Lévy processes. Finally, we explore the issue of model calibration for the proposed setting and illustrate its robustness on a number of numerical examples.     

Computing the Gerber–Shiu function by frame duality projection

Inspired by some works of Kirkby, J. L. [2015. Efficient option pricing by frame duality with the fast Fourier transform. SIAM Journal on Financial Mathematics 6(1), 713–747; 2016. An efficient transform method for Asian option pricing. SIAM Journal on Financial Mathematics 7(1), 845–892], we present a systematic study on effectively computing the Gerber–Shiu function in the Lévy risk model, where the frame duality projection is used for approximation. By introducing an auxiliary function, we provide a smooth extension of the Gerber–Shiu function, which has closed-form Fourier transform and is differentiable over the whole real line under some conditions. The objective function is approximated by its frame duality projection onto a Riesz basis, and the projection coefficients are readily computed by the fast Fourier transform algorithm. Error analysis is made and the effectiveness of our results will be further illustrated in the numerical experiments.     

Multivariate models for operational risk

Böcker and Klüppelberg [Risk Mag., 2005, December, 90–93] presented a simple approximation of OpVaR of a single operational risk cell. The present paper derives approximations of similar quality and simplicity for the multivariate problem. Our approach is based on the modelling of the dependence structure of different cells via the new concept of a Lévy copula.     

COS method for option pricing under a regime-switching model with time-changed Lévy processes

We extend the regime-switching model to the rich class of time-changed Lévy processes and use the Fourier cosine expansion (COS) method to price several options under the resulting models. The extension of the COS method to price under the regime-switching model is not straightforward because it requires the evaluation of the characteristic function which is based on a matrix exponentiation which is not an easy task. For a two-state economy, we give an analytical expression for computing this matrix exponential, and for more than two states, we use the Carathéodory–Fejér approximation to find the option prices efficiently. In the new framework developed here, it is possible to allow switches not only in the model parameters as is commonly done in literature, but we can also completely switch among various popular financial models under different regimes without any additional computational cost. Calibration of the different regime-switching models with real market data shows that the best models are the regime-switching time-changed Lévy models. As expected by the error analysis, the COS method converges exponentially and thus outperforms all other numerical methods that have been proposed so far.     

Time-changed Lévy LIBOR market model: Pricing and joint estimation of the cap surface and swaption cube

We propose a novel time-changed Lévy LIBOR (London Interbank Offered Rate) market model for jointly pricing of caps and swaptions. The time changes are split into three components. The first component allows matching the volatility term structure, the second generates stochastic volatility, and the third accommodates for stochastic skew. The parsimonious model is flexible enough to accommodate the behavior of both caps and swaptions. For the joint estimation we use a comprehensive data set spanning the financial crisis of 2007–2010. We find that, even during this period, neither market is as fragmented as suggested by the previous literature.     

What is the Natural Scale for a Lévy Process in Modelling Term Structure of Interest Rates?

This paper gives examples of explicit arbitrage-free term structure models with Lévy jumps via the state price density approach. By generalizing quadratic Gaussian models, it is found that the probability density function of a Lévy process is a “natural” scale for the process to be the state variable of a market.      

VaR limits for pension funds: an evaluation

In a series of papers during the last ten years an interest rate theory with models which are driven by Lévy or more general processes has been developed. In this paper we derive explicit formulas for the correlations of interest rates as well as zero coupon bonds with different maturities. The models considered in this general setting are the forward rate (HJM), the forward process and the LIBOR model as well as the multicurrency extension of the latter. Specific subclasses of the class of generalized hyperbolic Lévy motions are studied as driving processes. Based on a data set of parametrized yield curves derived from German government bond prices we estimate correlations. In a second step the empirical correlations are used to calibrate the Lévy forward rate model. The superior performance of the Lévy driven models becomes obvious from the graphs.     

Equilibrium asset pricing: with non-Gaussian factors and exponential utilities

We analyse the equilibrium asset pricing implications for an economy with single period return exposures to explicit non-Gaussian systematic factors, that may be both skewed and long-tailed, and Gaussian idiosyncratic components. Investors maximize expected exponential utility and equilibrium factor prices are shown to reflect exponentially tilted prices for non-Gaussian factor risk exposures. It is shown that these prices may be directly estimated from the univariate probability law of the factor exposure, given an estimate of average risk aversion in the economy. In addition, a residual form of the capital asset pricing model continues to hold and prices the idiosyncratic or Gaussian risks. The theory is illustrated on data for the US economy using independent components analysis to identify the factors and the variance gamma model to describe the probability law of the non-Gaussian factors. It is shown that the residual CAPM accounts for no more than 1% of the pricing of risky assets, while the exponentially tilted systematic factor risk exposures account for the bulk of risky asset pricing.     

Optimal constrained investment in the Cramer-Lundberg model

We consider an insurance company whose surplus is represented by the classical Cramer-Lundberg process. The company can invest its surplus in a risk-free asset and in a risky asset, governed by the Black-Scholes equation. There is a constraint that the insurance company can only invest in the risky asset at a limited leveraging level; more precisely, when purchasing, the ratio of the investment amount in the risky asset to the surplus level is no more than a; and when short-selling, the proportion of the proceeds from the short-selling to the surplus level is no more than b. The objective is to find an optimal investment policy that minimizes the probability of ruin. The minimal ruin probability as a function of the initial surplus is characterized by a classical solution to the corresponding Hamilton-Jacobi-Bellman (HJB) equation. We study the optimal control policy and its properties. The interrelation between the parameters of the model plays a crucial role in the qualitative behavior of the optimal policy. For example, for some ratios between a and b, quite unusual and at first ostensibly counterintuitive policies may appear, like short-selling a stock with a higher rate of return to earn lower interest, or borrowing at a higher rate to invest in a stock with lower rate of return. This is in sharp contrast with the unrestricted case, first studied in Hipp and Plum, or with the case of no short-selling and no borrowing studied in Azcue and Muler.     

Application of Heston’s Model to the Chinese Stock Market

This article applies Heston’s (1993) stochastic volatility model to the Chinese stock market indices and subsequently assesses its pricing performance. A two-step estimation procedure is adopted to calibrate Heston’s model. First, we find that the option price is affected by both the moneyness and the maturity. Second, Heston’s model is more likely to overprice options, whereas the BS model tends to underestimate options. Finally, Heston’s model, by employing volatility as a random process, significantly improves the pricing accuracy compared to the BS model. Therefore, Heston’s model is tractable to analyze the Chinese stock market indices, and there is volatility risk that must not be overlooked in the Chinese stock market.     

Lévy processes driven by stochastic volatility

In this paper we extend option pricing under Lévy dynamics, by assuming that the volatility of the Lévy process is stochastic. We, therefore, develop the analog of the standard stochastic volatility models, when the underlying process is not a standard (unit variance) Brownian motion, but rather a standardized Lévy process. We present a methodology that allows one to compute option prices, under virtually any set of diffusive dynamics for the parameters of the volatility process. First, we use ‘local consistency’ arguments to approximate the volatility process with a finite, but sufficiently dense Markov chain; we then use this regime switching approximation to efficiently compute option prices using Fourier inversion. A detailed example, based on a generalization of the popular stochastic volatility model of Heston (Rev Financial Stud 6 (1993) 327), is used to illustrate the implementation of the algorithms. Computer code is available at www.theponytail.net/