# set up 4-dimensional D-vine in R-Vine notation
d <- 4
dd <- d * (d-1) / 2
order <- 1:d
family <- c(1, 2, 3, 4, 7, 3)
par <- c(0.5, 0.4, 2, 1.5, 1.2, 1.5)
par2 <- c(0, 5, 0, 0, 2, 0)
RVM <- D2RVine(order, family, par, par2)
# simulate from C-vine model
simdata <- RVineSim(500, RVM)
# calculate likelihood
out2 <- RVineLogLik(simdata, RVM)
out2$loglik
Run the code above in your browser using DataLab