The Effect of Non-Stationarity on Extreme Sea-Level Estimation

The sea-level is the composition of astronomical tidal and meteorological surge processes. It exhibits temporal non-stationarity due to a combination of long-term trend in the mean level, the deterministic tidal component, surge seasonality and interactions between the tide and surge. We assess the effect of these non-stationarities on the estimation of the distribution of extreme sea-levels. This is important for coastal flood assessment as the traditional method of analysis assumes that, once the trend has been removed, extreme sea-levels are from a stationary sequence. We compare the traditional approach with a recently proposed alternative that incorporates the knowledge of the tidal component and its associated interactions, by applying them to 22 UK data sites and through a simulation study. Our main finding is that if the tidal non-stationarity is ignored then a substantial underestimation of extreme sea-levels results for most sites. In contrast, if surge seasonality and the tide–surge interaction are not modelled the traditional approach produces little additional bias. The alternative method is found to perform well but requires substantially more statistical modelling and better data quality.     

Maximum term of a particular autoregressivesequence with discrete margins

In this paper we consider a stationary sequence of discrete random variables with marginal distribution H(x), obtained by a simple transformation from the max-AR(1) sequence considered by Alpuim (1989). Because discrete distributions impose severe restrictions on the convergence of the normalized maxima to an extreme value distribution, it is seen that in this particular case, whenever H(x) belongs to the domain of attraction of any max-stable distribution, the sequence possesses an extremal index 0 = 0. Nevertheless, it, is possible to obtain a nondegenerate limiting distribution for the linearized maxima by choosing other sets of normalizing constants. Whenever H(x) does not belong to the domain of attraction of any max-stable distribution, but, satisfies adequate conditions, the maxima nearly possess an asymptotic stability with the presence of an extremal index 0 <θ<1.

Motivated by the behaviour of these sequences we obtained a more general result extending the results of Anderson (1970) and Me (Jon nick and Park (1992) over the mixing conditionsD ^(k)(u_n), defined by Chermck et al (1991).

Several examples, obtained after simulation, are presented in order to illustrate the different situations that may occur.     

A location-invariant non-positive moment-type estimator of the extreme value index

This paper investigates a class of location invariant non-positive moment-type estimators of extreme value index, which is highly flexible due to the tuning parameter involved. Its asymptotic expansions and its optimal sample fraction in terms of minimal asymptotic mean square error are derived. A small scale Monte Carlo simulation turns out that the new estimators, with a suitable choice of the tuning parameter driven by the data itself, perform well compared to the known ones. Finally, the proposed estimators with a bootstrap optimal sample fraction are applied to an environmental data set.     

Extremal analysis of short series with outliers: sea-levels and athletics records

The analysis of extreme values is often required from short series which are biasedly sampled or contain outliers. Data for sea-levels at two UK east coast sites and data on athletics records for women's 3000 m track races are shown to exhibit such characteristics. Univariate extreme value methods provide a poor quantification of the extreme values for these data. By using bivariate extreme value methods we analyse jointly these data with related observations, from neighbouring coastal sites and 1500 m races respectively. We show that using bivariate methods provides substantial benefits, both in these applications and more generally with the amount of information gained being determined by the degree of dependence, the lengths and the amount of overlap of the two series, the homogeneity of the marginal characteristics of the variables and the presence and type of the outlier.     

Too Connected to Fail? Inferring Network Ties From Price Co-Movements

We use extreme value theory methods to infer conventionally unobservable connections between financial institutions from joint extreme movements in credit default swap spreads and equity returns. Estimated pairwise co-crash probabilities identify significant connections among up to 186 financial institutions prior to the crisis of 2007/2008. Financial institutions that were very central prior to the crisis were more likely to be bailed out during the crisis or receive the status of systemically important institutions. This result remains intact also after controlling for indicators of too-big-to-fail concerns, systemic, systematic, and idiosyncratic risks. Both credit default swap (CDS)-based and equity-based connections are significant predictors of bailouts. Supplementary materials for this article are available online.     

多元极值的参数建模方法及其金融应用：最新进展述评

 由于现实中的极值事件往往倾向于同时或相继发生,因此多元极值研究正成为极值统计学的理论前沿和研究热点。本文对该领域中参数建模方法的最新进展做了系统性述评,包括经典多元极值理论、Ledford-Tawn-Ramos方法和Heffernan和Tawn条件法等,并指出了这些建模方法的优缺点以及未来可能的理论突破点。本文还全面分析了近年来多元极值分析方法在金融领域的国内外应用现状,并探讨其未来的应用前景,可能是在金融传染、组合问题和系统性风险管理等方面。     

Diagnostics for dependence within time series extremes

Summary. The analysis of extreme values within a stationary time series entails various assumptions concerning its long- and short-range dependence. We present a range of new diagnostic tools for assessing whether these assumptions are appropriate and for identifying structure within extreme events. These tools are based on tail characteristics of joint survivor functions but can be implemented by using existing estimation methods for extremes of univariate independent and identically distributed variables. Our diagnostic aids are illustrated through theoretical examples, simulation studies and by application to rainfall and exchange rate data. On the basis of these diagnostics we can explain characteristics that are found in the observed extreme events of these series and also gain insight into the properties of events that are more extreme than those observed.     

Anticipating catastrophes through extreme value modelling

Summary. When catastrophes strike it is easy to be wise after the event. It is also often argued that such catastrophic events are unforeseeable, or at least so implausible as to be negligible for planning purposes. We consider these issues in the context of daily rainfall measurements recorded in Venezuela. Before 1999 simple extreme value techniques were used to assess likely future levels of extreme rainfall, and these gave no particular cause for concern. In December 1999 a daily precipitation event of more than 410 mm, almost three times the magnitude of the previously recorded maximum, caused devastation and an estimated 30000 deaths. We look carefully at the previous history of the process and offer an extreme value analysis of the data—with some methodological novelty—that suggests that the 1999 event was much more plausible than the previous analyses had claimed. Deriving design parameters from the results of such an analysis may have had some mitigating effects on the consequences of the subsequent disaster. The themes of the new analysis are simple: the full exploitation of available data, proper accounting of uncertainty, careful interpretation of asymptotic limit laws and allowance for non-stationarity. The effect on the Venezuelan data analysis is dramatic. The broader implications are equally dramatic; that a naïve use of extreme value techniques is likely to lead to a false sense of security that might have devastating consequences in practice.     

A semi-parametric model for multivariate extreme values

Threshold methods for multivariate extreme values are based on the use of asymptotically justified approximations of both the marginal distributions and the dependence structure in the joint tail. Models derived from these approximations are fitted to a region of the observed joint tail which is determined by suitably chosen high thresholds. A drawback of the existing methods is the necessity for the same thresholds to be taken for the convergence of both marginal and dependence aspects, which can result in inefficient estimation. In this paper an extension of the existing models, which removes this constraint, is proposed. The resulting model is semi-parametric and requires computationally intensive techniques for likelihood evaluation. The methods are illustrated using a coastal engineering application.     

On the estimation and application of max-stable processes

The theory of max-stable processes generalizes traditional univariate and multivariate extreme value theory by allowing for processes indexed by a time or space variable. We consider a particular class of max-stable processes, known as M4 processes, that are particularly well adapted to modeling the extreme behavior of multiple time series. We develop procedures for determining the order of an M4 process and for estimating the parameters. To illustrate the methods, some examples are given for modeling jumps in returns in multivariate financial time series. We introduce a new measure to quantify and predict the extreme co-movements in price returns.     

Extremal analysis of processes sampled at different frequencies

The observed extremes of a discrete time process depend on the process itself and the sampling frequency. We develop theoretical results which show how to account for the effect of sampling frequency on extreme values, thus enabling us to analyse systematically extremal data from series with different sampling rates. We present statistical methodology based on these results which we illustrate though simulations and by applications to sea-waves and rainfall data.     

An extended Gaussian max-stable process model for spatial extremes

The extremes of environmental processes are often of interest due to the damage that can be caused by extreme levels of the processes. These processes are often spatial in nature and modelling the extremes jointly at many locations can be important. In this paper, an extension of the Gaussian max-stable process is developed, enabling data from a number of locations to be modelled under a more flexible framework than in previous applications. The model is applied to annual maximum rainfall data from five sites in South-West England. For estimation we employ a pairwise likelihood within a Bayesian analysis, incorporating informative prior information.     

Modeling of maximum precipitation using maximal generalized extreme value distribution

Distribution of maximum or minimum values (extreme values) of a dataset is especially used in natural phenomena including sea waves, flow discharge, wind speeds, and precipitation and it is also used in many other applied sciences such as reliability studies and analysis of environmental extreme events. So if we can explain the extremal behavior via statistical formulas, we can estimate how their behavior would be in the future. In this paper, we study extreme values of maximum precipitation in Zahedan using maximal generalized extreme value distribution, which all maxima of a data set are modeled using it. Also, we apply four methods to estimate distribution parameters including maximum likelihood estimation, probability weighted moments, elemental percentile and quantile least squares then compare estimates by average scaled absolute error criterion and obtain quantiles estimates and confidence intervals. In addition, goodness-of-fit tests are described. As a part of result, the return period of maximum precipitation is computed.     

Asymptotic distributions for the intervals estimators of the extremal index and the cluster-size probabilities

Under appropriate long range dependence conditions, the point process of exceedances of a stationary sequence weakly converges to a homogeneous compound Poisson point process. This limiting point process can be characterized by the extremal index and the cluster-size probabilities. In this paper we address the problem of estimating these quantities and we consider the intervals estimators introduced in Ferro and Segers [2003. Inference for clusters of extreme values. J. Roy. Statist. Soc. Ser. B 545–556] and in Ferro [2004. Statistical methods for clusters of extreme values. Ph.D. Thesis, Lancaster University]. We establish asymptotic weak convergence to Gaussian random variables and we give their asymptotic variance.     

Comparison of approaches for estimating the probability of coastal flooding

Coastal flooding is typically caused by combinations of extreme water-levels and large waves. Two extreme value methods, one univariate and the other multivariate, have been used for estimating the probability of coastal flooding at an existing flood defence structure and for aiding the design of a new structure. The properties of these two methods are compared in terms of extrapolation, sophistication and use of information for a range of extremal dependence structures. We find that, when applied to the assessment of the safety offered by an existing Dutch dike, the multivariate approach provides the more useful and accurate design information and has the substantial benefits of consistency and reduced statistical analysis when applied to several sites along a Dutch coastline.     

A Bayesian approach to zero-inflated data in extremes

Abstract

The generalized extreme value (GEV) distribution is known as the limiting result for the modeling of maxima blocks of size n, which is used in the modeling of extreme events. However, it is possible for the data to present an excessive number of zeros when dealing with extreme data, making it difficult to analyze and estimate these events by using the usual GEV distribution. The Zero-Inflated Distribution (ZID) is widely known in literature for modeling data with inflated zeros, where the inflator parameter w is inserted. The present work aims to create a new approach to analyze zero-inflated extreme values, that will be applied in data of monthly maximum precipitation, that can occur during months where there was no precipitation, being these computed as zero. An inference was made on the Bayesian paradigm, and the parameter estimation was made by numerical approximations of the posterior distribution using Markov Chain Monte Carlo (MCMC) methods. Time series of some cities in the northeastern region of Brazil were analyzed, some of them with predominance of non-rainy months. The results of these applications showed the need to use this approach to obtain more accurate and with better adjustment measures results when compared to the standard distribution of extreme value analysis.     

Bias and size corrections in extreme value modeling

Extreme value theory models have found applications in myriad fields. Maximum likelihood (ML) is attractive for fitting the models because it is statistically efficient and flexible. However, in small samples, ML is biased to O(N^?1) and some classical hypothesis tests suffer from size distortions. This paper derives the analytical Cox–Snell bias correction for the generalized extreme value (GEV) model, and for the model's extension to multiple order statistics (GEVr). Using simulations, the paper compares this correction to bootstrap-based bias corrections, for the generalized Pareto, GEV, and GEVr. It then compares eight approaches to inference with respect to primary parameters and extreme quantiles, some including corrections. The Cox–Snell correction is not markedly superior to bootstrap-based correction. The likelihood ratio test appears most accurately sized. The methods are applied to the distribution of geomagnetic storms.     

Nonparametric estimation of extreme conditional quantiles

The estimation of extreme conditional quantiles is an important issue in different scientific disciplines. Up to now, the extreme value literature focused mainly on estimation procedures based on independent and identically distributed samples. Our contribution is a two-step procedure for estimating extreme conditional quantiles. In a first step nonextreme conditional quantiles are estimated nonparametrically using a local version of [Koenker, R. and Bassett, G. (1978). Regression quantiles. Econometrica, 46, 33–50.] regression quantile methodology. Next, these nonparametric quantile estimates are used as analogues of univariate order statistics in procedures for extreme quantile estimation. The performance of the method is evaluated for both heavy tailed distributions and distributions with a finite right endpoint using a small sample simulation study. A bootstrap procedure is developed to guide in the selection of an optimal local bandwidth. Finally the procedure is illustrated in two case studies.