Modelling Dependence Uncertainty in the Extremes of Markov Chains

Abstract

General theory on the extremes of stationary processes leads only to a limited representation for extreme-state behaviour, usually summarised by the extremal index. In practice this means that other quantities such as the duration of extreme episodes or aggregate of threshold exceedances within a cluster require stronger model assumptions. In this paper we propose a model based on a Markov assumption for the underlying process, with high-level transitions determined by an asymptotically motivated distribution. This idea is not new: Smith et al. (1997) first developed the statistical basis for such a procedure, which was subsequently extended by Bortot and Tawn (1998) to better handle the case of weak extremal temporal dependence for which the extremal index is unity. We adopt similar procedures to each of these earlier works, but suggest a different model for the Markov transitions. The model we use was developed by Coles and Pauli (2002) to enable a Bayesian inference of multivariate extremes that provides a posterior distribution on the status of asymptotic independence. By adopting this model in the Markov framework, we show here that the model has all the flexibility of the model developed by Bortot and Tawn (1998), but with the additional advantage of providing a posterior probability on the extremal index and inferences that take full account of the uncertainty in the extremal index. We demonstrate the methodology on both simulated data and a time series of daily rainfall that exhibit weak temporal dependence at extreme levels.