norm.cop <- normalCopula(c(0.5, 0.6, 0.7), dim = 3, dispstr = "un")
t.cop <- tCopula(c(0.5, 0.3), dim = 3, dispstr = "toep",
df = 2, df.fixed = TRUE)
## from the wrapper
norm.cop <- ellipCopula("normal", param = c(0.5, 0.6, 0.7),
dim = 3, dispstr = "un")
if(require("scatterplot3d")) {
## 3d scatter plot of 1000 random observations
scatterplot3d(rcopula(norm.cop, 1000))
scatterplot3d(rcopula(t.cop, 1000))
}
Run the code above in your browser using DataLab