n = 800
Z = stats::runif(n = n)
CKT = 0.2 * as.numeric(Z <= 0.3) +
0.5 * as.numeric(Z > 0.3 & Z <= 0.5) +
- 0.8 * as.numeric(Z > 0.5)
simCopula = VineCopula::BiCopSim(N = n,
par = VineCopula::BiCopTau2Par(CKT, family = 1), family = 1)
X1 = simCopula[,1]
X2 = simCopula[,2]
partition = cbind(Z <= 0.3, Z > 0.3 & Z <= 0.5, Z > 0.5)
condPseudoObs = bCond.pobs(X = cbind(X1, X2),
partition = partition)
Run the code above in your browser using DataLab