Truncated pair-wise likelihood for the Brown-Resnick process with applications to maximum temperature data

Max-stable processes are a natural extension of multivariate extreme value theory important for modeling the spatial dependence of environmental extremes. Inference for max-stable processes observed at several spatial locations is challenging due to the intractability of the full likelihood function. Composite likelihood methods avoid these difficulties by combining a number of low-dimensional likelihood objects, typically defined on pairs or triplets of spatial locations. This work develops a new truncation procedure based on ℓ1-penalty to reduce the complexity and computational cost associated with the composite likelihood function. The new method is shown to achieve a favorable trade-off between computational burden and statistical efficiency by dropping a number of noisy sub-likelihoods. The properties of the new method are illustrated through numerical simulations and an application to real extreme temperature data.

[1]  S. Padoan,et al.  Likelihood-Based Inference for Max-Stable Processes , 2009, 0902.3060.

[2]  Bruce G. Lindsay,et al.  ISSUES AND STRATEGIES IN THE SELECTION OF COMPOSITE LIKELIHOODS , 2011 .

[3]  Marc G. Genton,et al.  Tapered composite likelihood for spatial max-stable models , 2014 .

[4]  Harry Joe,et al.  Composite Likelihood Methods , 2012 .

[5]  C. Heyde Quasi-likelihood and its application : a general approach to optimal parameter estimation , 1998 .

[6]  A. Davison,et al.  Composite likelihood estimation for the Brown–Resnick process , 2013 .

[7]  Anthony C. Davison,et al.  Bayesian inference for the Brown-Resnick process, with an application to extreme low temperatures , 2015, 1506.07836.

[8]  Stefano Castruccio,et al.  High-Order Composite Likelihood Inference for Max-Stable Distributions and Processes , 2014, 1411.0086.

[9]  Marco Oesting,et al.  Bayesian inference for multivariate extreme value distributions , 2016, 1611.05602.

[10]  Marco Oesting,et al.  Asymptotic properties of the maximum likelihood estimator for multivariate extreme value distributions , 2016, 1612.05178.

[11]  Marc G. Genton,et al.  On the likelihood function of Gaussian max-stable processes , 2011 .

[12]  D. Cox,et al.  A note on pseudolikelihood constructed from marginal densities , 2004 .

[13]  Sidney I. Resnick,et al.  Extreme values of independent stochastic processes , 1977 .

[14]  Laurens de Haan,et al.  Stationary max-stable fields associated to negative definite functions. , 2008, 0806.2780.

[15]  Anthony C. Davison,et al.  Statistics of Extremes , 2015 .

[16]  P. Naveau,et al.  Variograms for spatial max-stable random fields , 2006 .

[17]  N. Reid,et al.  AN OVERVIEW OF COMPOSITE LIKELIHOOD METHODS , 2011 .

[18]  Marc G. Genton,et al.  Full likelihood inference for max‐stable data , 2017, Stat.

[19]  A. Davison,et al.  Statistical Modeling of Spatial Extremes , 2012, 1208.3378.

[20]  Jonathan A. Tawn,et al.  A dependence measure for multivariate and spatial extreme values: Properties and inference , 2003 .

[21]  Jianqing Fan,et al.  A Selective Overview of Variable Selection in High Dimensional Feature Space. , 2009, Statistica Sinica.

[22]  R. Tibshirani,et al.  REJOINDER TO "LEAST ANGLE REGRESSION" BY EFRON ET AL. , 2004, math/0406474.