## For 'matrix' objects
cop <- gumbelCopula(2, dim = 3)
n <- 1000
set.seed(271)
U <- rCopula(n, copula = cop)
cloud2(U, xlab = expression(U[1]), ylab = expression(U[2]), zlab = expression(U[3]))
## For 'Copula' objects
set.seed(271)
cloud2(cop, n = n) # same as above
## For 'mvdc' objects
mvNN <- mvdc(cop, c("norm", "norm", "exp"),
list(list(mean = 0, sd = 1), list(mean = 1), list(rate = 2)))
cloud2(mvNN, n = n)
Run the code above in your browser using DataLab