copstress function

Calculates copstress for given MDS object

copstress(
  obj,
  stressweight = 1,
  cordweight = 5,
  q = 1,
  minpts = 2,
  epsilon = 10,
  rang = NULL,
  verbose = 0,
  normed = TRUE,
  scale = c("std", "sd", "proc", "none"),
  init,
  ...
)

Arguments

• obj: MDS object (supported are sammon, cmdscale, smacof, rstress, powermds)
• stressweight: weight to be used for the fit measure; defaults to 1
• cordweight: weight to be used for the cordillera; defaults to 0.5
• q: the norm of the cordillera; defaults to 1
• minpts: the minimum points to make up a cluster in OPTICS; defaults to 2
• epsilon: the epsilon parameter of OPTICS, the neighbourhood that is checked; defaults to 10
• rang: 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.
• verbose: numeric value hat prints information on the fitting process; >2 is very verbose (copstress level), >3 is extremely (up to MDS optimization level)
• normed: should the cordillera be normed; defaults to TRUE
• scale: should the configuration be scale adjusted.
• init: a reference configuration when doing procrustes adjustment
• ...: additional arguments to be passed to the cordillera function

Returns

A list with the components

• copstress: the weighted loss value
• OC: the Optics cordillera value
• parameters: the parameters used for fitting (kappa, lambda)
• cordillera: the cordillera object