Cláudia Neves

This paper addresses important weaknesses in current methodology for the estimation of multivariate extreme event distributions. The estimation of the residual dependence index η∈(0,1] is notoriously problematic. This motivates us to introduce a flexible class of gradient estimators for this parameter, designed to ameliorate the usual problems of threshold selection while furnishing a unifying approach to the usual marginal transforms. The efficacy of this semi-parametric gradient estimation of η originates from a hitherto neglected exponentially decaying term in the hidden regular variation characterisation. Numerical simulations aimed at assessing the performance for finite samples over a range of familiar copulas indicate an improved performance relatively to the existing estimators, notably the Hill’s. Paired with a reduced-bias procedure that we also construct, the gradient estimation offers much desired statistical guarantees as it further pinpoints the optimal number of intermediate data points to enter the estimation (in terms of minimising asymptotic mean squared error) in an automatic and demonstrable way. Our leading application illustrates how asymptotic independence can be discerned from monsoon-related rainfall occurrences at different locations in Ghana. The considerations involved in extending this framework to inference on the extreme value index drawing on max-domains of attraction are briefly discussed.

Unlike univariate extreme value theory, multivariate extreme value distributions cannot be specified through a finite-dimensional parameter family of distributions. Instead, the many facets of multivariate extremes are mirrored in the inherent dependence structure of component-wise maxima which must be dissociated from the limiting extreme behaviour of its marginal distribution functions before a probabilistic characterisation of an extreme value quality can be determined. Mechanisms applied to elicit extremal dependence typically rely on standardisation of the unknown marginal distribution functions from which pseudo-observations for either Pareto or Fr\'echet marginals result. The relative merits of both of these choices for transformation of marginals have been discussed in the literature, particularly in the context of domains of attraction of an extreme value distribution. This paper is set within this context of modelling penultimate dependence as it proposes a unifying class of estimators for the residual dependence index that eschews consideration of choice of marginals. In addition, a reduced bias variant of the new class of estimators is introduced and their asymptotic properties are developed. The pivotal role of the unifying marginal transform in effectively removing bias is borne by a comprehensive simulation study. The leading application in this paper comprises an analysis of asymptotic independence between rainfall occurrences originating from monsoon-related events at several locations in Ghana.

This article compares the modeling of nonstationary extreme events using parametric models with local parametric and semiparametric approaches also motivated by extreme value theory. Specifically, three estimators are compared based on (a) (local) semiparametric moment estimation, (b) (local) maximum likelihood estimation, and (c) spline‐based maximum likelihood estimation. Inference is performed in a sequential manner, highlighting the synergies between the different approaches to estimating extreme quantiles, including the T‐year level and right endpoint when finite. We present a novel heuristic to estimate nonstationary extreme value threshold with exceedances varying on a circular domain, and hypothesis‐testing procedures for identifying max‐domain of attraction in the nonstationary setting. Bootstrapping is used to estimate nonstationary confidence bounds throughout. We provide step‐by‐step guides for estimation, and explore the different inference strategies in application to directional modeling of hindcast storm peak significant wave heights recorded in the North Sea.

A new approach for evaluating time-trends in extreme values accounting also for spatial dependence is proposed. Based on exceedances over a space-time threshold, estimators for a trend function and for extreme value parameters are given, leading to a homogenization procedure for then applying stationary extreme value processes. Extremal dependence over space is further evaluated through variogram analysis including anisotropy. We detect significant inhomogeneities and trends in the extremal behaviour of daily precipitation data over a time period of 84 years and from 68 observational weather stations in North-West Germany. We observe that the trend is not monotonous over time in general. Asymptotic normality of the estimators under maximum domain of attraction conditions are proven.