The extremal index \(\theta\) is a measure of the degree of local dependence in the extremes of a stationary process. The exdex package performs frequentist inference about \(\theta\) using two types of methodology.
One type (Northrop, 2015) is based on a model that relates the distribution of block maxima to the marginal distribution of the data, leading to a semiparametric maxima estimator. Two versions of this type of estimator are provided, following Northrop, 2015 and Berghaus and Bücher, 2018. A slightly modified version of the latter is also provided. Estimates are produced using both disjoint and sliding block maxima, the latter providing greater precision of estimation. A graphical block size diagnostic is provided.
The other type of methodology uses a model for the distribution of threshold inter-exceedance times (Ferro and Segers, 2003). Three versions of this type of approach are provided: the iterated weight least squares approach of Süveges (2007), the \(K\)-gaps model of Süveges and Davison (2010) and a similar approach of Holesovsky and Fusek (2020) that we refer to as D-gaps. For the \(K\)-gaps and \(D\)-gaps models the exdex package allows missing values in the data, can accommodate independent subsets of data, such as monthly or seasonal time series from different years, and can incorporate information from censored inter-exceedance times. Graphical diagnostics for the threshold level and the respective tuning parameters \(K\) and \(D\) are provided.
The following code estimates the extremal index using the semiparametric maxima estimators, for an example dataset containing a time series of sea surges measured at Newlyn, Cornwall, UK over the period 1971-1976. The block size of 20 was chosen using a graphical diagnostic provided by choose_b().
Now we estimate \(\theta\) using the \(K\)-gaps model. The threshold \(u\) and runs parameter \(K\) were chosen using the graphical diagnostic provided by choose_uk().
To get the current released version from CRAN:
VignetteSee vignette("exdex-vignette", package = "exdex") for an overview of the package.
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4