This article compares the modeling of nonstationary extreme events using parametric models with local parametric and semiparametric approaches also motivated by extreme value theory. Specifically, three estimators are compared based on (a) (local) semiparametric moment estimation, (b) (local) maximum likelihood estimation, and (c) spline‐based maximum likelihood estimation. Inference is performed in a sequential manner, highlighting the synergies between the different approaches to estimating extreme quantiles, including the T‐year level and right endpoint when finite. We present a novel heuristic to estimate nonstationary extreme value threshold with exceedances varying on a circular domain, and hypothesis‐testing procedures for identifying max‐domain of attraction in the nonstationary setting. Bootstrapping is used to estimate nonstationary confidence bounds throughout. We provide step‐by‐step guides for estimation, and explore the different inference strategies in application to directional modeling of hindcast storm peak significant wave heights recorded in the North Sea.
We present and describe the GPFDA package for R. The package provides flexible functionalities for dealing with Gaussian process regression (GPR) models for functional data. Multivariate functional data, functional data with multidimensional inputs, and nonseparable and/or nonstationary covariance structures can be modeled. In addition, the package fits functional regression models where the mean function depends on scalar and/or functional covariates and the covariance structure is modeled by a GPR model. In this paper, we present the versatility of GPFDA with respect to mean function and covariance function specifications and illustrate the implementation of estimation and prediction of some models through reproducible numerical examples.
Models for extreme values accommodating non-stationarity have been amply studied and evaluated from a parametric perspective. Whilst these models are flexible, in the sense that many parametrizations can be explored, they assume an asymptotic distribution as the proper fit to observations from the tail. This paper provides a holistic approach to the modelling of nonstationary extreme events by iterating between parametric and semi-parametric approaches, thus providing an automatic procedure to estimate a moving threshold with respect to a periodic covariate in circular data. By exploiting advantages and mitigating pitfalls of each approach, a unified framework is provided as the backbone for estimating extreme quantiles, including that of the T -year level and finite right endpoint, which seeks to optimize bias-variance trade-off. To this end, two tuning parameters related to the spread of peaks over threshold are introduced. We provide guidance for applying the methodology to the directional modelling of hindcast storm peak significant wave heights recorded in the North Sea. Although the theoretical underpinning for adaptation of well-known estimators in statistics of extremes to circular data is given in some detail, the derivation of their asymptotic properties lays beyond the scope of this paper. A bootstrap technique is implemented for obtaining direction-driven confidence bounds in such a way as to account for the relevant boundary restrictions with minimal sensitivity to initial point. This provides a template for other applications where the analysis of directional extremes is of importance.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.