level.curvesCOP(cop=P, para=NULL, delt=.02) # Product copula
level.curvesCOP(cop=PSP, para=NULL) # PSP copula
# Produce the level curves for the Plackett copula for
# a theta parameter < 0 (positive association)
para <- c(0.0001)
level.curvesCOP(cop=PLACKETTcop, para=para)
# Produce the level curves for the Plackett copula for
# a theta parameter > 0 (positive association)
para <- c(10)
level.curvesCOP(cop=PLACKETTcop, para=para)
# Now create a composite copula
para <- list(cop1=PLACKETTcop, cop2=PSP,
para1=c(0.4), para2=NULL,
alpha=0.6, beta=0.09)
level.curvesCOP(cop=composite2COP, para=para, delt=.1)
Run the code above in your browser using DataLab