The score is computed via linear interpolation for the shape parameter in a neighborhood of zero
rlarg.score(par, dat)
par
vector of loc, scale and shape
loc
scale
shape
dat
an n by r sample matrix, ordered from largest to smallest in each row
n
r
score vector of size 3
rlarg