A mathematical analysis of the Gumbel test for jumps in stochastic volatility models

This article gives an exhaustive mathematical analysis of the Gumbel test for additive jump components based on extreme value theory. The Gumbel test was first introduced by Lee and Mykland in 2008 from an economical point of view. They consider a continuous-time stochastic volatility model with a general continuous volatility process and observe it under a high-frequency sampling scheme. The test statistics based on the maximum of increments converges to the Gumbel distribution under the null hypothesis of no additive jump component and to infinity otherwise. Our article presents a moment method based technique that provides some deeper mathematical insights into the convergence and divergence case of the test statistics. In the non-jump case we are able to prove the convergence to the Gumbel distribution under greatly weak assumptions: The volatility process has to be merely pathwise Hölder continuous with an arbitrary random Hölder exponent and we have no restrictions concerning an additional drift term. Therefore, for example, we are allowing for long and short-range dependence. In the case of existing additive jumps, we give divergence results in a general semimartingale setting and investigate the speed of divergence depending on the jump activity. As a by-product of our analysis we also deduce an optimal pathwise estimator for the spot volatility process. Moreover, we provide a detailed simulation study that compares the power of the Gumbel test with the power of the jump test proposed by Barndorff–Nielsen and Shephard in 2006 for Hölder exponents close to zero. Finally, both tests are applied to a real dataset.     

FFT based option pricing under a mean reverting process with stochastic volatility and jumps

Numerous studies present strong empirical evidence that certain financial assets may exhibit mean reversion, stochastic volatility or jumps. This paper explores the valuation of European options when the underlying asset follows a mean reverting log-normal process with stochastic volatility and jumps. A closed form representation of the characteristic function of the process is derived for the computation of European option prices via the fast Fourier transform.     

Estimating spot volatility in the presence of infinite variation jumps

We propose a kernel estimator for the spot volatility of a semi-martingale at a given time point by using high frequency data, where the underlying process accommodates a jump part of infinite variation. The estimator is based on the representation of the characteristic function of Lévy processes. The consistency of the proposed estimator is established under some mild assumptions. By assuming that the jump part of the underlying process behaves like a symmetric stable Lévy process around 0, we establish the asymptotic normality of the proposed estimator. In particular, with a specific kernel function, the estimator is variance efficient. We conduct Monte Carlo simulation studies to assess our theoretical results and compare our estimator with existing ones.     

On the explicit evaluation of the Geometric Asian options in stochastic volatility models with jumps

In the present paper we provide a semiexplicit valuation formula for Geometric Asian options, with fixed and floating strike under continuous monitoring, when the underlying stock price process exhibits both stochastic volatility and jumps. More precisely, we shall work in the Barndorff-Nielsen and Shephard (BNS) model framework. We shall provide some numerical illustrations of the results obtained.     

The long memory HEAVY process: modeling and forecasting financial volatility

This paper studies the bivariate HEAVY system of volatility regression equations and its various extensions that are directly applicable to the day-to-day business treasury operations of trading in foreign exchange and commodities, investing in bond and stock markets, hedging out market risk, and capital budgeting. We enrich the HEAVY framework with powers, asymmetries, and long memory that improve its forecasting accuracy significantly. Other findings are as follows. First, hyperbolic memory fits the realized measure better, whereas fractional integration is more suitable for the powered returns. Second, the structural breaks applied to the bivariate system capture the time-varying behavior of the parameters, in particular during and after the global financial crisis of 2007/2008.

     

On the Esscher transforms and other equivalent martingale measures for Barndorff-Nielsen and Shephard stochastic volatility models with jumps

We compute and then discuss the Esscher martingale transform for exponential processes, the Esscher martingale transform for linear processes, the minimal martingale measure, the class of structure preserving martingale measures, and the minimum entropy martingale measure for stochastic volatility models of the Ornstein–Uhlenbeck type as introduced by Barndorff-Nielsen and Shephard. We show that in the model with leverage, with jumps both in the volatility and in the returns, all those measures are different, whereas in the model without leverage, with jumps in the volatility only and a continuous return process, several measures coincide, some simplifications can be made and the results are more explicit. We illustrate our results with parametric examples used in the literature.     

The tail empirical process for long memory stochastic volatility sequences

This paper describes the limiting behaviour of tail empirical processes associated with long memory stochastic volatility models. We show that such a process has dichotomous behaviour, according to an interplay between the Hurst parameter and the tail index. On the other hand, the tail empirical process with random levels never suffers from long memory. This is very desirable from a practical point of view, since such a process may be used to construct the Hill estimator of the tail index. To prove our results we need to establish new results for regularly varying distributions, which may be of independent interest.     

Testing the Gumbel Hypothesis Via the Pot-Method

We consider an i.i.d. sample, generated by some distribution function, which belongs to the domain of attraction of an extreme value distribution with unknown shape and scale parameters. We treat the scale parameter as a nuisance parameter and establish for the hypothesis of Gumbel domain of attraction an asymptotically optimal test based on those observations among the sample, which exceed a given threshold sequence. Asymptotic optimality is achieved along certain contiguous extreme value alternatives within the concept of local asymptotic normality (LAN). Adaptive test procedures exist under restrictive assumptions. The finite sample size behavior of the proposed test is studied by simulations and it is compared to that of a test based on the sample coefficient of variation.     

A fast numerical approach to option pricing with stochastic interest rate,stochastic volatility and double jumps

This study proposes a pricing model through allowing for stochastic interest rate and stochastic volatility in the double exponential jump-diffusion setting. The characteristic function of the proposed model is then derived. Fast numerical solutions for European call and put options pricing based on characteristic function and fast Fourier transform (FFT) technique are developed. Simulations show that our numerical technique is accurate, fast and easy to implement, the proposed model is suitable for modeling long-time real-market changes. The model and the proposed option pricing method are useful for empirical analysis of asset returns and risk management in firms.     

Gumbel fluctuations for cover times in the discrete torus

This work proves that the fluctuations of the cover time of simple random walk in the discrete torus of dimension at least three with large side-length are governed by the Gumbel extreme value distribution. This result was conjectured for example in Aldous and Fill (Reversible Markov chains and random walks on graphs, in preparation). We also derive some corollaries which qualitatively describe “how” covering happens. In addition, we develop a new and stronger coupling of the model of random interlacements, introduced by Sznitman (Ann Math (2) 171(3):2039–2087, 2010), and random walk in the torus. This coupling is used to prove the cover time result and is also of independent interest.     

Testing the Gumbel hypothesis by Galton's ratio

We study a generalization of a ratio of spacings introduced by Galton in 1902. The ratio proves to be an important building block in the construction of a large sample test for the hypothesis that a distribution from an extremal domain of attraction belongs to the domain of attraction of the Gumbel law.     

A test for a parametric form of the volatility in second-order diffusion models

Second-order diffusion models have been found to be promising in analyzing financial market data. Based on nonparametric fitting, Nicolau (Stat Probabil Lett 78(16):2700–2704, 2008) suggested that the quadratic function may be an appropriate specification of the volatility when a second-order diffusion model is used to analyze some European and American financial market data sets, which motivates us to develop a formal statistical test for this finding. To achieve the task, a generalized likelihood ratio test is proposed in this paper and a residual-based bootstrap is suggested to compute the p value of the test. The analysis of many real-world financial market data sets demonstrates that the quadratic specification of the volatility function is in general reasonable.     

The distribution of extremes in the degree sequence: A Gumbel distribution approach

In this work, an approximation of the asymptotics of the distribution for the maximum and minimum of the degree sequence defined over a random graph is determined using a new approach in terms of the Gumbel distribution for extremes.     

Singularities of the matrix exponent of a Markov additive process with one-sided jumps

We analyze the number of zeros of det(F(α))$det\left(F\left(\alpha \right)\right)$, where F(α)$F\left(\alpha \right)$ is the matrix exponent of a Markov Additive Process (MAP) with one-sided jumps. The focus is on the number of zeros in the right half of the complex plane, where F(α)$F\left(\alpha \right)$ is analytic. In addition, we also consider the case of a MAP killed at an independent exponential time. The corresponding zeros can be seen as the roots of a generalized Cramér–Lundberg equation. We argue that our results are particularly useful in fluctuation theory for MAPs, which leads to numerous applications in queueing theory and finance.     

Gumbel分布的简单贝叶斯估计

在实际应用中,两参数Gumbel分布的贝叶斯估计往往需要预先知道Gumbel参数的二维联合先验分布。由于获取先验分布的主观性和统计推断的复杂性,目前有关Gumbel分布贝叶斯估计理论及其性质的讨论还比较少,更不要说获得较为简单的Gumbel分布的贝叶斯估计。本文基于Kaminskiy和Vasiliy提出的简单贝叶斯估计过程,利用可靠度函数估计的区间形式表示先验信息,从而得到两个参数Gumbel分布的简单贝叶斯估计。基于此先验信息,该估计过程构造了Gumbel参数的连续联合先验分布,给出了在给定任意时点的可靠度(或累积密度)及其标准差的后验估计,为可靠性与风险评估中简单快速的使用贝叶斯估计刻画极端事件提供了可能.     

Turing jumps in the Ershov hierarchy

We look at infinite levels of the Ershov hierarchy in the natural system of notation, which are proper for jumps of sets. It is proved that proper infinite levels for jumps are confined to D_a ^{- 1} Delta_a^{ - 1} -levels, where a stands for an ordinal ωⁿ > 1.     

基于多维Gumbel Copula函数的投资组合VaR分析

基于目前国内有关Copula函数的实证研究主要是研究二种资产的相关性为主,文章根据Copula函数在构建反映随机变量实际分布与相关性的联合分布函数上具有的优势,首先利用GJR模型构建资产的边缘分布,接着利用多元阿基米德Copula函数族中的Gumbel Copula函数构建了反映多个资产收益实际分布和相关性的联合分布函数,并使用蒙特卡罗模拟技术,分析在不同置信度下的投资组合的最小风险价值(VaR)及其资产组成,实证说明根据文章提出的模型度量资产的风险,可以使投资者选择的资产更加稳健,同时也有利于投资者对投资组合整体风险进行分散和监管。     

Estimation of limiting conditional distributions for the heavy tailed long memory stochastic volatility process

We consider Stochastic Volatility processes with heavy tails and possible long memory in volatility. We study the limiting conditional distribution of future events given that some present or past event was extreme (i.e. above a level which tends to infinity). Even though extremes of stochastic volatility processes are asymptotically independent (in the sense of extreme value theory), these limiting conditional distributions differ from the i.i.d. case. We introduce estimators of these limiting conditional distributions and study their asymptotic properties. If volatility has long memory, then the rate of convergence and the limiting distribution of the centered estimators can depend on the long memory parameter (Hurst index).