Refine
Language
- English (34)
Keywords
- Monte Carlo simulation (3)
- optimal stopping (2)
- American options (1)
- Arbitrage (1)
- BSDE (1)
- Bermudan options (1)
- Finance (1)
- Financial derivatives (1)
- Forward backward SDE (1)
- LIBOR market model (1)
Project
- E5 (34) (remove)
Application Area
- E (34) (remove)
The CreditRisk model launched by CSFB in 1997 is widely used by practitioners in the banking sector as a simple means for the quantification of credit
risk, primarily of the loan book. We present an alternative numerical recursion scheme for CreditRisk, equivalent to an algorithm recently proposed by
Giese, based on well-known expansions of the logarithm and the exponential
of a power series. We show that it is advantageous to the Panjer recursion
advocated in the original CreditRisk
document, in that it is numerically stable. The crucial stability arguments are explained in detail. Furthermore, the
computational complexity of the resulting algorithm is stated.
We introduce a new Monte Carlo method for constructing the exercise
boundary of an American option in a generalized Black-Scholes framework.
Based on a known exercise boundary, it is shown how to price and hedge the
American option by Monte Carlo simulation of suitable probabilistic represen-
tations in connection with the respective parabolic boundary value problem.
The methods presented are supported by numerical simulation experiments.
We present a new iterative procedure for solving the multiple stopping
problem in discrete time and discuss the stability of the algorithm.
The algorithm produces monotonically increasing approximations of the
Snell envelope, which coincide with the Snell envelope after finitely many
steps. Contrary to backward dynamic programming, the algorithm allows
to calculate approximative solutions with only a few nestings of conditional
expectations and is, therefore, tailor-made for a plain Monte-Carlo
implementation.
In this project we propose the use of some widespread prediction techniques in the last few years for modeling derivatives. In order to do that, we have reviewed the state-of-the-art of the prediction models dealing with stochastic processes. In the oil futures sector, Schwartz suggested a model in which the oil futures price was split in two factors: the long-term equilibrium price and the short-term variations. As a result, we propose a Hull-White discrete-time two-factor interest rate model, whose factors are the short and the long term.
In Kolodko & Schoenmakers (2004) and Bender & Schoenmakers (2004) a policy iteration was introduced which allows to achieve tight lower approximations of the price for early exercise options via a nested Monte-Carlo simulation in a Markovian setting. In this paper we enhance the algorithm by a scenario selection method. It is demonstrated by numerical examples that the scenario selection can significantly reduce the number of actually performed inner simulations, and thus can heavily speed up the method (up to factor 10 in some examples). Moreover, it is shown that the modified algorithm retains the desirable properties of the original one such as the monotone improvement property, termination after a finite number of iteration steps, and numerical stability.
We propose a valuation method for callable structures in a multi-factor Libor model which are path-dependent in the sense that, after calling, one receives a sequence of cash-flows in the future, instead of a well specified cash-flow at the calling date. The method is based on a Monte Carlo procedure for standard Bermudans recently developed in Kolodko & Schoenmakers (2004), and is applied to the cancelable snowball interest rate swap. The proposed procedure is quite generic, straightforward to implement, and can be easily adapted to other related path-dependent products.
We apply theoretical results of S. Peng on supersolutions for BS-DEs
to the problem of finding optimal superhedging strategies in a
Black-Scholes market under constraints. Constraints may be imposed
simultaneously on wealth process and portfolio. They may be nonconvex,
time-dependent, and random. Constraints on the portfolio may
e.g. be formulated in terms of the amount of money invested, the portfolio
proportion, or the number of shares held.
We introduce a forward scheme to simulate backward SDEs. Compared
to existing schemes, we avoid high order nestings of conditional
expectations backwards in time. In this way the error, when
approximating the conditional expectation, in dependence of the
time partition is significantly reduced. Besides this generic
result, we present an implementable algorithm and provide an error
analysis for it. Finally, we demonstrate the strength of the new
algorithm by solving some financial problems numerically.
We show that pricing a big class of relevant options by hedging
and no-arbitrage can be extended beyond semimartingale models. To
this end we construct a subclass of self-financing portfolios that
contains hedges for these options, but does not contain arbitrage
opportunities, even if the stock price process is a
non-semimartingale of some special type.
Moreover, we show that the option prices depend
essentially only on a path property of the stock price process,
viz. on the quadratic variation. As a consequence, we can
incorporate many stylized facts to a pricing model without
changing the option prices.
In this paper we lay the foundation for a numerical algorithm to
simulate high-dimensional coupled FBSDEs under weak coupling or
monotonicity conditions. In particular we prove convergence of a
time discretization and a Markovian iteration. The iteration
differs from standard Picard iterations for FBSDEs in that the
dimension of the underlying Markovian process does not increase
with the number of iterations. This feature seems to be
indispensable for an efficient iterative scheme from a numerical
point of view. We finally suggest a fully explicit numerical
algorithm and present some numerical examples with up to
10-dimensional state space.
In a rather general setting of multivariate stochastic volatility market models we derive global iterative probabilistic schemes for computing the free boundary and its Greeks for a generic class of American derivative models using front-fixing methods. Establishment of convergence is closely linked to a proof of global regularity of the free boundary surface.
In this paper we carry over the concept of reverse probabilistic representa-
tions developed in Milstein, Schoenmakers, Spokoiny (2004) for diffusion pro-
cesses, to discrete time Markov chains. We outline the construction of reverse
chains in several situations and apply this to processes which are connected
with jump-diffusion models and finite state Markov chains. By combining
forward and reverse representations we then construct transition density esti-
mators for chains which have root-N accuracy in any dimension and consider
some applications.
In this paper we introduce efficient Monte Carlo estimators for the valuation
of high-dimensional derivatives and their sensitivities (”Greeks”).
These estimators are based on an analytical, usually approximative representation
of the underlying density. We study approximative densities
obtained by the WKB method. The results are applied in the context of
a Libor market model.
We consider regular polynomial interpolation algorithms on recursively defined sets of interpolation points which approximate global solutions of arbitrary well-posed systems of linear partial differential equations. Convergence of the "limit" of the recursively constructed family of polynomials to the solution and error estimates are obtained from a priori estimates for some standard classes of linear partial differential equations, i.e. elliptic and hyperbolic equations. Another variation of the algorithm allows to construct polynomial interpolations which preserve systems of linear partial differential equations at the interpolation points. We show how this can be applied in order to compute higher order terms of WKB-approximations of fundamental solutions of a large class of linear parabolic equations. The error estimates are sensitive to the regularity of the solution. Our method is compatible with recent developments for solution of higher dimensional partial differential equations, i.e. (adaptive) sparse grids, and weighted Monte-Carlo, and has obvious applications to mathematical finance and physics.
In this paper we propose a Libor model with a high-dimensional specially structured system of
driving CIR volatility processes. A stable calibration procedure which takes into account
a given local correlation structure is presented. The calibration algorithm is FFT based, so fast and easy
to implement.
We present a generic non-nested Monte Carlo procedure for computing true upper bounds for Bermudan products, given an approximation of the Snell envelope. The pleonastic ``true'' stresses that, by construction, the estimator is biased above the Snell envelope. The key idea is a regression estimator for the Doob martingale part of the approximative Snell envelope, which preserves the martingale property. The so constructed martingale may be employed for computing dual upper bounds without nested simulation. In general, this martingale can also be used as a control variate for simulation of conditional expectations. In this context, we develop a variance reduced version of the nested primal-dual estimator (Anderson & Broadie (2004)) and nested consumption based (Belomestny & Milstein (2006)) methods . Numerical experiments indicate the efficiency of the non-nested Monte Carlo algorithm and the variance reduced nested one.
We present two approximation methods for pricing of CMS spread options in Libor market models. Both approaches are based on approximating the underlying swap rates with lognormal processes under suitable measures. The first method is derived straightforwardly from the Libor market model. The second one uses a convexity adjustment technique under a linear swap model assumption. A numerical study demonstrates that both methods provide satisfactory approximations of spread option prices and can be used for calibration of a Libor market model to the CMS spread option market.
The Real Multiple Dual
(2009)
In this paper we present a dual representation for the multiple stopping
problem, hence multiple exercise options. As such it is a natural generalization of the
method in Rogers (2002) and Haugh and Kogan (2004) for the standard stopping
problem for American options. We consider this representation as the real dual as it is
solely expressed in terms of an infimum over martingales rather than an infimum over
martingales and stopping times as in Meinshausen and Hambly (2004). For the multiple
dual representation we present three Monte Carlo simulation algorithms which require
only one degree of nesting.