Bayesian Dirichlet mixture model for multivariate extremes: a re-parametrization

Anne Sabourin (1, 2) , Philippe Naveau (2, 3)

1 TSI - Département Traitement du Signal et des Images
2 LSCE - Laboratoire des Sciences du Climat et de l'Environnement [Gif-sur-Yvette]
3 ESTIMR - Extrèmes : Statistiques, Impacts et Régionalisation

The probabilistic framework of extreme value theory is well-known: the dependence structure of large events is characterized by an angular measure on the positive orthant of the unit sphere. The family of these angular measures is non-parametric by nature. Nonetheless, any angular measure may be approached arbitrarily well by a mixture of Dirichlet distributions. The semi-parametric Dirichlet mixture model for angular measures is theoretically valid in arbitrary dimension, but the original parametrization is subject to a moment constraint making Bayesian inference very challenging in dimension greater than three. A new unconstrained parametrization is proposed. This allows for a natural prior specification as well as a simple implementation of a reversible-jump MCMC. Posterior consistency and ergodicity of the Markov chain are verified and the algorithm is tested up to dimension five. In this non identifiable setting, convergence monitoring is performed by integrating the sampled angular densities against Dirichlet test functions.

Domaines: Applications [stat.AP]; Calcul [stat.CO]; Méthodologie [stat.ME]; Statistiques [math.ST]; Théorie [stat.TH]

Anne Sabourin, Philippe Naveau. Bayesian Dirichlet mixture model for multivariate extremes: a re-parametrization. Computational Statistics and Data Analysis, 2013, ⟨10.1016/j.csda.2013.04.021⟩. ⟨hal-00880879⟩