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:
'liouv' denote cores, prefixes
'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
liouv.maxim.mm. 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.
K.plot by Pr. Christian Genest.
Maintainer: Leo Belzile <email@example.com>
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.