Learn R Programming

ExtremalDep: Extremal Dependence Models

A set of procedures for parametric and non-parametric modelling of the dependence structure of multivariate extreme-values is provided. The statistical inference is performed with non-parametric estimators, likelihood-based estimators and Bayesian techniques. It adapts the methodologies of Beranger and Padoan (2015) doi:10.48550/arXiv.1508.05561, Marcon et al. (2016) doi:10.1214/16-EJS1162, Marcon et al. (2017) doi:10.1002/sta4.145, Marcon et al. (2017) doi:10.1016/j.jspi.2016.10.004 and Beranger et al. (2021) doi:10.1007/s10687-019-00364-0. This package also allows for the modelling of spatial extremes using flexible max-stable processes. It provides simulation algorithms and fitting procedures relying on the Stephenson-Tawn likelihood as per Beranger at al. (2021) doi:10.1007/s10687-020-00376-1.

Version 0.0.4-2 updates:

  • Fixes a typo in dExtDep() for in the Asymetric Logistic model;
  • Replaces the Calloc() and Free() calls in the .C files by the R_* prefixed counterparts since STRICT_R_HEADERS=1 becomes the default with R 4.5.0;
  • Improves some entries in the manual.

Version 0.0.4-3 updates:

  • Inclusion of the PAMfmado() function written by Philippe Naveau.

Version 0.0.4-4 updates:

  • Fixes resetting graphical parameters in plot_ExtDep.np() when type = "Qsets";
  • Replaces closeAllConnections() by stopCluster() in fExtDepSpat().

Version 0.0.4-5 updates:

  • Improvement to dExtDep() function when model="HR" and "ET";
  • New lambda.hr() function that can be used to define the parameters of the trivariate Husler-Reiss model. Given two parameters, a range for the third parameter is provided to ensure positive definite matrices in the exponent function;
  • Improve efficiency when manipulating matrices:
    • t(A) %*% B is replaced by crossprod(A,B) (and vice versa);
    • solve(A) is replaced by chol2inv(chol(A));
    • t(x) %*% solve(A) %*% x is replaced by sum(forwardsolve(t(chol(A)),x));
  • Define outputs of fExtDep(), fExtDep.np() and fExtDepSpat() to be of class ExtDep_Freq, ExtDep_Bayes, ExtDep_npFreq, ExtDep_npBayes, ExtDep_npEmp and ExtDep_Spat;
  • Make use of class methods for plot() and summary() functions replacing the previous plot.ExtDep(), etc;
  • New est(), StdErr(), logLik(), bic(), tic() to extract outputs from objects of class ExtDep_Freq, ExtDep_Bayes and ExtDep_Spat.

Copy Link

Version

Install

install.packages('ExtremalDep')

Monthly Downloads

301

Version

0.0.4-5

License

GPL (>= 2)

Maintainer

Simone Padoan

Last Published

June 29th, 2025

Functions in ExtremalDep (0.0.4-5)

beed.boot

Bootstrap Resampling and Bernstein Estimation of Extremal Dependence
dim_ExtDep

Dimensions calculations for parametric extremal dependence models
desn

Univariate extended skew-normal distribution
bic

Extract the Bayesian Information Criterion
beed.confband

Nonparametric Bootstrap Confidence Intervals
dest

Univariate extended skew-t distribution
dmesn

Bivariate and trivariate extended skew-normal distribution
dExtDep

Parametric and non-parametric density of Extremal Dependence
diagnostics

Diagnostics plots for MCMC algorithm.
dGEV

The Generalized Extreme Value Distribution
fExtDep.np

Non-parametric extremal dependence estimation
est

Extract the estimated parameter
dmest

Bivariate and trivariate extended skew-t distribution
index.ExtDep

Index of extremal dependence
fGEV

Fitting of the Generalized Extreme Value Distribution
ellipse

Level sets for bivariate normal, student-t and skew-normal distributions probability densities.
heat

Summer temperature maxima in Melbourne, Australia between 1961 and 2010.
fExtDep

Extremal dependence estimation
lambda.HR

Valid set of parameters for the 3-dimensional Husler-Reiss model.
fExtDepSpat

Fitting of a max-stable process
model

Extract the model attribute
method

Extract the method attribute
logReturns

Monthly maxima of log-return exchange rates of the Pound Sterling (GBP) against the US dollar (USD) and the Japanese yen (JPY), between March 1991 and December 2014.
rExtDepSpat

Random generation of max-stable processes
pExtDep

Parametric and non-parametric distribution function of Extremal Dependence
pFailure

Probability of falling into a failure region
tic

Extract the Takeuchi Information Criterion
simplex

Definition of a multivariate simplex
trans2GEV

Transformation to GEV distribution
pollution

Air quality datasets containing daily maxima of air pollutants (PM10, NO, NO2, 03 and S02) recorded in Leeds (U.K.), during five winter seasons (November-Februrary) between 1994 and 1998.
pickands.plot

Plot for the Pickands dependence function.
madogram

Madogram-based estimation of the Pickands Dependence Function
returns

Compute return values
rExtDep

Parametric and semi-parametric random generator of extreme events
trans2UFrechet

Transformation to unit Frechet distribution
returns.plot

Plot some return levels.
angular.plot

Angular density plot.
angular

Estimation of the angular density, angular measure and random generation from the angular distribution.
PAMfmado

Clustering of maxima
Wind

Weekly maximum wind speed data collected over 4 stations across Oklahoma, USA, over the March-May preiod between 1996 and 2012.
PrecipFrance

Weekly maxima of hourly rainfall in France
WindSpeedGust

Hourly wind gust, wind speed and air pressure at Lingen (GER), Ossendorf (GER) and Parcay-Meslay (FRA).
beed

Bernstein Polynomials Based Estimation of Extremal Dependence
MilanPollution

Pollution data for summer and winter months in Milan, Italy.
ExtQ

Univariate Extreme Quantile
StdErr

Extract the standard errors of estimated parameters