ExtremeRisks (version 0.0.4-1)

extMultiQuantile: Multidimensional Value-at-Risk (VaR) or Extreme Quantile (EQ) Estimation

Description

Computes point estimates and \((1-\alpha)100\%\) confidence regions for d-dimensional VaR based on the Weissman estimator.

Usage

extMultiQuantile(data, tau, tau1, var=FALSE, varType="asym-Ind-Log", bias=FALSE,
                 k=NULL, alpha=0.05, plot=FALSE)

Value

A list with elements:

• ExtQHat: an estimate of the d-dimensional VaR or \(\tau'_n\)-th d-dimensional quantile;

• VarCovExQHat: an estimate of the asymptotic variance-covariance of the d-dimensional VaR estimator;

• EstConReg: an estimate of the approximate \((1-\alpha)100\%\) confidence regions for the d-dimensional VaR.

Arguments

data

A matrix of \((n \times d)\) observations.

tau

A real in \((0,1)\) specifying the intermediate level \(\tau_n\). See Details.

tau1

A real in \((0,1)\) specifying the extreme level \(\tau'_n\). See Details.

var

If var=TRUE then an estimate of the asymptotic variance-covariance matrix of the d-dimensional VaR estimator is computed.

varType

A string specifying the type of asymptotic variance-covariance matrix to compute. By default varType="asym-Ind-Log" specifies that the variance-covariance matrix is obtained assuming dependent variables and exploiting the logarithmic scale. See Details.

bias

A logical value. By default biast=FALSE specifies that no bias correction is computed. See Details.

k

An integer specifying the value of the intermediate sequence \(k_n\). See Details.

alpha

A real in \((0,1)\) specifying the confidence level \((1-\alpha)100\%\) of the approximate confidence region for the d-dimensional VaR.

plot

A logical value. By default plot=FALSE specifies that no graphical representation of the estimates is not provided. See Details.

Author

Simone Padoan, , https://faculty.unibocconi.it/simonepadoan/; Gilles Stupfler, , https://math.univ-angers.fr/~stupfler/

Details

For a dataset data of d-dimensional observations and sample size \(n\), the VaR or EQ, correspoding to the extreme level tau1, is computed by applying the d-dimensional Weissman estimator. The definition of the Weissman estimator depends on the estimation of the d-dimensional tail index \(\gamma\). Here, \(\gamma\) is estimated using the Hill estimation (see MultiHTailIndex). The data are regarded as d-dimensional temporal independent observations coming from dependent variables. See Padoan and Stupfler (2020) for details.

• The so-called intermediate level tau or \(\tau_n\) is a sequence of positive reals such that \(\tau_n \to 1\) as \(n \to \infty\). Practically, for each variable, \((1-\tau_n)\in (0,1)\) is a small proportion of observations in the observed data sample that exceed the \(tau_n\)-th empirical quantile. Such proportion of observations is used to estimate the individual \(tau_n\)-th quantile and tail index \(\gamma\).

• The so-called extreme level tau1 or \(\tau'_n\) is a sequence of positive reals such that \(\tau'_n \to 1\) as \(n \to \infty\). For each variable, the value \((1-tau'_n) \in (0,1)\) is meant to be a small tail probability such that \((1-\tau'_n)=1/n\) or \((1-\tau'_n) < 1/n\). It is also assumed that \(n(1-\tau'_n) \to C\) as \(n \to \infty\), where \(C\) is a positive finite constant. The value \(C\) is the expected number of exceedances of the individual \(\tau'_n\)-th quantile. Typically, \(C \in (0,1)\) which means that it is expected that there are no observations in a data sample exceeding the individual quantile of level \((1-\tau_n')\).

• If var=TRUE then an estimate of the asymptotic variance-covariance matrix of the \(tau'_n\)-th d-dimensional quantile is computed. The data are regarded as temporal independent observations coming from dependent variables. The asymptotic variance-covariance matrix is estimated exploiting the formula in Section 5 of Padoan and Stupfler (2020). In particular, the variance-covariance matrix is computed exploiting the asymptotic behaviour of the normalized quantile estimator which is expressed in logarithmic scale. This is achieved through varType="asym-Ind-Log". If varType="asym-Ind" then the variance-covariance matrix is computed exploiting the asymptotic behaviour of the d-dimensional relative quantile estimator appropriately normalized (see formula in Section 5 of Padoan and Stupfler (2020)).

• If bias=TRUE then an estimate of each individual \(\tau'_n\)--th quantile is estimated using the formula in page 330 of de Haan et al. (2016), which provides a bias corrected version of the Weissman estimator. However, in this case the asymptotic variance is not estimated using the formula in Haan et al. (2016) Theorem 4.2. For simplicity standard the variance-covariance matrix is still computed using formula in Section 5 of Padoan and Stupfler (2020).

• k or \(k_n\) is the value of the so-called intermediate sequence \(k_n\), \(n=1,2,\ldots\). Its represents a sequence of positive integers such that \(k_n \to \infty\) and \(k_n/n \to 0\) as \(n \to \infty\). Practically, for each marginal distribution, the value \(k_n\) specifies the number of k\(+1\) larger order statistics to be used to estimate the individual \(\tau_n\)-th empirical quantile and individual tail index \(\gamma_j\) for \(j=1,\ldots,d\). The intermediate level \(\tau_n\) can be seen defined as \(\tau_n=1-k_n/n\).

• Given a small value \(\alpha\in (0,1)\) then an estimate of an asymptotic confidence region for \(tau'_n\)-th d-dimensional quantile, with approximate nominal confidence level \((1-\alpha)100\%\), is computed. The confidence regions are computed exploiting the asymptotic behaviour of the normalized quantile estimator in logarithmic scale. This is an "asymmetric" region and it is achieved through varType="asym-Ind-Log". A "symmetric" region is obtained exploiting the asymptotic behaviour of the relative quantile estimator appropriately normalized, see formula in Section 5 of Padoan and Stupfler (2020). This is achieved through varType="asym-Ind".

• If plot=TRUE then a graphical representation of the estimates is not provided.

References

Simone A. Padoan and Gilles Stupfler (2022). Joint inference on extreme expectiles for multivariate heavy-tailed distributions, Bernoulli 28(2), 1021-1048.

de Haan, L., Mercadier, C. and Zhou, C. (2016). Adapting extreme value statistics to financial time series: dealing with bias and serial dependence. Finance and Stochastics, 20, 321-354.

de Haan, L. and Ferreira, A. (2006). Extreme Value Theory: An Introduction. Springer-Verlag, New York.

See Also

MultiHTailIndex, estMultiExpectiles, predMultiExpectiles

Examples

# Extreme quantile estimation at the extreme level tau1 obtained with
# d-dimensional observations simulated from a joint distribution with
# a Gumbel copula and equal Frechet marginal distributions.
library(plot3D)
library(copula)
library(evd)

# distributional setting
copula <- "Gumbel"
dist <- "Frechet"

# parameter setting
dep <- 3
dim <- 3
scale <- rep(1, dim)
shape <- rep(3, dim)
par <- list(dep=dep, scale=scale, shape=shape, dim=dim)

# Intermediate level (or sample tail probability 1-tau)
tau <- 0.95
# Extreme level (or tail probability 1-tau1 of unobserved quantile)
tau1 <- 0.9995

# sample size
ndata <- 1000

# Simulates a sample from a multivariate distribution with equal Frechet
# marginals distributions and a Gumbel copula
data <- rmdata(ndata, dist, copula, par)
scatter3D(data[,1], data[,2], data[,3])

# High d-dimensional expectile (intermediate level) estimation
extQHat <- extMultiQuantile(data, tau, tau1, TRUE)

extQHat$ExtQHat
extQHat$VarCovExQHat
# run the following command to see the graphical representation
# \donttest{
 extQHat <- extMultiQuantile(data, tau, tau1, TRUE, plot=TRUE)
# }