amvevd | R Documentation |
Calculate the dependence function A
for the multivariate
logistic and multivariate asymmetric logistic models; plot the
estimated function in the trivariate case.
amvevd(x = rep(1/d,d), dep, asy, model = c("log", "alog"), d = 3, plot =
FALSE, col = heat.colors(12), blty = 0, grid = if(blty) 150 else 50,
lower = 1/3, ord = 1:3, lab = as.character(1:3), lcex = 1)
x |
A vector of length |
dep |
The dependence parameter(s). For the logistic model,
should be a single value. For the asymmetric logistic model,
should be a vector of length |
asy |
The asymmetry parameters for the asymmetric logistic
model. Should be a list with |
model |
The specified model; a character string. Must be
either |
d |
The dimension; an integer greater than or equal to two.
The trivariate case |
plot |
Logical; if |
col |
A list of colours (see |
blty |
The border line type, for the border that surrounds
the triangular image. By default |
grid |
For plotting, the function is evaluated at |
lower |
The minimum value for which colours are plotted. By
defualt |
ord |
A vector of length three, which should be a permutation
of the set |
lab |
A character vector of length three, in which case the
|
lcex |
A numerical value giving the amount by which the
labels should be scaled relative to the default. Ignored
if |
Let z = (z_1,\dots,z_d)
and
w = (w_1,\dots,w_d)
.
Any multivariate extreme value distribution can be written as
G(z) = \exp\left\{- \left\{\sum\nolimits_{j=1}^{d} y_j
\right\} A\left(\frac{y_1}{\sum\nolimits_{j=1}^{d} y_j}, \dots,
\frac{y_d}{\sum\nolimits_{j=1}^{d} y_j}\right)\right\}
for some function A
defined on the simplex
S_d = \{w \in R^d_+ : \sum\nolimits_{j=1}^{d}
w_j = 1\}
,
where
y_i = \{1+s_i(z_i-a_i)/b_i\}^{-1/s_i}
for 1+s_i(z_i-a_i)/b_i > 0
and
i = 1,\dots,d
, and where the (generalized extreme value)
marginal parameters are given by
(a_i,b_i,s_i)
, b_i > 0
.
If s_i = 0
then y_i
is defined by
continuity.
A
is called (by some authors) the dependence function.
It follows that A(w) = 1
when w
is one of the d
vertices of S_d
, and that A
is a convex function with
\max(w_1,\dots,w_d) \leq A(w)\leq 1
for
all w
in S_d
.
The lower and upper limits of A
are obtained under complete
dependence and mutual independence respectively.
A
does not depend on the marginal parameters.
A numeric vector of values. If plotting, the smallest evaluated function value is returned invisibly.
amvnonpar
, abvevd
,
rmvevd
, image
amvevd(dep = 0.5, model = "log")
s3pts <- matrix(rexp(30), nrow = 10, ncol = 3)
s3pts <- s3pts/rowSums(s3pts)
amvevd(s3pts, dep = 0.5, model = "log")
## Not run: amvevd(dep = 0.05, model = "log", plot = TRUE, blty = 1)
amvevd(dep = 0.95, model = "log", plot = TRUE, lower = 0.94)
asy <- list(.4, .1, .6, c(.3,.2), c(.1,.1), c(.4,.1), c(.2,.3,.2))
amvevd(s3pts, dep = 0.15, asy = asy, model = "alog")
amvevd(dep = 0.15, asy = asy, model = "al", plot = TRUE, lower = 0.7)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.