This model was introduced in Coles and Tawn (1991); the present method uses the simulation algorithm of Boldi (2009) for the extremal Dirichlet model
.rdirspec(n, d, alpha, irv = FALSE)an n by d sample from the spectral distribution
sample size
dimension of sample
vector of Dirichlet parameters of dimension d, or \(d+1\) vector with the d Dirichlet parameters and an index of regular variation in \([0, 1]\)
should the usual model (FALSE) or the general scaled version (TRUE) be used
Boldi (2009). A note on the representation of parametric models for multivariate extremes. Extremes 12, 211--218.