PAMK fits PAM (or CLARA) for every candidate in krange and keeps the one
scoring best under criterion. There is therefore no k to set, and the
number of clusters chosen is reported by the resulting Clustering.
Argument names are rtemis' own. They map to fpc::pamk() as:
use_pam -> usepam; the criterion's type -> criterion and its
n_subsets -> ns; the rest keep their names.
fpc::pamk() also accepts a per-variable numeric vector for scaling; only
the logical form is exposed here. Scale by column with preprocess instead,
where the decision is recorded and replayed with the data.