## Notice that, for a given degree of dependence,
## these copulas are strikingly similar.
tau <- 0.33
gumbel.cop <- evCopula("gumbel", param=2)
gumbel.cop@parameters <- calibKendallsTau(gumbel.cop, tau)
galambos.cop <- galambosCopula(2)
galambos.cop@parameters <- calibKendallsTau(galambos.cop, tau)
huslerReiss.cop <- huslerReissCopula(2)
huslerReiss.cop@parameters <- calibKendallsTau(huslerReiss.cop, tau)
tawn.cop <- tawnCopula(0.3)
tawn.cop@parameters <- calibKendallsTau(tawn.cop, tau)
tev.cop <- tevCopula(0.7)
tev.cop@parameters[1] <- calibKendallsTau(tev.cop, tau)
curve(Afun(gumbel.cop, x), 0, 1)
curve(Afun(galambos.cop, x), lty=2, add=TRUE)
curve(Afun(huslerReiss.cop, x), lty=3, add=TRUE)
curve(Afun(tawn.cop, x), lty=4, add=TRUE)
curve(Afun(tev.cop, x), lty=5, add=TRUE)
Run the code above in your browser using DataLab