The function is used internally for optimization.
.pliouv.opt_old(theta, data, family, alphavec, MC.approx = TRUE)value of marginal density
parameter of the corresponding Archimedean copula
sample matrix from a Liouville copula
family of the Liouville copula. Either "clayton", "gumbel", "frank", "AMH" or "joe"
vector of Dirichlet allocations (must be a vector of integers)
whether to use Monte-Carlo approximation for the inverse survival function (default is TRUE)