Tail index estimation for heavy tails: accommodation of bias in the excesses over a high threshold

In statistics of extremes, inference is often based on the excesses over a high random threshold. Those excesses are approximately distributed as the set of order statistics associated to a sample from a generalized Pareto model. We then get the so-called “maximum likelihood” estimators of the tail index γ. In this paper, we are interested in the derivation of the asymptotic distributional properties of a similar “maximum likelihood” estimator of a positive tail index γ, based also on the excesses over a high random threshold, but with a trial of accommodation of bias in the Pareto model underlying those excesses. We next proceed to an asymptotic comparison of the two estimators at their optimal levels. An illustration of the finite sample behaviour of the estimators is provided through a small-scale Monte Carlo simulation study. Research partially supported by FCT/POCTI and POCI/FEDER.     

On the block thresholding wavelet estimators with censored data

We consider block thresholding wavelet-based density estimators with randomly right-censored data and investigate their asymptotic convergence rates. Unlike for the complete data case, the empirical wavelet coefficients are constructed through the Kaplan-Meier estimators of the distribution functions in the censored data case. On the basis of a result of Stute [W. Stute, The central limit theorem under random censorship, Ann. Statist. 23 (1995) 422-439] that approximates the Kaplan-Meier integrals as averages of i.i.d. random variables with a certain rate in probability, we can show that these wavelet empirical coefficients can be approximated by averages of i.i.d. random variables with a certain error rate in L². Therefore we can show that these estimators, based on block thresholding of empirical wavelet coefficients, achieve optimal convergence rates over a large range of Besov function classes , p≥2, q≥1 and nearly optimal convergence rates when 1≤p<2. We also show that these estimators achieve optimal convergence rates over a large class of functions that involve many irregularities of a wide variety of types, including chirp and Doppler functions, and jump discontinuities. Therefore, in the presence of random censoring, wavelet estimators still provide extensive adaptivity to many irregularities of large function classes. The performance of the estimators is tested via a modest simulation study.     

Local asymptotic normality in a stationary model for spatial extremes

De Haan and Pereira (2006) [6] provided models for spatial extremes in the case of stationarity, which depend on just one parameter β>0 measuring tail dependence, and they proposed different estimators for this parameter. We supplement this framework by establishing local asymptotic normality (LAN) of a corresponding point process of exceedances above a high multivariate threshold. Standard arguments from LAN theory then provide the asymptotic minimum variance within the class of regular estimators of β. It turns out that the relative frequency of exceedances is a regular estimator sequence with asymptotic minimum variance, if the underlying observations follow a multivariate extreme value distribution or a multivariate generalized Pareto distribution.     

Bias Free Threshold Estimation for Jump Intensity Function

In this paper, combining the threshold technique, we reconstruct Nadaraya-Watson estimation using Gamma asymmetric kernels for the unknown jump intensity function of a diffusion process with finite activity jumps. Under mild conditions, we obtain the asymptotic normality for the proposed estimator. Moreover, we have verified the better finite-sampling properties such as bias correction and effciency gains of the underlying estimator compared with other nonparametric estimators through a Monte Carlo experiment.     

Optimal kernel estimation of densities

Precise asymptotic behavior for mean integrated squared error (MISE) is determined for sequences of kernel estimators of a density in a broad class, including discontinuous and possibly unbounded densities. The paper shows that the sequence using the kernel optimal at each fixed sample size is asymptotically more efficient than a sequence generated by changing the bandwidth of a fixed kernel shape, regardless of the kernel shape. The class of densities considered are those whose characteristic functions behave at large arguments like the product of a Fourier series and a regularly varying function. This condition may be related to the smoothness of an m-th derivative of the density.Partially supported by National Science Foundation Grant DMS-8711924.     

Controlled jump processes

Finite and infinite planning horizon Markov decision problems are formulated for a class of jump processes with general state and action spaces and controls which are measurable functions on the time axis taking values in an appropriate metrizable vector space. For the finite horizon problem, the maximum expected reward is the unique solution, which exists, of a certain differential equation and is a strongly continuous function in the space of upper semi-continuous functions. A necessary and sufficient condition is provided for an admissible control to be optimal, and a sufficient condition is provided for the existence of a measurable optimal policy. For the infinite horizon problem, the maximum expected total reward is the fixed point of a certain operator on the space of upper semi-continuous functions. A stationary policy is optimal over all measurable policies in the transient and discounted cases as well as, with certain added conditions, in the positive and negative cases.     

��������������ɢģ���л��ֲ����ʵķDz�������

??This paper studies nonparametric estimation of the integrated volatility of Poisson jump-diffusion processes with noisy high-frequency data. We propose jump-robust two-scale and multi-scale estimators. The estimators are based on a combination of the multi-scale method and threshold technique, which serves to remove microstructure noise and jumps, respectively. Furthermore, asymptotic properties of the proposed estimators, such as consistency, are established.     

On the link between infinite horizon control and quasi-stationary distributions

We study infinite horizon control of continuous-time non-linear branching processes with almost sure extinction for general (positive or negative) discount. Our main goal is to study the link between infinite horizon control of these processes and an optimization problem involving their quasi-stationary distributions and the corresponding extinction rates. More precisely, we obtain an equivalent of the value function when the discount parameter is close to the threshold where the value function becomes infinite, and we characterize the optimal Markov control in this limit. To achieve this, we present a new proof of the dynamic programming principle based upon a pseudo-Markov property for controlled jump processes. We also prove the convergence to a unique quasi-stationary distribution of non-linear branching processes controlled by a Markov control conditioned on non-extinction.     

Mean-Variance Hedging on Uncertain Time Horizon in a Market with a Jump

In this work, we study the problem of mean-variance hedging with a random horizon T∧τ, where T is a deterministic constant and τ is a jump time of the underlying asset price process. We first formulate this problem as a stochastic control problem and relate it to a system of BSDEs with a jump. We then provide a verification theorem which gives the optimal strategy for the mean-variance hedging using the solution of the previous system of BSDEs. Finally, we prove that this system of BSDEs admits a solution via a decomposition approach coming from filtration enlargement theory.     

Stopping of functionals with discontinuity at the boundary of an open set

We explore properties of the value function and existence of optimal stopping times for functionals with discontinuities related to the boundary of an open (possibly unbounded) set O. The stopping horizon is either random, equal to the first exit from the set O, or fixed (finite or infinite). The payoff function is continuous with a possible jump at the boundary of O. Using a generalization of the penalty method, we derive a numerical algorithm for approximation of the value function for general Feller-Markov processes and show existence of optimal or ε-optimal stopping times.     

Reduced‐bias kernel estimators of a positive extreme value index

In this paper, we deal with the semi‐parametric estimation of the extreme value index, an important parameter in extreme value analysis. It is well known that many classic estimators, such as the Hill estimator, reveal a strong bias. This problem motivated the study of two classes of kernel estimators. Those classes generalize the classical Hill estimator and have a tuning parameter that enables us to modify the asymptotic mean squared error and eventually to improve their efficiency. Since the improvement in efficiency is not very expressive, we also study new reduced bias estimators based on the two classes of kernel statistics. Under suitable conditions, we prove their asymptotic normality. Moreover, an asymptotic comparison, at optimal levels, shows that the new classes of reduced bias estimators are more efficient than other reduced bias estimator from the literature. An illustration of the finite sample behaviour of the kernel reduced‐bias estimators is also provided through the analysis of a data set in the field of insurance.     

Recursive Estimation of Regression Functions by Local Polynomial Fitting

The recursive estimation of the regression function m(x) = E(Y/X = x) and its derivatives is studied under dependence conditions. The examined method of nonparametric estimation is a recursive version of the estimator based on locally weighted polynomial fitting, that in recent articles has proved to be an attractive technique and has advantages over other popular estimation techniques. For strongly mixing processes, expressions for the bias and variance of these estimators are given and asymptotic normality is established. Finally, a simulation study illustrates the proposed estimation method.     

时间序列中方差的结构变点的小波识别(英文)

本文给出了时间序列中方差的小波系数的两种估计:连续估计和离散估计.这两种估计可以用来检测时间序列中方差的结构变点.利用这两种估计我们给出了方差变点的位置和跳跃幅度的估计,并且显示出这些估计可达到最佳收敛速度.同时,我们还给出了这些估计的收敛速度以及检验统计量的渐进分布!     

Optimal importance sampling for Lévy processes

We develop importance sampling estimators for Monte Carlo pricing of European and path-dependent options in models driven by Lévy processes. Using results from the theory of large deviations for processes with independent increments, we compute an explicit asymptotic approximation for the variance of the pay-off under a time-dependent Esscher-style change of measure. Minimizing this asymptotic variance using convex duality, we then obtain an importance sampling estimator of the option price. We show that our estimator is logarithmically optimal among all importance sampling estimators. Numerical tests in the variance gamma model show consistent variance reduction with a small computational overhead.     

Stochastic stability for Markovian jump linear systems associated with a finite number of jump times

This paper deals with a stochastic stability concept for discrete-time Markovian jump linear systems. The random jump parameter is associated to changes between the system operation modes due to failures or repairs, which can be well described by an underlying finite-state Markov chain. In the model studied, a fixed number of failures or repairs is allowed, after which, the system is brought to a halt for maintenance or for replacement. The usual concepts of stochastic stability are related to pure infinite horizon problems, and are not appropriate in this scenario. A new stability concept is introduced, named stochastic τ-stability that is tailored to the present setting. Necessary and sufficient conditions to ensure the stochastic τ-stability are provided, and the almost sure stability concept associated with this class of processes is also addressed. The paper also develops equivalences among second order concepts that parallels the results for infinite horizon problems.     

Robust Estimation of the Generalized Pareto Distribution

One approach used for analyzing extremes is to fit the excesses over a high threshold by a generalized Pareto distribution. For the estimation of the shape and scale parameters in the generalized Pareto distribution, under some restrictions on the value of the scale parameter, maximum likelihood, method of moments and probability weighted moments' estimators are available. However, these are not robust estimators. In this paper we implement a robust estimation procedure known as the method of medians (He and Fung, 1999) to estimate the parameters in the generalized Pareto distribution. The asymptotic distribution of our estimator is normal for any value of the shape parameter except –1.     

Set Intersection Theorems and Existence of Optimal Solutions

The question of nonemptiness of the intersection of a nested sequence of closed sets is fundamental in a number of important optimization topics, including the existence of optimal solutions, the validity of the minimax inequality in zero sum games, and the absence of a duality gap in constrained optimization. We consider asymptotic directions of a sequence of closed sets, and introduce associated notions of retractive, horizon, and critical directions, based on which we provide new conditions that guarantee the nonemptiness of the corresponding intersection. We show how these conditions can be used to obtain simple and unified proofs of some known results on existence of optimal solutions, and to derive some new results, including a new extension of the Frank–Wolfe Theorem for (nonconvex) quadratic programming.     

Nonparametric estimation of the spectral measure,and associated dependence measures,for multivariate extreme values using a limiting conditional representation

The traditional approach to multivariate extreme values has been through the multivariate extreme value distribution G, characterised by its spectral measure H and associated Pickands’ dependence function A. More generally, for all asymptotically dependent variables, H determines the probability of all multivariate extreme events. When the variables are asymptotically dependent and under the assumption of unit Fréchet margins, several methods exist for the estimation of G, H and A which use variables with radial component exceeding some high threshold. For each of these characteristics, we propose new asymptotically consistent nonparametric estimators which arise from Heffernan and Tawn’s approach to multivariate extremes that conditions on variables with marginal values exceeding some high marginal threshold. The proposed estimators improve on existing estimators in three ways. First, under asymptotic dependence, they give self-consistent estimators of G, H and A; existing estimators are not self-consistent. Second, these existing estimators focus on the bivariate case, whereas our estimators extend easily to describe dependence in the multivariate case. Finally, for asymptotically independent cases, our estimators can model the level of asymptotic independence; whereas existing estimators for the spectral measure treat the variables as either being independent, or asymptotically dependent. For asymptotically dependent bivariate random variables, the new estimators are found to compare favourably with existing estimators, particularly for weak dependence. The method is illustrated with an application to finance data.     

Nonparametric regression on random fields with random design using wavelet method

We consider non-linear wavelet-based estimators of spatial regression functions with (known) random design on strictly stationary random fields, which are indexed by the integer lattice points in the \(N\)-dimensional Euclidean space and are assumed to satisfy some mixing conditions. We investigate their asymptotic rates of convergence based on thresholding of empirical wavelet coefficients and show that these estimators achieve nearly optimal convergence rates within a logarithmic term over a large range of Besov function classes \(B^{s}_{p,q}\). Therefore, wavelet estimators still achieve nearly optimal convergence rates for random fields and provide explicitly the extraordinary local adaptability.