Simulate the r
-largest observations from a Poisson point process with intensity
$$\Lambda(x) = (1+\xi(x-\mu)/\sigma)^{-1/\xi}$$.
rrlarg(n, r, loc = 0, scale = 1, shape = 0)
an n
by r
matrix of samples from the point process, ordered from largest to smallest in each row.
sample size
number of observations per block
location parameter
scale parameter
shape parameter