Extremes of Markov chains with tail switching potential

Summary.  A recent advance in the utility of extreme value techniques has been the characteri- zation of the extremal behaviour of Markov chains. This has enabled the application of extreme value models to series whose temporal dependence is Markovian, subject to a limitation that prevents switching between extremely high and extremely low levels. For many applications this is sufficient, but for others, most notably in the field of finance, it is common to find series in which successive values switch between high and low levels. We term such series Markov chains with tail switching potential, and the scope of this paper is to generalize the previous theory to enable the characterization of the extremal properties of series displaying this type of behaviour. In addition to theoretical developments, a modelling procedure is proposed. A simulation study is made to assess the utility of the model in inferring the extremal dependence structure of autoregressive conditional heteroscedastic processes, which fall within the tail switching Markov family, and generalized autoregressive conditional heteroscedastic processes which do not, being non-Markov in general. Finally, the procedure is applied to model extremal aspects of a financial index extracted from the New York Stock Exchange compendium.     

A Latent Process Model for Temporal Extremes

This paper presents a hierarchical approach to modelling extremes of a stationary time series. The procedure comprises two stages. In the first stage, exceedances over a high threshold are modelled through a generalized Pareto distribution, which is represented as a mixture of an exponential variable with a Gamma distributed rate parameter. In the second stage, a latent Gamma process is embedded inside the exponential distribution in order to induce temporal dependence among exceedances. Unlike other hierarchical extreme‐value models, this version has marginal distributions that belong to the generalized Pareto family, so that the classical extreme‐value paradigm is respected. In addition, analytical developments show that different choices of the underlying Gamma process can lead to different degrees of temporal dependence of extremes, including asymptotic independence. The model is tested through a simulation study in a Markov chain setting and used for the analysis of two datasets, one environmental and one financial. In both cases, a good flexibility in capturing different types of tail behaviour is obtained.     

A Bayesian kriged Kalman model for short-term forecasting of air pollution levels

Summary.  Short-term forecasts of air pollution levels in big cities are now reported in news-papers and other media outlets. Studies indicate that even short-term exposure to high levels of an air pollutant called atmospheric particulate matter can lead to long-term health effects. Data are typically observed at fixed monitoring stations throughout a study region of interest at different time points. Statistical spatiotemporal models are appropriate for modelling these data. We consider short-term forecasting of these spatiotemporal processes by using a Bayesian kriged Kalman filtering model. The spatial prediction surface of the model is built by using the well-known method of kriging for optimum spatial prediction and the temporal effects are analysed by using the models underlying the Kalman filtering method. The full Bayesian model is implemented by using Markov chain Monte Carlo techniques which enable us to obtain the optimal Bayesian forecasts in time and space. A new cross-validation method based on the Mahalanobis distance between the forecasts and observed data is also developed to assess the forecasting performance of the model implemented.     

A conditional approach for multivariate extreme values (with discussion)

Summary.  Multivariate extreme value theory and methods concern the characterization, estimation and extrapolation of the joint tail of the distribution of a d -dimensional random variable. Existing approaches are based on limiting arguments in which all components of the variable become large at the same rate. This limit approach is inappropriate when the extreme values of all the variables are unlikely to occur together or when interest is in regions of the support of the joint distribution where only a subset of components is extreme. In practice this restricts existing methods to applications where d is typically 2 or 3. Under an assumption about the asymptotic form of the joint distribution of a d -dimensional random variable conditional on its having an extreme component, we develop an entirely new semiparametric approach which overcomes these existing restrictions and can be applied to problems of any dimension. We demonstrate the performance of our approach and its advantages over existing methods by using theoretical examples and simulation studies. The approach is used to analyse air pollution data and reveals complex extremal dependence behaviour that is consistent with scientific understanding of the process. We find that the dependence structure exhibits marked seasonality, with ex- tremal dependence between some pollutants being significantly greater than the dependence at non-extreme levels.     

Estimating the probability of simultaneous rainfall extremes within a region: a spatial approach

In this paper we investigate the impact of model mis-specification, in terms of the dependence structure in the extremes of a spatial process, on the estimation of key quantities that are of interest to hydrologists and engineers. For example, it is often the case that severe flooding occurs as a result of the observation of rainfall extremes at several locations in a region simultaneously. Thus, practitioners might be interested in estimates of the joint exceedance probability of some high levels across these locations. It is likely that there will be spatial dependence present between the extremes, and this should be properly accounted for when estimating such probabilities. We compare the use of standard models from the geostatistics literature with max-stables models from extreme value theory. We find that, in some situations, using an incorrect spatial model for our extremes results in a significant under-estimation of these probabilities which – in flood defence terms – could lead to substantial under-protection.     

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.     

Statistics for Extreme Sea Currents

Estimates of various characteristics of extreme sea currents, such as speeds and their directions, are required when designing offshore structures. This paper extends standard statistical methods for extreme values to handle the directionality, temporal dependence and tidal non-stationarity that are present in sea current extremes. The methods are applied to a short period of data from the Inner Dowsing Light Tower in the North Sea. Substantial benefits, over existing methods, are obtained from our analysis of the sea current by decomposing it into tide and surge currents. In particular, we find that at the Inner Dowsing the strong directionality in extreme sea current speeds is completely explained by the tidal current and directionality in the non-extreme surge currents. This finding aids model fitting and extrapolation.     

Semi‐Parametric Models for the Multivariate Tail Dependence Function – the Asymptotically Dependent Case

Abstract. In general, the risk of joint extreme outcomes in financial markets can be expressed as a function of the tail dependence function of a high‐dimensional vector after standardizing marginals. Hence, it is of importance to model and estimate tail dependence functions. Even for moderate dimension, non‐parametrically estimating a tail dependence function is very inefficient and fitting a parametric model to tail dependence functions is not robust. In this paper, we propose a semi‐parametric model for (asymptotically dependent) tail dependence functions via an elliptical copula. Under this model assumption, we propose a novel estimator for the tail dependence function, which proves favourable compared to the empirical tail dependence function estimator, both theoretically and empirically.     

Dependence between oil and commodities markets using time-varying Archimedean copulas and effectiveness of hedging strategies

The aim of this work is to study in a first step the dependence between oil and some commodity prices (cotton, rice, wheat, sucre, coffee, and silver) using copula theory, and then in a second step to determine the optimal hedging strategy for oil–commodity portfolio against the risk of negative variation in commodity markets prices. The model is implemented with an AR-GARCH model with innovations that follow t distribution for the marginal distribution and the extreme value copula for the joint distribution and parameters and dependence indices are re-estimated in each new day which allow taking into account nonlinear dependence, tails behavior, and their development over time. Various copula functions are used to model the dependence structure between oil and commodity markets. Empirical results show an increase in the dependence during the last 6 years. Volatility for commodity prices registered record levels in the same time with the increase in uncertainty. Optimal hedging ratio varies over time as a consequence of the change in the dependence structure.     

A new class of models for bivariate joint tails

Summary.  A fundamental issue in applied multivariate extreme value analysis is modelling dependence within joint tail regions. The primary focus of this work is to extend the classical pseudopolar treatment of multivariate extremes to develop an asymptotically motivated representation of extremal dependence that also encompasses asymptotic independence. Starting with the usual mild bivariate regular variation assumptions that underpin the coefficient of tail dependence as a measure of extremal dependence, our main result is a characterization of the limiting structure of the joint survivor function in terms of an essentially arbitrary non-negative measure that must satisfy some mild constraints. We then construct parametric models from this new class and study in detail one example that accommodates asymptotic dependence, asymptotic independence and asymmetry within a straightforward parsimonious parameterization. We provide a fast simulation algorithm for this example and detail likelihood-based inference including tests for asymptotic dependence and symmetry which are useful for submodel selection. We illustrate this model by application to both simulated and real data. In contrast with the classical multivariate extreme value approach, which concentrates on the limiting distribution of normalized componentwise maxima, our framework focuses directly on the structure of the limiting joint survivor function and provides significant extensions of both the theoretical and the practical tools that are available for joint tail modelling.     

Testing asymptotic independence in bivariate extremes

Bivariate extreme value condition (see (1.1) below) includes the marginal extreme value conditions and the existence of the (extreme) dependence function. Two cases are of interest: asymptotic independence and asymptotic dependence. In this paper, we investigate testing the existence of the dependence function under the null hypothesis of asymptotic independence and present two suitable test statistics. Small simulations are studied and the application for a real data is shown. The other case with the null hypothesis of asymptotic dependence is already investigated.     

Optimal design of air quality networks detecting warning and alert conditions

A statistical method is presented to determine the optima design of air quality networks detecting warning and alert levels. A simulation model is used to describe temporal and spatial variations of atmospheric pollutants; air quality patterns serve as the database of the procedure to design the network. Only the sites exceeding warning and alert levels, at different meteorological scenarios, are considered as potential monitoring stations. For the selection of the optima set, spatial and temporal representativity criteria are introduced; accordingly, the optima set provides a complete representativity of the space and time considered. The method is applied to the Mestre urban area, in Venice district, for the carbon monoxide pollutant.     

Trend and Return Level of Extreme Snow Events in New York City

Abstract

A major winter storm brought up to 42 inches of snow in parts of the Mid-Atlantic and Northeast states for January 22–24, 2016. The blizzard of January 2016 impacted about 102.8 million people, claiming at least 55 lives and $500 million to $3 billion in economic losses. This article studies two important aspects of extreme snowfall events: 1. trends in annual maxima and threshold exceedances and 2. return levels for extreme snowfall. Applying extreme value methods to the extreme snow data in the New York City area, we quantify linear trends in extreme snowfall and assess how severe the 2016 blizzard is in terms of return levels. To find a more realistic standard error for the extreme value methods, we extend Smith’s method to adapt to both spatial and temporal correlations in the snow data. Our results show increasing, but insignificant trends in the annual maximum snowfall series. However, we find that the 87.5th percentile snowfall has significantly increased by 0.564 inches per decade, suggesting that, while the maximum snowfall is not significantly increasing, there have been increases in the snowfall among the larger storms. We also find that the 2016 blizzard is indeed an extreme snow event equivalent to about a 40-year return level in the New York City area. The extreme value methods used in this study are thoroughly illustrated for general readers. Data and modularized programming codes are to be available online to aid practitioners in using extreme value methods in applications.     

Non-parametric Estimation of Tail Dependence

Abstract.  Dependencies between extreme events (extremal dependencies) are attracting an increasing attention in modern risk management. In practice, the concept of tail dependence represents the current standard to describe the amount of extremal dependence. In theory, multi-variate extreme-value theory turns out to be the natural choice to model the latter dependencies. The present paper embeds tail dependence into the concept of tail copulae which describes the dependence structure in the tail of multivariate distributions but works more generally. Various non-parametric estimators for tail copulae and tail dependence are discussed, and weak convergence, asymptotic normality, and strong consistency of these estimators are shown by means of a functional delta method. Further, weak convergence of a general upper-order rank-statistics for extreme events is investigated and the relationship to tail dependence is provided. A simulation study compares the introduced estimators and two financial data sets were analysed by our methods.     

Likelihood Inference for Multivariate Extreme Value Distributions Whose Spectral Vectors have known Conditional Distributions

Multivariate extreme value statistical analysis is concerned with observations on several variables which are thought to possess some degree of tail dependence. The main approaches to inference for multivariate extremes consist in approximating either the distribution of block component‐wise maxima or the distribution of the exceedances over a high threshold. Although the expressions of the asymptotic density functions of these distributions may be characterized, they cannot be computed in general. In this paper, we study the case where the spectral random vector of the multivariate max‐stable distribution has known conditional distributions. The asymptotic density functions of the multivariate extreme value distributions may then be written through univariate integrals that are easily computed or simulated. The asymptotic properties of two likelihood estimators are presented, and the utility of the method is examined via simulation.     

Spatial modeling of extreme rainfall in northeast Thailand

It is well recognized that the generalized extreme value (GEV) distribution is widely used for any extreme events. This notion is based on the study of discrete choice behavior; however, there is a limit for predicting the distribution at ungauged sites. Hence, there have been studies on spatial dependence within extreme events in continuous space using recorded observations. We model the annual maximum daily rainfall data consisting of 25 locations for the period from 1982 to 2013. The spatial GEV model that is established under observations is assumed to be mutually independent because there is no spatial dependency between the stations. Furthermore, we divide the region into two regions for a better model fit and identify the best model for each region. We show that the regional spatial GEV model reflects the spatial pattern well compared with the spatial GEV model over the entire region as the local GEV distribution. The advantage of spatial extreme modeling is that more robust return levels and some indices of extreme rainfall can be obtained for observed stations as well as for locations without observed data. Thus, the model helps to determine the effects and assessment of vulnerability due to heavy rainfall in northeast Thailand.     

A Case Study for Modelling Cancer Incidence Using Bayesian Spatio‐Temporal Models

Researchers familiar with spatial models are aware of the challenge of choosing the level of spatial aggregation. Few studies have been published on the investigation of temporal aggregation and its impact on inferences regarding disease outcome in space–time analyses. We perform a case study for modelling individual disease outcomes using several Bayesian hierarchical spatio‐temporal models, while taking into account the possible impact of spatial and temporal aggregation. Using longitudinal breast cancer data from South East Queensland, Australia, we consider both parametric and non‐parametric formulations for temporal effects at various levels of aggregation. Two temporal smoothness priors are considered separately; each is modelled with fixed effects for the covariates and an intrinsic conditional autoregressive prior for the spatial random effects. Our case study reveals that different model formulations produce considerably different model performances. For this particular dataset, a classical parametric formulation that assumes a linear time trend produces the best fit among the five models considered. Different aggregation levels of temporal random effects were found to have little impact on model goodness‐of‐fit and estimation of fixed effects.     

Modeling Dependence in High Dimensions With Factor Copulas

This article presents flexible new models for the dependence structure, or copula, of economic variables based on a latent factor structure. The proposed models are particularly attractive for relatively high-dimensional applications, involving 50 or more variables, and can be combined with semiparametric marginal distributions to obtain flexible multivariate distributions. Factor copulas generally lack a closed-form density, but we obtain analytical results for the implied tail dependence using extreme value theory, and we verify that simulation-based estimation using rank statistics is reliable even in high dimensions. We consider “scree” plots to aid the choice of the number of factors in the model. The model is applied to daily returns on all 100 constituents of the S&P 100 index, and we find significant evidence of tail dependence, heterogeneous dependence, and asymmetric dependence, with dependence being stronger in crashes than in booms. We also show that factor copula models provide superior estimates of some measures of systemic risk. Supplementary materials for this article are available online.