Tail persistence of extremes at finite and asymptotic levels with Bayesian hierarchical inference
Abstract
Tail persistence characterizes the relative escalation of observations beyond an extreme threshold and complements conventional summaries such as return levels and tail indices. Existing literature is mainly distribution-based and does not account for uncertainty, nonstationarity, or covariate effects. We develop a Bayesian framework for modeling tail persistence under generalized extreme value and generalized Pareto models. We first make explicit the connection between tail persistence and the extreme value index under regular variation, which identifies the asymptotic persistence with 1/(1−ξ) and links persistence to tail heaviness. We then derive closed-form expressions for persistence at finite quantile levels under both models and show that, for any location--scale extreme value family with a common shape parameter, finite-level persistence depends on the location and scale parameters only through their ratio. Motivated by these results and by the weak identifiability of observation-specific shape parameters, we formulate a Bayesian hierarchical extreme value model in which the shape parameter is constant over time, while the location and scale parameters vary with covariates and latent effects. The asymptotic persistence is then a single, well-identified quantity, whereas finite-level persistence evolves over time through the location--scale structure. The framework provides posterior inference for both quantities, with credible intervals and probabilistic assessment of temporal evolution and covariate effects. Through simulation studies and real-data applications, we show that persistence provides an interpretable characterization of the escalation of extremes beyond high thresholds, complementing classical extreme value summaries with comprehensive uncertainty quantification in nonstationary settings. The methodology is implemented in an open-source \texttt{R} package.
Related
- A Bayesian hierarchical model for preferentially sampled spatial extremes data: an application to air quality
- Behavioral Mapping of Hierarchical Dynamical Models Accelerates Bayesian Parameter Inference
- A Hierarchical Bayesian Cox Model with Variational Inference for High-Dimensional Survival Analysis
- Bayesian model averaging of risk set probabilities using a geometric representation of multivariate extremes
- A Geometric Framework for the Bayesian Estimation of Joint Environmental Extremes
BibTeX
@misc{rawat2026tail,
title = {{Tail persistence of extremes at finite and asymptotic levels with Bayesian hierarchical inference}},
author = {Nikita Rawat and Shiv Kumar Yadav and Dreamlee Sharma and Rishikesh Yadav},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.03288},
eprint = {2611.03288},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.03288},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Tail persistence of extremes at finite and asymptotic levels with Bayesian hierarchical inference
AU - Nikita Rawat
AU - Shiv Kumar Yadav
AU - Dreamlee Sharma
AU - Rishikesh Yadav
PY - 2026
DA - 2026/11/02
PB - arXiv
DO - 10.48550/arxiv.2611.03288
UR - https://arxiv.org/abs/2611.03288
AB - Tail persistence characterizes the relative escalation of observations beyond an extreme threshold and complements conventional summaries such as return levels and tail indices. Existing literature is mainly distribution-based and does not account for uncertainty, nonstationarity, or covariate effects. We develop a Bayesian framework for modeling tail persistence under generalized extreme value and generalized Pareto models. We first make explicit the connection between tail persistence and the extreme value index under regular variation, which identifies the asymptotic persistence with 1/(1−ξ) and links persistence to tail heaviness. We then derive closed-form expressions for persistence at finite quantile levels under both models and show that, for any location--scale extreme value family with a common shape parameter, finite-level persistence depends on the location and scale parameters only through their ratio. Motivated by these results and by the weak identifiability of observation-specific shape parameters, we formulate a Bayesian hierarchical extreme value model in which the shape parameter is constant over time, while the location and scale parameters vary with covariates and latent effects. The asymptotic persistence is then a single, well-identified quantity, whereas finite-level persistence evolves over time through the location--scale structure. The framework provides posterior inference for both quantities, with credible intervals and probabilistic assessment of temporal evolution and covariate effects. Through simulation studies and real-data applications, we show that persistence provides an interpretable characterization of the escalation of extremes beyond high thresholds, complementing classical extreme value summaries with comprehensive uncertainty quantification in nonstationary settings. The methodology is implemented in an open-source \texttt{R} package.
ER -