Numerical analysis and multi-precision computational methods applied to the extant problems of Asian option pricing and simulating stable distributions and unit root densities
MetadataShow full item record
This thesis considers new methods that exploit recent developments in computer technology to address three extant problems in the area of Finance and Econometrics. The problem of Asian option pricing has endured for the last two decades in spite of many attempts to find a robust solution across all parameter values. All recently proposed methods are shown to fail when computations are conducted using standard machine precision because as more and more accuracy is forced upon the problem, round-off error begins to propagate. Using recent methods from numerical analysis based on multi-precision arithmetic, we show using the Mathematica platform that all extant methods have efficacy when computations use sufficient arithmetic precision. This creates the proper framework to compare and contrast the methods based on criteria such as computational speed for a given accuracy. Numerical methods based on a deformation of the Bromwich contour in the Geman-Yor Laplace transform are found to perform best provided the normalized strike price is above a given threshold; otherwise methods based on Euler approximation are preferred. The same methods are applied in two other contexts: the simulation of stable distributions and the computation of unit root densities in Econometrics. The stable densities are all nested in a general function called a Fox H function. The same computational difficulties as above apply when using only double-precision arithmetic but are again solved using higher arithmetic precision. We also consider simulating the densities of infinitely divisible distributions associated with hyperbolic functions. Finally, our methods are applied to unit root densities. Focusing on the two fundamental densities, we show our methods perform favorably against the extant methods of Monte Carlo simulation, the Imhof algorithm and some analytical expressions derived principally by Abadir. Using Mathematica, the main two-dimensional Laplace transform in this context is reduced to a one-dimensional problem.
Thesis, PhD Doctor of Philosophy
Creative Commons Attribution 4.0 Internationalhttp://creativecommons.org/licenses/by/4.0/
Except where otherwise noted within the work, this item's license for re-use is described as Creative Commons Attribution 4.0 International
Items in the St Andrews Research Repository are protected by copyright, with all rights reserved, unless otherwise indicated.
Showing items related by title, author, creator and subject.
Probing the average distribution of water in organic hydrate crystal structures with radial distribution functions (RDFs) Skyner, Rachael Elaine; Mitchell, John B. O.; Groom, Colin (2017-01-28) - Journal articleThe abundance of crystal structures of solvated organic molecules reflects the common role of solvent in the crystallisation process. An understanding of solvation is therefore important for crystal engineering, with solvent ...
Bates, Dominic J.; Tojeiro, Rita; Newman, Jeffrey A.; Gonzalez-Perez, Violeta; Comparat, Johan; Schneider, Donald P.; Lima, Marcos; Streblyanska, Alina (2019-07) - Journal articleThis paper presents stellar mass functions and i-band luminosity functions for Sloan Digital Sky Survey (SDSS) galaxies with i < 21 using clustering redshifts. From these measurements, we also compute targeting completeness ...
Monitoring the assembly-disassembly-organisation-reassembly process of germanosilicate UTL through in situ pair distribution function analysis Henkelis, Susan E.; Morris, Sam A.; Mazur, Michal; Wheatley, Paul S.; McHugh, Lauren N.; Morris, Russell (2018-08-06) - Journal articleA study into the disassembly and organisation steps of the ADOR process has been undertaken through in situ Pair Distribution Function (PDF) analysis. Three aqueous systems (water, 6 M HCl and 12 M HCl) were introduced to ...