Function returning the density of an n
sample from the GP distribution.
gpd.ll.optim
returns the negative log likelihood parametrized in terms of log(scale)
and shape
in order to perform unconstrained optimization. The function is coded in such a way that the log likelihood is infinite when \(\xi < -1\).