vgislw(lw, wcp = 0.2, thresh = 100, kmax = 2, wtrunc = 3/4,
cores = parallel::detectCores())lw =
-log_lik (see extract_log_lik) and is an $S$ by $N$
matrix where $S$ is the number of simulations and $N$ is the number
of 100*wcp% largest weights are used as the sample
from which to estimate the parameters of the generalized Pareto
distribution.100*wcp% of the importance weights are modified according to
pmax(x, max(x) - thresh) for numerical stability.wtrunc. Set
to zero for no truncation.lw_smooth (modified log
weights) and pareto_k (estimated generalized Pareto shape parameters
$k$).loo-package.loo-package, loo_and_waic