A general approach to generate random variates for multivariate copulae

N. Tajvidi, B. A. Turlach

Research output: Contribution to journalArticlepeer-review

4 Citations (Scopus)


We suggest two methods for simulating from a multivariate copula in an arbitrary dimension. Although our main emphasis in this paper is on multivariate extreme value distributions, the proposed methods can be applied to any copula. The basic idea is to approximate the (unknown) density of the copula by a distribution that has a piece-wise constant (histogram) density. This is achieved by partitioning the support of a given copula C into a large number of hyper-rectangles and using them to generate random variates from an approximation of the copula. We suggest two methods for finding this approximation which correspond to either finding hyper-rectangles which have equal probability mass with respect to C, or determining a partition using hyper-squares of equal volume and finding the corresponding probability mass of each hyper-square. We also discuss how the generated random variates can be used as proposals in a Metropolis–Hastings algorithm, when C is an absolutely continuous distribution function, to generate a sequence of random variates from C. An implementation of the proposed methodologies is provided for the statistical computing and graphics environment R in our package called SimCop.

Original languageEnglish
Pages (from-to)140-155
Number of pages16
JournalAustralian and New Zealand Journal of Statistics
Issue number1
Publication statusPublished - 1 Mar 2018


Dive into the research topics of 'A general approach to generate random variates for multivariate copulae'. Together they form a unique fingerprint.

Cite this