madogram | R Documentation |
Computes a non-parametric estimate Pickands dependence function, A(w) for multivariate data, based on the madogram estimator.
madogram(w, data, margin = c("emp","est","exp","frechet","gumbel"))
w |
(m x d) design matrix (see Details). |
data |
(n x d) matrix of data or data frame with |
margin |
string, denoting the type marginal distributions ( |
The estimation procedure is based on the madogram as proposed in Marcon et al. (2017). The madogram is defined by
ν(w) = E[max_{i=1,...,d}(F^{1/w_i}_{i}(X_{i})) - 1/d ∑_{i=1,...,d}F^{1/w_i}_{i}(X_{i})],
where 0<w_i<1 and w_d=1-(w_1+...+w_{d-1}).
Each row of the design matrix w
is a point in the unit
d
-dimensional simplex.
If X is a d
-dimensional max-stable distributed random vector, with exponent measure function V(x) and Pickands dependence function A(w), then
ν(w) = V(1/w_1,...,1/w_d)/(1+V(1/w_1,...,1/w_d))-c(w), where c(w)=d^{-1}∑_{i=1}^{d}{w_i/(1+w_i)}.
From this, it follows that
V(1/w_1,...,1/w_d)=(ν(w)+c(w))/(1-ν(w)-c(w)),
and
A(w)=(ν(w)+c(w))/(1-ν(w)-c(w)).
An empirical transformation of the marginals is performed when margin="emp"
.
A max-likelihood fitting of the GEV distributions is implemented when margin="est"
.
Otherwise it refers to marginal parametric GEV theorethical distributions (margin="exp", "frechet", "gumbel"
).
A numeric vector of estimates.
Simone Padoan, simone.padoan@unibocconi.it, https://faculty.unibocconi.it/simonepadoan/; Boris Beranger, borisberanger@gmail.com https://www.borisberanger.com; Giulia Marcon, giuliamarcongm@gmail.com
Marcon, G., Padoan, S.A., Naveau, P., Muliere, P. and Segers, J. (2017) Multivariate Nonparametric Estimation of the Pickands Dependence Function using Bernstein Polynomials. Journal of Statistical Planning and Inference, 183, 1-17.
Naveau, P., Guillou, A., Cooley, D., Diebolt, J. (2009) Modelling pairwise dependence of maxima in space, Biometrika, 96(1), 1-17.
beed
, beed.confband
x <- simplex(2) data <- evd::rbvevd(50, dep = 0.4, model = "log", mar1 = c(1,1,1)) Amd <- madogram(x, data, "emp") Amd.bp <- beed(data, x, 2, "md", "emp", 20, plot=TRUE) lines(x[,1], Amd, lty = 1, col = 2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.