Asymptotic theory for Bayesian inference and prediction: from the ordinary to a conditional Peaks-Over-Threshold method
Clément Dombry, Simone A. Padoan, Stefano Rizzelli
TL;DR
This work develops a rigorous Bayesian asymptotic theory for Peaks Over Threshold (POT) methods in univariate extremes and extends it to conditional extremes via a proportional tail model. It establishes posterior consistency, $\sqrt{k}$-contraction rates, a Bernstein–von Mises limit, and Wasserstein-consistent posterior predictive distributions, addressing the intrinsic GP misspecification and threshold-dependence challenges. For extremes conditional on covariates, it introduces a Dirichlet Process prior on the covariate distribution to model the scedasis function $c_0(x)$, proving contraction and BvM results and enabling predictive inference for conditional extreme quantiles through the posterior predictive distribution. Simulation studies and a financial-crisis application demonstrate good finite-sample performance, with accurate quantile forecasts and the ability to test covariate effects on extremes. Overall, the paper provides a principled, predictive Bayesian framework for both unconditional and covariate-adjusted extreme events, with practical implications for risk assessment and forecasting of future extreme occurrences.
Abstract
The Peaks Over Threshold (POT) method is the most popular statistical method for the analysis of univariate extremes. Even though there is a rich applied literature on Bayesian inference for the POT, the asymptotic theory for such proposals is missing. Even more importantly, the ambitious and challenging problem of predicting future extreme events according to a proper predictive statistical approach has received no attention to date. In this paper we fill this gap by developing the asymptotic theory of posterior distributions (consistency, contraction rates, asymptotic normality and asymptotic coverage of credible intervals) and prediction within the Bayesian framework in the POT context. We extend this asymptotic theory to account for cases where the focus is on the tail properties of the conditional distribution of a response variable given a vector of random covariates. To enable accurate predictions of extreme events more severe than those previously observed, we derive the posterior predictive distribution as an estimator of the conditional distribution of an out-of-sample random variable, given that it exceeds a sufficiently high threshold. We establish Wasserstein consistency of the posterior predictive distribution under both the unconditional and covariate-conditional approaches and derive its contraction rates. Simulations show the good performances of the proposed Bayesian inferential methods. The analysis of the change in the frequency of financial crises over time shows the utility of our methodology.
