# a hac obejct is defined
M = matrix(c(2, 0, 0, 0, 3, 4, 0, 0, ~X1 + X5, ~X2, ~X3, ~X4 + X6),
ncol = 4, byrow = TRUE)
g.model = hac(HAC_GUMBEL, M)
# sample from copula g.model
sample = rHAC(100, g.model)
# cdf at each vector of the sample
values = pHAC(sample, g.model)
Run the code above in your browser using DataLab