smacof (version 1.8-13)

jackknife: SMACOF Jackknife

Description

These methods perform a SMACOF Jackknife and plot the corresponding solution.

Usage

"jackknife"(object, eps = 1e-6, itmax = 100, verbose = FALSE)
"plot"(x, plot.dim = c(1,2), hclpar = list(c = 50, l = 70), 
plot.lines = TRUE, main, xlab, ylab, xlim, ylim, ...)

Arguments

object

Object of class "smacofB", i.e., an MDS solution from smacofSym()

itmax

Maximum number of iterations

eps

Convergence criterion

verbose

If TRUE, intermediate stress is printed out

x

Object of class "smacofJK"

plot.dim

Vector with dimensions to be plotted.

hclpar

Chroma and luminance to be used for HCL colors (further details see rainbow_hcl)

plot.lines

If TRUE, the Jackknife configurations are plotted are connected with their centroid

main

Plot title.

xlab

Label of x-axis.

ylab

Label of y-axis.

xlim

Scale x-axis.

ylim

Scale y-axis.

...

Further plot arguments passed: see plot for detailed information.

Value

smacof.conf

SMACOF configurations

jackknife.conf

An array of n-1 configuration matrices for each Jackknife MDS solution

comparison.conf

Centroid Jackknife configurations (comparison matrix)

stab

Stability measure

cross

Cross validity

disp

Dispersion

loss

Value of the loss function

ndim

Number of dimensions

call

Model call

niter

Number of iterations

nobj

Number of objects

Details

In order to examine the stability solution of an MDS, a Jackknife on the configurations can be performed (see de Leeuw & Meulman, 1986) and plotted. The plot shows the jackknife configurations which are connected to their centroid. In addition, the original smacof configuration (transformed through Procrustes) is plotted. The Jackknife function itself returns also a stability measure (as ratio of between and total variance), a measure for cross validity, and the dispersion around the original smacof solution.

References

De Leeuw, J., & Meulman, J. (1986). A special jackknife for multidimensional scaling. Journal of Classification, 3, 97-112.

See Also

smacofSym, smacofConstraint, plot.smacof

Examples

## symmetric smacof
data(kinshipdelta)
res.smacof <- smacofSym(kinshipdelta)
res.jk <- jackknife(res.smacof)
res.jk

plot(res.jk)
plot(res.jk, hclpar = list(c = 80, l = 40))
plot(res.jk, hclpar = list(c = 80, l = 40), plot.lines = FALSE)