Abstract:The generalized extreme value (GEV) distribution is a popular model for analyzing and forecasting extreme weather data. To increase prediction accuracy, spatial information is often pooled via a latent Gaussian process (GP) on the GEV parameters. Inference for GEV‐GP models is typically carried out using Markov Chain Monte Carlo (MCMC) methods, or using approximate inference methods such as the integrated nested Laplace approximation (INLA). However, MCMC becomes prohibitively slow as the number of spatial locations increases, whereas INLA is applicable in practice only to a limited subset of GEV‐GP models. In this article, we revisit the original Laplace approximation for fitting spatial GEV models. In combination with a popular sparsity‐inducing spatial covariance approximation technique, we show through simulations that our approach accurately estimates the Bayesian predictive distribution of extreme weather events, is scalable to several thousand spatial locations, and is several orders of magnitude faster than MCMC. A case study in forecasting extreme snowfall across Canada is presented. Résumé La distribution de valeurs extrêmes généralisée (GEV) est couramment utilisée pour l'analyse et la prévision des phénomènes météorologiques extrêmes. Pour améliorer la précision prédictive, les informations spatiales sont souvent regroupées via un processus gaussien latent (GP) appliqué aux paramètres GEV. L'inférence des modèles GEV‐GP utilise généralement les méthodes de Monte‐Carlo par chaînes de Markov (MCMC) ou des approximations comme l'approximation de Laplace imbriquée intégrée (INLA). Toutefois, les méthodes MCMC deviennent très lentes quand le nombre de sites spaciaux augmente, et l'INLA ne s'applique qu'à un nombre restreint de modèles GEV‐GP. Dans cet article, les auteurs réexaminent l'approximation de Laplace classique pour ajuster les modèles GEV spatiaux. En la combinant à une technique d'approximation de la covariance spatiale induisant la parcimonie, ils montrent par simulation que leur approche estime avec précision la distribution prédictive bayésienne des événements météorologiques extrêmes, qu'elle s'adapte à plusieurs milliers de sites spatiaux, et qu'elle est nettement plus rapide que les méthodes MCMC. L'article présente également une étude de cas portant sur la prévision des chutes de neige extrêmes au Canada.

Clustering spatial functional data using a geographically weighted Dirichlet process

Bayesian statistics in spatial epidemiology]

A Bayesian clustering ensemble Gaussian process model for network-wide traffic flow clustering and prediction

Bayesian cluster geographically weighted regression for spatial heterogeneous data

Clustering Spatially Correlated Functional Data with Multiple Scalar Covariates.

Bayesian Function-on-Function Regression for Spatial Functional Data

Spatial non-parametric Bayesian clustered coefficients

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Spatially Clustered Regression

Bayesian cluster analysis for registration and clustering homogeneous subgroups in multidimensional functional data

A Bayesian approach to disease clustering using restricted Chinese restaurant processes

Powered Dirichlet Process for Controlling the Importance of "Rich-Get-Richer" Prior Assumptions in Bayesian Clustering

Geographically Weighted Regression Analysis for Spatial Economics Data: a Bayesian Recourse

Fast and scalable inference for spatial extreme value models

Bayesian modeling and clustering for spatio-temporal areal data: An application to Italian unemployment

Generalized Bayesian nonparametric clustering framework for high-dimensional spatial omics data

Combining heterogeneous spatial datasets with process-based spatial fusion models: A unifying framework

A Regularized Spatial Market Segmentation Method with Dirichlet Process Gaussian Mixture Prior

On clustering uncertain and structured data with Wasserstein barycenters and a geodesic criterion for the number of clusters

Bayesian Spatial Homogeneity Pursuit for Survival Data with an Application to the SEER Respiratory Cancer Data

Spatiotemporal clustering using Gaussian processes embedded in a mixture model