composite1COP {copBasic} | R Documentation |
Composition of a Single Symmetric Copula with Two Compositing Parameters
Description
The composition of a single copula (Salvadori et al., 2006, p. 266, prop. C.3) is created by the following result related to “composition of copulas” in that reference. Suppose is a symmetric copula (see
COP
) with parameters and
(for
see
P
), then a family of generally asymmetric copulas with two compositing parameters
, and
, which also includes just the copula
as a limiting case for
and is given by
The composite1COP
function provides the means for inserting permutation asymmetry from a permutation symmetric copula as described by Joe (2017, p. 124), but do so in a more general way through the provision of two and not just one parameter. Joe's description is supported herein if one of the or
is held at zero. Very loosely, the
kicks probability density down towards the lower right corner, whereas
kicks density up towards the upper left corner. Finally, the
composite2COP
function is based on a slighty more general result (see composite2COP
for further details of copula composition).
Usage
composite1COP(u, v, para, ...)
Arguments
u |
Nonexceedance probability |
v |
Nonexceedance probability |
para |
A special parameter |
... |
Additional arguments to pass to the copula. |
Value
Value(s) for the composited copula are returned.
Note
The following descriptions list in detail the structure and content of the para
argument:
alpha
— The
compositing parameter;
beta
— The
compositing parameter;
cop1
— Function of the copula
; and
para1
— Vector of parameters
for
.
For the para
argument, the same nomenclature as used for composite2COP
is used with obviously cop2
and para2
dropped for composite1COP
. The cop1
and para1
names remain enumerated for composite1COP
so that the para
argument of the more general composite2COP
function could be used directly in composite1COP
. Albeit, the second copula and its parameters would not be used. A more complex (extended) composition in composite3COP
extends this basic parameter structure.
Author(s)
W.H. Asquith
References
Joe, H., 2017, Parametric copula families for statistical models (chap. 8) in Copulas and dependence models with applications—Contributions in honor of Roger B. Nelsen, eds. Flores, U.M., Amo Artero, E., Durante, F., Sánchez, J.F.: Springer, Cham, Switzerland, ISBN 978–3–319–64220–9, doi:10.1007/978-3-319-64221-5.
Salvadori, G., De Michele, C., Kottegoda, N.T., and Rosso, R., 2007, Extremes in Nature—An approach using copulas: Springer, 289 p.
See Also
COP
, breveCOP
, composite2COP
, composite3COP
,
convexCOP
, glueCOP
Examples
## Not run:
alpha <- 0.24; beta <- 0.23; Theta1 <- NA;
# W() does not use a parameter, but show how a parameter would be set if needed.
para <- list(alpha=alpha, beta=beta, cop1=W, para1=Theta1)
t <- composite1COP(0.4, 0.6, para)
if( t != W(0.4, 0.6)) message("Not equal as expected")
# Next use this as a chance to check logic flow through the various
# "compositing" operators and their as needed dispatch to COP().
my.para <- list(cop1=GHcop, para1=exp(+1.098612) + 1,
cop2=PLcop, para2=exp(-1.203973),
alpha=0.5, beta=0.25,
kappa=0.1, gamma=0.1,
weights=c(0.95, 0.05))
# uses cop1/2, para1/2, only weights
nustarCOP(cop=convexCOP, para=my.para) # 0.8570434
# uses cop1/2, para1/2, only alpha
nustarCOP(cop=convex2COP, para=my.para) # 0.2697063
# uses cop1, para1, only alpha / beta
nustarCOP(cop=composite1COP, para=my.para) # 0.5103119
# uses cop1/2, para1/2, only alpha / beta
nustarCOP(cop=composite2COP, para=my.para) # 0.0714571
# uses cop1/2, para1/2, only alpha, beta, kappa, gamma
nustarCOP(cop=composite3COP, para=my.para) # 0.0792634
## End(Not run)