Sub‐asymptotic motivation for new conditional multivariate extreme models

Anthony Christopher Davison, Thomas Lugrin, Jonathan A. Tawn
2021
Journal paper

Abstract

Statistical models for extreme values are generally derived from non-degenerate probabilistic limits that can be used to approximate the distribution of events that exceed a selected high threshold. If convergence to the limit distribution is slow, then the approximation may describe observed extremes poorly, and bias can only be reduced by choosing a very high threshold at the cost of unacceptably large variance in any subsequent tail inference. An alternative is to use sub-asymptotic extremal models, which introduce more parameters but can provide better fits for lower thresholds. We consider this problem in the context of the Heffernan–Tawn conditional tail model for multivariate extremes, which has found wide use due to its flexible handling of dependence in high-dimensional applications. Recent extensions of this model appear to improve joint tail inference. We seek a sub-asymptotic justification for why these extensions work and show that they can improve convergence rates by an order of magnitude for certain copulas. We also propose a class of extensions of them that may have wider value for statistical inference in multivariate extremes.

Official source

https://infoscience.epfl.ch/record/288683?ln=en

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Related publications

Contributions to Likelihood-Based Modelling of Extreme Values

Léo Raymond-Belzile

xtreme value analysis is concerned with the modelling of extreme events such as floods and heatwaves, which can have large impacts. Statistical modelling can be useful to better assess risks even if, due to scarcity of measurements, there is inherently very large residual uncertainty in any analysis. Driven by the increase in environmental databases, spatial modelling of extremes has expanded rapidly in the last decade. This thesis presents contributions to such analysis. The first chapter is about likelihood-based inference in the univariate setting and investigates the use of bias-correction and higher-order asymptotic methods for extremes, highlighting through examples and illustrations the unique challenge posed by data scarcity. We focus on parametric modelling of extreme values, which relies on limiting distributional results and for which, as a result, uncertainty quantification is complicated. We find that, in certain cases, small-sample asymptotic methods can give improved inference by reducing the error rate of confidence intervals. Two data illustrations, linked to assessment of the frequency of extreme rainfall episodes in Venezuela and the analysis of survival of supercentenarians, illustrate the methods developed. In the second chapter, we review the major methods for the analysis of spatial extremes models. We highlight the similarities and provide a thorough literature review along with novel simulation algorithms. The methods described therein are made available through a statistical software package. The last chapter focuses on estimation for a Bayesian hierarchical model derived from a multivariate generalized Pareto process. We review approaches for the estimation of censored components in models derived from (log)-elliptical distributions, paying particular attention to the estimation of a high-dimensional Gaussian distribution function via Monte Carlo methods. The impacts of model misspecification and of censoring are explored through extensive simulations and we conclude with a case study of rainfall extremes in Eastern Switzerland.

EPFL2019

Inference on the Angular Distribution of Extremes

Claudio Andri Semadeni

The spectral distribution plays a key role in the statistical modelling of multivariate extremes, as it defines the dependence structure of multivariate extreme-value distributions and characterizes the limiting distribution of the relative sizes of the components of large multivariate observations. No parametric family captures all possible types of multivariate dependence, and numerous parametric models have been proposed. Inference on the spectral distribution is typically based on the pseudo-angles of

large' observations under the assumption that they follow the spectral distribution. There has been little attention on studying the impact of this approximation on inference, and it turns out that it can yield significantly biased estimates. We provide a characterization of the angular distribution of excesses corresponding to the distribution of pseudo-angles of

large'observations that improves direct inference on the spectral distribution in the bivariate setting. Extremal dependence is at the heart of extreme value modelling and numerous measures to quantify it have been proposed in the literature. In many applications, datasets seem to exhibit asymmetry in the dependence between the variables. Many parametric multivariate extreme-value models can accommodate asymmetry in the sense that the spectral density can be asymmetric, resulting in a non-exchangeable dependence structure. There has been little attention paid to quantifying asymmetry at extreme levels, which can be useful for diagnosis and model checking. We propose a coefficient of extremal asymmetry that quantifies the asymmetry at extreme levels for pairs of variables. We also propose two non-parametric estimators of the coefficient of extremal asymmetry and compare their properties through numerical simulation. The two estimators have diametrically opposed bias-variance trade-offs. The estimator based on maximum empirical likelihood performs well and is nearly unbiased.

EPFL2020

The thesis is a contribution to extreme-value statistics, more precisely to the estimation of clustering characteristics of extreme values. One summary measure of the tendency to form groups is the inverse average cluster size. In extreme-value context, this parameter is called the extremal index, and apart from its relation with the size of groups, it appears as an important parameter measuring the effects of serial dependence on extreme levels in time series. Although several methods exist for its estimation in univariate sequences, these methods are only applicable for strictly stationary series satisfying a long-range asymptotic independence condition on extreme levels, cannot take covariates into consideration, and yield only crude estimates for the corresponding multivariate quantity. These are strong restrictions and great drawbacks. In climatic time series, both stationarity and asymptotic independence can be broken, due to climate change and possible long memory of the data, and not including information from simultaneously measured linked variables may lead to inefficient estimation. The thesis addresses these issues. First, we extend the theorem of Ferro and Segers (2003) concerning the distribution of inter-exceedance times: we introduce truncated inter-exceedance times, called K-gaps, and show that they follow the same exponential-point mass mixture distribution as the inter-exceedance times. The maximization of the likelihood built on this distribution yields a simple closed-form estimator for the extremal index. The method can admit covariates and can be applied with smoothing techniques, which allows its use in a nonstationary setting. Simulated and real data examples demonstrate the smooth estimation of the extremal index. The likelihood, based on an assumption of independence of the K-gaps, is misspecified whenever K is too small. This motivates another contribution of the thesis, the introduction into extreme-value statistics of misspecification tests based on the information matrix. For our likelihood, they are able to detect misspecification from any source, not only those due to a bad choice of the truncation parameter. They provide help also in threshold selection, and show whether the fundamental assumptions of stationarity or asymptotic independence are broken. Moreover, these diagnostic tests are of general use, and could be adapted to many kinds of extreme-value models, which are always approximate. Simulated examples demonstrate the performance of the misspecification tests in the context of extremal index estimation. Two data examples with complex behaviour, one univariate and the other bivariate, offer insight into their power in discovering situations where the fundamental assumptions of the likelihood model are not valid. In the multivariate case, the parameter corresponding to the univariate extremal index is the multivariate extremal index function. As in the univariate case, its appearance is linked to serial dependence in the observed processes. Univariate estimation methods can be applied, but are likely to give crude, unreasonably varying, estimates, and the constraints on the extremal index function implied by the characteristics of the stable tail dependence function are not automatically satisfied. The third contribution of the thesis is the development of methodology based on the M4 approximation of Smith and Weissman (1996), which can be used to estimate the multivariate extremal index, as well as other cluster characteristics. For this purpose, we give a preliminary cluster selection procedure, and approximate the noise on finite levels with a flexible semiparametric model, the Dirichlet mixtures used widely in Bayesian analysis. The model is fitted by the EM algorithm. Advantages and drawbacks of the method are discussed using the same univariate and bivariate examples as the likelihood methods.

EPFL2009