Estimation of jump-diffusion processes via empirical characteristic functions

Preface The starting point for this work and eventually the subject of the whole thesis was the question: how to estimate parameters of the affine stochastic volatility jump-diffusion models. These models are very important for contingent claim pricing. Their major advantage, availability T of analytical solutions for characteristic functions, made them the models of choice for many theoretical constructions and practical applications. At the same time, estimation of parameters of stochastic volatility jump-diffusion models is not a straightforward task. The problem is coming from the variance process, which is non-observable. There are several estimation methodologies that deal with estimation problems of latent variables. One appeared to be particularly interesting. It proposes the estimator that in contrast to the other methods requires neither discretization nor simulation of the process: the Continuous Empirical Characteristic function estimator (EGF) based on the unconditional characteristic function. However, the procedure was derived only for the stochastic volatility models without jumps. Thus, it has become the subject of my research. This thesis consists of three parts. Each one is written as independent and self contained article. At the same time, questions that are answered by the second and third parts of this Work arise naturally from the issues investigated and results obtained in the first one. The first chapter is the theoretical foundation of the thesis. It proposes an estimation procedure for the stochastic volatility models with jumps both in the asset price and variance processes. The estimation procedure is based on the joint unconditional characteristic function for the stochastic process. The major analytical result of this part as well as of the whole thesis is the closed form expression for the joint unconditional characteristic function for the stochastic volatility jump-diffusion models. The empirical part of the chapter suggests that besides a stochastic volatility, jumps both in the mean and the volatility equation are relevant for modelling returns of the S&P500 index, which has been chosen as a general representative of the stock asset class. Hence, the next question is: what jump process to use to model returns of the S&P500. The decision about the jump process in the framework of the affine jump- diffusion models boils down to defining the intensity of the compound Poisson process, a constant or some function of state variables, and to choosing the distribution of the jump size. While the jump in the variance process is usually assumed to be exponential, there are at least three distributions of the jump size which are currently used for the asset log-prices: normal, exponential and double exponential. The second part of this thesis shows that normal jumps in the asset log-returns should be used if we are to model S&P500 index by a stochastic volatility jump-diffusion model. This is a surprising result. Exponential distribution has fatter tails and for this reason either exponential or double exponential jump size was expected to provide the best it of the stochastic volatility jump-diffusion models to the data. The idea of testing the efficiency of the Continuous ECF estimator on the simulated data has already appeared when the first estimation results of the first chapter were obtained. In the absence of a benchmark or any ground for comparison it is unreasonable to be sure that our parameter estimates and the true parameters of the models coincide. The conclusion of the second chapter provides one more reason to do that kind of test. Thus, the third part of this thesis concentrates on the estimation of parameters of stochastic volatility jump- diffusion models on the basis of the asset price time-series simulated from various "true" parameter sets. The goal is to show that the Continuous ECF estimator based on the joint unconditional characteristic function is capable of finding the true parameters. And, the third chapter proves that our estimator indeed has the ability to do so. Once it is clear that the Continuous ECF estimator based on the unconditional characteristic function is working, the next question does not wait to appear. The question is whether the computation effort can be reduced without affecting the efficiency of the estimator, or whether the efficiency of the estimator can be improved without dramatically increasing the computational burden. The efficiency of the Continuous ECF estimator depends on the number of dimensions of the joint unconditional characteristic function which is used for its construction. Theoretically, the more dimensions there are, the more efficient is the estimation procedure. In practice, however, this relationship is not so straightforward due to the increasing computational difficulties. The second chapter, for example, in addition to the choice of the jump process, discusses the possibility of using the marginal, i.e. one-dimensional, unconditional characteristic function in the estimation instead of the joint, bi-dimensional, unconditional characteristic function. As result, the preference for one or the other depends on the model to be estimated. Thus, the computational effort can be reduced in some cases without affecting the efficiency of the estimator. The improvement of the estimator s efficiency by increasing its dimensionality faces more difficulties. The third chapter of this thesis, in addition to what was discussed above, compares the performance of the estimators with bi- and three-dimensional unconditional characteristic functions on the simulated data. It shows that the theoretical efficiency of the Continuous ECF estimator based on the three-dimensional unconditional characteristic function is not attainable in practice, at least for the moment, due to the limitations on the computer power and optimization toolboxes available to the general public. Thus, the Continuous ECF estimator based on the joint, bi-dimensional, unconditional characteristic function has all the reasons to exist and to be used for the estimation of parameters of the stochastic volatility jump-diffusion models.

Special Sentence Revocations High, Data Download, Iowa Department of Corrections, May 2010

In addition to their original sentence, persons convicted of sexual abuse, incest or sexual exploitation of a minor also receive a “special sentence” of ten years, or in some cases, life. In its prison population forecast, the Iowa Division of Criminal and Juvenile Justice Planning noted “an unexpectedly high rate of revocation among those released to the special sentence, particularly given past research that has shown Iowa sex offenders having very low rates of re-arrest and/or return to prison.”

Special Sentence Sex Offender Caseloads Adding Up, Data Downloads, Iowa Department of Corrections, December 2008

Several years ago the General Assembly increased the penalties for certain sex offenses by an additional ten years of community-based supervision, and in some case lifetime supervision. The Division of Criminal and Juvenile Justice Planning (CJJP) studied the effect the new law would have on CBC supervision caseloads:

Special investigation of the City of Scarville for the period January 1, 2006 through March 31, 2011

Special investigation of the City of Scarville for the period January 1, 2006 through March 31, 2011

Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Report on a special investigation of the City of Bloomfield Volunteer Fire Department for the period July 1, 2006 through June 30, 2010

Report on a special investigation of the City of Bloomfield Volunteer Fire Department for the period July 1, 2006 through June 30, 2010

Report on a special investigation of the Mills County Treasurer’s Office for the period January 1, 2005 through March 31, 2011

Report on a special investigation of the Mills County Treasurer’s Office for the period January 1, 2005 through March 31, 2011

Report on a special investigation of the City of Davis City for the period October 1, 2006 through November 26, 2011

Report on a special investigation of the City of Davis City for the period October 1, 2006 through November 26, 2011

Report on a special investigation of the Amana Colonies Land Use District for the period May 21, 2007 through March 31, 2011

Report on a special investigation of the Amana Colonies Land Use District for the period May 21, 2007 through March 31, 2011

Report on a special investigation of the City of Callender for the period January 1, 2009 through January 3, 2012

Report on a special investigation of the City of Callender for the period January 1, 2009 through January 3, 2012

Report on a special investigation of the Benton County Conservation Department for the period June 1, 2010 through November 7, 2011

Report on a special investigation of the Benton County Conservation Department for the period June 1, 2010 through November 7, 2011

On the configuration of Herman rings of meromorphic functions

We prove some results concerning the possible configuration s of Herman rings for transcendental meromorphic functions. We show that one pole is enough to obtain cycles of Herman rings of arbitrary period a nd give a sufficient condition for a configuration to be realizable.

Option valuation as an expectation in the complex domain: The Black-Scholes case

It is very well known that the first succesful valuation of a stock option was done by solving a deterministic partial differential equation (PDE) of the parabolic type with some complementary conditions specific for the option. In this approach, the randomness in the option value process is eliminated through a no-arbitrage argument. An alternative approach is to construct a replicating portfolio for the option. From this viewpoint the payoff function for the option is a random process which, under a new probabilistic measure, turns out to be of a special type, a martingale. Accordingly, the value of the replicating portfolio (equivalently, of the option) is calculated as an expectation, with respect to this new measure, of the discounted value of the payoff function. Since the expectation is, by definition, an integral, its calculation can be made simpler by resorting to powerful methods already available in the theory of analytic functions. In this paper we use precisely two of those techniques to find the well-known value of a European call

Report on a special investigation of the City of Blencoe for the period January 1, 2008 through September 30, 2011

Report on a special investigation of the City of Blencoe for the period January 1, 2008 through September 30, 2011