EQUILIBRIUM ASSET AND OPTION PRICING UNDER JUMP DIFFUSION

This paper develops an equilibrium asset and option pricing model in a production economy under jump diffusion. The model provides analytical formulas for an equity premium and a more general pricing kernel that links the physical and risk‐neutral densities. The model explains the two empirical phenomena of the negative variance risk premium and implied volatility smirk if market crashes are expected. Model estimation with the S&P 500 index from 1985 to 2005 shows that jump size is indeed negative and the risk aversion coefficient has a reasonable value when taking the jump into account.     

A STATE-SPACE PARTITIONING METHOD FOR PRICING HIGH-DIMENSIONAL AMERICAN-STYLE OPTIONS

The pricing of American-style options by simulation-based methods is an important but difficult task primarily due to the feature of early exercise, particularly for high-dimensional derivatives. In this paper, a bundling method based on quasi-Monte Carlo sequences is proposed to price high-dimensional American-style options. The proposed method substantially extends Tilley's bundling algorithm to higher-dimensional situations. By using low-discrepancy points, this approach partitions the state space and forms bundles. A dynamic programming algorithm is then applied to the bundles to estimate the continuation value of an American-style option. A convergence proof of the algorithm is provided. A variety of examples with up to 15 dimensions are investigated numerically and the algorithm is able to produce computationally efficient results with good accuracy.     

ON THE ERROR IN THE MONTE CARLO PRICING OF SOME FAMILIAR EUROPEAN PATH-DEPENDENT OPTIONS

This paper studies the relative error in the crude Monte Carlo pricing of some familiar European path-dependent multiasset options. For the crude Monte Carlo method it is well known that the convergence rate   O ( n ^−1/2)  , where n is the number of simulations, is independent of the dimension of the integral. This paper also shows that for a large class of pricing problems in the multiasset Black-Scholes market the constant in   O ( n ^−1/2)  is independent of the dimension. To be more specific, the constant is only dependent on the highest volatility among the underlying assets, time to maturity, and degree of confidence interval.     

THE DOTHAN PRICING MODEL REVISITED

We compute zero‐coupon bond prices in the Dothan model by solving the associated PDE using integral representations of heat kernels and Hartman–Watson distributions. We obtain several integral formulas for the price P(t, T) at time t > 0 of a bond with maturity T > 0 that complete those of the original paper of Dothan, which are shown not to always satisfy the boundary condition P(T, T) = 1 .     

MODEL UNCERTAINTY AND ITS IMPACT ON THE PRICING OF DERIVATIVE INSTRUMENTS

Uncertainty on the choice of an option pricing model can lead to "model risk" in the valuation of portfolios of options. After discussing some properties which a quantitative measure of model uncertainty should verify in order to be useful and relevant in the context of risk management of derivative instruments, we introduce a quantitative framework for measuring model uncertainty in the context of derivative pricing. Two methods are proposed: the first method is based on a coherent risk measure compatible with market prices of derivatives, while the second method is based on a convex risk measure. Our measures of model risk lead to a premium for model uncertainty which is comparable to other risk measures and compatible with observations of market prices of a set of benchmark derivatives. Finally, we discuss some implications for the management of "model risk."     

MCMC ESTIMATION OF LÉVY JUMP MODELS USING STOCK AND OPTION PRICES

We examine the performances of several popular Lévy jump models and some of the most sophisticated affine jump‐diffusion models in capturing the joint dynamics of stock and option prices. We develop efficient Markov chain Monte Carlo methods for estimating parameters and latent volatility/jump variables of the Lévy jump models using stock and option prices. We show that models with infinite‐activity Lévy jumps in returns significantly outperform affine jump‐diffusion models with compound Poisson jumps in returns and volatility in capturing both the physical and risk‐neutral dynamics of the S&P 500 index. We also find that the variance gamma model of Madan, Carr, and Chang with stochastic volatility has the best performance among all the models we consider.     

A STRUCTURAL RISK‐NEUTRAL MODEL FOR PRICING AND HEDGING POWER DERIVATIVES

We develop a structural risk‐neutral model for energy market modifying along several directions the approach introduced in Aïd et al. In particular, a scarcity function is introduced to allow important deviations of the spot price from the marginal fuel price, producing price spikes. We focus on pricing and hedging electricity derivatives. The hedging instruments are forward contracts on fuels and electricity. The presence of production capacities and electricity demand makes such a market incomplete. We follow a local risk minimization approach to price and hedge energy derivatives. Despite the richness of information included in the spot model, we obtain closed‐form formulae for futures prices and semiexplicit formulae for spread options and European options on electricity forward contracts. An analysis of the electricity price risk premium is provided showing the contribution of demand and capacity to the futures prices. We show that when far from delivery, electricity futures behave like a basket of futures on fuels.     

FAST MONTE CARLO GREEKS FOR FINANCIAL PRODUCTS WITH DISCONTINUOUS PAY‐OFFS

We introduce a new class of numerical schemes for discretizing processes driven by Brownian motions. These allow the rapid computation of sensitivities of discontinuous integrals using pathwise methods even when the underlying densities postdiscretization are singular. The two new methods presented in this paper allow Greeks for financial products with trigger features to be computed in the LIBOR market model with similar speed to that obtained by using the adjoint method for continuous pay‐offs. The methods are generic with the main constraint being that the discontinuities at each step must be determined by a one‐dimensional function: the proxy constraint. They are also generic with the sole interaction between the integrand and the scheme being the specification of this constraint.