Additive and multiplicative duals for American option pricing

We investigate and compare two dual formulations of the American option pricing problem based on two decompositions of supermartingales: the additive dual of Haugh and Kogan (Oper. Res. 52:258–270, 2004) and Rogers (Math. Finance 12:271–286, 2002) and the multiplicative dual of Jamshidian (Minimax optimality of Bermudan and American claims and their Monte- Carlo upper bound approximation. NIB Capital, The Hague, 2003). Both provide upper bounds on American option prices; we show how to improve these bounds iteratively and use this to show that any multiplicative dual can be improved by an additive dual and vice versa. This iterative improvement converges to the optimal value function. We also compare bias and variance under the two dual formulations as the time horizon grows; either method may have smaller bias, but the variance of the multiplicative method typically grows much faster than that of the additive method. We show that in the case of a discrete state space, the additive dual coincides with the dual of the optimal stopping problem in the sense of linear programming duality and the multiplicative method arises through a nonlinear duality.      

Interacting particle systems for the computation of rare credit portfolio losses

In this paper, we introduce the use of interacting particle systems in the computation of probabilities of simultaneous defaults in large credit portfolios. The method can be applied to compute small historical as well as risk-neutral probabilities. It only requires that the model be based on a background Markov chain for which a simulation algorithm is available. We use the strategy developed by Del Moral and Garnier in (Ann. Appl. Probab. 15:2496–2534, 2005) for the estimation of random walk rare events probabilities. For the purpose of illustration, we consider a discrete-time version of a first passage model for default. We use a structural model with stochastic volatility, and we demonstrate the efficiency of our method in situations where importance sampling is not possible or numerically unstable.      

Credit default swap calibration and derivatives pricing with the SSRD stochastic intensity model

We introduce the two-dimensional shifted square-root diffusion (SSRD) model for interest-rate and credit derivatives with (positive) stochastic intensity. The SSRD is the unique explicit diffusion model allowing an automatic and separated calibration of the term structure of interest rates and of credit default swaps (CDSs), and retaining free dynamics parameters that can be used to calibrate option data. We propose a new positivity preserving implicit Euler scheme for Monte Carlo simulation. We discuss the impact of interest-rate and default-intensity correlation and develop an analytical approximation to price some basic credit derivatives terms involving correlated CIR processes. We hint at a formula for CDS options under CIR + + CDS-calibrated stochastic intensity.Received: March 2004, Mathematics Subject Classification (2000): 60H10, 60J60, 60J75, 91B70JEL Classification: G13     

Simulation-based Value-at-Risk for nonlinear portfolios

Value-at-risk (VaR) has been playing the role of a standard risk measure since its introduction. In practice, the delta-normal approach is usually adopted to approximate the VaR of portfolios with option positions. Its effectiveness, however, substantially diminishes when the portfolios concerned involve a high dimension of derivative positions with nonlinear payoffs; lack of closed form pricing solution for these potentially highly correlated, American-style derivatives further complicate the problem. This paper proposes a generic simulation-based algorithm for VaR estimation that can be easily applied to any existing procedures. Our proposal leverages cross-sectional information and applies variable selection techniques to simplify the existing simulation framework. Asymptotic properties of the new approach demonstrate faster convergence due to the additional model selection component introduced. We have also performed sets of numerical results that verify the effectiveness of our approach in comparison with some existing strategies.     

证券衍生品种法制构建与《证券法》修改

本文提出，《证券法》修改时应明确证券衍生品种的法律地位，并应考虑证券衍生品种的特殊性，相关发行、上市、交易、产品认定等环节的具体制度设计应为证券衍生品种的发展留下应有的法律空间。本文认为，证券衍生品种交易场所定位于证券交易所是合理的。     

A class of Lévy process models with almost exact calibration to both barrier and vanilla FX options

Vanilla (standard European) options are actively traded on many underlying asset classes, such as equities, commodities and foreign exchange (FX). The market quotes for these options are typically used by exotic options traders to calibrate the parameters of the (risk-neutral) stochastic process for the underlying asset. Barrier options, of many different types, are also widely traded in all these markets but one important feature of the FX options markets is that barrier options, especially double-no-touch (DNT) options, are now so actively traded that they are no longer considered, in any way, exotic options. Instead, traders would, in principle, like to use them as instruments to which they can calibrate their model. The desirability of doing this has been highlighted by talks at practitioner conferences but, to our best knowledge (at least within the realm of the published literature), there have been no models which are specifically designed to cater for this. In this paper, we introduce such a model. It allows for calibration in a two-stage process. The first stage fits to DNT options (or other types of double barrier options). The second stage fits to vanilla options. The key to this is to assume that the dynamics of the spot FX rate are of one type before the first exit time from a ‘corridor’ region but are allowed to be of a different type after the first exit time. The model allows for jumps (either finite activity or infinite activity) and also for stochastic volatility. Hence, not only can it give a good fit to the market prices of options, it can also allow for realistic dynamics of the underlying FX rate and realistic future volatility smiles and skews. En route, we significantly extend existing results in the literature by providing closed-form (up to Laplace inversion) expressions for the prices of several types of barrier options as well as results related to the distribution of first passage times and of the ‘overshoot’.     

Smoothing the payoff for efficient computation of Basket option prices

We consider the problem of pricing basket options in a multivariate Black–Scholes or Variance-Gamma model. From a numerical point of view, pricing such options corresponds to moderate and high-dimensional numerical integration problems with non-smooth integrands. Due to this lack of regularity, higher order numerical integration techniques may not be directly available, requiring the use of methods like Monte Carlo specifically designed to work for non-regular problems. We propose to use the inherent smoothing property of the density of the underlying in the above models to mollify the payoff function by means of an exact conditional expectation. The resulting conditional expectation is unbiased and yields a smooth integrand, which is amenable to the efficient use of adaptive sparse-grid cubature. Numerical examples indicate that the high-order method may perform orders of magnitude faster than Monte Carlo or Quasi Monte Carlo methods in dimensions up to 35.     

Sensitivity estimates for portfolio credit derivatives using Monte Carlo

Portfolio credit derivatives are contracts that are tied to an underlying portfolio of defaultable reference assets and have payoffs that depend on the default times of these assets. The hedging of credit derivatives involves the calculation of the sensitivity of the contract value with respect to changes in the credit spreads of the underlying assets, or, more generally, with respect to parameters of the default-time distributions. We derive and analyze Monte Carlo estimators of these sensitivities. The payoff of a credit derivative is often discontinuous in the underlying default times, and this complicates the accurate estimation of sensitivities. Discontinuities introduced by changes in one default time can be smoothed by taking conditional expectations given all other default times. We use this to derive estimators and to give conditions under which they are unbiased. We also give conditions under which an alternative likelihood ratio method estimator is unbiased. We illustrate the application and verification of these conditions and estimators in the particular case of the multifactor Gaussian copula model, but the methods are more generally applicable.      

Optimal importance sampling with explicit formulas in continuous time

In the Black–Scholes model, consider the problem of selecting a change of drift which minimizes the variance of Monte Carlo estimators for prices of path-dependent options. Employing large deviations techniques, the asymptotically optimal change of drift is identified as the solution to a one-dimensional variational problem, which may be reduced to the associated Euler–Lagrange differential equation. Closed-form solutions for geometric and arithmetic average Asian options are provided. The authors acknowledge the support of the National Science Foundation under grants DMS-0532390 (Guasoni) and DGE-0221680 (Robertson) at Boston University.     

蒙特卡罗方差减小技术以及在金融中的应用

本文主要介绍MonteCarlo方差减小技术及其在金融中衍生证券定价的一些的应用。对近几年来国际上取得的主要研究成果作了简要的介绍和比较。并提出一些目前待解决的问题。     

Primal–dual quasi-Monte Carlo simulation with dimension reduction for pricing American options

The pricing of American options is one of the most challenging problems in financial engineering due to the involved optimal stopping time problem, which can be solved by using dynamic programming (DP). But applying DP is not always practical, especially when the state space is high dimensional. However, the curse of dimensionality can be overcome by Monte Carlo (MC) simulation. We can get lower and upper bounds by MC to ensure that the true price falls into a valid confidence interval. During the recent decades, progress has been made in using MC simulation to obtain both the lower bound by least-squares Monte Carlo method (LSM) and the upper bound by duality approach. However, there are few works on pricing American options using quasi-Monte Carlo (QMC) methods, especially to compute the upper bound. For comparing the sample variances and standard errors in the numerical experiments, randomized QMC (RQMC) methods are usually used. In this paper, we propose to use RQMC to replace MC simulation to compute both the lower bound (by the LSM) and the upper bound (by the duality approach). Moreover, we propose to use dimension reduction techniques, such as the Brownian bridge, principal component analysis, linear transformation and the gradients based principle component analysis. We perform numerical experiments on American–Asian options and American max-call options under the Black–Scholes model and the variance gamma model, in which the options have the path-dependent feature or are written on multiple underlying assets. We find that RQMC in combination with dimension reduction techniques can significantly increase the efficiency in computing both the lower and upper bounds, resulting in better estimates and tighter confidence intervals of the true price than pure MC simulation.