lcopula-package: Liouville Copulas

Description Author(s) References See Also


lcopula provides functions to evaluate and sample Liouville copulas. The code adapted from GumbelLiouville.R used in Genest and Neslehova (2013) and from the routine LiouvilleFunction.R used by McNeil and Neslehova (2010) and coded by Alexander McNeil. New implementation draws heavily on functions implemented in the copula package and derived from work by Hofert, Machler and McNeil (2012).

Warning: a different definition is used by the latter and in the copula package for the Clayton copula. The implementation below uses the generator definition found in the first edition of the book Quantitative Risk Management.

The naming conventions are as follows: 'archi' and 'liouv' denote cores, prefixes 'd', 'p', 'r', 's', 'is' denote respectively density, distribution function, random number generation, survival function and inverse survival function. Suffixes are 'm' for marginal, '_m' for simultaneous multiple margins, '_p' for vectorized versions of the functions.

The functions mostly draws from the theory laid out in McNeil and Neslehova (2010) and Belzile (2014), which are the references for definitions and expressions.

The Archimedean families implemented are Clayton, Gumbel, Frank, Ali-Mikhail-Haq (abbreviated AMH) and Joe. Random number generation from the copula is available via the function rliouv. Method-of-moment estimators can be used for fitting for the Clayton and Gumbel families for bivariate samples, using the function More general models can be fitted using liouv.maxim, but optimization may be slow. The function theta.bci allows to construct bootstrap confidence interval for the parameter theta. Unless the method-of-moment estimator can be used, it is very slow since it relies on full optimization of each replicated sample under the model.

The copula and survival copula domains of attraction are implemented and the corresponding spectral density and Pickands dependence function can be computed and plotted in the bivariate case.


Leo Belzile, with adapted routines from Pr. Alexander J McNeil, Pr. Johanna G Neslehova and Pr. Christian Genest. Function K.plot by Pr. Christian Genest.

Maintainer: Leo Belzile <>


McNeil A.J. and Neslehova, J.G. (2010) From Archimedean to Liouville Copulas. J. Multivar. Anal., 101(8): 1772–1790.

Belzile L. (2014) Extremal and inferential properties of Liouville copulas. Master thesis, McGill.

Hofert, M., Machler, M., and McNeil, A. J. (2012) Likelihood inference for Archimedean copulas in high dimensions under known margins. J. Multivar. Anal., 110, 133–150.

Genest, C. and Neslehova, J. G. (2013) Assessing and Modeling Asymmetry in Bivariate Continuous Data. In P. Jaworski, F. Durante, and W. K. Hardle (Eds.), Copulae in Mathematical and Quantitative Finance, Lecture Notes in Statistics, 91–114, Springer: Berlin Heidelberg.

See Also

copula package

Search within the lcopula package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.