Empirical Bayes inference for the block maxima method

Abstract

The block maxima method is one of the most popular approaches for extreme value analysis with independent and identically distributed observations in the domain of attraction of an extreme value distribution. The lack of a rigorous study on the Bayesian inference in this context has limited its use for statistical analysis of extremes. In this paper we propose an empirical Bayes procedure for inference on the block maxima law and its related quantities. We show that the posterior distributions of the tail index of the data distribution and of the return levels (representative of future extreme episodes) are consistent and asymptotically normal. These properties guarantee the reliability of posterior-based inference. We also establish contraction rates of the posterior predictive distribution, the key tool in Bayesian probabilistic forecasting. Posterior computations are readily obtained via an efficient adaptive Metropolis-Hasting type of algorithm. Simulations show its excellent inferential performances already with modest sample sizes. The utility of our proposal is showcased analysing extreme winds generated by hurricanes in Southeastern US.