Estimation of Jump-Diffusion Processes via Empirical Characteristic Functions

In this paper, we propose a fractional stochastic volatility jump-diffusion model which extends the Bates(1996) model. Where we model the volatility as a fractional process. Extensive empirical studies show that the distributions of the logarithmic returns of financial asset usually exhibit properties of self-similarity and long-range dependence and since the fractional Brownian motion has these two important properties, it has the ability to capture the behavior of underlying asset price. Further incorporating jumps into the stochastic volatility framework gives further freedom to financial mathematicians to fit both the short and long end of the implied volatility surface. We propose a stochastic model which contains both fractional and jump process. Then we price options using Monte Carlo simulations along with a variance reduction technique(antithetic variates). We use market data from the S&P 500 index and we compare our results with the Heston and Bates model using error measures. The results show our model greatly outperforms previous models in terms of estimation accuracy.

The Stochastic Volatility Jump Diffusion (SVJD) model is a common tool used for option pricing. Existing literature has employed many techniques in order to estimate the parameters of this model, including both Bayesian and frequentist approaches. In 2010, Duan and Yeh [Jump and volatility risk premiums implied by VIX. J Econ Dyn Control. 2010;34(11):2232–2244] used a frequentist approach with maximum likelihood estimation that included the VIX, the Chicago Board Options Exchange (CBOE) Volatility Index, which carries real time S&P 500 option price information in addition to the S&P 500 data. In this paper, we combine the previous Bayesian approach in Cape et al. [Estimating Heston's and Bates' models parameters using Markov Chain Monte Carlo simulation. J Stat Comput Simul. 2015;85(11):2295–2314] with the inclusion of VIX by Duan and Yeh. We provide the derivations of the conditional posterior distributions of the SVJD model's parameters and a guide to a Markov Chain Monte Carlo (MCMC) algorithm used to estimate the parameters. We apply the algorithm to the historical S&P 500 and VIX data and provide the subsequent parameter estimates. We offer these tools for others performing similar research with our R code available upon request.

What Jump Process to Use to Model S&P 500 Returns?

Estimating volatility and model parameters of stochastic volatility models with jumps using particle filter

This thesis presents our study on using the hybrid stochastic-local volatility model for option pricing. Many researchers have demonstrated that stochastic volatility models cannot capture the whole volatility surface accurately, although the model parameters have been calibrated to replicate the market implied volatility data for near at-the-money strikes. On the other hand, the local volatility model can reproduce the implied volatility surface, whereas it does not consider the stochastic behaviour of the volatility. To combine the advantages of stochastic volatility (SV) and local volatility (LV) models, a class of stochastic-local volatility (SLV) models has been developed. The SLV model contains a stochastic volatility component represented by a volatility process and a local volatility component represented by a so-called leverage function. The leverage function can be roughly seen as a ratio between local volatility and conditional expectation of stochastic volatility. The difficulty of implementing the SLV model lies in the calibration of the leverage function. In the thesis, we first review the fundamental theories of stochastic differential equations and the classic option pricing models, and study the behaviour of the volatility in the context of FX market. We then introduce the SLV model and illustrate our implementation of the calibration and pricing procedure. We apply the SLV model to exotic option pricing in the FX market and compare pricing results of the SLV model with pure local volatility and pure stochastic volatility models. Numerical results show that the SLV model can match the implied volatility surface very well as well as improve the pricing performance for barrier options. In addition, we further discuss some extensions of the SLV project, such as parallelization potential for accelerating option pricing and pricing techniques for window barrier options. Although the SLV model we use in the thesis is not entirely new, we contribute to the research in the following aspects: 1) we investigate the hybrid volatility modeling thoroughly from theoretical backgrounds to practical implementations; 2) we resolve some critical issues in implementing the SLV model such as developing a fast and stable numerical method to derive the leverage function; and 3) we build a robust calibration and pricing platform under the SLV model, which can be extended for practical uses.

This thesis presents our study on using the hybrid stochastic-local volatility model for option pricing. Many researchers have demonstrated that stochastic volatility models cannot capture the whole volatility surface accurately, although the model parameters have been calibrated to replicate the market implied volatility data for near at-the-money strikes. On the other hand, the local volatility model can reproduce the implied volatility surface, whereas it does not consider the stochastic behaviour of the volatility. To combine the advantages of stochastic volatility (SV) and local volatility (LV) models, a class of stochastic-local volatility (SLV) models has been developed. The SLV model contains a stochastic volatility component represented by a volatility process and a local volatility component represented by a so-called leverage function. The leverage function can be roughly seen as a ratio between local volatility and conditional expectation of stochastic volatility. The difficulty of implementing the SLV model lies in the calibration of the leverage function. In the thesis, we first review the fundamental theories of stochastic differential equations and the classic option pricing models, and study the behaviour of the volatility in the context of FX market. We then introduce the SLV model and illustrate our implementation of the calibration and pricing procedure. We apply the SLV model to exotic option pricing in the FX market and compare pricing results of the SLV model with pure local volatility and pure stochastic volatility models. Numerical results show that the SLV model can match the implied volatility surface very well as well as improve the pricing performance for barrier options. In addition, we further discuss some extensions of the SLV project, such as parallelization potential for accelerating option pricing and pricing techniques for window barrier options. Although the SLV model we use in the thesis is not entirely new, we contribute to the research in the following aspects: 1) we investigate the hybrid volatility modeling thoroughly from theoretical backgrounds to practical implementations; 2) we resolve some critical issues in implementing the SLV model such as developing a fast and stable numerical method to derive the leverage function; and 3) we build a robust calibration and pricing platform under the SLV model, which can be extended for practical uses.

Variational Bayes approximation of factor stochastic volatility models

We examine the effect of varying the temperature points on MEMS inertial sensors’ noise models using Allan variance and least-squares spectral analysis (LSSA). Allan variance is a method of representing root-mean-square random drift error as a function of averaging times. LSSA is an alternative to the classical Fourier methods and has been applied successfully by a number of researchers in the study of the noise characteristics of experimental series. Static data sets are collected at different temperature points using two MEMS-based IMUs, namely MotionPakII and Crossbow AHRS300CC. The performance of the two MEMS inertial sensors is predicted from the Allan variance estimation results at different temperature points and the LSSA is used to study the noise characteristics and define the sensors’ stochastic model parameters. It is shown that the stochastic characteristics of MEMS-based inertial sensors can be identified using Allan variance estimation and LSSA and the sensors’ stochastic model parameters are temperature dependent. Also, the Kaiser window FIR low-pass filter is used to investigate the effect of de-noising stage on the stochastic model. It is shown that the stochastic model is also dependent on the chosen cut-off frequency.

We consider a class of stochastic path-dependent volatility models where the stochastic volatility, whose square follows the Cox--Ingersoll--Ross model, is multiplied by a (leverage) function of th...

A better estimation of surface reaction efficiency of semiconductor-grade graphite with atomic nitrogen, as well as the calibration error are calculated using Bayesian updating based on experimental data. Compared with a conventional deterministic model, the stochastic model approach is a powerful tool in the sense that the model is capable of taking into account underlying error correlations among the data quantities. In this paper, we investigate four different stochastic models (called “stochastic system model classes” herein) corresponding to different descriptions of modeling and measurement error structures, given one deterministic physical model. These stochastic system model classes differ in the covariance matrix structure that is used in the uncertainty model to represent uncertainties associated with the physical model and experimental measurements. For each model class, Bayesian inference is used to estimate the posterior probabilities of the physical model parameters as well as of the stochastic model parameters. Model comparison and selection are then applied based on two measures including Bayesian evidence and Bayesian information criterion, as well as the deviance information criterion. Both measures suggest the stochastic model class, which considers that a correlation between errors in two data quantities among different data points is the most plausible. With the stochastic model class, the range of uncertainty in surface reaction efficiency is estimated to be about two orders of magnitude at .

This paper introduces a general class of stochastic volatility models that can serve as a basis for modeling and estimating simultaneous equation systems involving multivariate volatility processes. The class consists of processes obtained as monotone polynomial transformations of so-called stochastic autoregressive volatility (SARV) models. The class permits a flexible modeling of volatility and avoids strong distributional assumptions. Most of the standard stochastic volatility models are incorporated in the framework, including the lognormal autoregressive model and stochastic volatility generalizations of GARCH and EGARCH. General conditions for strict statinarity, ergodicity, and the existence of finite lower-order unconditional moments for SARV are established. Closed-form expressions for the lowerorder moments constitute the basis for the proposed generalized method of moments estimation procedure. As an illustration we consider an informationdriven model of the return volatility–trading volume system for financial assets. Parameters with structural interpretation are estimated and an extended Kalmanfilter procedure allows for the construction of volatility forecasts based on the bivariate return-volume series. Empirical results are based on daily NYSE data for IBM over 1973–1991.

Stochastic models specified by statistical source parameters have been proposed to describe the details of fault plane irregularities. The statisical parameters have been estimated through the analysis of strong-motion records by using these models. In the study herein, a new semi-empirical method is developed to estimate strong ground-motions in the short period range, taking into account the fault plane irregularities.The fault plane of large earthquake is assumed to consist of some small segments with a specific size, which is called fault patch size. The displacement Dij of each segment is divided into two parts. One is the average displacement D0 over the fault plane and the other is the deviation ΔDij from D0. Seismic wave SL(t) due to D0 and the seismic wave SS(t) due to ΔDij are synthesized separately as follows:SL(t)=nΣi=1nΣj=1nΣk=1(Re/Rij)Se(t-τijk)SS(t)=nΣi=1nΣj=1κij(Re/Rij)Se(t-τij1)where Se(t) shows an observed seismic wave of element earthquake, of which fault area should be consistent with the fault patch size, Re is the source to site distance of the element earthquake, Rij is a distance from each segment on the fault plane to the site, and τijk is a delay time in relation to rupture propagation and travel time of seismic waves.SL(t) can be derived based on the similarity conditions of source parameters and n is determined from the cubic root of the seismic-moment ratio between the large and the element earthquakes. κij is introduced as a ratio between the deviation ΔDij of each segment on the fault plane and the displacement of element earthquake, and is assumed to be a probabiristic parameter of normal distribution with average of 0 and with standard deviation of SD. SD shows the degree of inhomogeneity of displacement on the fault plane. The final result of synthesized strong ground-motion can be expressed as a sum of SL(t)and SS(t).The new method is applied to synthesis of strong ground motions at Akita during the 1983 Nihonkai-chubu earthquake (M=7.7) and the relations among SD and some statistical parameters of other stochastic fault models are discussed theoretically.

The penalized simulated maximum likelihood (PSML) approach can be used to estimate parameters for a stochastic differential equation model based on completely or partially observed discrete-time observations. The PSML uses an auxiliary variable importance sampler and parameters are estimated in a penalized maximum likelihood framework. In this paper, we extend the PSML to allow for measurement error, including unknown initial conditions. Simulation studies for two stochastic models and a real world example aimed at understanding the dynamics of chronic wasting disease illustrate that our method has favorable performance in the presence of measurement error. PSML reduces both the bias and root mean squared error as compared to existing methods. Lastly, we establish consistency and asymptotic normality for the proposed estimators.

Wastewater production, like many other engineered and environmental processes, is inherent stochastic in nature and requires the use of complex stochastic models, for example, to predict realistic patterns of down-the-drain chemicals or pharmaceuticals and personal care products. Up until now, a formal method of statistical inference has been lacking for many of those models, where explicit likelihood functions were intractable. In this Article, we investigate Approximate Bayesian Computation (ABC) methods to infer important parameters of stochastic environmental models. ABC methods have been recently suggested to perform model-based inference in a Bayesian setting when model likelihoods are analytically or computationally intractable and have not been applied to environmental systems analysis or water quality modeling before. In a case study, we investigate the performance of three different algorithms to infer the number of wastewater pulses contained in three high-resolution data series of benzotriazole and total nitrogen loads in sewers. We find that all algorithms perform well and that the uncertainty in the inferred number of corresponding wastewater pulses varies between 6% and 28%. In our case, the results are more sensitive to substance characteristics than to catchment properties. Although the application of ABC methods requires careful tuning and attention to detail, they have a great general potential to update stochastic model parameters with monitoring data and improve their predictive capabilities.

Increasing attention has been focused on the analysis of the realized volatility, which can be treated as a proxy for the true volatility. In this paper, we study the potential use of the realized volatility as a proxy in a stochastic volatility model estimation. We estimate the leveraged stochastic volatility model using the realized volatility computed from five popular methods across six sampling-frequency transaction data (from 1-min to 60- min) based on the trust region method. Availability of the realized volatility allows us to estimate the model parameters via the MLE and thus avoids computational challenge in the high dimensional integration. Six stock indices are considered in the empirical investigation. We discover some consistent findings and interesting patterns from the empirical results. In general, the significant leverage effect is consistently detected at each sampling frequency and the volatility persistence becomes weaker at the lower sampling frequency.