19 research outputs found
Recommended from our members
High-order compact finite difference schemes for option pricing in stochastic volatility jump-diffusion models
This thesis focusses on the derivation and implementation of high-order compact finite difference schemes to price a variety of options under various stochastic volatility and jump models, with the inclusion of further studies relating to the derivatives of options and the practice of hedging.
First, we derive a implicit-explicit high-order compact finite difference scheme for pricing European options under the Bates model. The resulting scheme is fourth order accurate in space and second order accurate in time. In the numerical study this scheme is compared to both a second order finite difference scheme and high-order finite element methods, where it outperforms both in terms of convergence, computational speed and required memory allocation. A numerical stability study is conducted which indicates unconditional stability of the scheme.
Second, we introduce the practice of hedging and give examples of hedging strategies created from a combination of option payoffs, we show the important role the derivatives of the option price play in forming profitable strategies. We go on to complete a study of the convergence of derivatives of the option price, the so-called Greeks. We conduct studies into Delta, vega and gamma hedging, where the derived high-order compact scheme outperforms a second-order finite difference method. Examples are provided to display how this increase in computational efficiency may assist financial practictoners.
Third, we extend the high-order compact scheme to price European options under the stochastic volatility with comtemporaneous jumps model. The derived scheme is fourth order accurate in space and second order accurate in time. We conduct numerical studies to test the new high-order compact schemes convergence, computational speed and required memory allocation against a second-order finite difference scheme, where the results show improvements in convergence at the expense of computational time. Further studies of numerical stability indicate unconditional stability of the high-order compact scheme
Generalized averaged Gaussian quadrature and applications
A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal
MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications
Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described
Properties and advances of probabilistic and statistical algorithms with applications in finance
This thesis is concerned with the construction and enhancement of algorithms involving probability
and statistics. The main motivation for these are problems that appear in finance and
more generally in applied science. We consider three distinct areas, namely, credit risk modelling,
numerics for McKean Vlasov stochastic differential equations and stochastic representations of
Partial Differential Equations (PDEs), therefore the thesis is split into three parts.
Firstly, we consider the problem of estimating a continuous time Markov chain (CTMC)
generator from discrete time observations, which is essentially a missing data problem in
statistics. These generators give rise to transition probabilities (in particular probabilities of
default) over any time horizon, hence the estimation of such generators is a key problem in
the world of banking, where the regulator requires banks to calculate risk over different time
horizons. For this particular problem several algorithms have been proposed, however, through
a combination of theoretical and numerical results we show the Expectation Maximisation (EM)
algorithm to be the superior choice. Furthermore we derive closed form expressions for the
associated Wald confidence intervals (error) estimated by the EM algorithm. Previous attempts
to calculate such intervals relied on numerical schemes which were slower and less stable. We
further provide a closed form expression (via the Delta method) to transfer these errors to the
level of the transition probabilities, which are more intuitive. Although one can establish more
precise mathematical results with the Markov assumption, there is empirical evidence suggesting
this assumption is not valid. We finish this part by carrying out empirical research on non-Markov
phenomena and propose a model to capture the so-called rating momentum. This model has
many appealing features and is a natural extension to the Markov set up.
The second part is based on McKean Vlasov Stochastic Differential Equations (MV-SDEs),
these Stochastic Differential Equations (SDEs) arise from looking at the limit, as the number of
weakly interacting particles (e.g. gas particles) tends to infinity. The resulting SDE has coefficients
which can depend on its own law, making them theoretically more involved. Although MV-SDEs
arise from statistical physics, there has been an explosion in interest recently to use MV-SDEs
in models for economics. We firstly derive an explicit approximation scheme for MV-SDEs with
one-sided Lipschitz growth in the drift. Such a condition was observed to be an issue for standard
SDEs and required more sophisticated schemes. There are implicit and explicit schemes one
can use and we develop both types in the setting of MV-SDEs. Another main issue for MVSDEs
is, due to the dependency on their own law they are extremely expensive to simulate
compared to standard SDEs, hence techniques to improve computational cost are in demand.
The final result in this part is to develop an importance sampling algorithm for MV-SDEs, where
our measure change is obtained through the theory of large deviation principles. Although
importance sampling results for standard SDEs are reasonably well understood, there are several
difficulties one must overcome to apply a good importance sampling change of measure in this
setting. The importance sampling is used here as a variance reduction technique although our
results hint that one may be able to use it to reduce propagation of chaos error as well.
Finally we consider stochastic algorithms to solve PDEs. It is known one can achieve numerical
advantages by using probabilistic methods to solve PDEs, through the so-called probabilistic
domain decomposition method. The main result of this part is to present an unbiased stochastic
representation for a first order PDE, based on the theory of branching diffusions and regime
switching. This is a very interesting result since previously (Itô based) stochastic representations
only applied to second order PDEs. There are multiple issues one must overcome in order to
obtain an algorithm that is numerically stable and solves such a PDE. We conclude by showing
the algorithm’s potential on a more general first order PDE
Essays in financial asset pricing
Three essays in financial asset pricing are given; one concerning the partial differential equation (PDE) pricing and hedging of a class of continuous/generalized power mean Asian options, via their (optimal) Lie point symmetry groups, leading to practical pricing formulas. The second presents high-frequency predictions of S&P 500 returns via several machine learning models, statistically significantly demonstrating short-horizon market predictability and economically significantly profitable (beyond transaction costs) trading strategies. The third compares profitability between these [(mean) ensemble] strategies and Asian option Δ-hedging, using results of the first. Interpreting bounds on arithmetic Asian option prices as ask and bid values, hedging profitability depends largely on securing prices closer to the bid, and settling midway between the bid and ask, significant profits are consistently accumulated during the years 2004-2016. Ensemble predictive trading the S&P 500 yields comparatively very small returns, despite trading much more frequently. The pricing and hedging of (arithmetic) Asian options are difficult and have spurred several solution approaches, differing in theoretical insight and practicality. Multiple families of exact solutions to relaxed power mean Asian option pricing boundary-value problems are explicitly established, which approximately satisfy the full pricing problem, and in one case, converge to exact solutions under certain parametric restrictions. Corresponding hedging parameters/ Greeks are derived. This family consists of (optimal) invariant solutions, constructed for the corresponding pricing PDEs. Numerical experiments explore this family behaviorally, achieving reliably accurate pricing. The second chapter studies intraday market return predictability. Regularized linear and nonlinear tree-based models enjoy significant predictability. Ensemble models perform best across time and their return predictability realizes economically significant profits with Sharpe ratios after transaction costs of 0.98. These results strongly evidence that intraday market returns are predictable during short time horizons, beyond that explainable by transaction costs. The lagged constituent returns are shown to hold significant predictive information not contained in lagged market returns or price trend and liquidity characteristics. Consistent with the hypothesis that predictability is driven by slow-moving trader capital, predictability decreased post-decimalization, and market returns are more predictable midday, on days with high volatility or illiquidity, and during financial crises