copBasic (version 2.1.4)

composite1COP: 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 \(\mathbf{C}(u,v)\) is a symmetric copula (see COP) with parameters \(\Theta\) and \(\mathbf{C} \ne \mathbf{\Pi}\) (for \(\mathbf{\Pi}\) see P), then a family of generally asymmetric copulas \(\mathbf{C}_{\alpha,\beta; \Theta}\) with two compositing parameters \(0 < \alpha,\beta < 1, \mbox{and\ } \alpha \ne \beta\), which also includes just the copula \(\mathbf{C}(u,v)\) as a limiting case for \(\alpha = \beta = 0\) and is given by

$$\mathbf{C}_{\alpha,\beta}(u,v) = u^\alpha v^\beta \cdot \mathbf{C}(u^{1-\alpha},v^{1-\beta})\mbox{.}$$

The composite1COP function provides the means for permutation asymmetry from a permutation symmetric copula as also 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 \(\alpha\) or \(\beta\) is held at zero. Very loosely, the \(\alpha > 0\) kicks probability density down towards the lower right corner, whereas \(\beta > 0\) 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 \(u\) in the \(X\) direction;

v

Nonexceedance probability \(v\) in the \(Y\) direction;

para

A special parameter list (see Note); and

...

Additional arguments to pass to the copula.

Value

Value(s) for the composited copula are returned.

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<U+00E1>nchez, J.F.: Springer, Cham, Switzerland, ISBN 978--3--319--64220--9.

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, composite2COP, composite3COP

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=.5, beta=0.25,
                kappa=.1, gamma=.1,
                weights=c(0.95, 0.05))
# uses cop1/2, para1/2, only weights
nustarCOP(cop=convexCOP,     para=my.para)

# uses cop1/2, para1/2, only alpha
nustarCOP(cop=convex2COP,    para=my.para)

# uses cop1,   para1,   only alpha/beta
nustarCOP(cop=composite1COP, para=my.para)

# uses cop1/2, para1/2, only alpha/beta
nustarCOP(cop=composite2COP, para=my.para)

# uses cop1/2, para1/2, only alpha,beta,kappa,gamma
nustarCOP(cop=composite3COP, para=my.para)
# }