BAYEX: Spatiotemporal Bayesian hierarchical modeling of storm surge extremes with max-stable processes

BAYEX offers spatiotemporal Bayesian hierarchical modeling of storm surge extremes using max-stable and latent processes. As a key feature, BAYEX makes estimates of both the GEV parameters (location, scale and shape) and the annual maxima at any arbitrary location, either gauged or ungauged, while providing realistic uncertainty estimates. Inference in BAYEX is performed using Hamiltonian Monte Carlo as implemented by the Stan probabilistic programming language.This version 1.0 of the code was released to accompany the paper: Calafat, F. M., and M. Marcos (2020), Probabilistic reanalysis of storm surge extremes in Europe. Proc. Natl. Acad. Sci. U. S. A. Please cite that paper when using this code