copBasic (version 1.5.1)

derCOPinv2: The Inverse of a Numeric Derivative for U with respect to V of a Copula

Description

Compute the inverse of a numeric derivative for $U$ with respect to $V$ of a copula or $$t = c_v(u) = \frac{\delta \mathbf{C}(u,v)}{\delta v}\mbox{,}$$ and solving for $u$.

Usage

derCOPinv2(cop=NULL, v, t, delv=.Machine$double.eps^0.50,
           para=NULL, ...)

Arguments

cop

A copula function,

v

A nonexceedance probability in $Y$ direction,

t

A nonexceedance probability in $t$ direction,

delv

The $\delta v$ interval for the derivative,

para

Vector of parameters or other data structure, if needed, to pass to the copula, and

...

Additional arguments to pass.

Value

• The value for the copula is returned.

References

Nelson, R.B., 2006, An introduction to copulas: New York, Springer, 269 p.

See Also

derCOP2, W, P, M, PSP

Examples

v <- runif(1)
t <- runif(1)
derCOPinv2(cop=W,v,t)
derCOPinv2(cop=P,v,t)
derCOPinv2(cop=M,v,t)
derCOPinv2(cop=PSP,v,t)

# Simulate 1000 values from product (independent)
# copula
plot(c(0,1),c(0,1),type="n")
for(i in 1:1000) {
  v <- runif(1); t <- runif(1)
  u <- derCOPinv2(cop=P,v,t)
  points(u,v,cex=0.5,pch=16)
}
# Now simulate 1000 from the PSP copula and note
# the lower tail dependence, but absence of upper
# tail dependence.
for(i in 1:1000) {
  v <- runif(1); t <- runif(1)
  u <- derCOPinv2(cop=PSP,v,t)
  points(u,v,cex=0.5,pch=16,col=2)
}