library(copula)
data("SEAex")
cop <- normalCopula(param = 0.5, dim = 4)
distList <- c("RG", "SN1", "RG", "RG")
allDistrs <- list(list(mu = 31.07, sigma = 0.28),
list(mu = 34.4, sigma = 0.98, nu = 1.7),
list(mu = 31.4, sigma = 0.34),
list(mu = 25.6, sigma = 0.24))
copSEA <- mvdc(cop, distList, allDistrs)
pbx <- make_pbox(data = SEAex, cop = copSEA)
print(class(pbx))
Run the code above in your browser using DataLab