Conditional sampling for barrier option pricing under the LT method

We develop a conditional sampling scheme for pricing knock-out barrier options under the Linear Transformations (LT) algorithm from Imai and Tan (2006). We compare our new method to an existing conditional Monte Carlo scheme from Glasserman and Staum (2001), and show that a substantial variance reduction is achieved. We extend the method to allow pricing knock-in barrier options and introduce a root-finding method to obtain a further variance reduction. The effectiveness of the new method is supported by numerical results

Conditional Quasi-Monte Carlo with Constrained Active Subspaces

Conditional Monte Carlo or pre-integration is a useful tool for reducing variance and improving regularity of integrands when applying Monte Carlo and quasi-Monte Carlo (QMC) methods. To choose the variable to pre-integrate with, one need to consider both the variable importance and the tractability of the conditional expectation. For integrals over a Gaussian distribution, one can pre-integrate over any linear combination of variables. Liu and Owen (2022) propose to choose the linear combination based on an active subspace decomposition of the integrand. However, pre-integrating over such selected direction might be intractable. In this work, we address this issue by finding the active subspaces subject to the constraints such that pre-integration can be easily carried out. The proposed method is applied to some examples in derivative pricing under stochastic volatility models and is shown to outperform previous methods

New Brownian bridge construction in quasi-Monte Carlo methods for computational finance

AbstractQuasi-Monte Carlo (QMC) methods have been playing an important role for high-dimensional problems in computational finance. Several techniques, such as the Brownian bridge (BB) and the principal component analysis, are often used in QMC as possible ways to improve the performance of QMC. This paper proposes a new BB construction, which enjoys some interesting properties that appear useful in QMC methods. The basic idea is to choose the new step of a Brownian path in a certain criterion such that it maximizes the variance explained by the new variable while holding all previously chosen steps fixed. It turns out that using this new construction, the first few variables are more “important” (in the sense of explained variance) than those in the ordinary BB construction, while the cost of the generation is still linear in dimension. We present empirical studies of the proposed algorithm for pricing high-dimensional Asian options and American options, and demonstrate the usefulness of the new BB

Modeling and quasi-Monte Carlo simulation of risk in credit portfolios

Credit risk is the risk of losing contractually obligated cash flows promised by a counterparty such as a corporation, financial institution, or government due to default on its debt obligations. The need for accurate pricing and hedging of complex credit derivatives and for active management of large credit portfolios calls for an accurate assessment of the risk inherent in the underlying credit portfolios. An important challenge for modeling a credit portfolio is to capture the correlations within the credit portfolio. For very large and homogeneous portfolios, analytic and semi-analytic approaches can be used to derive limiting distributions. However, for portfolios of inhomogeneous default probabilities, default correlations, recovery values, or position sizes, Monte Carlo methods are necessary to capture their underlying dynamic evolutions. Since the feasibility of the Monte Carlo methods is limited by their relatively slow convergence rate, methods to improve the efficiency of simulations for credit portfolios are highly desired. In this dissertation, a comparison of the commonly employed single step models for credit portfolios, referred to as the copula-based default time approach, with our novel applications of multi-step models was made at first. Comparison of simulation results indicates that the dependency structure may be better incorporated by the multi-step models, since the default time models can introduce substantially skewed correlations within credit portfolios, a shortcoming which has become more evident in the recent subprime crisis. Next, to improve the efficiency of simulations, quasi- random sequences were introduced into both the single step and multi-step models by devising several new algorithms involving the Brownian bridge construction and principal component analysis. The simulation results from tests under various scenarios suggest that quasi-Monte Carlo methods can substantially improve simulation effectiveness not only for the problems of computing integrals but also for those of order statistics, indicating significant advantage when calculating a number of risk quantities such as Value at Risk (VaR). Finally, the performance of the simulations based on the above credit portfolio models and the quasi-Monte Carlo methods was examined in the context of modeling and valuation of credit portfolio derivatives. The results suggest that these methods can considerably improve the simulation of complex financial instruments involving portfolio credit risk

The LIBOR Market Model

Student Number : 0003819T - MSc dissertation - School of Computational and Applied Mathematics - Faculty of ScienceThe over-the-counter (OTC) interest rate derivative market is large and rapidly developing. In March 2005, the Bank for International Settlements published its “Triennial Central Bank Survey” which examined the derivative market activity in 2004 (http://www.bis.org/publ/rpfx05.htm). The reported total gross market value of OTC derivatives stood at $6.4 trillion at the end of June 2004. The gross market value of interest rate derivatives comprised a massive 71.7% of the total, followed by foreign exchange derivatives (17.5%) and equity derivatives (5%). Further, the daily turnover in interest rate option trading increased from 5.9% (of the total daily turnover in the interest rate derivative market) in April 2001 to 16.7% in April 2004. This growth and success of the interest rate derivative market has resulted in the introduction of exotic interest rate products and the ongoing search for accurate and efficient pricing and hedging techniques for them. Interest rate caps and (European) swaptions form the largest and the most liquid part of the interest rate option market. These vanilla instruments depend only on the level of the yield curve. The market standard for pricing them is the Black (1976) model. Caps and swaptions are typically used by traders of interest rate derivatives to gamma and vega hedge complex products. Thus an important feature of an interest rate model is not only its ability to recover an arbitrary input yield curve, but also an ability to calibrate to the implied at-the-money cap and swaption volatilities. The LIBOR market model developed out of the market’s need to price and hedge exotic interest rate derivatives consistently with the Black (1976) caplet formula. The focus of this dissertation is this popular class of interest rate models. The fundamental traded assets in an interest rate model are zero-coupon bonds. The evolution of their values, assuming that the underlying movements are continuous, is driven by a finite number of Brownian motions. The traditional approach to modelling the term structure of interest rates is to postulate the evolution of the instantaneous short or forward rates. Contrastingly, in the LIBOR market model, the discrete forward rates are modelled directly. The additional assumption imposed is that the volatility function of the discrete forward rates is a deterministic function of time. In Chapter 2 we provide a brief overview of the history of interest rate modelling which led to the LIBOR market model. The general theory of derivative pricing is presented, followed by a exposition and derivation of the stochastic differential equations governing the forward LIBOR rates. The LIBOR market model framework only truly becomes a model once the volatility functions of the discrete forward rates are specified. The information provided by the yield curve, the cap and the swaption markets does not imply a unique form for these functions. In Chapter 3, we examine various specifications of the LIBOR market model. Once the model is specified, it is calibrated to the above mentioned market data. An advantage of the LIBOR market model is the ability to calibrate to a large set of liquid market instruments while generating a realistic evolution of the forward rate volatility structure (Piterbarg 2004). We examine some of the practical problems that arise when calibrating the market model and present an example calibration in the UK market. The necessity, in general, of pricing derivatives in the LIBOR market model using Monte Carlo simulation is explained in Chapter 4. Both the Monte Carlo and quasi-Monte Carlo simulation approaches are presented, together with an examination of the various discretizations of the forward rate stochastic differential equations. The chapter concludes with some numerical results comparing the performance of Monte Carlo estimates with quasi-Monte Carlo estimates and the performance of the discretization approaches. In the final chapter we discuss numerical techniques based on Monte Carlo simulation for pricing American derivatives. We present the primal and dual American option pricing problem formulations, followed by an overview of the two main numerical techniques for pricing American options using Monte Carlo simulation. Callable LIBOR exotics is a name given to a class of interest rate derivatives that have early exercise provisions (Bermudan style) to exercise into various underlying interest rate products. A popular approach for valuing these instruments in the LIBOR market model is to estimate the continuation value of the option using parametric regression and, subsequently, to estimate the option value using backward induction. This approach relies on the choice of relevant, i.e. problem specific predictor variables and also on the functional form of the regression function. It is certainly not a “black-box” type of approach. Instead of choosing the relevant predictor variables, we present the sliced inverse regression technique. Sliced inverse regression is a statistical technique that aims to capture the main features of the data with a few low-dimensional projections. In particular, we use the sliced inverse regression technique to identify the low-dimensional projections of the forward LIBOR rates and then we estimate the continuation value of the option using nonparametric regression techniques. The results for a Bermudan swaption in a two-factor LIBOR market model are compared to those in Andersen (2000)