PCOPS version of elastic scaling with powers
cop_powerelastic(
dis,
theta = c(1, 1),
weightmat = 1 - diag(nrow(dis)),
init = NULL,
ndim = 2,
itmaxi = 10000,
...,
stressweight = 1,
cordweight = 0.5,
q = 1,
minpts = ndim + 1,
epsilon = 10,
rang = NULL,
verbose = 0,
scale = "sd",
normed = TRUE,
stresstype = c("default", "stress1", "rawstress", "normstress", "enormstress",
"enormstress1")
)A list with the components
stress: the stress
stress.m: default normalized stress
copstress: the weighted loss value
OC: the Optics cordillera value
parameters: the parameters used for fitting (kappa, lambda)
fit: the returned object of the fitting procedure
cordillera: the cordillera object
numeric matrix or dist object of a matrix of proximities
the theta vector of powers; a vector of length two where the first element is kappa (for the fitted distances), the second lambda (for the observed proximities). If a scalar for the free parameters is given it is recycled. Defaults to 1 1.
(optional) a matrix of nonnegative weights
(optional) initial configuration
number of dimensions of the target space
number of iterations. default is 10000.
additional arguments to be passed to the fitting procedure
weight to be used for the fit measure; defaults to 1
weight to be used for the cordillera; defaults to 0.5
the norm of the cordillera; defaults to 1
the minimum points to make up a cluster in OPTICS; defaults to ndim+1
the epsilon parameter of OPTICS, the neighbourhood that is checked; defaults to 10
range of the distances (min distance minus max distance). If NULL (default) the cordillera will be normed to each configuration's maximum distance, so an absolute value of goodness-of-clusteredness.
numeric value hat prints information on the fitting process; >2 is extremely verbose
should the configuration be scale adjusted
should the cordillera be normed; defaults to TRUE
which stress to report? Defaults to explicitly normed stress