MNM (version 0.95-0)

plotShape: Pairwise Scatterplot Matrix of Shape Matrices

Description

Function for visual comparisons for up to three shape matrices.

Usage

plotShape(est1, est2 = NULL, est3 = NULL, X = NULL, alim = NULL, 
          color.ell = 2:4, color.points = grey(0.5), lty.ell = rep(1, 3), 
          pch.ell = rep(16, 3), pch.points = 1, level = 0.5, npoints = 100, 
          x.legend, y.legend, cex.legend = 1, pty = "s", gap = 1, 
          oma.bottom, labels, cex.labels = 2, main, ...)

Value

• A scatter plot matrix.

Details

All scatter matrices are standardized to have determinant 1. If X is given, the Mahalanobis distances based on the location and shape estimates are computed, and t in the function ellipse is the level quantile of the Mahalanobis distances. If no X is provided t equals level. The location of the legend is currently problematic and it is recommended that the user should provide the coordinates for the legend.

See Also

ellipse, plotMvloc

Examples

X <- rmvt(100, diag(3), df=3)

EST1 <- list(location=colMeans(X), scatter=cov(X), est.name="COV")
HR <- HR.Mest(X)
EST2 <- list(location=HR$center, scatter=HR$scatter, est.name="Tyler")
plotShape(EST1,EST2, X=X)